我正在使用带有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()进行家庭评估，logit-link是评估您对x变量DWL的预测y响应的正确方法吗？
谢谢，
P Farrell