Я пытаюсь использовать ols в rms пакете, чтобы соответствовать ограничительному кубическому сплайну c линейной регрессии. Но я не знаю, как использовать ggplot, чтобы нарисовать ограничительную кубическую c сплайн-диаграмму коэффициента регрессии? Как на картинке:
введите описание изображения здесь
А кто-нибудь может мне помочь и поделится кодом? Спасибо большое!
library(ggplot2)
library(rms)
set.seed(100)
x1 <- rnorm(30, 100, 20)
x2 <- rnorm(30, 200, 20)
x3 <- rnorm(30, 300, 20)
y <- rnorm(30, 400, 20)
ddist <- datadist(x1,x2,x3)
options(datadist='ddist')
f <- ols(y ~ rcs(x1,3) + x2 + x3 )
anova(f)
ggplot()