1

我是 R 的新手。我想添加一列，将其解释为当前数据帧的截距列。任何有用的评论将不胜感激。

谢谢！伊萨克

4

2 回答 2

2

由于回收，您只需要在数据框中d$one <- 1创建一列具有名称的列。但是@AntoineVernet 是完全正确的，即使您有一个复杂的混合模型，这几乎肯定是不必要的。oned

试试这些模型（当然比你的更简单）并检查输出以了解我的意思。

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)
于 2013-02-12T19:21:07.703 回答
0

这应该有效：

incpt <- 代表（1，nrows（df））

df <- as.data.frame(cbind(incpt,df))

于 2013-02-12T18:24:17.250 回答