Как закодировать Фиксированные эффекты модели Пуассона в R?

Question

Я пытаюсь подогнать модель Пуассона с фиксированными эффектами в R, используя функцию pglm. Мне нужно использовать как индивидуальные, так и фиксированные по времени эффекты в модели. Как я могу включить их обоих? Мой пример данных:

    library(pglm)
    library(readstata13)
    library(lmtest)
    library(MASS)

    ships<-readstata13::read.dta13("http://www.stata- 
    press.com/data/r13/ships.dta")

    ships$lnservice=log(ships$service)

    ships$time <- rep(1:8, 5)

Я пытался использовать следующий код для оценки модели по времени и отдельным фиксированным эффектам:

    res <- pglm(accident ~ 
    op_75_79+co_65_69+co_70_74+co_75_79+lnservice,family = poisson, data = 
    ships, effect = "twoways", model="within", index = c("ship", "time"))
    summary(res)

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

    res <- pglm(accident ~ 
    op_75_79+co_65_69+co_70_74+co_75_79+lnservice,family = poisson, data = 
    ships, effect = "individual", model="within", index = c("ship"))
    summary(res)

Знаете ли вы, где может быть проблема? Я не ожидаю, что две модели будут давать одинаковые оценки.