В настоящее время я анализирую занятость ящиков с летучими мышами и факторы, влияющие на занятость. Чтобы определить наиболее значимые факторы, я бегу GLM. Моя переменная ответа - двоичные данные (0 = занято / 1 = не занято), а мои объяснительные переменные являются числовыми, за исключением одной категориальной переменной с 4 уровнями. Бат-бокс установлен на дереве / столбе / балконе / фасаде).
Итак, япоместил все мои факторы в glm:
modelb <- glm(occupation ~ TreeCov500 + distance + Water2000 + greenvol + height + age + shadow + exposition + number_of_boxes + mounted_on)
итоговая доходность (не могу добавить изображение из-за низкой репутации):
Коэффициент p-значения
(Перехват) 0.2704
TreeCov500 0,2859
расстояние 0,3412
Вода2000 0,3165
greenvol 0,2194
высота 0,2305
возраст 0,3940
тень 0,2939
экспозиция 0,6039
число_фо_стор 0,0154
смонтировано_onTREES 0,767
mount_onFACADES 0.4319
mount_onPOLES 0.0770
Моя проблема: я хочу отбросить переменную с самым высоким p-значением. mount_onTREES имеет самое высокое значение p. Как я могу отбросить этот единственный аспект? Я не хочу отбрасывать всю переменную mount_on (так как mount_onPOLES имеет очень низкое значение p).
Я пытался:
modelc<-update(modelb,~. -mounted_onTREES)
Но это не сработало. Так есть ли смысл отбрасывать отдельные аспекты одной переменной? Если нет, то как мне обрабатывать переменную mount_on? Должен ли я хранить его, пока один из аспектов имеет очень низкое значение p?
Я очень благодарен за любую помощь!