lme4 - error message when performing Gamma glmer in R- PIRLS step-halvings failed to reduce deviance in pwrssUpdate -
i trying perform glmer in r using gamma error family. error message: "error: (maxstephalfit) pirls step-halvings failed reduce deviance in pwrssupdate"
my response variable flower mass. fixed effects base mass, f1 treatment, , fertilisation method. random effects line , maternal id nested within line.
when perform same analysis using integer response (ie. flower number) error not occur.
here sample of data:
line maternal_id f1treat flwr_mass base_mass 17 81 stress s 2.7514 9.488 5 41 control o 0.3042 1.809 37 89 control o 2.3749 6.694 5 41 stress s 3.6140 9.729 9 5 control s 0.5020 7.929 13 7 stress s 0.4914 0.969 35 88 stress s 0.4418 1.840 1 57 control o 2.1531 6.673 13 7 stress s 3.0191 7.131
here code using:
library(lme4) m <- glmer(data=mydata, flwr_mass~base_mass*f1treat*so+(1 |line/maternal_id),family=gamma)
(i using r 3.0.3 windows)
@hongooi answered question in comments, repeat here else having issue. suggested changing
family=gamma
to
family=gamma(link=log)
Comments
Post a Comment