search

pitakakariki / simr

Power Analysis of Generalised Linear Mixed Models by Simulation
68 stars 19 forks source link

Mislabelled default LR test #224

Open pitakakariki opened 2 years ago

pitakakariki commented 2 years ago 
> simdata$g <- as.character(simdata$g)
> fm <- lmer(y ~ g + (1|x), data=simdata)
> doTest(fm, fixed("g"))
boundary (singular) fit: see ?isSingular
p-value for predictor 'g': 1.434264e-13
          --------------------
Test: Kenward Roger (package pbkrtest)
> doTest(fm, fixed("g", "lr"))
boundary (singular) fit: see ?isSingular
p-value for predictor 'g': 1.434264e-13
          --------------------
Test: Likelihood ratio
> doTest(fm, fixed("g", "kr"))
boundary (singular) fit: see ?isSingular
p-value for predictor 'g': 1.035445e-10
          --------------------
Test: Kenward Roger (package pbkrtest)
pitakakariki commented 2 years ago

Error seems to be in defaultdesc, which checks for factor but not character.

pitakakariki commented 2 years ago

Might apply to lm too - t-test vs LR.