error in gamma link specification for lmer

Asked Mar 28 '17 at 17:51

Active Mar 28 '17 at 17:55

Viewed 1,256 times

I'm trying to fit a mixed-effects model with a gamma distribution. The most basic model has one fixed predictor and 1 random effect. No matter which link I specify (I've tried log, identity and inverse), I obtain the following error. My real data has zeros in Y, but even when I apply simulated data with only positive Y as below, it throws the same error.

mockdf = data.frame(y =  rnorm(100,77,6.5), x1 = sample(letters,100,replace = T), x2 = seq(1900,1999,1))
mod = lmer(y ~ (1|x1) + x2, family = gamma(link = 'identity'), na.action = na.exclude, data = mockdf)

Error in gamma(link = "identity") : 
  supplied argument name 'link' does not match 'x'

I searched through SO and couldn't find another person who ran into this error. Is my syntax incorrect?

Thanks for your help.

edited Mar 28 '17 at 17:55

Roman Luštrik

69,533
24
154
197

asked Mar 28 '17 at 17:51

makai

2

Hum, lmer function from lme4 doesn't support a `family` argument... – Roman Luštrik Mar 28 '17 at 18:01
2

using `glmer`, as Roman mentioned... the family is `Gamma` (big G) – user20650 Mar 28 '17 at 18:18
Ah, I see. Thanks for the quick response. Is there any way to incorporate random effects in a gamma-distributed model? – makai Mar 28 '17 at 18:37
Sorry, I didn't sees @user20650 response until just now. Will attempt with edits. – makai Mar 28 '17 at 18:38

error in gamma link specification for lmer

0 Answers0