Contexte Métier
Une usine analyse la durée de vie de composants mécaniques provenant de différents fournisseurs. L'information 'Fournisseur' est parfois manquante dans les logs de maintenance, mais ces pièces doivent être incluses dans l'analyse stratifiée. De plus, le 'Type' de pièce est une variable catégorielle.
Préparation des Données
Données de fiabilité avec une variable 'Fournisseur' contenant des valeurs manquantes et une variable catégorielle 'Type'.
Copié !
| 1 | |
| 2 | DATA casuser.parts_reliability; |
| 3 | call streaminit(789); |
| 4 | DO i=1 to 300; |
| 5 | lifetime = rand('EXPONENTIAL', 500); |
| 6 | fail = 1; |
| 7 | IF rand('UNIFORM') < 0.1 THEN supplier=' '; |
| 8 | ELSE IF rand('UNIFORM') < 0.6 THEN supplier='A'; |
| 9 | ELSE supplier='B'; |
| 10 | IF rand('UNIFORM') > 0.5 THEN type='Rotatif'; |
| 11 | ELSE type='Fixe'; |
| 12 | OUTPUT; |
| 13 | END; |
| 14 | |
| 15 | RUN; |
| 16 |
Étapes de réalisation
1
Modèle de Cox stratifié autorisant les valeurs manquantes dans la strate.
Copié !
| 1 | |
| 2 | PROC CAS; |
| 3 | phreg.cox / TABLE="parts_reliability" class={{vars={"type"}}} model={depVars={{name="lifetime"}}, censor="fail", censVals={0}, effects={{vars={"type"}}}} strata="supplier" strataMissing=TRUE; |
| 4 | |
| 5 | RUN; |
| 6 | |
| 7 | QUIT; |
| 8 |
Résultat Attendu
L'analyse s'exécute sans erreur. Les observations où 'supplier' est manquant sont traitées comme une strate valide (grâce à strataMissing=TRUE). La variable 'type' est correctement gérée comme catégorielle.