R-asreml : effets maternels
Ajouté par Herve Chapuis il y a environ 7 ans
Bonjour
Quelqu'un aurait-il la moindre idée de comment ajuster un modèle avec effet maternel génétique (corrélé à l'effet génétique direct) avec asreml-R ?
Je trouve R bien pratique pour simuler des données et, en théorie, asreml est bien pratique pour analyser ces données.
Or je ne trouve pas la syntaxe pour effectuer ce calcul (que j'imaginais "basique"). Les maigres résultats trouvés sur le web me conduisent à écrire :
reml.ainv<-asreml.Ainverse(reml.gen)$ginv #calcul matrice inv(A)
effmat<-asreml(perfo~sexe+generation,random=str(~ped(animal)+ped(mere),us(2):ped(animal)),ginverse=list(animal=reml.ainv,mere=reml.ainv),data=reml.perfo,maxiter=25)
mais le résultat est :
Class 'formula'
Error in eval(expr, envir, enclos) :
impossible de trouver la fonction "us"
Et c'est pareil si j'utilise id(2)
Merci !
Hervé