我是R的新手,我想知道我的研究发生在不同领域的随机效应是否会影响我的结果,并对我的条件产生误报 . 我的数据基于在9个区域的四种不同条件下进行的自然观测 . 它是一个主题设计之间 - 每一行是对不同主题的观察 .

条件=因素(4级)

面积=随机效应(9级)

GeneralDisplayingBehaviours = DV

ObsTime =观察发生了多少小步

这就是我的模型

data2$Condition<-as.factor(data2$Condition)
data2$Area<-as.factor(data2$Area)
data2$Condition<-relevel(data2$Condition,ref="22")
data2$Area<-relevel(data2$Area,ref="1")

mod<-glmmadmb(GeneralDisplayingBehaviours~Condition+ObsTime+(1|Area), data=data2, family="nbinom", zeroInflation=TRUE)

这是输出:

Call:
glmmadmb(formula = GeneralDisplayingBehaviours ~ Condition +
ObsTime + (1 | Area), data = data2, family = "nbinom", zeroInflation = TRUE)

AIC: 2990.1 

Coefficients:
            Estimate Std. Error z value Pr(>|z|)  
(Intercept)   1.6233     0.5019    3.23   0.0012 ** 
Condition12   1.3291     0.1330    9.99   <2e-16 ***
Condition11   1.2965     0.1294   10.02   <2e-16 ***
Condition21   0.0715     0.1351    0.53   0.5966    
ObsTime       0.0829     0.0341    2.43   0.0151 *  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Number of observations: total=360, Area=9 
Random effect variance(s):
Group=Area

             Variance    StdDev
(Intercept) 8.901e-09 9.434e-05

Negative binomial dispersion parameter: 1.7376 (std. err.: 0.16112)
Zero-inflation: 0.16855  (std. err.:  0.02051 )

Log-likelihood: -1487.06

然后我会继续改变每个级别的条件参考,我想我必须对区域做同样的事情 .

我如何解释这些结果并报告它们?这告诉我有关区域随机效应的信息是什么?它会影响DV吗? DV的条件是否明显不同,哪些条件明显不同?

对于最后一个问题,我有一种感觉,我需要进行多次比较,那么我将如何在GLMMABMD中执行此操作?

谢谢