2017-10-19 7 views

답변

0

당신은 자신을 쉽게 테스트 할 수 있습니다 :

DF <- data.frame(x1 = 1:2, x2 = 3:4, x3 = 5:6) 
with(DF, poly(x1, x2, x3, degree = 2, raw = TRUE)) 
#  1.0.0 2.0.0 0.1.0 1.1.0 0.2.0 0.0.1 1.0.1 0.1.1 0.0.2 
#[1,]  1  1  3  3  9  5  5 15 25 
#[2,]  2  4  4  8 16  6 12 24 36 
#attr(,"degree") 
#[1] 1 2 1 2 2 1 2 2 2 
#attr(,"class") 
#[1] "poly" "matrix" 

열 이름은 세 가지 변수의 제품과 각 변수는이 제품에이 정도를 보여줍니다. 예 : 1.1.0x1^1 + x2^1 + x3^0을 의미합니다.

물론 회귀 모델 출력에서도이를 볼 수 있습니다.

계수가 원시 다항식에 해당하도록하려면 raw = TRUE이 필요합니다 (즉, alpha0 + alpha11 * x1^1 + alpha12 * x1^2 + ...). 필요하지 않으면 이 아닌raw = TRUE으로 설정되어야합니다. 왜냐하면 회귀 분석에 바람직한 특성을 가진 orthogonal polynomials이 있기 때문입니다.