| Author |
Message |
< ASReml ~ restrict on fixed effects for mixed model in R-Asreml
|
| nicobardol |
Posted: Mon Jun 25, 2012 9:49 am |
|
|
|
Joined: 13 Jan 2010
Posts: 7
Location: INRA Gif sur Yvette
|
Dear all,
I am doing a mixed model, and I would like put some restrict for the estimation of my fixed effects.
For example, I have a fixed population effect. There is 7 populations, so I give to asreml a matrix with 7 columns and use the function grp() to consider one factor.
I would like to add a restriction for their estimation effect as :
pop1 + pop2 + pop3 + pop4 + pop5 + pop6 + pop7 = 0
How is possible to do?
My model is the following one :
model2.asr = asreml(fixed = Pheno ~ 1+grp(pop), group = QTL, random=~ ped(Ind, var=T), ginverse=list(Ind=ColInvSimil), data=dataASREML,sed=T)
Thanks for your answers,
Regards,
Nicolas. |
|
|
| Back to top |
|
| luansheng |
Posted: Tue Jun 26, 2012 2:29 am |
|
|
|
Joined: 05 Aug 2008
Posts: 32
|
nicobardol wrote: Dear all,
I am doing a mixed model, and I would like put some restrict for the estimation of my fixed effects.
For example, I have a fixed population effect. There is 7 populations, so I give to asreml a matrix with 7 columns and use the function grp() to consider one factor.
I would like to add a restriction for their estimation effect as :
pop1 + pop2 + pop3 + pop4 + pop5 + pop6 + pop7 = 0
How is possible to do?
My model is the following one :
model2.asr = asreml(fixed = Pheno ~ 1+grp(pop), group = QTL, random=~ ped(Ind, var=T), ginverse=list(Ind=ColInvSimil), data=dataASREML,sed=T)
Thanks for your answers,
Regards,
Nicolas.
You can try setting population as random effect, then the sum of population is to zero. |
|
|
| Back to top |
|
|
|
All times are GMT
|
|
You cannot post new topics in this forum You cannot reply to topics in this forum You cannot edit your posts in this forum You cannot delete your posts in this forum You cannot vote in polls in this forum You can attach files in this forum You can download files in this forum
|
|
|