I'm adapting the BIVARIATE Boulder-skript atm and I'm trying to freely estimate the paths connecting the As of the two phenotypes for the DOS-twins. So it should be 0.5 for the DZ twins, but freely estimated for the DOS twins.
With this skrip, I get ONE estimate for both paths:
pathRados <- mxMatrix (type="Full", nrow=1, ncol=1, free=TRUE, values=0.5, label="rados", lbound=-1, ubound=1, name="ra")
expCovDOS <- mxAlgebra(name = "expCovDOS",
expression = rbind (cbind(Am+Cm+Em, (ra%x%(am%%t(af)))+cm%%t(cf)),
But I'd like to have two estimates:
Estimate 1: The path connecting the A of phenotype 1, first twin, with the A of phenotype 1, second twin.
Estimate 2: The path connecting the A of phenotype 2, first twin, with the A of phenotype 2, second twin.
I couldn't find an elegant solution yet - do you have any ideas or perhaps even a working skript?
Any help would be very much appreciated!