Exemple de la fonction Etude de stabilité avec un facteur de lot fixe

Un ingénieur qualité travaillant pour un laboratoire pharmaceutique souhaite déterminer la durée de stockage maximale de pilules contenant un nouveau médicament. La concentration du médicament dans les pilules diminue au fil du temps. L'ingénieur souhaite déterminer le moment où la concentration des pilules tombe à 90 % de la concentration souhaitée. Comme il s'agit d'un nouveau médicament, l'entreprise ne dispose que de 5 lots pilotes à utiliser pour estimer la durée de stockage. L'ingénieur teste une pilule pour chaque lot, à neuf moments différents.

Pour estimer la durée de stockage, l'ingénieur effectue une étude de stabilité. L'ingénieur échantillonnant tous les lots, le lot est un facteur fixe et non un facteur aléatoire.

  1. Ouvrez les données échantillons, DuréeStockage.MTW.
  2. Sélectionnez Stat > Régression > Etude de stabilité > Etude de stabilité.
  3. Dans la zone Réponse, saisissez %Médicament.
  4. Dans la zone Durée, saisissez Mois.
  5. Dans la zone Lot, saisissez Lot.
  6. Dans la zone Spécification inférieure, saisissez 90.
  7. Cliquez sur Graphiques.
  8. Sous Diagramme de durée de stockage, dans la deuxième liste déroulante, sélectionnez Pas de graphique pour les lots individuels.
  9. Sous Graphiques des valeurs résiduelles, sélectionnez Quatre en un.
  10. Cliquez sur OK dans chaque boîte de dialogue

Interpréter des résultats

Pour se conformer aux recommandations de 2003 de la conférence internationale sur l'harmonisation des exigences techniques pour l'enregistrement des médicaments à usage humain (CIH), l'ingénieur sélectionne une valeur de p de 0,25 pour les termes à inclure dans le modèle. La valeur de p pour l'interaction Mois*Lot est de 0,048. Cette valeur de p étant inférieure au seuil de signification de 0,25, l'ingénieur conclut que les pentes sont différentes pour les équations de régression de chaque lot. Le lot 3 est celui qui présente la plus forte pente, -0,1630, ce qui indique que la concentration décroît plus rapidement dans ce lot. Le lot 2 étant celui qui présente la plus courte durée de stockage, 54,79, la durée de stockage globale est alignée sur celle du lot 2.

Les valeurs résiduelles présentent une normalité satisfaisante et sont dispersées de façon aléatoire autour de 0. Sur le diagramme des valeurs résiduelles en fonction des valeurs ajustées, on observe un plus petit nombre de points sur le côté gauche que sur le côté droit. Cette disposition est due au fait que l'ingénieur qualité a collecté davantage de données au début de l'étude, quand les concentrations étaient élevées. Elle ne contredit pas les hypothèses de l'analyse.

Méthode

Lignes non utilisées5

Informations sur les facteurs

FacteurTypeNombre de
niveaux
Niveaux
LotFixe51; 2; 3; 4; 5

Sélection du modèle avec α = 0,25

SourceDLSomCar séqCM séqValeur FValeur de p
Mois1122,460122,460345,930,000
Lot42,5870,6471,830,150
Mois*Lot43,8500,9622,720,048
Erreur3010,6200,354   
Total39139,516     
Termes dans le modèle sélectionné : Mois; Lot; Mois*Lot

Récapitulatif du modèle

SR carréR carré
(ajust)
R carré
(prév)
0,59498392,39%90,10%85,22%

Coefficients

TermeCoeffCoef ErTValeur de TValeur de pFIV
Constante100,0850,143701,820,000 
Mois-0,136330,00769-17,740,0001,07
Lot         
  1-0,2320,292-0,800,4323,85
  20,0680,2920,230,8183,85
  30,3940,2751,430,1623,41
  4-0,3170,292-1,080,2873,85
  50,0880,2750,320,752*
Mois*Lot         
  10,04540,01642,760,0104,52
  2-0,02410,0164-1,470,1524,52
  3-0,02670,0136-1,960,0603,65
  40,00140,01640,080,9354,52
  50,00400,01360,300,769*

Equation de régression

Lot
1%Médicament=99,853 - 0,0909 Mois
       
2%Médicament=100,153 - 0,1605 Mois
       
3%Médicament=100,479 - 0,1630 Mois
       
4%Médicament=99,769 - 0,1350 Mois
       
5%Médicament=100,173 - 0,1323 Mois

Ajustements et diagnostics pour les observations aberrantes

Observation%MédicamentValeur
ajustée
RésiduelleVal.
résid.
norm.
1198,00199,190-1,189-2,21R 
4392,24292,655-0,413-1,47  X
4494,06993,8230,2460,87  X
R : Valeur résiduelle élevée
X : Valeur de X aberrante

Estimation de la durée de stockage

Limite de spéc. inférieure = 90
Durée de stockage = durée pendant laquelle vous pouvez être sûr à 95 % qu'au moins 50 % de la
     réponse se trouvent au-dessus de la limite de spécification inférieure
LotDurée de
stockage
183,552
254,790
357,492
460,898
566,854
Global54,790