Hi all,

Is it possible to use mcmcsamp() in lme4 after fitting a model with a gamma
response variable and a log link function?

eg.. model fitted like this..

best_fm_bio <- lmer(Bio ~ P_pca +(1|Site), family = Gamma(link = "log"))

# then

mcmcsamp(best_fm_bio, n = 1000)

# I am getting this error ...

# Error in .local(object, n, verbose, ...) : Update not yet written

Is there another package I could use to run the mcmc for this?

Alternatively, I'm thinking that I could transform the response variable
and use a linear mixed model but I would prefer to avoid that.

Many thanks,

Ruth

p.s. I can supply a dataset if necessary.

	[[alternative HTML version deleted]]


