Блог

Когда я использую proc logistic в SAS, на выходе он возвращает доверительное значение интервала и p-значение отношения шансов, как я могу вывести стандартную ошибку отношения шансов?

 proc logistic data=edu;                                                                                  
model school = age sex income/ clodds=wald orpvalue;
oddsratio age;
run; 
  

Вывод нравится

                  Odds Ratio Estimates and Wald confidence interval

                                      Point          95% Wald
                    Effect         Estimate      Confidence Limits  p-value

                    age              1.21         0.74      2.001   < 0.01
  

• Совет: На странице документации Proc Logistic Details -> ODS Table Names перечислены все таблицы, которые процедура создаст для ODS. ODDSRATIO ... /CL=WALD ...; Оператор создает выходную таблицу с именем OddsRatiosWald .
  • В ODS TRACE ON инструкции также будут записаны имена таблиц, которые создает шаг процесса для вывода ODS.
• Сохраните таблицу как набор выходных данных, используя ODS OUTPUT инструкцию.

Пример:

Код из образцов SAS изменен для сохранения ODS OUTPUT .

 * Example 76.4 Nominal Response Data: Generalized Logits Model;

data school;
   length Program $ 9;
   input School Program $ Style $ Count @@;
   datalines;
1 regular   self 10  1 regular   team 17  1 regular   class 26
1 afternoon self  5  1 afternoon team 12  1 afternoon class 50
2 regular   self 21  2 regular   team 17  2 regular   class 26
2 afternoon self 16  2 afternoon team 12  2 afternoon class 36
3 regular   self 15  3 regular   team 15  3 regular   class 16
3 afternoon self 12  3 afternoon team 12  3 afternoon class 20
;

ods trace on;
ods graphics on;

ods html file='logistic.html';

proc logistic data=school;
   freq Count;
   class School Program(ref=first);
   model Style(order=data)=School Program School*Program / link=glogit;
   oddsratio program / cl=wald;

   ods output OddsRatiosWald=or_program;
run;

proc print data=or_program;
  title "Logistic Odds Ratios CL=Wald output data";
run;

ods html close;
ods trace off;
title;
  

Выходные данные, проверенные viewtable в базовом SAS