我正在通过带有 family = cox.ph() 的mgcv包运行GAM模型，并按层（ strata = id ）对我的数据进行分组 . 该数据对应于单个动物的一个使用位置和与该个体相关的20个随机位置可供使用 .
require(mgcv)
require(survival)
require(smoothHR)
gam1 = gam(time1~s(DWL)+strata(id),family=cox.ph(),method = "REML",data=dataset, weight = event1)
模型运行平稳，但我不确定如何绘制与x变量的关系 . DWL 是一个连续变量 . 我使用以下方法来预测图表：
x = seq(0,120) #extent of DWL values
plot(gam1,residuals=T,trans=function(x)exp(x)/(1+exp(x)),shade=T)
我对绘图语法中使用 trans 参数感到有点困惑 . 使用 cox.ph() 为您的家庭agrument， logit-link 是评估您对x变量 DWL 的预测y响应的正确方法吗？
谢谢，
P法瑞尔