Attachment | Size |
---|---|
SingleFactor24.R [6] | 15.98 KB |
SingleFactorCohortManifestFixed24.R [7] | 15.99 KB |
SingleFactorNoCohort24.R [8] | 14.57 KB |
SingleFactorCohortManifest24.R [9] | 15.97 KB |
CohortACE.csv [10] | 1.48 KB |
Hi all,
I've been playing around with various ways to include a cohort covariate in a single-factor model. I've tried four approaches, three of which have pretty consistent results and the fourth, which is fairly different. I've attached the four scripts here, as well as a table of selected output (ACE estimates with CIs, factor loadings, and expected means for observed variables).
For all models, the factor is identified by fixing the variance to 1 and there are four observed drug use variables (AlcFreq, AlcQuant, CPD_TOB, and MJ_FREQ). Models are fit separately by sex and by zygosity (MZ males, MZ females, DZ males, DZ females). I have three intake cohorts, indexed by two dummy variables OC and YC.
The first script, SingleFactorNoCohort24.R evaluates a factor model for which cohort is not regressed out at any point.
The second script, SingleFactor24.R, evaluates a factor model for which the mean of the factor is fixed to 0 and cohort is regressed onto the factor.
The third script, SingleFactorCohortManifest24.R, evaluates a factor model for which cohort is regressed onto the manifest variables and the effect of cohort is freely estimated for each of the four variables.
The fourth script, SingleFactorCohortManifestFixed24.R, evaluates a factor model for which cohort is regressed onto the manifest variables and the effect of cohort is fixed to be the same across drugs.
The part of the output that I am finding confusing is how different the ACE estimates are for the SingleFactor24.R script where cohort is regressed on the factor compared to the other 3. The C effect is large (.62 for males, .68 for females) in this version, whereas none of the other 3 show much of a C effect (0-.09) and I am generally surprised when I see a C effect of that magnitude at that developmental point (~24 years old). There also seem to be some differences in the factor loadings.
Does anyone have any insight as to why the ACE estimates would be so different when cohort is regressed on the factor rather than on the manifest variables or not at all?
Thanks!
-Stephanie