我是R的新手。我想添加一列1,将其解释为当前数据帧的intercept列。任何有用的意见都将不胜感激。
谢谢!Issac
发布于 2013-02-13 03:21:08
由于要进行回收,您只需要使用d$one <- 1在数据框d中创建一列名为one的数据。但@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)发布于 2013-02-13 02:24:17
这应该是可行的:
incpt <- rep(1,nrow(Df))
df <- as.data.frame(cbind(incpt,df))
https://stackoverflow.com/questions/14838390
复制相似问题