[R-es] pedigree list

[Alejandro Martinez Meier]

Estimados
Quisiera saber si alguno tiene experiencia en la utilización de la 
función pedigreemm para la estimación de los blups individuales en 
función de la información de pedigree
Dispongo del siguiente ejemplo

rep flia DAP
1 1 15.8
1 1 11.1
1 2 19.3
1 2 4.8
1 3 19.1
1 3 19.3
2 1 12.3
2 1 13.8
2 2 16.6
2 2 22.6
2 3 16.2
2 3 18.1
3 1 9.7
3 1 6.4
3 2 23.1
3 2 20.8
3 3 12.4
3 3 9.5

En lmer mi modelo es el siguiente
model<-lmer(data[,"variable"]~rep+(1|flia),data=data)
Esto me permite obtener los blups familiars pero no los individuals. 
Estuve explorando la función pedigreemm para estimar los blups 
individuales pero no tengo experiencia en el armado del pedigree list 
objects necesario para la estimación de los blups individuales.
Si alguno de ustedes tiene experiencia, desde ya muchas gracias.
Saludos, Alejandro




------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20100409/534ec4f1/attachment.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: almarti.vcf
Type: text/x-vcard
Size: 361 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20100409/534ec4f1/attachment.vcf>