R - 创建虚拟对象来评估一周中的几天效果

如何解决R - 创建虚拟对象来评估一周中的几天效果

我试图创建工作日假人，开始日期是星期日，总观察期数是 271。我使用了这些代码：

mon <- rep(c(0,1,0),T/7) 
tue <- rep(c(0,T/7)
wed <- rep(c(0,T/7)
thu <- rep(c(0,T/7)
fri <- rep(c(0,T/7)
sat <- rep(c(0,1),T/7)
sun <- rep(c(1,T/7) 

linreg4 <- lm(formula = hits ~ 0 + mon + tue + wed + thu + fri + sat + + sun) 

summary(linreg4)

并得到错误提示：

“model.frame.default 中的错误（公式 = hits ~ 0 + mon + tue + wed + thu + ： 可变长度不同（发现'mon'）”

请帮我修复这个错误，当我重新排列顺序时，错误仍然存​​在。

谢谢！

解决方法

你得到这个错误是因为 271 /7 不是一个整数，所以当你使用 rep 时它实际上是向下取整的：

total = 271
hits = rnorm(271,10,1)

mon <- rep(c(0,1,0),total/7) 

length(mon)
[1] 266

lm(hits ~ mon)

Error in model.frame.default(formula = hits ~ mon,drop.unused.levels = TRUE) : 
  variable lengths differ (found for 'mon')

您基本上可以重复从星期日开始的天数向量，并且只取前 271 天：

days = rep(c("sun","mon","tue","wed","thu","fri","sat"),ceiling(total/7))
days = days[1:total]

linreg4 <- lm(formula = hits ~ 0 + days)