요인 설계를 위한 이항 반응 분석에 대한 주요 결과 해석

요인 설계를 분석하려면 다음 단계를 수행하십시오. 주요 결과에는 Pareto 차트, p-값, 계수, 모형 요약 통계량 및 잔차 그림이 포함됩니다.

1단계: 반응에 가장 큰 영향을 미치는 항 확인

표준화된 효과의 Pareto 차트를 사용하여 주효과와 교호작용 효과의 상대적 크기 및 통계적 유의성을 비교할 수 있습니다.

Minitab은 각 표준화된 효과의 절대값을 내림차순으로 표시합니다. 이때 차트의 기준선은 유의한 효과를 나타냅니다. 기본적으로 Minitab에서는 0.05의 유의 수준을 사용하여 기준선을 그립니다.

주요 결과: Pareto 차트

이 결과에서는 방부제 유형(A), 진공 밀봉 압력(B), 오염 수준(C) 등 세 가지 주효과가 통계적으로 유의합니다(α = 0.05).

또한 방부제 유형(A) 막대의 길이가 가장 길기 때문에 이 효과가 가장 크다는 것을 알 수 있습니다. 마찬가지로 방부제*냉각 온도 교호작용(AD) 막대의 길이가 가장 짧기 때문에 이 효과가 가장 작습니다.

2단계: 반응에 통계적으로 유의한 영향을 미치는 항 확인

반응과 모형의 각 항 간의 연관성이 통계적으로 유의한지 여부를 확인하려면 항에 대한 p-값을 유의 수준과 비교하여 귀무 가설을 평가합니다. 귀무 가설은 항의 계수가 0으로, 항과 반응 간에 연관성이 없다는 것을 나타냅니다. 일반적으로 0.05의 유의 수준(α 또는 알파로 표시됨)이 적절합니다. 0.05의 유의 수준은 실제로 연관성이 없는데 연관성이 존재한다는 결론을 내릴 위험이 5%라는 것을 나타냅니다.
p-값 ≤ α: 연관성이 통계적으로 유의함
p-값이 유의 수준보다 작거나 같으면 반응 변수와 항 간에 통계적으로 유의한 연관성이 있다는 결론을 내릴 수 있습니다.
p-값 > α: 연관성이 통계적으로 유의하지 않음
p-값이 유의 수준보다 크면 반응 변수와 항 간에 통계적으로 유의한 연관성이 있다는 결론을 내릴 수 없습니다. 항 없이 모형을 다시 적합시킬 수도 있습니다.
반응과 통계적으로 유의한 연관성이 없는 예측 변수가 여러 개 있는 경우 한 번에 하나씩 항을 줄여 모형을 축소할 수 있습니다. 모형에서 항을 제거하는 방법은 모형 축소에서 확인하십시오.
모형 항이 통계적으로 유의하면 해석은 항의 유형에 따라 다릅니다. 해석은 다음과 같습니다.
  • 요인에 대한 계수가 유의하면 사건의 확률이 요인의 모든 수준에 대해 동일하지는 않다는 결론을 내릴 수 있습니다.
  • 공변량에 대한 계수가 유의하면 변수 값의 변화가 확률의 변화와 연관성이 있습니다.
  • 교호작용 항에 대한 계수가 유의하면 요인과 반응의 관계가 항의 요인에 따라 다릅니다. 이 경우에는 교호작용 효과를 고려하지 않고 주효과를 해석해서는 안 됩니다.
  • 블럭에 대한 계수가 통계적으로 유의하면 블럭에 대한 연결 함수가 평균 값과 다르다는 결론을 내릴 수 있습니다.

코드화된 계수

