Предсказать и использовать график взаимодействия без единого эффекта
У меня есть проблема в построении взаимодействия и в использовании функции предсказания. Это моя последняя модель:
`mymodel <- glm(y ~ SITE + TREAT*HET, family= quasibinomial, data = PRESU)`.
Я хочу построить взаимодействие между Treat (категориальный, 2 уровня) и Het (Непрерывное). Следуя главе 16 Кроули (более или менее, но в примере нет другой переменной), я попробовал этот код:
x <-seq(0,1,length=100)
v <- rep("treat1", length(x))
y <- predict(mymodel, list(TREAT=factor(v),HET=x), type="response")
plot(HET[TREAT=="treat1"], fitted(mymodel)[TREAT=="treat1"],col="black")
lines(x,y,col="Black",lwd=1.5,lty=3)
t <- rep("treat2", length(x))
s <- predict(mymodel, list(TREAT=factor(t),HET=x), type="response")
points(HET[TREAT=="treat2"],fitted(mymodel)[TREAT=="treat2"],col="red")
lines(x,s,col="red",lwd=1.5)
Я получил этот вывод по прогнозу
Ошибка в eval(expr, envir, enclos): объект 'SITE' не найден.
Я не могу понять, как я могу построить свое взаимодействие с моделью. Если я удалю SITE код работает, но я хочу сохранить его в окончательной модели. Я надеюсь, что я объяснил ясно. Спасибо за то, кто может помочь.
Набор данных
SITE TREAT1 HET Tot Prey
1 treat1 0.3 30 9
1 treat2 0.5 30 16
2 treat1 0.8 30 19
2 treat2 0.76 30 25
3 treat2 0.97 30 1
3 treat1 0.75 30 23
4 treat1 0.67 30 18
4 treat2 0.89 30 29
У меня есть больше копий для каждого сайта, но допустим, что это что-то похожее на мой набор данных. Переменная ответа: у = (добыча, добыча)