Attachment | Size |
---|---|
Test Data ACE Model.csv [6] | 964 bytes |
ACE Skript tested with own Data.R [7] | 6.16 KB |
Hey guys. I am trying to do an ACE Model with my own twin data for my master thesis.
So far I had not a lot to do with statistics and also the Program R was new for me, so I spent a lot of time to get into the topic.
As basic skript I used https://openmx.ssri.psu.edu/docs/OpenMx/2.6.7/GeneticEpi_Matrix.html
I have some questions for you and hope someone can help me with them.
-
I tested the skript with the standard twinData set and tried to achieve the same output as they did in the paper: Estimating Heritability from Twin Studies (from Karin Verweij). I was able to reconstruct the same results expect the c^2 value was different. Could someone tell me how to get the correct c^2 value (compare skript)?
-
In the attachment you will find a file with some variables from my twin study. How do I determine the starting values? Can I use the Falconer's formula to set the start values?
PS: Skript and dataset are attached! Duration 1 stands for sleep duration in twin 1 and duration 2 for the sleep duration twin 2.
I would be greatful for your help.
Best
Andrea