using label="mean"
Posted on
wanlin
Joined: 10/28/2009
Forums
Dear all
When set
mxMatrix(type="Full", nrow=1, ncol=2, free=TRUE, values=c(0,0), label= "mean", name="expMean")
in twin ACE model, i found if i set label="mean", the
twinACEFit@output$estimate return 4 values with label (mean, a,c,e)
if the label item not set
mxMatrix(type="Full", nrow=1, ncol=2, free=TRUE, values=c(0,0), name="expMean")
When set
mxMatrix(type="Full", nrow=1, ncol=2, free=TRUE, values=c(0,0), label= "mean", name="expMean")
in twin ACE model, i found if i set label="mean", the
twinACEFit@output$estimate return 4 values with label (mean, a,c,e)
if the label item not set
mxMatrix(type="Full", nrow=1, ncol=2, free=TRUE, values=c(0,0), name="expMean")
twinACEFit@output$estimate return 5 values with label (NA,NA, a,c,e) and the first two mean values are not exactly the same.
Which one is correct?
Wanlin Which is correct
Which is correct depends on which model you wanted to fit. The first constrains the means to be equal. Only one label was provided ("mean") which was automatically re-used by R, so both elements of expMean were set to the same free parameter. This is the usual way of specifying ACE models (constrain means of twin1 and twin2 to be equal). The second approach is unconventional, allowing different means for twin1 and twin2. Sometimes - typically when using raw data - this approach is used as part of a saturated model so that the fit of the ACE model can be compared to it.
Log in or register to post comments
In reply to Wanlin Which is correct by neale
1. Thanks for clarification
2. if the variance of twin1 and twin2 is not equal, could i standardize both twin1 and twin2 to make them have equal mean=0 and equal variance = 1?
3. As i don't know how to combine the regression into the ACE model, all my data were adjusted age and sex using glm before input to SEM, is it appropriate?
Log in or register to post comments
In reply to 1. Thanks for clarification by wanlin
Ah, brain structure data,
Unequal Twin 1 and Twin 2 variances should only happen by chance, if their ordering is random. If it is not (birth order perhaps) then it may be worth seeing if they cause loss of fit. Any loss of fit above expectations (if you run 20 ROI's then you would expect 1 test to be signficant at .05 after all) may indicate birth order effects which seem unlikely to me. Much more likely is failure of the assumption of multivariate normality. It may be best to hunt for outliers.
I would not standardize twin 1 and twin 2 just to make the model fit. Model failure is worth knowing about and investigating. It should not be swept under the carpet.
Log in or register to post comments
In reply to Ah, brain structure data, by neale
Thanks neale. If the
If the means and variances of Twin1 and Twin2 for both MZ and DZ are not significantly different. Can i standardize them.
Some times, i found if the data were standardized, the ACE model will fit the saturated model very well with a fixed initial value.
Another question, sometimes i find the -2 log likelihood is minus, but results from examples or papers always report positive likelihood, should the value always positive?
Log in or register to post comments
In reply to Thanks neale. If the by wanlin
No, it is not a good idea to
-2lnL can be either a positive or a negative value. Usually, most of the data points have a likelihood less than 1, so the log-likelihoods are negative, and the -2lnL's are positive. However, with normal theory (FIML) analysis with variables with small variances, the likelihoods can often exceed 1.0, so the log-likelihoods are positive, and the -2lnL's are negative.
While we're on the topic, it's an idea to avoid variables with very small variances, simply because during estimation trial values for a variance of zero necessarily cannot be evaluated (likelihood is zero, -2lnL is negative infinity).
Log in or register to post comments