Use OpenMx to quickly fit a univariate ACE model
fit_uniACE(
data_1,
data_2,
GroupRel = c(1, 0.5),
GroupR_c = c(1, 1),
lbound = FALSE
)Returns a list with the following:
A data.frame displaying the nested comparison model between ACE, AE, CE, E models
A list of all model fit information generated from OpenMx
A n by 2 data.frame consisting of the group1 kin pairs
A n by 2 data.frame consisting of the group2 kin pairs
A numeric vector specifying two genetic relatedness values of two groups of kin pairs
A numeric vector specifying two common environment correlation coefficients of two groups of kin pairs
A logical value indicating if a lower boundary of .0001 will be imposed to the estimated A, C and E components