Internal function to perform martini test
epi(y, Z, G = NULL)
y
Z
genomic information matrix
G
kinship matrix
Estimated breeding values