【发布时间】:2018-04-15 18:35:11
【问题描述】:
我在运行 lmer 时收到此警告消息。我知道为什么会这样。有什么办法可以隐藏吗?
fixed-effect model matrix is rank deficient so dropping 1 column / coefficient
【问题讨论】:
我在运行 lmer 时收到此警告消息。我知道为什么会这样。有什么办法可以隐藏吗?
fixed-effect model matrix is rank deficient so dropping 1 column / coefficient
【问题讨论】:
这是一个消息(a milder form of informational message than a warning:警告使用前缀Warning message:),所以你可以使用suppressMessages():
library(lme4)
ss <- transform(sleepstudy,Days2=Days) ## create duplicate variable
m1 <- lmer(Reaction~Days+Days2+(1|Subject),ss)
## fixed-effect model matrix is rank deficient so dropping 1 column / coefficient
m2 <- suppressMessages(lmer(Reaction~Days+Days2+(1|Subject),ss))
但是,您可以通过包含control=lmerControl(check.rankX="silent.drop.cols") 来告诉lmer 禁止仅显示此消息:
m3 <- update(m1, control=lmerControl(check.rankX="silent.drop.cols"))
一般来说,如果可能的话,最好首先在模型中避免这些共线项(尽管它们是无害的)。 (我看到 OP 说他们知道发生了什么,所以这些信息更适合未来的读者。)您可以通过查看模型矩阵的属性来查看正在删除哪些变量:
attr(getME(m2,"X"),"col.dropped")
## Days2
#3 3
【讨论】: