Hi,

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)),

cbind((ra%x%(af%*%t(am)))+cf%*%t(cm), Af+Cf+Ef))

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!

Regards

Charlotte