[R-sig-ME] Is there a leverage diagnostic test for glmmTMB models?
bbo|ker @end|ng |rom gm@||@com
Mon Jul 27 22:25:09 CEST 2020
As I mentioned in my CrossValidated answer
, this does work with the add-on influence_mixed function in the glmmTMB
package ; you just need to specify group=".cases" (or leave the group
blank, since ".cases" is the default value).
This doesn't use hatvalues, which as I commented there aren't
available (yet) for glmmTMB; instead, it does brute-force leave-one-out
computations (although it can do them in parallel).
It would be a good idea making this functionality more easily
available, but I've messed up in the past when trying to export 'car'
functionality, so I'd like to be careful.
On 7/27/20 3:47 PM, John Fox wrote:
> Dear April,
> Unfortunately, there's no hatvalues() method for "glmmTMB" objects,
> nor do functions in the car package, where the influenceIndexPlot()
> function that you tried to use originates, have "glmmTMB" methods. car
> does have methods for the "merMod" objects produced by functions in
> the lme4 package and for the "lme" objects produced by the lme()
> function in the nlme package, but those probably aren't useful to you.
> I'm willing to look into providing the necessary infrastructure in the
> car package to support glmmTMB() but can't promise to do that before
> the Fall.
> John Fox
> Professor Emeritus
> McMaster University
> Hamilton, Ontario, Canada
> web: https://socialsciences.mcmaster.ca/jfox/
> On 7/27/2020 1:40 PM, April Martinig wrote:
>> I've posted my question to Cross Validated, but I am also linking to
>> it here
>> because it might have been missed.
>> I am trying to use the DHARMa package to test for how much leverage
>> specific points (on one extreme) for my response variable has on my
>> In my question I have posted the data, what I have tried, and linked
>> to other sources. Thank you in advance.
>> Take care,
>> April Martinig, M.Sc.
>> Ph.D. Candidate in Ecology
>> University of Alberta
>> “I do not accept the conventional order of things as a given."
>> [[alternative HTML version deleted]]
>> R-sig-mixed-models using r-project.org mailing list
> R-sig-mixed-models using r-project.org mailing list
More information about the R-sig-mixed-models