Beispiel für Stabilitätsuntersuchung mit Charge als Zufallsfaktor

Auswählen des Modells

Ein Qualitätstechniker eines Pharmaunternehmens möchte die Haltbarkeit eines Medikaments bestimmen. Die Konzentration des Wirkstoffs im Medikament nimmt mit der Zeit ab. Der Techniker möchte feststellen, zu welchem Zeitpunkt die Konzentration ein Niveau von 90 % der beabsichtigten Konzentration erreicht. Der Techniker wählt nach dem Zufallsprinzip 8 Chargen von Medikamenten aus einer größeren Population möglicher Chargen aus und testet eine Probe aus jeder Charge zu neun verschiedenen Zeitpunkten.

Zum Schätzen der Haltbarkeit führt der Techniker eine Stabilitätsuntersuchung durch. Da es sich bei den Chargen um eine Zufallsstichprobe aus einer größeren Grundgesamtheit möglicher Chargen handelt, ist die Charge ein Zufallsfaktor und kein fester Faktor.

  1. Öffnen Sie die Beispieldaten HaltbarkeitZufallschargen.MWX.
  2. Wählen Sie Statistik > Regression > Stabilitätsuntersuchung > Stabilitätsuntersuchung aus.
  3. Wählen Sie Charge ist ein Zufallsfaktor (gemischtes Modell) aus.
  4. Geben Sie im Feld Antwort die Spalte Wirkstoff% ein.
  5. Geben Sie im Feld Zeit die Spalte Monat ein.
  6. Geben Sie im Feld Charge die Spalte Charge ein.
  7. Geben Sie im Feld Untere Spezifikationsgrenze den Wert 90 ein.
  8. Klicken Sie auf Grafiken.
  9. Wählen Sie unter Residuendiagramme die Option Vier-in-Eins aus.
  10. Klicken Sie in den einzelnen Dialogfeldern auf OK.

Interpretieren der Ergebnisse

Der p-Wert, der die Modelle mit und ohne die Interaktion "Monat nach Batch" vergleicht, beträgt 0,059. Da der p-Wert kleiner als das Signifikanzniveau von 0,25 ist, wird in der Analyse das Modell mit der Interaktion Monat für Charge verwendet. Die Haltbarkeit, die etwa 53 Monate beträgt, ist eine Schätzung der Zeit, in der der Ingenieur zu 95 % sicher sein kann, dass 95 % des Arzneimittels über der unteren Spezifikationsgrenze liegen. Der Schätzwert gilt für alle Chargen, die der Techniker zufällig aus dem Prozess auswählt.

Die marginalen Residuen folgen möglicherweise keiner Normalverteilung mit konstanter Varianz. Die Punkte im Wahrscheinlichkeitsnetz (Normal) folgen der Linie nicht gut. Ein Grund für das nicht normale Verhalten der marginalen Residuen besteht darin, dass die Varianz der marginalen Residuen von der Zeitvariablen abhängt und möglicherweise nicht konstant ist, wenn das endgültige Modell die Interaktion zwischen Charge und Zeit enthält. Sie können die bedingten Residuen verwenden, um die Normalität des Fehlerterms im Modell zu überprüfen.

Faktorinformationen

FaktorTypAnzahl
der
Stufen
Stufen
ChargeZufällig81; 2; 3; 4; 5; 6; 7; 8

Modellauswahl mit α = 0,25

Modell-2 Log-LikelihoodDifferenzp-Wert
Monat Charge Monat*Charge128,599   
Monat Charge133,4244,824760,059
Terme im ausgewählten Modell: Monat; Charge; Monat*Charge

Varianzkomponenten

QuelleVar% von GesamtSE Varz-Wertp-Wert
Charge0,52740972,91%0,3038531,7357390,041
Monat*Charge0,0001740,02%0,0001421,2241020,110
Fehler0,19573927,06%0,0367525,3259320,000
Gesamt0,723322       

Zusammenfassung des Modells

SR-QdR-Qd(kor)
0,44242496,91%96,87%

Koeffizienten

TermKoefSE KoefDFt-Wertp-Wert
Konstante100,0602470,2687067,22372,3783470,000
Monat-0,1387660,0057947,22-23,9501960,000

Prognosen für Zufallseffekte

TermBLUPStdAbwDFt-Wertp-Wert
Charge         
  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
Monat*Charge         
  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

Randanpassungen und Bewertung für ungewöhnliche Beobachtungen

BeobWirkstoff%AnpassungDFResidStd. Resid
10101,56400099,6439507,043681,9200502,375254R
31100,61800098,8113547,052731,8066462,213787R
5598,48100096,7298668,873831,7511342,033482R
R  Großes Residuum

Schätzung der Haltbarkeit

Untere Spez.-Grenze = 90
Haltbarkeit = Zeitraum, während dessen Sie zu 95 % sicher sein können, dass mindestens 95 %
     der Antwortwerte über der unteren Spezifikationsgrenze liegen
Haltbarkeit für alle Chargen = 53,1818

Prüfen der bedingten Residuen

  1. Wählen Sie Statistik > Regression > Stabilitätsuntersuchung > Stabilitätsuntersuchung aus.
  2. Klicken Sie auf Grafiken.
  3. Wählen Sie im Feld Residuen für Diagramme die Option Bedingt regulär aus.
  4. Klicken Sie in den einzelnen Dialogfeldern auf OK.

Interpretieren der Ergebnisse

In diesen Ergebnissen folgen die bedingten Residuen anscheinend einer Normalverteilung. Das vollständige Modell scheint gut an die Daten angepasst zu sein.