R Under development (unstable) (2023-11-03 r85470 ucrt) -- "Unsuffered Consequences" Copyright (C) 2023 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > ## drop1 may not work right with contrasts: make up an example something like this ... > ## options(contrasts=c("contr.sum","contr.poly")) > ## drop1(fecpoiss_lm3,test="Chisq",scope=.~.) > > > if (.Platform$OS.type != "windows") withAutoprint({ + library(lme4) + oldopts <- options(contrasts=c("contr.sum","contr.poly")) + fm1 <- lmer(Reaction~Days+(Days|Subject),data=sleepstudy) + drop1(fm1,test="Chisq") + ## debug(lme4:::drop1.merMod) + drop1(fm1,test="Chisq",scope=.~.) + + fm0 <- lm(Reaction~Days+Subject,data=sleepstudy) + drop1(fm0,test="Chisq",scope=.~.) + options(oldopts) ## restore original contrasts + + ff <- function() { + lmer(Reaction~Days+(Days|Subject),data=sleepstudy) + } + drop1(ff()) ## OK because sleepstudy is accessible! + + }) ## skip on windows (for speed) > > proc.time() user system elapsed 0.10 0.03 0.14