我正在使用网站的数据来检验不同因素在几个晚上对网站退出的影响.我用日落后的几分钟来研究蝙蝠离开现场的情况.

我想看看骚乱前的5个晚上和骚乱后的5个晚上.我删除了我的分析中的干扰之夜(Night0).

我的问题是:我是否可以取前5天(自然变率)的变量(日落后的分钟数)的平均值,然后将其与夜间+1、+2、+3、+4和+5的输出进行比较?Is this statistically valid?

我不愿用日落后的分钟数的平均值来影响"之前"的因素,或者叠加前5个晚上的所有行数并影响它到"之前"的因素.

我希望我的问题是清楚的.

非常感谢您的回复

推荐答案

我不会接受骚乱前的夜晚,不会的.我会将原始数据集中在一个"干扰前"因素下,然后将它们与"夜间1"、"夜间2"等汇集的数据进行比较.如果你有多个站点,你需要加入一个随机效果.此外,考虑到事件发生前的响应时间为几分钟,您需要使用伽马分布.下面是在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)

R相关问答推荐

高质量地将R格式的图表从Word中输出

使用ggcorrplot在相关性矩阵上标注supertitle和index标签

无法在我的情节中表现出显着的差异

使用tidyverse / Mutate的存款账户余额

更改默认系列1以更改名称

derrr mutate case_when grepl不能在R中正确返回值

如果第一个列表中的元素等于第二个列表的元素,则替换为第三个列表的元素

R—将各种CSV数字列转换为日期

矩阵的堆叠条形图,条形图上有数字作为标签

如何提取所有完美匹配的10个核苷酸在一个成对的匹配与生物字符串在R?>

在不丢失空值的情况下取消列出嵌套列表

根据1个变量绘制 colored颜色 发散的 map ,由另一个变量绘制饱和度,ggplot2不工作

按多列统计频次

调换行/列并将第一行(原始数据帧的第一列)提升为标题的Tidyr类似功能?

减go R中列表的所有唯一元素对

如何计算增加10米(0.01公里)的行?

在ggplot2上从多个数据框创建复杂的自定义图形

如何在使用Alpha时让geom_curve在箭头中显示恒定透明度

如何将EC50值绘制在R中的剂量-react 曲线上?

从多行中 Select 最小值