二项式GLM的先验对比

时间:2019-01-17 14:13:25

标签: glm contrast

经过大量阅读后,我决定写信,因为我找不到解决问题的方法。

我已经对具有正态分布的连续变量进行了先验对比。现在,我有另一个变量,该变量是二项式的,我可以为其进行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)

测试对比30与所有对比(有4个类别要比较)

k3010R1<-matrix(c(3,-1,-1,-1),1)
k3010R1
t3010<-glht(model4B3.2,linfct=k3010R1)
summary(t3010)

但是不起作用,我确定它应该起作用。

是因为我的解释变量是算术的吗? 还是在某些类别中全为0时无法对二项式进行对比?

先谢谢您

0 个答案:

没有答案