[R-sig-ME] using glmer with user-defined link function

[Ben Bolker]

  I tried to deal with this at one point to put in a negative
binomial model.  I, too, was hopeful that just
adding the name of my link function to the list
of allowed names would be all that was
needed. Unfortunately, the variance functions are
hard coded within the glmer code, so it will not be completely
trivial to do this ...

   Ben


vito muggeo wrote:
> dear all,
> I am trying to use an user-specified link function in glmer() (via the
> argument family), but the error is
> 
> Error in famType(glmFit$family) : unknown link: ‘my_link’
> 
> my_link is a proper link function (of class "link-glm") which works
> correctly with glm() and also with MASS::glmmPQL()
> 
> It appears that the function lme4::famType() checks for the name by
> comparing it with the names in "famNms" which includes "gaussian",
> "binomial", "Gamma", ... however I do not know whether it is sufficient
> to include "my_link" in famNms, and possibly I do not how to do it.. :-)
> 
> many thanks,
> vito
> 
> 
>