SAS中文论坛

标题: 有关1：1配对logistic程序问题 [打印本页]

作者: shiyiming 时间: 2004-10-3 23:36
标题: 有关1：1配对logistic程序问题
各位大哥大姐：
<img src="{SMILIES_PATH}/icon_cry.gif" alt=":cry:" title="Crying or Very sad" /> 要做一个1：1配对的logistic 分析，假如有x1 x2 x3 x4共4个自变量，一个因变量y, 问其sas程序如何？如果 x1 为等级变量，共3个等级，在程序中如何体现其亚变量？如果还想分析x1与x2有无相互作用，在程序中又如何体现？
帮帮忙啊！好人一生平安！

作者: shiyiming 时间: 2004-10-4 22:22
标题: try it.
配对资料可用条件logistic回归模型。

1. SAS程序框架：

proc logistic data=yourdata;
class x1 x2;
model events/trials=x1 x2 x3 x4 x1*x2;
run;
quit;

2. SAS中无需考虑哑变量问题，SAS会自动按照分类变量的水平建立相应的哑变量。

3. 交互作用的分析正如以上程序中显示的那样，将欲分析交互作用的变量用星号连结起来放入模型表达式中即可。

4. 配对资料的logistic回归应当有两个应变量，一个代表事件发生数，即程序中的events，另一个代表试验次数，即程序中的trials。

希望对你有所帮助。

不妥之处也请各位批评指正。

作者: shiyiming 时间: 2004-10-5 11:59
标题: 1：1配对logistic程序问题
<img src="{SMILIES_PATH}/icon_biggrin.gif" alt=":D" title="Very Happy" /> [i:c4aaa][b:c4aaa]首先感谢student给予的指点；[/b:c4aaa][/i:c4aaa]
我有一个1:1配对logistic回归程序，如下
data dataname;
drop id1 x11 x21 x31;
retain id1 x11 x21 x31 0;
input id outcome x1 x2 x3 @@;
if(id=id1) then do;
x1=x11-x1;
x2=x21-x2;
x3=x31-x3;
output;
end;
else do;
id1=id;
x11=x1;
x21=x2;
x31=x3;
end;
cards;
1 1 1 3 0 1 0 1 0 1
2 1 0 3 1 2 0 1 3 0
3 1 0 1 2 3 0 0 2 0
4 1 1 2 0 4 0 1 0 0
5 1 1 1 1 5 0 1 2 1
6 1 0 2 2 6 0 2 0 0
7 1 1 1 1 7 0 0 0 0
8 1 1 1 2 8 0 0 0 0
9 1 3 3 2 9 0 2 2 0
10 1 2 2 2 10 0 0 0 0
;
proc logistic;
model outcome=x1 x2 x3/noint selection=stepwise sle=0.25 sls=0.25;
run;
有些问题仍要求助各位！上面例子只有3个自变量，而我要处理的数据有25个自变量，共120对！因此在cards;语句后输入自变量赋值将是很麻烦的事情！1:1配对的资料又不能象非条件的logistic那样在SAS viewtable窗口方便直接输入，我想问还有什么办法能更好的输入那些自变量赋值吗？
<img src="{SMILIES_PATH}/icon_cry.gif" alt=":cry:" title="Crying or Very sad" /> <img src="{SMILIES_PATH}/icon_cry.gif" alt=":cry:" title="Crying or Very sad" /> <img src="{SMILIES_PATH}/icon_cry.gif" alt=":cry:" title="Crying or Very sad" />
希望得到您的帮助！

作者: shiyiming 时间: 2004-10-5 17:53
标题: infile statement.
如果数据已经在外部文件中，用infile语句或import过程。

欢迎光临 SAS中文论坛 (http://mysas.net/forum/)