Hello,
This is my first time using this forum...any help is greatly appreciated!
Are there any OpenMx scripts for determining power in a twin sample using bivariate ACE analyses?
In a twin sample, I'm conducting bivariate ACE analysis of a personality trait and depressive symptoms. I found 4.6% in shared additive genetic variance, which is smaller than expected (compared to previous reports in the literature)
My univariate ACE results: A^2 for personality trait = 62%; A^2 for depressive symptoms = 32%
My sample size is on the small side (MZ=112, DZ=100; total N=212), so I'm trying to determine if the relatively small effect size is due to power issues. I'd like to assess power of my sample.
Are there any OpenMx R scripts for assessing power in a twin sample using bivariate ACE ?
OpenMx version: 2.0.1.4157
R version: R version 3.1.3 (2015-03-09)
Platform: x86_64-apple-darwin10.8.0
Default optimiser: NPSOL