Я пытаюсь создать модель в R, которая учитывает мой опрос.
В исследовании участвуют многие объясняющие переменные, ответная переменная называется HUNTED, а данные о 3 разных времени PERIODS.
Мой набор данных слишком большой, я только что добавил его часть, чтобы вы могли понять, как он распространяется.
Мой вопрос: если я оставлю переменную PERIOD в скрипте, как я это сделал, программа "поймет", что на уровне PERIOD, называемом "до", есть было меньше охотящихся видов, чем на уровне PERIOD, называемом «В ТЕЧЕНИЕ»? И поэтому из-за того, что уровень DURING кажется более значимым, или я неправильно понимаю результаты?
Obs .: Уровень PERIOD, называемый "AFTER", не отображается в консоли, поскольку он используется как базовый уровень.
Результаты, которые я получил в консоли:
Random effects:
Groups Name Variance Std.Dev.
Specie (Intercept) 1.728 1.314
Number of obs: 5622, groups: Specie, 25
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -1.611e+01 2.510e+00 -6.420 1.36e-10 *
eco 9.337e-02 6.522e-02 1.432 0.152254
TrophicLevel 4.403e-01 3.561e-01 1.236 0.216324
age 8.676e-04 2.616e-03 0.332 0.740175
log(body_mass) 8.717e+00 1.495e+00 5.833 5.46e-09
habitat -5.238e-01 8.269e-02 -6.335 2.37e-10
ABUND 4.432e-01 4.938e-02 8.976 < 2e-16 *
BEFORE -2.850e-01 1.090e-01 -2.615 0.008932 **
DURING 3.478e-01 9.103e-02 3.821 0.000133 ***
taboo -1.057e+00 1.016e+00 -1.041 0.297926
Мой код
allspecies4<-read.csv(file= "dataGLMM.csv", header= TRUE, sep= ";" )
attach(allspecies4)
library(lme4)
allspecies4<- glmer(HUNTED~ eco+TrophicLevel+age+ log(body_mass)+
habitat+ABUND + PERIOD + tabu +(1|Specie),
data=dataGLMM, family=binomial)
summary(allspecies4)
Пример данных
Specie body_mass TrophicLevel age eco habitat PERIOD HUNTED ABUND tabu
Cerc_mit 3.70 2 57 2 F DURING 0 2 1
Cerc_mit 3.70 2 57 2 F AFTER 0 3 1
Cerc_mit 3.70 2 57 2 F BEFORE 0 4 1
Cerc_mit 3.70 2 67 2 F DURING 1 2 1
Cerc_mit 3.70 2 67 2 F AFTER 0 2 1
Chlor_cyn 3.70 2 53 2 S DURING 0 3 0
Chlor_cyn 3.70 2 74 2 S DURING 0 3 0
Chlor_cyn 3.70 2 30 2 S DURING 0 4 0
Chlor_cyn 3.70 2 63 2 S DURING 0 4 0
Chlor_cyn 3.70 2 54 2 S DURING 0 3 0
Chlor_cyn 3.70 2 30 2 S DURING 0 4 0
Chlor_cyn 3.70 2 30 2 S DURING 0 3 0
Phil_mont 3.69 2 24 3 F DURING 0 3 0
Phil_mont 3.69 2 24 3 F BEFORE 0 4 0
Phil_mont 3.69 2 33 3 F AFTER 1 4 0
Phil_mont 3.69 2 33 3 F BEFORE 0 4 0
Phil_mont 3.69 2 33 3 F DURING 0 4 0
Phil_mont 3.69 2 43 3 F BEFORE 0 4 0
Phil_mont 3.69 2 43 3 F DURING 0 4 0