변량 배치 요인을 사용한 안정성 연구의 예

모형 선택

한 제약 회사의 품질 엔지니어가 약물의 저장 수명을 확인하려고 합니다. 약물에 포함된 활성 성분의 농도는 시간이 지남에 따라 감소합니다. 엔지니어는 농도가 언제 90%에 도달하는지 확인하려고 합니다. 엔지니어는 가능한 배치의 모집단에서 8개의 약물 배치를 랜덤하게 선택하고 각 배치에서 추출한 표본 하나를 9번 검사합니다.

저장 수명을 추정하기 위해 엔지니어는 안정성 연구를 수행합니다. 배치는 가능한 배치의 모집단에서 추출한 랜덤 표본이기 때문에 고정 요인이 아니라 변량 요인입니다.

  1. 표본 데이터 저장수명랜덤배치.MTW을 엽니다.
  2. 통계분석 > 회귀 분석 > 안정성 연구 > 안정성 연구을 선택합니다.
  3. 반응 변수약의 농도을 입력합니다.
  4. 시간을 입력합니다.
  5. 배치배치을 입력합니다.
  6. 규격 하한90을 입력합니다.
  7. 옵션을 클릭합니다.
  8. 드롭다운 리스트에서 배치가 변량 인자임을 선택합니다.
  9. 확인을 클릭한 다음 그래프을 클릭합니다.
  10. 잔차 그림에서 네 개 모두을 선택합니다.
  11. 각 대화 상자에서 확인을 클릭합니다.

결과 해석

월 x 배치 교호작용이 포함된 모형과 포함되지 않은 모형을 비교하는 p-값은 0.059입니다. p-값이 유의 수준 0.25보다 작기 때문에 분석에서는 월 x 배치 교호작용이 포함된 모형을 사용합니다. 저장 수명(약 53개월)은 약의 95%가 규격 하한 위에 있다고 엔지니어가 95% 확신할 수 있는 기간의 추정치입니다. 추정치는 엔지니어가 공정에서 랜덤하게 선택하는 모든 배치에 적용됩니다.

주변 잔차는 분산이 일정한 정규 분포를 따르지 않을 수도 있습니다. 정규 확률도의 점들이 선을 잘 따르지 않습니다. 주변 잔차의 비정규 동작에 대한 한 가지 이유는 최종 모형에 배치*시간 교호작용이 포함되어 있는 경우 주변 잔차의 분산이 시간 변수에 의존하고 일정하지 않을 수 있기 때문입니다. 조건부 잔차를 사용하여 모형 내 오차 항의 정규성을 확인할 수 있습니다.

요인 정보

요인유형수준 개수수준
배치랜덤81, 2, 3, 4, 5, 6, 7, 8

α = 0.25인 상태에서 모형 선택

모형-2 로그 우도차이P-값
월 배치 월*배치128.599   
월 배치133.4244.824760.059
선택된 모형 안의 항: 월, 배치, 월*배치

분산 성분

출처분산총계의 %SE 분산Z-값P-값
배치0.52740972.91%0.3038531.7357390.041
월*배치0.0001740.02%0.0001421.2241020.110
오차0.19573927.06%0.0367525.3259320.000
총계0.723322       

모형 요약

SR-제곱R-제곱(수정)
0.44242496.91%96.87%

계수

계수SE 계수DFT-값P-값
상수100.0602470.2687067.22372.3783470.000
-0.1387660.0057947.22-23.9501960.000

랜덤 효과 예측

BLUP표준 편차DFT-값P-값
배치         
  11.3594330.31398812.454.3295670.001
  20.3953750.31398812.451.2592030.231
  30.1091510.31398812.450.3476290.734
  4-0.4093220.31398812.45-1.3036230.216
  5-0.1356430.31398812.45-0.4320010.673
  6-1.0647360.31398812.45-3.3910060.005
  70.0494200.31398812.450.1573940.877
  8-0.3036780.31398812.45-0.9671640.352
월*배치         
  10.0062810.00858110.490.7319250.480
  20.0199050.00858110.492.3195370.042
  3-0.0138310.00858110.49-1.6117420.137
  40.0034680.00858110.490.4041730.694
  50.0012400.00858110.490.1444550.888
  60.0002760.00858110.490.0321440.975
  7-0.0109610.00858110.49-1.2772720.229
  8-0.0063780.00858110.49-0.7432200.474

비정상적 관측치에 대한 주변 적합치 및 진단

관측약의 농도적합치DF잔차표준화 잔차
10101.56400099.6439507.043681.9200502.375254R
31100.61800098.8113547.052731.8066462.213787R
5598.48100096.7298668.873831.7511342.033482R
R  큰 잔차

저장 수명 추정치

규격 하한 = 90
저장 수명 = 반응 변수의 95% 이상이 규격 하한 위에 있다고 95% 신뢰할 수 있는 기간
모든 배치에 대한 저장 수명 = 53.1818

조건부 잔차 확인

  1. 통계분석 > 회귀 분석 > 안정성 연구 > 안정성 연구을 선택합니다.
  2. 그래프을 클릭합니다.
  3. 사용할 잔차에서 조건부 정규 잔차을 선택합니다.
  4. 각 대화 상자에서 확인을 클릭합니다.

결과 해석

이 결과에서 조건부 잔차는 정규 분포를 따르는 것으로 보입니다. 전체 모형이 데이터에 적합한 것으로 보입니다.