tbates / umx

Making Structural Equation Modeling (SEM) in R quick & powerful
https://tbates.github.io/
44 stars 17 forks source link

umxSummaryACE not handling second covariate #201

Closed tbates closed 1 year ago

tbates commented 1 year ago

h/t Eren Kafadar for reporting!

require(umx)
data(twinData)
twinData= umx_scale_wide_twin_data(data=twinData, varsToScale= "ht", sep="")
twinData$cohort1 = twinData$cohort2 =twinData$part
mzData = twinData[twinData$zygosity %in% "MZFF",]
dzData = twinData[twinData$zygosity %in% "DZFF",]

m1 = umxACE(selDVs = 'ht', selCovs = c("age","cohort"), sep = "", dzData = dzData, mzData = mzData)

Table: Standardized parameter estimates from a 1-factor Cholesky ACE model. A: additive genetic; C: common environment; E: unique environment.

a1 c1 e1
ht 0.93 0.08 0.36

Error in dimnames(x) <- dn : length of 'dimnames' [1] not equal to array extent

tbates commented 1 year ago

fixed

require(umx)
data(twinData)
twinData= umx_scale_wide_twin_data(data=twinData, varsToScale= "ht", sep="")
twinData$cohort1 = twinData$cohort2 =twinData$part
mzData = twinData[twinData$zygosity %in% "MZFF",]
dzData = twinData[twinData$zygosity %in% "DZFF",]

m1 = umxACE(selDVs = 'ht', selCovs = c("age","cohort"), sep = "", dzData = dzData, mzData = mzData)
1 row(s) dropped from 'data' due to missing definition variable(s)

Table: Standardized parameter estimates from a 1-factor Cholesky ACE model. A: additive genetic; C: common environment; E: unique environment.

a1 c1 e1
ht 0.933 . 0.36

Table: Means and (raw) betas from model$top$intercept and model$top$meansBetas

ht1 ht2
intercept -0.256 -0.256
age -0.006 -0.006
cohort -0.048 -0.048