Attachment | Size |
---|---|
bivariate ACE binary.R [6] | 20.76 KB |
data.csv [7] | 217.86 KB |
Hi everyone, I am fitting a bivariate ACE model and I am very confused about the Ra (genetic correlation of two traits) that is larger than 1. Since Ra is a correlation coefficient, it should range from -1 to 1, am I right? I don't know how to explain this result.
In the ACE model, Ra=1.37 ,Rc=-0.04, Re=0.22.
Then I fitted ADE model, Ra=0.07, Rd=NA, Re=0.22.
The results of the goodness-of-fit test between saturated model and ACE or ADE are as follows.
base comparison ep minus2LL df AIC diffLL diffdf p
1 Sat 24 7734.750 7728 -7721.250 NA NA NA
2 Sat ACE 15 7749.792 7737 -7724.208 15.04198 9 0.08978923
3 Sat ADE 15 7749.741 7737 -7724.259 14.88488 9 0.09414712
Thus, I wonder why Ra can be larger than 1?
And in this case, should I report Ra of ACE model (if so, how to explain it), or report Ra of ADE model?
I have enclosed the a sample data and my code for your review.
Thank you for your time.