R: как я могу построить кривые отклика переменной GLM успеха в сравнении с пробной моделью, используя biomod2 и ggplot2? - PullRequest
0 голосов
/ 31 января 2019

Я использую R и GLM для подгонки некоторых данных:

glm_out <- glm(presence ~ 1+ data1 + data2,family=binomial (logit), data=dataFrame)

Мне удается это сделать успешно, а затем я хочу получить кривые отклика кофакторов, используя biomod2:

rp <- response.plot2(model = "glm_out", Data = dataFrame[, c("data1","data2")],
show.variables = c("data1","data2"), fixed.var.metric = "mean", 
                 plot = TRUE, use.formal.names = TRUE)

, который также работает.Теперь я хотел бы получить такие же кривые отклика, но для данных испытаний и успеха.GLM может быть получен следующим образом:

glm_out <- glm(cbind(success,failure) ~ 1+ data1 + data2,family=binomial (logit), data=dataFrame)

Когда я пытаюсь построить кривые отклика, используя те же методы, я получаю эту ошибку:

Warning message in assign(mod_tmp@model_name, mod_tmp, envir = parent.frame(n = 1)):
“only the first element is used as variable name”

Error in get(model): object 'success_AllData_57.892408_GLM' not found

Любая идея, как решить эту проблему

...