Skip to content

incorporate pedigree into lmer

2 messages · Doran, Harold, Jon Hallander

#
Jon

I think you will find the members of this list helpful, but you have given zero information that we can use to help you with your problem. Please provide a representation of your model and a description of your data. Or better yet, some sample data. I don't have any idea what an "animal model" is (maybe others do), but if you can provide a statistical representation, maybe we can walk you through this.

Harold
#
Oh, sorry for being unclear. I am using a mixed linear model including both random and fixed effects:
y = Xb + Za + e
where the vector y are observed individual data, b is a vector of fixed effects, a is an individual additive random effect that are normally distributed vectors with (co)variance A*var(a). X and Z are incidence matrices relating the fixed and random effects (b and a) with the observations y. Furthermore, e is the error term, w normally distributed with mean zero and variance var(e). A is the additive genetic relationship matrix which is computed using the pedigree information (given below). A typical data might look like this:
Id Sire Dam y (i.e. weight in kg)
1  Na   Na  4.3 	
2  Na   Na  4.7
3  1    2   4.4
4  1    Na  4.1 
5  4    2   4.3
6  4    2   4.2

Thanks for the help. 

Best regards,
Jon

-----Original Message-----
From: r-sig-mixed-models-bounces at r-project.org [mailto:r-sig-mixed-models-bounces at r-project.org] On Behalf Of Doran, Harold
Sent: den 6 februari 2007 13:51
To: Jon Hallander; r-sig-mixed-models at r-project.org
Subject: Re: [R-sig-ME] incorporate pedigree into lmer

Jon

I think you will find the members of this list helpful, but you have given zero information that we can use to help you with your problem. Please provide a representation of your model and a description of your data. Or better yet, some sample data. I don't have any idea what an "animal model" is (maybe others do), but if you can provide a statistical representation, maybe we can walk you through this.

Harold
_______________________________________________
R-sig-mixed-models at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
###########################################

This message has been scanned by F-Secure Anti-Virus for Mic...{{dropped}}