将每晚的平均值与每晚的值进行比较,统计有效性?

我不会取扰动前夜晚的平均值。我会将原始数据汇集在“扰动前”因子下,然后将其与“夜晚1”、“夜晚2”等汇集的数据进行比较。如果你有多个站点,你需要加入一个随机效果。此外,假设您的响应是在事件发生前几分钟,则需要使用gamma分布。以下是如何在R和lme4中实现这一点的代码:

library(lme4)

my.data$pool <- relevel(my.data$pool, ref="pre-disturbance")
#This is setting your model's reference level to the pre-disturbance pooled data. 
#When you run summary() on the model object, it will compare your daily disturbance 
#pools directly to the 'pre-disturbance' pool. 

model <- glmer(num.minutes ~ data.pool + (1|site), family=Gamma(link='log'), data=my.data)
#The model

summary(model)