更新＆＃34; merModLmerTest＆＃34; object返回class＆＃34; lmerMod＆＃34;的对象

Question

只是想知道这是否是众所周知/可取的行为......？

可重复的例子：

## import lmerTest package
library(lmerTest)

## an object of class merModLmerTest
m <- lmer(Informed.liking ~ Gender+Information+Product +(1|Consumer), data=ham)

class(m)
class(update(m,REML=FALSE))

Answer 1

至少对我来说并不明显。您可以使用as(.,"merModLmerTest") ...

解决此问题

library(lmerTest)
m <- lmer(Reaction~Days+(Days|Subject),sleepstudy)
m2 <- update(m,REML=FALSE)
m3 <- as(m2,"merModLmerTest")
class(m3)  ## "merModLmerTest"

或者您可以使用此辅助函数：

update2 <- function(x,...) {
   z <- update(x,...)
   return(as(z,"merModLmerTest"))
}
update2(m,REML=FALSE)