我是R.的新手。我想添加一列,将解释为当前数据框的拦截列。任何有用的评论都会很有意义。
谢谢! 伊萨克
答案 0 :(得分:2)
由于回收,您只需要d$one <- 1
在数据框one
中创建一个名为d
的列。但@AntoineVernet完全正确,即使你有一个复杂的混合模型,这几乎肯定是不必要的。
尝试这些模型(这比你的模型更简单)并检查输出以了解我的意思。
set.seed(5)
d <- data.frame(x = rep(1:3, 5),
y = round(rnorm(15, mean=10), 1),
g = factor(rep(1:5, each=3)) )
d$one <- 1
library(lme4)
lmer(y ~ x + (1|g), data=d)
lmer(y ~ (1|g), data=d)
lmer(y ~ 1 + (1|g), data=d)
lmer(y ~ 0 + one + (1|g), data=d)
答案 1 :(得分:0)
这应该有效:
incpt&lt; - rep(1,nrows(df))
df&lt; - as.data.frame(cbind(incpt,df))