Proc logistic et valeurs manquantes [Résolu]

[dachboden]

Bonjour,

Je fais une proc logistic sur une table avec plusieurs variables explicatives dont quelques unes ont des valeurs manquantes.

J'ai lu dans la doc online SAS (version 8) que les enregistrements avec valeurs manquantes sont automatiquement exclus de l'analyse (*) .

Pourtant mes valeurs manquantes apparaissent comme des modalites dans l'output de ma proc. Un beta est calcule pour la dummy variable qui leur correspond.

Est ce normal ? Comment puis je effectivement eviter d'avoir mes valeurs manquantes prises en comptes dans l'analyse (je ne veux pas mettre en input une table ou elles sont deja elimines parce que si j'utilise l'option stepwise ou tout autre procedure iterative d'inclusion/exclusion des variables dans le modele je veux pouvoir garder le maximum d'enregistrements pour les variables qui n'ont pas de valeurs manquantes...)

Merci bp !

(*)

"
Any observation with missing values for the response, offset, or explanatory variables
is excluded from the analysis. The estimated linear predictor and its standard error
estimate, the fitted probabilities and confidence limits, and the regression diagnostic
statistics are not computed for any observation with missing offset or explanatory
variable values. However, if only the response value is missing, the linear predictor,
its standard error, the fitted individual and cumulative probabilities, and confidence
limits for the cumulative probabilities can be computed and output to a data set using
the OUTPUT statement."

[MEGAMIND2]

Je constate que les valeurs manquantes n'ont pas été estimées.

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
DATA test;
input x y;
cards;
10 1
58 1
36 0
.  1
87 0
46 0
.  1
.  0
38 0
;
run;
 
proc logistic DATA=test outest=bb ;
model y=x;
output out=pred p=phat lower=lcl upper=ucl;
run;

[dachboden]

Effectivement,
Vous avez raison.
J'ai confondu dans ma table input la categorie '.' avec la categorie valeur manquante...
Il ne s'agissait donc pas de vraies valeurs manquantes...
Merci pour votre reponse.
Tout s'explique

[MEGAMIND2]

Si t'es pas convaincu, tu peux rajouter une clause where afin d'éviter de modifier ta table

Code :

1
2
3
4
5
proc logistic DATA=test outest=bb ;
model y=x;
output out=pred p=phat lower=lcl upper=ucl;
WHERE x ne .;
run;