после долгих чтений я решил написать, потому что не могу найти решение своего вопроса.
Я уже делал априорные контрасты для непрерывной переменной с нормальным распределением. Теперь у меня есть другая переменная, которая является биномиальной, и я могу сделать GLM для нее. Но когда я делаю априорные контрасты, это не дает результата в случаях, когда все данные равны 0 (не то, что нет данных, они просто все 0 в категории, и я хочу сравнить это с другими, которые имеют значения ).
Структура данных выглядит следующим образом:
голова (нора)
дата день угощение
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 (burrow ~ Treat, family = binomial (link = "logit"), data = D4B) </p>
И я сделал контраст следующим образом:
require(multcomp)
Тест сравнивает 30 против всех (есть 4 категории для сравнения)
k3010R1<-matrix(c(3,-1,-1,-1),1)
k3010R1
t3010<-glht(model4B3.2,linfct=k3010R1)
summary(t3010)
Но не работает, и я уверен, что это должно работать.
Может быть, потому что моя объясняющая переменная является категорической?
Или просто невозможно сделать контраст для биномиального, когда у вас есть все 0 в какой-либо категории?
Заранее спасибо,