R语言gamRR包中的gamRR()函数,在R中用help(‘gamRR’)查看,对ref解释是:a vector of the independent variables at referenced level。下面举例是
require("mgcv")
dat <- gamSim(1,100,dist="poisson",scale=.25)
fit <- gam(y~s(x0)+s(x1)+s(x2)+s(x3),family=poisson,dat,method="REML")
plot(fit,select=2)
gamRR(
fit=fit,
ref=c(x0=dat$x0[1],x1=dat$x1[1],x2=dat$x2[1],x3=dat$x3[1]),
est="x1",
data=dat,
n.points=10,
plot=TRUE,
ylim=NULL)
我用gam拟合,后想计算RR for gam,数据是fit <- gam(二分类结局变量ys(连续变量自变量x, k=7)+age+sex, data=dat, family = binomial())写到ref这里有点疑惑
ref=c(x=dat$x[1],age=dat$age[1],sex=dat$sex[1]) 是这样理解gamRR()中的ref吗?