0

我在 GAMLSS 中运行可能性测试时遇到了问题。

我正在尝试测量观察者对蜥蜴遭遇概率的影响。几位观察者在不同的日子调查了一系列地块并测量了蜥蜴遭遇的概率,因此我的数据集由一系列概率值组成(值从 0 到 1,强烈膨胀为 0)。因此,我采用 Beta 膨胀分布 (beinf) 进行建模。

运行模型选择后,我最终得到以下最佳模型(最低 AIC):model<-gamlss(observations~Date*Year+Observer+random(Plot),data=dat.tf,family=BEINF,trace= F)

为了评估观察者的重要性,我比较了这个模型和空模型 null.model<-gamlss(observations~Date*Year+random(Plot),data=dat.tf,family=BEINF,trace=F)

LR.test(null.model,model) 并返回错误消息:“空模型的偏差比替代模型小。模型应该嵌套”

我在哪里做错了?

4

0 回答 0