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. Er wählt nach dem Zufallsprinzip acht Chargen des Medikaments aus einer umfassenderen Grundgesamtheit möglicher Chargen aus und testet eine Stichprobe 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 umfassenderen Grundgesamtheit möglicher Chargen handelt, ist Charge ein Zufallsfaktor und kein fester Faktor.

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

Interpretieren der Ergebnisse

Der p-Wert für den Vergleich der Modelle mit und ohne Wechselwirkung Monat*Charge beträgt 0,059. Da der p-Wert kleiner als das Signifikanzniveau 0,25 ist, wird in der Analyse das Modell mit der Wechselwirkung Monat*Charge verwendet. Die Haltbarkeit, die ungefähr 53 Monate beträgt, ist ein Schätzwert der Zeit, während der der Techniker sich zu 95 % sicher sein kann, dass 95 % des Medikaments über der unteren Spezifikationsgrenze liegt. Der Schätzwert gilt für alle Chargen, die der Techniker zufällig aus dem Prozess auswählt.

Die Randresiduen folgen u. U. nicht der Normalverteilung mit konstanter Varianz. Die Punkte im Wahrscheinlichkeitsnetz für Normalverteilung folgen der Linie nicht gut. Ein Grund für das nichtnormale Verhalten der Randresiduen ist die Tatsache, dass die Varianz der Randresiduen von der Zeitvariablen abhängt und möglicherweise nicht konstant ist, wenn das endgültige Modell die Wechselwirkung Charge*Zeit enthält. Mit Hilfe der bedingten Residuen können Sie prüfen, ob der Fehlerterm im Modell normalverteilt ist.

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.