더 복잡한 모델링
R로 하는 설문 데이터 분석
Kelly McConville
Assistant Professor of Statistics
다중 선형 회귀
다중 선형 회귀
- 다중 선형 회귀 방정식은 다음과 같습니다:
$$E(y) = B_0 + B_1 x_1 + B_2x_2 + \ldots + B_p x_p$$
babies
# A tibble: 484 x 4
AgeMonths HeadCirc WTMEC4YR Gender
<int> <dbl> <dbl> <fct>
1 3 42.7 12915. male
2 4 42.8 12791. female
3 2 38.8 2359. female
4 0 36.0 4306. female
5 5 42.7 2922. female
6 2 41.9 5561. male
7 6 44.3 10416. female
# ... with 477 more rows
다중 선형 회귀
- 다중 선형 회귀 방정식은 다음과 같습니다:
$$E(y) = B_0 + B_1 x_1 + B_2x_2$$
babies
# A tibble: 484 x 4
AgeMonths HeadCirc WTMEC4YR Gender
<int> <dbl> <dbl> <fct>
1 3 42.7 12915. male
2 4 42.8 12791. female
3 2 38.8 2359. female
4 0 36.0 4306. female
5 5 42.7 2922. female
6 2 41.9 5561. male
7 6 44.3 10416. female
# ... with 477 more rows
다중 선형 회귀
babies <- mutate(babies, Gender2 = case_when(
Gender == "male" ~ 1,
Gender == "female" ~ 0))
babies
# A tibble: 484 x 5
AgeMonths HeadCirc WTMEC4YR Gender Gender2
<int> <dbl> <dbl> <fct> <dbl>
1 3 42.7 12915. male 1.
2 4 42.8 12791. female 0.
3 2 38.8 2359. female 0.
4 0 36.0 4306. female 0.
5 5 42.7 2922. female 0.
6 2 41.9 5561. male 1.
7 6 44.3 10416. female 0.
# ... with 477 more rows
다중 선형 회귀
- 다중 선형 회귀 방정식은 다음과 같습니다:
$$E(y) = B_0 + B_1 x_1 + B_2x_2$$
- 남성의 직선:
$$E(y) = (B_0 + B_2) + B_1 x_1$$
- 여성의 직선:
$$E(y) = B_0 + B_1 x_1$$
다중 선형 회귀
mod <- svyglm(HeadCirc ~ AgeMonths + Gender, design = NHANES_design)
summary(mod)
svyglm(formula = HeadCirc ~ AgeMonths + Gender, design = NHANES_design)
Survey design:
svydesign(data = NHANESraw, strata = ~SDMVSTRA, id = ~SDMVPSU,
nest = TRUE, weights = ~WTMEC4YR)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 37.48508 0.18320 204.613 < 2e-16 ***
AgeMonths 1.08658 0.05379 20.200 < 2e-16 ***
Gendermale 1.15034 0.16298 7.058 6.3e-08 ***
(Some output omitted)
다중 선형 회귀
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 37.48508 0.18320 204.613 < 2e-16 ***
AgeMonths 1.08658 0.05379 20.200 < 2e-16 ***
Gendermale 1.15034 0.16298 7.058 6.3e-08 ***
(Some output omitted)
귀무가설: 나이가 모델에 포함된 경우, 성별은 포함되지 않아야 합니다 ($B_2 = 0$).
대립가설: 나이가 모델에 포함된 경우, 성별도 포함되어야 합니다 ($B_2 \neq 0$).
검정 통계량: $t = \frac{b_2}{SE}$
다중 선형 회귀
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 37.48508 0.18320 204.613 < 2e-16 ***
AgeMonths 1.08658 0.05379 20.200 < 2e-16 ***
Gendermale 1.15034 0.16298 7.058 6.3e-08 ***
(Some output omitted)
귀무가설: 성별이 모델에 포함된 경우, 나이는 포함되지 않아야 합니다 ($B_1 = 0$).
대립가설: 성별이 모델에 포함된 경우, 나이도 포함되어야 합니다 ($B_1 \neq 0$).
검정 통계량: $t = \frac{b_1}{SE}$
다중 선형 회귀
$$E(y) = B_0 + B_1 x_1 + B_2x_2$$
연습해 봅시다!
R로 하는 설문 데이터 분석
