Skip to content

confint.merMod, bootstrap and weights

2 messages · Denis Haine, Ben Bolker

Original
Denis Haine
# 
That's what I thought, that the weighting was not taken into account.
The weights are used to estimate an underlying causal model. The Poisson
model is used as a Cox regression model with 2-level random effects.
I believe I have to write my own boostrap function, or be happy with the
Wald method for confint.

Denis
:

  
  


  


    

    

      
      

        


  

        

      Ben Bolker
    

    

      
      #
    

  


  

      
In general profile confidence intervals (the default method for 
confidence intervals) are a good intermediate between Wald intervals and 
bootstrap intervals, in both computational cost and accuracy.  If you 
can decide how you want weights to be used, you might be able to use 
assignInNamespace() to hack lme4:::simfunList$poisson, which is 
currently defined as

function (object, nsim, ftd = fitted(object), wts = weights(object))
{
     wts <- weights(object)
     if (any(wts != 1))
         warning("ignoring prior weights")
     rpois(nsim * length(ftd), ftd)
}

        
On 16-03-14 10:17 AM, Denis Haine wrote: