Não estou conseguindo reproduzir o pequeno exemplo abaixo, usando o modelo animal. Aparece o erro abaixo. Alguém poderia me dar uma luz?? É possível predizer valores ao nível de indivíduos (Animal) neste exemplo, usando o R (pedigreemm)?? O que esta errado??

Herd<-as.factor(c(1,2,1,2,1))

Animal<-as.factor(c(1,2,3,4,5,6))

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))

Erro em function (fr, FL, start, REML, verbose)  :

  Number of levels of a grouping factor for the random effects

must be less than the number of observations

Herd	Animal	Sire	Dam	Phenotyp
1	1	0	0	100
2	2	0	0	130
1	3	0	0	120
2	4	1	2	110
1	5	4	3	140