为什么加了交互作用项，得不到OR数据集

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

1. data Neuralgia;
   
2.    input Treatment $ Sex $ Age Duration Pain $ @@;
   
3.    datalines;
   
4. P  F  68   1  No   B  M  74  16  No  P  F  67  30  No
   
5. P  M  66  26  Yes  B  F  67  28  No  B  F  77  16  No
   
6. A  F  71  12  No   B  F  72  50  No  B  F  76   9  Yes
   
7. A  M  71  17  Yes  A  F  63  27  No  A  F  69  18  Yes
   
8. B  F  66  12  No   A  M  62  42  No  P  F  64   1  Yes
   
9. A  F  64  17  No   P  M  74   4  No  A  F  72  25  No
   
10. P  M  70   1  Yes  B  M  66  19  No  B  M  59  29  No
    
11. A  F  64  30  No   A  M  70  28  No  A  M  69   1  No
    
12. B  F  78   1  No   P  M  83   1  Yes B  F  69  42  No
    
13. B  M  75  30  Yes  P  M  77  29  Yes P  F  79  20  Yes
    
14. A  M  70  12  No   A  F  69  12  No  B  F  65  14  No
    
15. B  M  70   1  No   B  M  67  23  No  A  M  76  25  Yes
    
16. P  M  78  12  Yes  B  M  77   1  Yes B  F  69  24  No
    
17. P  M  66   4  Yes  P  F  65  29  No  P  M  60  26  Yes
    
18. A  M  78  15  Yes  B  M  75  21  Yes A  F  67  11  No
    
19. P  F  72  27  No   P  F  70  13  Yes A  M  75   6  Yes
    
20. B  F  65   7  No   P  F  68  27  Yes P  M  68  11  Yes
    
21. P  M  67  17  Yes  B  M  70  22  No  A  M  65  15  No
    
22. P  F  67   1  Yes  A  M  67  10  No  P  F  72  11  Yes
    
23. A  F  74   1  No   B  M  80  21  Yes A  F  69   3  No
    
24. ;
    
25. run;
    
26. 
    
27. ods output parameterEstimates=estimate OddsRatios=OR;
    
28. proc logistic data=Neuralgia;
    
29.         class Treatment Sex;
    
30.         model Pain= Treatment Sex treatment*sex;
    
31. run;
    
32. ods output close;

复制代码

为什么加了交互作用项，得不到OR数据集呢？谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：交互作用 数据集 得不到 treatment Parameter

用Enterprise Guide分析，得到了OR呢。不知对不对？

LIBNAME ECLIB000 "C:\";

/* -------------------------------------------------------------------
   由“SAS 任务”生成的代码

   生成时间: 2014年1月8日星期三 09:59:28
   任务: Logistic 回归

   输入数据: C:\www.sas7bdat
   服务器: Local
   ------------------------------------------------------------------- */
ODS GRAPHICS ON;

%_eg_conditional_dropds(WORK.SORTTempTableSorted);
/* -------------------------------------------------------------------
   对数据集 C:\www.sas7bdat 排序
   ------------------------------------------------------------------- */

PROC SQL;
        CREATE VIEW WORK.SORTTempTableSorted AS
                SELECT T.Pain, T.Sex, T.Treatment
        FROM ECLIB000.WWW as T
;
QUIT;
TITLE;
TITLE1 "Logistic 回归结果";
FOOTNOTE;
FOOTNOTE1 "由 SAS 系统 (&_SASSERVERNAME, &SYSSCPL) 于 %TRIM(%QSYSFUNC(DATE(), NLDATE20.))%TRIM(%SYSFUNC(TIME(), NLTIMAP16.)) 生成";
PROC LOGISTIC DATA=WORK.SORTTempTableSorted
                PLOTS(ONLY)=ALL
        ;
        CLASS Sex         (PARAM=EFFECT) Treatment         (PARAM=REF);
        MODEL Pain (Event = 'No')=Sex*Treatment                /
                SELECTION=NONE
                LINK=LOGIT
        ;
RUN;
QUIT;

/* -------------------------------------------------------------------
   任务代码的结尾。
   ------------------------------------------------------------------- */
RUN; QUIT;
%_eg_conditional_dropds(WORK.SORTTempTableSorted);
TITLE; FOOTNOTE;
ODS GRAPHICS OFF;

as cross effect was involved in the model, OR for one effect depends on other effects level, and versa. You can use ODDSRATIO statement to output OR for a specific effect. For example, you want OR for sex

proc ...; model ...; odsRatio sex; ods output OddsRatioWald =ORs; run;

However, you need to notice that under this situation, the OR needs to be interpreted cautiously.
JingJu

"contrast" will give what you want.

楼上两位说的很对，因为有了交互作用，计算OR值的时候就必须要考虑到交互变量是在什么水平，比如在算A vs P的OR值在sex = F 的水平上，SAS code如下：

1. ods output parameterEstimates=estimate  ContrastEstimate=OR;
   
2. proc logistic data=Neuralgia;
   
3.         class Treatment Sex;
   
4.         model Pain= Treatment sex treatment*sex;
   
5.                 contrast 'A vs P at gender = F' treatment 2  1  treatment*sex -2 -1 /estimate=exp;
   
6. run;
   
7. ods output close;

复制代码

heperwong 发表于 2014-1-8 13:41
楼上两位说的很对，因为有了交互作用，计算OR值的时候就必须要考虑到交互变量是在什么水平，比如在算A vs P ...

请问最后那行的treatment 2  1  treatment*sex -2 -1是什么意思呢，如果我想做改动应该怎么改呢