Определение функции-оболочки в glmulti - PullRequest
0 голосов
/ 30 апреля 2020

Я новичок в сценариях R и хотел бы запустить glmulti с функцией-оберткой для отрицательного биномиального семейства со смещением. Я хотел бы получить совет по следующему коду, так как я получаю сообщение об ошибке. Переменные a, b и c являются факторами и ковариатами. Переменная e используется в качестве смещения.

glmnb.glmulti <- function(formula, data, offset, ...)
  glm.nb(formula, offset=offset, data=data, ...)

selection<- glmulti(y~ a + b + c + d,level=2, 
fitfunction=glmnb.glmulti, offset=log(e), crit="aic")

Initialization...
 Error in model.frame.default(formula = formula, data = data, offset = offset,  : 
  invalid type (closure) for variable '(offset)' 

Любой совет будет высоко ценится.

1 Ответ

0 голосов
/ 01 мая 2020

Спасибо большое! Я также попробовал следующее, и это также работало

eqt <- as.formula (y ~ a + b + c + d) </p>

glmnb <- функция (формула, данные, веса) , ...) glm.nb (формула, данные = данные, веса = смещение (log (e)), ...) </p>

выбор <- glmulti (например, данные = данные, fitfunction = glmnb, крит = "ай c") </p>

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...