Cómo traducir la salida de un lm() encaja con un spline cúbico en una ecuación de regresión

Question

Tengo algún código y salida, y me gustaría construir un modelo. No sé cómo construir un modelo con esta salida:

 require("splines")
 x   <- c(0.2,   0.23,   0.26,   0.29,   0.33,   0.46,    0.53 )
 y   <- c(0.211, 0.2026, 0.2034, 0.2167, 0.2177, 0.19225, 0.182)
 fit <- lm(y ~ ns(x,3))
 summary(fit)

Tenga en cuenta que ns() genera la matriz de base de B-spline para un spline cúbico natural. Así este modelo regresa y contra una B-spline para x con tres grados de libertad. ¿Lo que le gustaría la ecuación para una mirada de modelo?

Answer 1

require(rms)
f <- ols(y ~ rcs(x, 3))  # 2 d.f. for x
Function(f)  # represent fitted function in simplest R form
latex(f)     # typeset algebraic representation of fit

RCS "spline cublic restringido" es otra representación de un spline natural.