我正在使用Tobit I
模型和GAM
包运行cenGAM
回归的示例。
但是,运行程序包手册中提供的示例会给出以下错误:
subset.default(data,Parameter == parm)出错: 对象'参数'未找到
这是我运行的代码:`
library(cenGAM)
# Generate random data
set.seed(1)
x <- matrix(2*rnorm(300), 100)
yn <- 2*x[,3] + 4*cos(x[,1]*2)
y <- yn + rnorm(100)
ycensored <- pmax(y, 0) # data left-censored at 0
ycensored <- pmin(ycensored, 4) # data right-censored at 4
par(mfrow = c(3,3))
# True model
plot(gam(y ~ s(x[,1]) + s(x[,2]) + s(x[, 3])), ylim=c(-5, 5), main = "True")
# Naive estimation
plot(gam(ycensored ~ s(x[,1]) + s(x[,2]) + s(x[, 3])), ylim=c(-5, 5), main = "Naive")
# Tobit I estimation
m <- gam(ycensored ~ s(x[,1]) + s(x[,2]) + s(x[, 3]), family = tobit1(left.threshold=0))
在我运行代码的最后一行以适应Tobit 1
GAM模型之前,一切正常。
任何人都知道如何解决它吗?
我正在使用R 3.4.1