어떤 반응값에 2가지의 설명변수만 넣었을때는 유의하지 않았던 인자가 ( PerfRatioMed,ADCepi)
같은 반응값에 4가지의 설명 변수를 넣었더니 유의하게 바뀌었습니다.
confounder effect로 생각해서 다중공선성을 살펴보았는데, 특별한 사항이 없었습니다.
어떻게 이해해야 할지 모르겠습니다.
> fit<-lm(DI0.3~PDR.in.central.ephysis+PerfRatioMed,data=LCP)
> summary(fit)
Call:
lm(formula = DI0.3 ~ PDR.in.central.ephysis + PerfRatioMed, data = LCP)
Residuals:
Min 1Q Median 3Q Max
-0.5291 -0.4136 -0.1998 0.3801 0.8964
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.2260148 0.1346602 1.678 0.1005
PDR.in.central.ephysis 0.0072108 0.0031079 2.320 0.0251 *
PerfRatioMed -0.0015097 0.0007537 -2.003 0.0515 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.453 on 43 degrees of freedom
Multiple R-squared: 0.1945, Adjusted R-squared: 0.1571
F-statistic: 5.192 on 2 and 43 DF, p-value: 0.009555
> fit<-lm(DI0.3~ADR.in.metaphysis+ADCepi,data=LCP)
> summary(fit)
Call:
lm(formula = DI0.3 ~ ADR.in.metaphysis + ADCepi, data = LCP)
Residuals:
Min 1Q Median 3Q Max
-1.1044 -0.2712 -0.1426 0.3596 0.7451
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.002722 0.120065 0.023 0.982018
ADR.in.metaphysis 0.007409 0.001860 3.983 0.000258 ***
ADCepi 0.000441 0.000871 0.506 0.615216
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.4209 on 43 degrees of freedom
Multiple R-squared: 0.3048, Adjusted R-squared: 0.2725
F-statistic: 9.429 on 2 and 43 DF, p-value: 0.0004024
> fit<-lm(DI0.3~PDR.in.central.ephysis+PerfRatioMed+ADR.in.metaphysis+ADCepi,data=LCP)
> summary(fit)
Call:
lm(formula = DI0.3 ~ PDR.in.central.ephysis + PerfRatioMed +
ADR.in.metaphysis + ADCepi, data = LCP)
Residuals:
Min 1Q Median 3Q Max
-0.83068 -0.20354 -0.08403 0.27185 0.64545
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.2357246 0.1359272 -1.734 0.09040 .
PDR.in.central.ephysis 0.0076449 0.0024446 3.127 0.00324 **
PerfRatioMed -0.0020127 0.0006146 -3.275 0.00215 **
ADR.in.metaphysis 0.0064848 0.0015599 4.157 0.00016 ***
ADCepi 0.0017448 0.0007808 2.235 0.03094 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.35 on 41 degrees of freedom
Multiple R-squared: 0.5416, Adjusted R-squared: 0.4969
F-statistic: 12.11 on 4 and 41 DF, p-value: 1.373e-06
> vif(fit)
PDR.in.central.ephysis PerfRatioMed ADR.in.metaphysis ADCepi
1.046292 1.124601 1.109561 1.267423
> sqrt(vif(fit))>2
PDR.in.central.ephysis PerfRatioMed ADR.in.metaphysis ADCepi
FALSE FALSE FALSE FALSE
Comment 3
-
cardiomoon
2015.04.17 18:06
-
cardiomoon
2015.04.17 18:33
제 생각에는 PDR.in.central.ephysis와 PerfRatioMed 중 하나가 ADCepi가 미치는 영향을 상쇄하지 않나 생각됩니다,. ADCepi와 PDR.in.central.ephysis, ADCepi와 PerfRatioMed 두개씩 해보실 것을 권합니다.
또한 다중회귀분석 전에 변수를 하나씩 넣은 회귀모형을 만들어 회귀모형진단을 해보실 것을 권합니다. 그리고 변수 4개가 있는 최종 모형도 회귀모형진단을 해보세요. web-r에 회귀모형 진단 부분은 아주 잘되어있읍니다. 회귀분석하실때 변수를 넣으시고 오른쪽위 체크박스를 모두 체크하시면 회귀모형진단을 할 수 있읍니다.
-
김승욱
2015.04.20 18:24
첫 번째와 두 번째 모델이 상당히 불안불안 하네요... 로버스트 회귀 정도 써보심이...
No. | Subject | Author | Date |
---|---|---|---|
1053 | 통계와 관련한 기초적인 질문 드립니다.. [2] | 감사합니다 | 2021.12.10 |
1052 | 메타분석 해석 문의 [2] | 루아로 | 2021.12.07 |
1051 | 초보라 질문드립니다 ㅜ [2] | 쵸쿄쵸쿄 | 2021.12.06 |
1050 | ㅠㅠㅠ 초보 질문입니다 [1] | 만고 | 2021.11.26 |
1049 | 라벨붙이기가 안됩니다 [1] | 찬미맘 | 2021.11.22 |
1048 | moonBook ANOVA 질문드립니다. [1] | 타쿠 | 2021.11.20 |
1047 | 3원 분할표? 에서 모델 만드는법 | 플랫슬립 | 2021.11.20 |
1046 | [회귀분석 초급 질문] [1] | 서빱스 | 2021.11.14 |
1045 | [ANOVA 분석 질문] [1] | 서빱스 | 2021.11.14 |
1044 | 대조군1과 치료군2,3,4 의 PSM 에 관해서 질문있습니다. [4] | yessul | 2021.11.11 |
1043 | 서버접속이 안됩니다. Reload [2] | silkywater | 2021.11.11 |
1042 | 메타분석 실행이 안됩니다. [1] | 애플K | 2021.11.08 |
1041 | pair plot 관련하여 질문드립니다. | ulysses11 | 2021.10.29 |
1040 | 메타분석 결과파일 다운로드 하려면 [1] | mono | 2021.10.23 |
1039 | 기본 baseline character table문의드립니다. [1] | ulysses11 | 2021.10.22 |
1038 | 이 error는 뭐가 문제인 걸까요 [1] | CSJU | 2021.10.19 |
1037 | logistic regression에서 p>0.05인 변수의 OR, p 값... | 브이콘65 | 2021.10.14 |
1036 | R 문제가 도저히 안풀려요.. [1] | 자하연미필 | 2021.10.09 |
1035 | 구조방정식에서 종속변수가 0,1로 코딩된 데이터 사용이 가능한가요? [1] | 순수엄마 | 2021.10.08 |
1034 | roc.test를 통해 AUC 비교 [2] | 나현둥 | 2021.10.07 |
재미있는 경우인데요. 네가지 설명변수가 모두 임상적으로는 의미가 있는 건가요? 서로 독립적이구요?
나온 결과로만 해석하면 이렇게 할수 있겠군요.
먼저 첫번쨰 모형입니다. PDR.in.central.ephysis와 PerfRatioMed는 모두 상당한 의미가 있어보입니다. PerfRatioMed의 p값이 0.0515인것은 0,05이하가 아니라고 의미가 없는 것이 아니라 상당히 의미가 있어보입니다. 또한 모형의 adjusted-R2값을 보면 0.1571로 이 모형이 반응변수(종속변수) DIO.3의 15.7%를 설명하고 있음을 알 수 있읍니다.
두번째 모형인 ADR.in.metaphysis + ADCepi를 넣은 모형의 경우 adjusted-R2값이 0.2725로 두개의 영향변수를 넣은 모형이 DIO.3의 27.3%를 를 설명하고 있읍니다. ADCepi의 p값은 0.61로 의미가 없어보입니다. 두개의 영향변수를 넣고 다중회귀분석을 할때 해석은 다음과 같이 할수 있읍니다. ADR.in.metaphysis가 일정할때 ADCepi의 변동은 종속변수 DIO.3에 의미있는 영향을 미치지 않는다고 할수 있읍니다.
세번쨰 모형은 네개의 영향변수를 넣었더니 adjusted-R2 값이 무려 0.4969입니다. 즉 네개의 영향변수로 만든 회귀모형이 반응변수(종속변수)인 DIO.3의 49.7%나 설명하고 있읍니다. 그리고 네개의 영향변수가 모두 유의하게 나옵니다. 해석을 어떻게 해야 할까요? 원래 의미가 있었던 세개의 변수 말고 ADCepi를 중심으로 생각해보면 세개의 영향변수 PDR.in.epiphysis, PerfRatioMed, ADR.in.metaphysis 가 일정하다면 ADCepi가 변하면 DIO.3이 의미있게 변한다는 것입니다.