1

我正在尝试在 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

其他人有同样的问题吗?甚至更好 - 解释?

谢谢

/基拉

4

0 回答 0