r - 用偏移项分段的 R 包

Question

我正在使用R 中分段glm的包。我首先使用包含偏移项的函数为 Poisson glm 安装了一个日志链接函数，例如：

M1=glm(Y~X1+X2+X3+offset(log(X)),data=dat.1,family=poisson)

M1 安装没有任何错误。接下来，我尝试glm使用 R 中分段的包来拟合分段：

library(segmented)
seg.1=segmented(M1,seg.Z=~X1,psi=list(X1=c(0.5)))

我有以下错误：

Error in offset(log(X)) : object 'X' not found

我的错误在哪里？非常感谢。

Answer 1

在调用中显式提供X（例如）的位置就可以了。sample_dat$Xglm

这是一个可重现的示例：

library(segmented)

# 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)