From: "silvano" <silv...@uel.br>
I fitted a genetic model using kinship2 and I compared it with MTDFReml program
output.
The residual variance of both are very close but the genetic variance are very
differents.
The output are:
MTDFReml:
genetic variance = 1.24015
residual variance = 5.93424
R (Kinship2):
genetic variance = 0.767187
residual variance = 5.6712
Both of them use REML method. Could someone tell why the difference?
Hi Silvano.
You haven't given us enough information, I'm afraid. We would need the
commands you used, the versions of the software, plus output from both
jobs, and your data, or a smaller subset that exhibits the same problem.
You may have made an error, there may be a problem in one program, or the
likelihood might be flat around those solutions (ie your data are not very
informative).
It might be better to continue this on the R-sig-mixed-models list.
Cheers, David Duffy.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
