重复测量的代码，R中的泊松回归？

Question

我正在尝试运行遵循Poisson分布和对数链接的重复测量模型。主题是州，每年我都记录一个结果，这是一个计数。让我们假装其制造的热狗数量。我还需要对州人口使用偏移量，因为我想要每百万人每年制造的热狗数量。

我到处都在寻找，但是找不到正确的代码。

我需要在R中做到这一点-我知道如何在SAS中做到-只需使用proc genmod，重复测量语句和offset语句，并注意相关结构即可。但我无法在这里解决。

我已经尝试过这些，但是不起作用：

fm1 <- gls(hotdogs ~ Unemploy + Ketchup, family=poisson(link = log), offset= lnpop, data=LSss,
           correlation = corAR1(form = ~ 1 | stateID))

glmmPQL(hotdogs ~ Unemploy + Ketchup, ~1 | StateID, family = poisson(link = "log"), offset=lnpop,
    data = LSss)

Answer 1

您可以使用Nz或lme4。例如：

gamlss

或

lme4::glmer(hotdogs ~ offset(log(pop)) + Unemploy + Ketchup + (1|stateID),
            family = poisson, data = LSss)

这假设gamlss::gamlss(hotdogs ~ offset(log(pop)) + Unemploy + Ketchup + random(stateID), family = PO(), data = LSss) 是数据pop中的一列，如响应和预测变量。