3 차 스플라인이있는 lm ()의 출력을 회귀 방정식으로 변환하는 방법

코드와 출력이 있으며 모델을 만들고 싶습니다. 이 출력을 사용하여 모델을 구성하는 방법을 모르겠습니다.

 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)

참고 ns()천연 큐빅 스플라인 대한 B 스플라인 기저 행렬을 생성한다. 따라서이 모형 은 3 자유도 y를 x사용 하기 위해 B- 스플라인에 대해 회귀 합니다. 그러한 모델의 방정식은 어떻게 생겼습니까?

r splines

— 아멜리아
소스

이 문제가 여기에 관심이 있지만,이 질문은 R당신이 무엇을 설명하지 않는 한 중심적인 방식으로 너무 많이 언급 되므로 SO에 속합니다 ns. ( R패키지의 일부는 아닙니다 .)

— whuber

@whuber, 여기 참조 : ? ns ; ns()스플라인 패키지의 일부입니다. 나는이 질문이 R 용어에 속한다는 것을 알고 있지만 여기에 주제가 있다는 의견이 있습니다.

— gung-모니 티 복원

@ gung 예, 패키지도 찾을 수 있었지만 요점은 아닙니다.이 질문을 여기에 유지하려면 R이 아닌 사용자도 이해할 수 있어야합니다.

— whuber

@ whuber 나는 약간의 설명을 추가했습니다. 이보다 더 찾고 계십니까?

— 복원 Monica Monica-G. Simpson

@Gavin 감사합니다. 나는 R이 아닌 사용자가 질문을 이해하고 (따라서 의미있는 답변을 제공 할 수 있도록) 한 줄을 더 추가 할 자유를 얻었습니다.

— whuber

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 "제한된 큐브 스플라인"은 자연 스플라인의 또 다른 표현입니다.

— 프랭크 하렐
소스

인상적입니다. 그러나 나는이 사이트에서 찾은 대답 (SO와 반대)이 일반적으로 스플라인의 방정식을 결정하는 방법을 설명 할 것이라고 생각합니다.

— whuber

감사합니다! f에 의해 생성 된 계수가 피팅에 의해 생성 된 계수와 다른 이유를 설명 할 수 있습니까?

— AmeliaBright

nsrcs

X

$X$

X^{2}

$X^{2}$

X - \bar{X}

$X-\bar{X}$

(X - \bar{X})^{2}

$(X-\bar{X})^{2}$