Displaying 1 result from an estimated 1 matches for "missid".
Did you mean:
missed
2012 Nov 09
0
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,
-------------...