我正在尝试在 R 中使用偏移量的泊松回归(版本 3.1.1)使用分段函数。我的分段包 i 版本 0.4-0.0。
在较早的帖子中,给出了一个可重现的示例:
# sample data
set.seed(101)
sample_dat <- data.frame(Y = rpois(100, 3), X1 = rnorm(100), X2 = rnorm(100), X = rpois(100, 100))
# fit model
M1 = glm(Y ~ X1 + X2, data = sample_dat, offset = log(sample_dat$X), family = poisson)
seg.1 = segmented(M1, seg.Z=~X1, psi=list(X1 = c(0.5)), data = sample_dat)
但是我无法重现该示例 - 我收到一条错误消息:
Error in model.frame.default(formula = Y ~ Y + sample_dat + Y ~ X1 + X2 + : object is not a matrix
其他人有同样的问题吗?甚至更好 - 解释?
谢谢
/基拉