95% CI of A,C,E estimates in Twin model

5 replies [Last post]
Sabha's picture
Offline
Joined: 05/18/2010

Can anyone help me with the function or syntex needed to calculate the 95% CI for A,C,E estimates in the existing
UnivariateTwinAnalysis_MatrixRaw.R script.

Thanks in advance.

neale's picture
Offline
Joined: 07/31/2009
See this part of the wiki:
Sabha's picture
Offline
Joined: 05/18/2010
Thanks for your help. I was

Thanks for your help. I was able to get the CI for A,C,E. I also would like to get the CI for standardized components ( e.g ACE.A/V,ACE.C/ACE.V,ACE.E/ACE.V). How do I compute it.

tbates's picture
Offline
Joined: 07/31/2009
To add to mikes reply: You

To add to mikes reply: You will need to create these standardizing algebras in your model, if they are not already present.

i.e., you can't just ask for mxCI("ACE.A/ACE.Vtot"), you have to create
mxAlgebra(ACE.A/ACE.Vtot, name="stdA")
mxCI(c('stdA')

On that note, wouldn't it be great if you could just include things in the mxCI statement and they would be automagically created!!

All the information is there.

Sabha's picture
Offline
Joined: 05/18/2010
You spelled my wish.

You spelled my wish. Initially, I did try that way but it didn't work. Thanks a lot to both of you. It worked. Great help.

neale's picture
Offline
Joined: 07/31/2009
Same way, just use the names

Same way, just use the names of the algebras instead of the names of the matrices in the list of elements for which you would like the CIs.