효과계수SE 계수VIF
상수  -2.73700.0479 
방부제0.44970.22490.04771.03
진공 압력0.25740.12870.04771.06
오염 수준0.29540.14770.04781.06
냉각 온도-0.1107-0.05540.04781.07
방부제*진공 압력-0.0233-0.01170.04731.05
방부제*오염 수준0.07220.03610.04741.06
방부제*냉각 온도0.00670.00340.04721.05
진공 압력*오염 수준-0.0430-0.02150.04691.04
진공 압력*냉각 온도-0.0115-0.00580.04651.02
오염 수준*냉각 온도0.15730.07860.04671.02
주요 결과: 계수

이 결과에서는 방부제, 진공 압력 및 오염 수준의 주효과에 대한 계수가 양수입니다. 냉각 온도의 주효과에 대한 계수는 음수입니다. 일반적으로 계수가 양이면 항의 값이 증가함에 따라 사건이 발생할 가능성이 더 높고 계수가 음이면 사건이 발생할 가능성이 낮습니다.

분산 분석

출처DF수정 분산수정 평균카이-제곱P-값
모형1046.21304.621346.210.000
  방부제122.683522.683522.680.000
  진공 압력17.33137.33137.330.007
  오염 수준19.62099.62099.620.002
  냉각 온도11.34411.34411.340.246
  방부제*진공 압력10.06080.06080.060.805
  방부제*오염 수준10.57800.57800.580.447
  방부제*냉각 온도10.00510.00510.010.943
  진공 압력*오염 수준10.21060.21060.210.646
  진공 압력*냉각 온도10.01530.01530.020.902
  오염 수준*냉각 온도12.84752.84752.850.092
오차50.96740.1935   
총계1547.1804     
주요 결과: p-값

이 결과에서는 방부제, 진공 압력 및 오염 수준의 주효과가 유의 수준 0.05에서 통계적으로 유의합니다. 이러한 변수의 변화가 반응 변수의 변화와 연관되어 있다는 결론을 내릴 수 있습니다.

교호작용 항은 통계적으로 유의하지 않습니다. 각 변수와 반응 간의 관계가 다른 변수의 값에 종속되지 않을 수도 있습니다.

3단계: 예측 변수의 효과 이해

예측 변수의 효과를 이해하려면 승산비를 사용합니다. 승산비에 대한 해석은 예측 변수가 범주형인지 아니면 계량형인 지에 따라 달라집니다. 모형에서 로짓 연결 함수를 사용하는 경우 Minitab에서는 승산비를 계산합니다.
계량형 예측 변수에 대한 승산비

승산비가 1보다 크면 예측 변수가 증가함에 따라 사건 발생 확률이 증가한다는 것을 나타냅니다. 승산비가 1보다 작으면 예측 변수가 증가함에 따라 사건 발생 확률이 감소한다는 것을 나타냅니다.

계량형 예측 변수에 대한 승산비

변경 단위승산비95% CI
투여량(mg)0.56.1279(1.7218, 21.8087)
주요 결과: 승산비

이 결과에서 모형은 한 약품의 투여량 수준을 사용하여 성인의 박테리아 존재 여부를 예측합니다. 박테리아가 존재하지 않는 것이 사건입니다. 각 약에는 0.5mg의 투여량이 포함되어 있으며, 따라서 연구자들은 0.5mg의 단위 변화를 사용합니다. 승산비는 약 6입니다. 성인이 약 하나를 추가로 복용할 때마다 환자에게 박테리아가 없을 확률은 약 6배 증가합니다.

범주형 예측 변수에 대한 승산비

범주형 예측 변수의 경우 승산비는 예측 변수의 서로 다른 2개 수준에서 사건이 발생할 확률을 비교합니다. Minitab에서는 수준 A와 수준 B, 두 열에 수준을 나열하여 비교를 설정합니다. 수준 B는 요인에 대한 기준 수준입니다. 승산비가 1보다 크면 수준 B에서 사건 발생 확률이 증가한다는 것을 나타냅니다. 승산비가 1보다 작으면 수준 B에서 사건 발생 확률이 감소한다는 것을 나타냅니다. 범주형 예측 변수의 코드화에 대한 자세한 내용은 범주형 예측 변수의 코드화 방법에서 확인하십시오.

