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