经过大量阅读后,我决定写信,因为我找不到解决问题的方法。
我已经对具有正态分布的连续变量进行了先验对比。现在,我有另一个变量,该变量是二项式的,我可以为其进行GLM。但是,当我进行先验对比时,在所有数据均为0的情况下没有结果(不是没有数据,它们在类别中全为0,因此我想将其与其他具有值的值进行比较) )。 数据结构如下:
头(穴) 日期天对待psu sp钻洞钻洞 1 3 0 30-30 36 B 1 0 2 3 0 30-30 36 B 1 0 3 3 0 15-30 36 B 0 1 4 3 0 15-30 36 B 0 1 5 3 0 15-30 36 B 0 1 6 3 0 10-25 36 B 0 1
我的模型是这样的: model4B2 <-glm(挖洞〜请款待,family = binomial(link =“ logit”),data = D4B)
我做了这样的对比:
require(multcomp)
k3010R1<-matrix(c(3,-1,-1,-1),1)
k3010R1
t3010<-glht(model4B3.2,linfct=k3010R1)
summary(t3010)
但是不起作用,我确定它应该起作用。
是因为我的解释变量是算术的吗? 还是在某些类别中全为0时无法对二项式进行对比?
先谢谢您