mgcv パッケージを使用して GAM モデルを実行してfamily = cox.ph()おり、データを階層 ( strata = id) ごとにグループ化しています。データは、個々の動物の 1 つの使用場所と、その個体に関連付けられた使用可能な 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 変数に対する予測された y 応答を評価する適切な方法はありますDWLか?
ありがとうございました、
P・ファレル