Message/Author 


In the presentation of ex5.21 for Mplus 4.1, was "h" meant to be heretability? Thanks 


Yes, it is. 


I'm trying to estimate a twin model that has two groups, but I need to estimate the effect of a covariate on y1 and y2. Specifically, I'm trying to estimate example 5.18 in the new user's manual. I'm trying to augment this model so that it is a MACE model. When I just add the covariate on y1 and y2, I can't seem to identify the model. That is, the model will not converge with a chisquare. Any thoughts as to why and how I can identify my model? 


To understand this, you would have to send the output and license number to support@statmodel.com. 


For CFA (perhaps also multiple group CFA), and when using a sample of twins, both MZs and DZs, is it possible to use the stratification option along with type = complex to handle nonindependence of twins in sample? 


The CLUSTER option deals with nonindependence of observations not STRATIFICATION. Nonindependence of twins is already handled in ACE twin models because the two twins are analyzed together in a multivariate model that models the nonindependence. 


Thanks Linda. I haven't considered the ACE approach because I've not seen examples with twin data using latent variables, especially more than two LVs. In my case I am also doing longitudinal CFA and SEM. 


what would be really helpful is if there was an example of a latent variable crosslag panel model (2 LVs across 6 year period, LV1 & LV2 at T1 predict LV1 & LV2 at T2) using only male MZs, and Dzs. 


You might find something helpful at: http://www.statmodel.com/geneticstopic.shtml 

Dan R. posted on Friday, September 18, 2015  9:06 am



I'm trying to extract factor scores (for individuals) from latent variables created from crosssectional data that includes family clusters with MZ & DZ twin pairs and siblings within those families. To make it a little more interesting, there are also some singletons. The intention is to create factor scores that control/adjust for the nonindependence of not only twins but the family structures that the observations are nested within. I'm not interested (I don't think) in factor structure differences on within and between levels, just a TYPE = COMPLEX model. I've create variables for 'family' and 'zygosity', but I'm not sure if CLUSTER = family zygosity; is what I'm looking for. Would saving the between level factor scores from a model using the following commands accomplish those goals? TYPE = TWOLEVEL COMPLEX; CLUSTER = family zygosity; Thank you for any clarification of these commands that you can provide, or relevant resources to learn more. 


I think you want to have Cluster=family, so not include zygosity because the usual twin approach has zygosity as a Grouping variable, that is, a fixed mode of variation. See our UG for twin examples. Now, if you use Type=Complex with this you will get the same factor scores as without Complex because Complex affects only the SEs of the model parameters. If you do a Type=Twolevel analysis you would get different withinlevel factor scores than in a Type=Complex or regular analysis. The betweenlevel would refer to the family component and not the twin component. 

Dan R. posted on Sunday, September 20, 2015  11:59 am



Thank you, Bengt. Would Cluster=family take each families' clustering/similarity into account (i.e., on an individual family level), or does it do this by findings an average degree of similarity within clusters (an ICC for example) and applying that to individual's data when they come from a clustered group? I guess I'm just looking for more information on how the cluster function works and what it does to the results. I couldn't find much beyond how to use the cluster function in the UG (though I could have missed a relevant section), so some direction to more technical details would be greatly appreciated. Thanks! 

Back to top 