[R] Kinship2 and GenABEL

[Silvano Cesar da Costa]

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