使用R包segmentet进行带偏移的泊松回归

时间:2014-07-31 14:47:13

标签: r regression glm piecewise

我试图在带有R(版本3.1.1)偏移的泊松回归上使用分段函数。我的分段包我的版本是0.4-0.0。

在早期post中,给出了一个可重现的例子:

# 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

还有其他人有同样的问题吗?或者甚至更好 - 一个解释?

谢谢

/基拉

0 个答案:

没有答案