如何在R中编码固定效应Poisson模型？

Question

我正在尝试使用pglm函数在R中拟合固定效应Poisson模型。我需要在模型中同时使用个人和时间固定效果。如何同时包含两者？ 我的示例数据是：

    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)

您知道问题可能出在哪里吗？我不希望两个模型产生相同的估计。