如何解决调整差异分析中的每月季节性影响,并在r中随机影响
我是第一次使用随机截距分析时序数据。我有前后位置指示器的不同地区的时间序列数据。也有一些位置属于一个干预组,而另一些则属于不同的干预组。这是一个玩具数据集,
data<-data.frame(Outcome=c(arima.sim(model=list(ar=c(.9,-.2)),n=36)+rnorm(36),arima.sim(model=list(ar=c(.5,-.1)),arima.sim(model=list(ar=c(.6,-.3)),n=36)+rnorm(36)),Time=rep(1:36,3),Treatment=c(rep("A",72),rep("B",36)),Location= c(rep("South",36),rep("North",rep("East",pre_post = factor(c(rep(0,20),rep(1,16),rep(0,16))))
这里Outcome
是3年每月报告,是从3个不同的Location
收集的,Treatment
是干预指标,pre_post
是事前或事后时间指示符。为了估算位置随机拦截的差异,我使用以下代码,
model <- lmer(Outcome~ pre_post+Treatment+ Treatment:pre_post+ (1|Location),data = data)
我的原始数据具有每月季节性。现在我的问题是我应该如何调整这里的每月季节性?任何建议都会很有帮助。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。