[R-br] Modelo Animal - Pedigreemm

[Marcelo Cardoso mello]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110520/33c3fca4/attachment.html>