[R] Differences between MTDFReml and kinship2

[silvano]

Hi,

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?

Thanks a lot,

Silvano.
	[[alternative HTML version deleted]]