Я хотел бы сгенерировать значения доверительного интервала для некоторых пояснительных переменных переменной ответа (называемой HUNTED). Данные о 3 различных периодах времени.
Тем не менее, когда я использую код:
библиотека (МАСС)
confint (allspecies4)
Я получаю только сообщения об ошибках (см. Ниже). Кто-нибудь знает, как улучшить мой сценарий, чтобы я мог получить доступ к доверительному интервалу?
Obs .: Мой набор данных слишком большой, я только что добавил его часть, чтобы вы могли понять, как он распространяется.
allspecies4<-read.csv(file= "dadosGLMM.csv", header= TRUE, sep= ";" )
attach(allspecies4)
library(lme4)
allspecies4<- glmer(HUNTED~ eco+TrophicLevel+age+ log(body_mass)+
habitat+ABUND + PERIOD + tabu +(1|Specie),
data=dadosGLMM, family=binomial)
library(MASS)
confint(allspecies4)
library(MASS)
> confint(allspecies4)
Computing profile confidence intervals ...
Error in zeta(shiftpar, start = opt[seqpar1][-w]) :
profiling detected new, lower deviance
In addition: Warning messages:
1: In nextpar(mat, cc, i, delta, lowcut, upcut) :
unexpected decrease in profile: using minstep
2: In nextpar(mat, cc, i, delta, lowcut, upcut) :
unexpected decrease in profile: using minstep
3: In nextpar(mat, cc, i, delta, lowcut, upcut) :
unexpected decrease in profile: using minstep
4: In nextpar(mat, cc, i, delta, lowcut, upcut) :
unexpected decrease in profile: using minstep
5: In (function (fn, par, lower = rep.int(-Inf, n), upper = rep.int(Inf, :
failure to converge in 10000 evaluations
6: In nextpar(mat, cc, i, delta, lowcut, upcut) :
unexpected decrease in profile: using minstep
7: In nextpar(mat, cc, i, delta, lowcut, upcut) :
unexpected decrease in profile: using minstep
8: In nextpar(mat, cc, i, delta, lowcut, upcut) :
unexpected decrease in profile: using minstep
9: In FUN(X[[i]], ...) : non-monotonic profile for (Intercept)
Пример данных
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