我目前正在尝试在负二项式模型上使用 lsmeans 来获得我数据集中每个参与者的最小二乘均值。我能够通过我感兴趣的曝光来确定 lsmeans,但不能为每个参与者调整均值。到目前为止,我的代码发布在下面:
library(MASS)
library(lsmeans)
m1 <- zeroinfl(exposure ~ BMI_Category,
data = data_final, dist = "negbin", EM = TRUE)
summary(m1)
lsBact <- lsmeans(m1, specs = "BMI_Category")
任何帮助将不胜感激!谢谢!