ř封裝，帶有偏移項

Question

我在R.使用包segmented分割我第一裝配有使用glm功能，其包括像偏移項日誌鏈接功能的泊松GLM：

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

M1裝配沒有任何錯誤。接下來，我嘗試用R中分割爲包裝，以適應分段glm：

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

在對glm的調用中明確地提供X（例如sample_dat$X）的位置是有用的。

這裏是一個重複的例子：

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)