[R-br] Modelo Animal - Pedigreemm
Fabio Mathias Corrêa
fabio.ufla em yahoo.com.br
Sexta Maio 20 13:08:37 BRT 2011
Eu só consegui reproduzir o exemplo depois de ler o artigo. Recomendo o mesmo.
São dois artigos citados na library. Veja na documentação!
Valeu!!!
Fábio Mathias Corrêa
Universidade Estadual de Santa Cruz
Departamento de Ciências Exatas e da Terra - DCET
Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna
CEP 45662-900. Ilhéus-Bahia
Tel.: 73-3680-5076
________________________________
De: Marcelo Cardoso mello <marcelo_cm32 em yahoo.com.br>
Para: r-br em listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman em yahoo.com.br>
Enviadas: Sexta-feira, 20 de Maio de 2011 13:04
Assunto: Re: [R-br] Modelo Animal - Pedigreemm
Certo Ivan. Corrigido o problema de digitação. Mas o erro persiste.
???????
Rode os dados abaixo:
Herd<-as.factor(c(1,2,1,2,1))
Animal<-as.factor(c(1,2,3,4,5))
Sire<-c(0,0,0,1,4)
Dam<-c(0,0,0,2,3)
Phenotyp<-c(100,130,120,110,140)
require(pedigreemm)
p<- pedigree(sire = Sire,
dam = Dam, label= 1:5)
model<-pedigreemm(Phenotyp~
(1|Animal)+ Herd, pedigree=list(Animal=p))
--- Em sex, 20/5/11, Ivan Bezerra Allaman <ivanalaman em yahoo.com.br> escreveu:
>De: Ivan Bezerra Allaman <ivanalaman em yahoo.com.br>
>Assunto: Re: [R-br] Modelo Animal - Pedigreemm
>Para: "R Brasil" <r-br em listas.c3sl.ufpr.br>
>Data: Sexta-feira, 20 de Maio de 2011, 11:59
>
>
>Olhe novamente para o seu exemplo. O fator 'Animal' está com seis observações!! Deveria estar como mesmo número de observações das demais fontes de variação e variáveis, ou seja, cinco observações.
>
>
>Allaman
>(S,f,P)
>
>M.Sc Ivan Bezerra Allaman
>Zootecnista
>Doutorando em Produção Animal/Aquicultura - UFLA
>email e msn - ivanalaman em yahoo.com.br
>Tel: (35)3826-6608/9900-2924
>-----Anexo incorporado-----
>
>
>_______________________________________________
>R-br mailing list
>R-br em listas.c3sl.ufpr.br
>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110520/4ef7a09d/attachment.html>
Mais detalhes sobre a lista de discussão R-br