ある製薬会社の品質エンジニアが、新しい薬品を含む錠剤の保存期間を特定したいと考えています。錠剤中の薬品の濃度は時間が経つにつれて減少します。エンジニアは、錠剤がいつ目的の濃度の90%に達するのかを特定したいと考えています。新しい薬品であるため、保存期間の推定に使用できる試験用のバッチは5つのみです。各バッチから1錠ずつの錠剤を9回の異なるタイミングで試験します。
保存期間を推定するため、エンジニアは安定性分析を行います。エンジニアはバッチのすべてをサンプリングしたため、バッチは変量因子ではなく固定因子です。
保存期間.MTW サンプルデータを開きます。
を選択します。
応答 に薬品濃度 を入力します。
時間 に月 を入力します。
バッチ にバッチ を入力します。
下側規格 に、90 と入力します。
グラフ をクリックします。
保存期間プロット の下の2番目のドロップダウンリストから、個別バッチのグラフなし を選択します。
残差プロット で、一覧表示 を選択します。
各ダイアログボックスのOK をクリックします。
結果を解釈する
ICH(医薬品規制調和国際会議;International Conference on Harmonisation of Technical Requirements for Registration of Pharmaceuticals for Human Use)の2003年のガイドラインに従うには、エンジニアは、モデルに含める項にp値0.25を選択します。バッチ交互作用別の「月」のp値は0.048です。p値は有意水準0.25未満のため、エンジニアは、バッチごとの回帰式の傾きは異なりますと結論付けます。バッチ 3 には最も急な勾配 -0.1630 があり、これは、濃度がバッチ 3 で最も速く減少することを示します。バッチ 2 は最短の保存期間である 54.79 を持つため、全体の保存期間は Batch 2 の保存期間です。
残差は正規に十分に近似しており、0付近にランダムに分布しています。残差対適合値プロットでは、右側よりも左側の方が点が少なくなっています。このパターンは、濃度が高い場合に品質エンジニアが研究の早期段階で収集するデータ量を増やしたため発生しました。このパターンは、分析の仮説を損なうものではありません。
安定性分析: 薬品濃度対月、バッチ
方法
未使用の行 5
因子情報
因子 タイプ 水準の数 水準
バッチ 固定 5 1, 2, 3, 4, 5
モデル選択 (α = 0.25)
逐次平
要因 自由度 逐次平方和 均平方 F値 p値
月 1 122.460 122.460 345.93 0.000
バッチ 4 2.587 0.647 1.83 0.150
月*バッチ 4 3.850 0.962 2.72 0.048
誤差 30 10.620 0.354
合計 39 139.516
選択したモデル内の項: 月, バッチ, 月*バッチ
モデル要約
R二乗 (調 R二乗
S R二乗 整済み) (予測)
0.594983 92.39% 90.10% 85.22%
係数
係数の標
項 係数 準誤差 t値 p値 VIF
定数 100.085 0.143 701.82 0.000
月 -0.13633 0.00769 -17.74 0.000 1.07
バッチ
1 -0.232 0.292 -0.80 0.432 3.85
2 0.068 0.292 0.23 0.818 3.85
3 0.394 0.275 1.43 0.162 3.41
4 -0.317 0.292 -1.08 0.287 3.85
5 0.088 0.275 0.32 0.752 *
月*バッチ
1 0.0454 0.0164 2.76 0.010 4.52
2 -0.0241 0.0164 -1.47 0.152 4.52
3 -0.0267 0.0136 -1.96 0.060 3.65
4 0.0014 0.0164 0.08 0.935 4.52
5 0.0040 0.0136 0.30 0.769 *
回帰式
バッチ
1 薬品濃度 = 99.853 - 0.0909 月
2 薬品濃度 = 100.153 - 0.1605 月
3 薬品濃度 = 100.479 - 0.1630 月
4 薬品濃度 = 99.769 - 0.1350 月
5 薬品濃度 = 100.173 - 0.1323 月
異常な観測値の適合値と診断
標準化
観測値 薬品濃度 適合値 残差 残差
11 98.001 99.190 -1.189 -2.21 R
43 92.242 92.655 -0.413 -1.47 X
44 94.069 93.823 0.246 0.87 X
R 大きな残差
X 異常なX
保存期間の推定
下側規格 = 90
保存期間 = 95%の信頼度で50%以上の応答が下側規格を上回る期間
バッチ 保存期間
1 83.552
2 54.790
3 57.492
4 60.898
5 66.854
全体 54.790