범주형 예측 변수에 대한 승산비

수준 A수준 B승산비95% CI
     
  211.1250(0.0600, 21.0834)
  313.3750(0.2897, 39.3165)
  417.7143(0.7461, 79.7592)
  512.2500(0.1107, 45.7172)
  616.0000(0.5322, 67.6397)
  323.0000(0.2547, 35.3325)
  426.8571(0.6556, 71.7169)
  522.0000(0.0976, 41.0019)
  625.3333(0.4679, 60.7946)
  432.2857(0.4103, 12.7323)
  530.6667(0.0514, 8.6389)
  631.7778(0.2842, 11.1200)
  540.2917(0.0252, 3.3719)
  640.7778(0.1464, 4.1326)
  652.6667(0.2124, 33.4861)
수준 B에 상대적인 수준 A에 대한 승산비
주요 결과: 승산비

이 결과에서 범주형 예측 변수는 호텔의 성수기가 시작된 이후의 달입니다. 반응은 투숙객이 예약을 취소하느냐 여부입니다. 이 예에서는 취소가 사건입니다. 최대 승산비는 수준 A가 월 4이고 수준 B가 월 1일 때 약 7.71입니다. 이는 월 4에 투숙객이 예약을 취소할 확률이 월 1에 예약을 취소할 확률보다 약 8배 높다는 것을 나타냅니다.

4단계: 모형이 데이터를 얼마나 잘 적합시키는지 확인

모형이 데이터를 얼마나 잘 적합시키는지 확인하려면 모형 요약 표의 적합도 통계량을 조사합니다.
참고

많은 모형 요약 및 적합도 통계량은 워크시트에서 데이터가 배열되는 방식 및 행당 시행 횟수에 영향을 받습니다. Hosmer-Lemeshow 검정은 데이터가 배열되는 방식에 영향을 받지 않으며 행당 시행 횟수에 관계없이 유사합니다. 자세한 내용은 이항 로지스틱 회귀 분석에서 데이터 형식이 적합도에 미치는 영향에서 확인하십시오.

이탈도 R-제곱

이탈도 R2이 클수록 모형이 데이터를 더 잘 적합시킵니다. 이탈도 R2은 항상 0%에서 100% 사이입니다.

모형에 항을 추가하면 이탈도 R2은 항상 증가합니다. 예를 들어, 최량 항이 5개인 모형은 최량 항이 4개인 모형보다 항상 R2 값이 큽니다. 따라서 이탈도 R2은 같은 크기의 모형을 비교할 때 가장 유용합니다.

데이터 배열은 이탈도 R2 값에 영향을 미칩니다. 일반적으로 행당 시행 횟수가 여러 번인 데이터에 대한 이탈도 R2이 행당 시행 횟수가 한 번인 데이터보다 큽니다. 이탈도 R2 값은 동일한 데이터 형식을 사용하는 모형 간에만 유사합니다.

적합도 통계량은 모형이 데이터를 얼마나 잘 적합시키는 지에 대한 하나의 측도에 지나지 않습니다. 모형에 바람직한 값이 있더라도 해당 모형이 데이터를 충족하는지 확인하려면 잔차 그림 및 적합도 검정을 확인해야 합니다.

이탈도 R-제곱(수정)

항 수가 다른 여러 모형을 비교하려면 수정 이탈도 R2을 사용하십시오. 모형에 항을 추가하면 이탈도 R2은 항상 증가합니다. 수정 이탈도 R2 값은 모형의 항 수에 통합되어 올바른 모형을 선택하는 데 유용합니다.

AIC, AICc 및 BIC

