Behavioral Genetics Models

hi,
just working on the 'univACEP.R' script in /failing

It gives this error.. which is obscure to me

> twinACEFit <- mxRun(twinACEModel)
Running twinACE
Error in mxRun(twinACEModel) :
REAL() can only be applied to a 'numeric', not a 'character'

Hi,
Is there any guidance on how people should be getting the fit of their models, esp multi group models?
running an ACE type model, the summary doesn't calculate any fit statistics, nor does it seem to know how many observations were being made etc so df is just -est parameters etc.

summary(fit)
         name      matrix row col parameter estimate error estimate
1             all.a_c   1   1       2.384409e-01     0.22324048
2             all.a_c   2   1       4.199091e-01     0.21043198
3             all.a_c   3   1       2.366656e-01     0.17050014

running a tri-bivariate ACE model, and the results can be all over the show (with no convergence warnings), depending on the start values: for instance C wandering from .998 to .054 for one variable.
Any thoughts?
Here are the results, and script fragment

Results starting at a, c, and e at .6

               A1     A2    A3    C1    C2    C3    E1    E2    E3
var1     0.284     NA    NA 0.282    NA    NA 0.434    NA    NA
var2    -0.292  0.027    NA 0.998 0.434    NA 0.294 0.539    NA
var3     0.224 -0.028 0.206 0.366 0.458 0.196 0.411 0.570 0.598
[1] 7550.711

Are there multivariate Twin scripts for pathic and/or matic scripts available to work with?

Anyone have a reference for the first twin model, using the distinction of MZ and DZ similarity?

In writing the wikipedia page, my current understanding is it was Weinberg, with the area being forgotten again or perhaps stifled by the thoughts of Fisher (1915) on (the lack of) MZ twins, to be reinvigorated by Seimens in 1924 with his poly symptom checklist. Thoughts?