R/johnson_neyman_lme4.R

In AlexHartmann00/lmmadd: Functions for Linear Mixed Models

#'@noRd
.johnson_neyman_lme4 <- function(object,X,M,modrange=NULL,resolution=10000,sig.thresh){
  B <- lme4::fixef(object)
  sigm <- vcov(object)
  mr <- range(model.matrix(object)[,colnames(model.matrix(object))==M])
  if(!is.null(modrange))mr <- modrange
  df <- nrow(model.matrix(object)) - length(B)
  return(.johnson_neyman(B,sigm,X,M,mr,df,resolution,sig.thresh))
}