我正在通过 mgcv 包运行 GAM 模型，family = cox.ph()并将我的数据按地层（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()您的家庭 agument，logit-link评估您对 x 变量的预测 y 响应的正确方法是DWL什么？
谢谢，
P法雷尔