Skip to content
Back to formatted view

Raw Message

Message-ID: <af2de572078cd2c65c766f4f927e47c9.squirrel@webmail.uel.br>
Date: 2012-11-09T17:27:24Z
From: Silvano
Subject: 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