Bonjour,
J'ai fait de l'imputation multiple sur une variable catégorielle de 4 catégories et analysé les résultats par proc mianalyzej'obtiens un HR et p-value pour chaque catégorie mais j'aimerais obtenir le p-value global, comment faire?
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25 Proc phreg data=pop1_imput1; format slevel slevel.; class slevel (param=ref ref='Primary'); model delriskc1*riskc1(0)=slevel/eventcode=1; by _Imputation_; ods output parameterEstimates = u_slevel; run; proc mianalyze parms(classvar=classval)=a_mvn; class slevel; modeleffects slevel; ods output parameterEstimates = u_slevel1; RUN; run; data u_slevel2; set u_slevel1; HR_est=EXP(ESTIMATE); LCI_HR=HR_est*EXP(-1.96*STDERR); UCI_HR=HR_est*EXP(+1.96*STDERR); RUN; proc print data=u_slevel2; var Parm slevel Probt HR_est LCI_HR UCI_HR; run;





Répondre avec citation
Partager