[R-sig-ME] HPDinterval MCMCglmm and using transformed variables in MCMCglmm

Jarrod Hadfield j.hadfield at ed.ac.uk
Mon Jul 12 12:30:33 CEST 2010


Dear Rebecca,

You probably need to load MCMCglmm (library(MCMCglmm)) or at least  
load coda. At the moment the response needs to be logged externally:

my.data$ylog<-log(my.data$y)

but the next version will "fix" this.

Cheers,

Jarrod


On 12 Jul 2010, at 10:31, Rebecca Ross wrote:

> Dear ME users,
>
> Apologies if this is a naive question, I am new to using MCMC glmm.  
> I have been working through the tutorial, using some of my own data,  
> and I thought I had followed the steps correctly. The model I have  
> been fitting is:
>
> AMstem3<-MCMCglmm(log(stem)~Dist_fr_A,random=~Block 
> +animal,pedigree=Ped,data=SapRO,prior=prior1.3,
> family="gaussian")
> Briefly, I want to see how the Distance variable (continuous)  
> affects stem weight, which I have log transformed as the raw data is  
> heteroscedastic. Block is field block, and I am fitting the animal  
> random effect to take into account different relatedness between  
> individuals (I'm not at the moment trying to estimate Va).
>
> When I tried this model last night, I could get the 95% CIs using  
> this function
> HPDinterval(AMstem3$Sol[,"Dist_fr_A"])
> HPDinterval(AMstem3$VCV)
>
> When I tried this morning, I get the following error message
> Error in function (classes, fdef, mtable)  :
>  unable to find an inherited method for function "HPDinterval", for  
> signature "mcmc"
> Can anyone explain what  has gone wrong? I've tried changing the  
> terms in the model, but that doesn't help, so I think I must have  
> changed something in R elsewhere? Or else missing something  
> completely obvious.
>
> Question 2: I dont think MCMCglmm is actually logging my stem  
> variable because the fixed effect parameter estimate from the model  
> above (AMstem3) are about a factor of 10 out from the fixed effect  
> estimate from
> a similar model, fit in lme4, without the animal effect and very  
> similar to the same model using the raw stem weight (see below). Do  
> I need to specify the log transform differently?
>
> AMstem3c Dist_fr_A  posterior mode: -0.004236659
> versus:
> stem5a<-lme(log(stem)~Dist_fr_A,random=~1| 
> Block,data=SapRO,na.action=na.omit)
> Fixed effect Estimate: Dist_fr_A   -0.0004442 (SE=0.00006136)
> and:
> stem5b<-lme(stem~Dist_fr_A,random=~1| 
> Block,data=SapRO,na.action=na.omit)
> Fixed effect Estimate: Dist_fr_A   -0.004342 (SE=0.0006714)
>
> Thank you very much for your help,
>
> Rebecca Ross
> DPhil Candidate,
> University of Oxford.
>
>
>
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>


-- 
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.




More information about the R-sig-mixed-models mailing list