Kinship2 and GenABEL
Hi, I'm using kinship2 to calculate heritabilty, but I would like calculate in GenABEL too. I trying the code:
require(kinship2) require(GenABEL)
pedig = with(Dados, pedigree(id=IID, dadid=PAT, momid=MAT, sex=SEX,
famid=FID, missid=0))
kmat = kinship(pedig)
(mod1 = polygenic(altura ~ SEX + idade, data=Dados, kin=kmat))
Erro em intI(i, n = d[1], dn[[1]], give.dn = FALSE) : invalid character indexing
How can I calculate heritabilty using GenABEL with kinship matrix of kinship2 package? Thanks, --------------------------------------------- Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ci?ncias Exatas Departamento de Estat?stica Fone: (43) 3371-4346