여러 모형을 비교하려면 AIC, AICc 및 BIC를 사용합니다. 각 통계량에 대해 작은 값을 사용하는 것이 바람직합니다. 그러나 예측 변수 집합에 대한 값이 가장 작은 모형이 반드시 데이터를 잘 적합시키는 것은 아닙니다. 모형이 데이터를 얼마나 잘 적합시키는지 평가하려면 적합도 검정과 잔차 그림도 사용하십시오.

모형 요약

이탈도 R-Sq이탈도 R-Sq(수정)AICAICcBIC
97.95%76.75%105.98171.98114.48
주요 결과: 이탈도 R-제곱, 이탈도 R-제곱(수정), AIC

이 결과에서 모형은 반응 변수에 있는 총 이탈도의 97.95%를 설명합니다. 이러한 데이터의 경우 이탈도 R2 값은 모형이 데이터를 잘 적합시킨다는 것을 나타냅니다. 다른 예측 변수를 사용하여 추가 모형이 적합되는 경우 수정된 이탈도 R2 값, AIC 값, AICc 값 및 BIC 값을 사용하여 모형들이 데이터를 얼마나 잘 적합시키는지 비교하십시오.

5단계: 모형이 데이터에 적합하지 않는지 여부 확인

적합도 검정은 예측 확률이 이항 분포에서 예측하지 않는 방식으로 관측된 확률에서 벗어나는지 확인하기 위해 사용합니다. 적합도 검정의 p-값이 선택한 유의 수준보다 낮으면 예측 확률이 이항 분포에서 예측하지 않는 방식으로 관측된 확률에서 벗어남을 의미합니다. 다음 리스트에는 이러한 편차의 일반적인 이유가 나와 있습니다.
  • 잘못된 연결 함수
  • 모형에 있는 변수에 대한 고차 항 제외됨
  • 모형에 없는 예측 변수 제외됨
  • 과대산포

편차가 통계적으로 유의한 경우 다른 연결 함수를 사용하거나 모형의 항을 변경할 수 있습니다.

다음 통계량은 적합도를 검정합니다. 이탈도 및 Pearson 통계량은 워크시트에서 데이터가 배열되는 방식 및 행당 시행 횟수의 영향을 받습니다.
  • 이탈도: 행당 단일 시행 형식으로 배열된 데이터에 대한 이탈도 검정의 p-값은 행당 시행 횟수가 여러 번인 데이터에 비해 작은 경향이 있으며, 일반적으로 행당 시행 횟수가 감소함에 따라 감소합니다. 행당 단일 시행 형식 데이터의 경우 Hosmer-Lemeshow 결과를 더 신뢰할 수 있습니다.
  • Pearson: Pearson 검정에서 사용하는 카이-제곱 분포에 대한 근사는 데이터의 행당 기대되는 사건 수가 작은 경우 부정확합니다. 그러므로, 데이터가 행당 단일 시행 형식인 경우 Pearson 적합도 검정이 부정확합니다.
  • Hosmer-Lemeshow: Hosmer-Lemeshow 검정은 다른 적합도 검정처럼 데이터의 행당 시행 횟수에 종속되지 않습니다.데이터의 행당 시행 횟수가 작은 경우 Hosmer-Lemeshow 검정이 모형이 데이터를 얼마나 잘 적합하는지 나타내는 더 신뢰할 수 있는 지표입니다.

반응 정보

변수카운트사건 이름
부패사건506Event
  비사건7482 
용기총계7988 

적합도 검정

검정DF카이-제곱P-값
이탈도50.970.965
Pearson50.970.965
Hosmer-Lemeshow60.101.000
사건/시행 형식에 대한 주요 결과: 반응 정보, 이탈도 검정, Pearson 검정, Hosmer-Lemeshow 검정

이 결과에서 모든 적합도 검정의 p-값은 일반적인 유의 수준인 0.05보다 높습니다. 검정은 예측 확률이 이항 분포에서 예측하지 않는 방식으로 관측된 확률에서 벗어난다는 증거를 제공하지 않습니다.