Я хочу запустить GLM по Гауссу со ссылкой на журнал и смещением. Возникают следующие проблемы:
y <- c(1,1,0,0)
t <- c(5,3,2,4)
Нет проблем:
exp(coef(glm(y~1 + offset(log(t)), family=poisson)))
с family = gaussian
, начальные значения необходимо указать, он работает здесь:
exp(coef(glm(y~1, family=gaussian(link=log), start=0)))
, но не работает здесь:
exp(coef(glm(y~1 + offset(log(t)), family=gaussian(link=log), start=0)))
Ошибка в eval (expr, envir, enclos): не удается найти допустимые начальные значения: укажите некоторые "
Кто-нибудь видит, что не так (надеюсь, только в моем коде)?