モデルの非線形性と標準誤差帯域による制限された 3 次スプラインを含む回帰の予測線をプロットしたいと思います。予測されたポイントを取得できますが、線とエラー バンドをプロットするかどうかはわかりません。ggplot が推奨されますが、ベース グラフィックスでも問題ありません。ありがとう。
ドキュメントの例を次に示します。
library(rms)
# Fit a complex model and approximate it with a simple one
x1 <- runif(200)
x2 <- runif(200)
x3 <- runif(200)
x4 <- runif(200)
y <- x1 + x2 + rnorm(200)
f <- ols(y ~ rcs(x1,4) + x2 + x3 + x4)
pred <- fitted(f) # or predict(f) or f$linear.predictors
f2 <- ols(pred ~ rcs(x1,4) + x2 + x3 + x4, sigma=1)
fastbw(f2, aics=100000)
options(datadist=NULL)
モデルの予測値のプロット:
plot(predict(f2))
