Attachment | Size |
---|---|
NTF.ASDE_.R [6] | 8.93 KB |
Hi all,
I'm running a nuclear twin family design model (an SEM model that contains a number of non-linear constraints). Two issues:
(1) The point estimate of the pathways are fine - and agree with those I simulated - but the standard errors are WAY off; e.g., the standard error for a path coefficient that cannot be dropped (e) is 30.3 and the point estimate is -.54! Here are the other coefficients from summary():
summary(ASDE.Fit)
name matrix row col Estimate Std.Error
1 AddGenPath MZNTF.a 1 1 0.63320641 1.9272561
2 DomPath MZNTF.d 1 1 -0.31596212 0.4459057
3 EnvPath MZNTF.e 1 1 -0.54266100 30.3050451
4 AMCopath MZNTF.mu 1 1 0.21215070 0.8786654
5 SibPath MZNTF.s 1 1 0.38618256 1.2540086
6 VarPhen MZNTF.Vp1 1 1 0.98593865 0.1272815
7 CovPhenGen MZNTF.delta1 1 1 0.69880636 0.3767022
8 LatentVarAddGen MZNTF.q1 1 1 1.10359962 2.9363220
9 mean MZNTF.expMeanMz 1 Tw1 -0.01938465 2.4590259
(2) Less critically, why do the observed statistics change when I drop a parameter? The difference in degrees of freedom between two models where a single parameter was fixed is "4" according to OpenMx.
sessionInfo()
R version 2.10.1 (2009-12-14)
x86_64-apple-darwin9.8.0
locale:
[1] C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] OpenMx_0.2.5-1050
loaded via a namespace (and not attached):
[1] tools_2.10.1
I've attached the script I'm working on. The code is self-contained. Any help with the two questions above would be appreciated.
Matt