Question

我是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中做到这一点？

谢谢

来自glmmadmb

0 个答案: