青椒 做1：2配对的条件Logistic回归的sas程序？

shiyiming

做1：2配对病例对照研究的条件Logistic回归的sas程序。有目标变量lung_cancer(=1患病即病例组有225例，=0正常即对照组450例)，还有自变量x1 x2 x3 x4 x5等，均为分类变量
欲研究x1与lung_cancer（即肺癌）的关系，且x1有3个等级，同时调整x2、x3、x3、x4等混杂因素变量（各有不同等级，有的有序，有的无序）。

请问如何实现呢？
data a;
input lung_cancer x1 x2 x3 x4 x5@@;
cards;
0  1   2  3  5  4
1  1   3  2  1  2
0  2   3  2  5  3
1  2   4  1  3  2
0  3   2  1  4  1
1  3   1  1  1  2
……/*共675条记录*/
；
proc phreg data=a;

model lung_cancer=x1;
run;

以上程序是做非条件logistic的，做条件的话中间少了一些内容，请教高手最精华的内容。
还有且不论事非条件，还是条件logistic回归，如何调整其余的混杂变量呢？
或者还需要补充什么资料？
是直接分析好了，还是做个r*2频数表分析。

shiyiming

小弟不才，你的这个是条件logistic回归？看上去怎么是生存分析的子程序？抛砖引玉吧：
1.首先，可以用proc logistic子程序。
2.由于X1是有序变量，所以直接用其值。但是对于无序变量，有两种方法解决：一种是定义class语句；一种是生成哑变量，注意哑变量个数应该是无序变量水平数-1，感觉你是学统计的，应该知道个中原因。

shiyiming

phreg的确可以做条件logistic回归。只是程序的写法有点问题吧。

shiyiming

to jingju11
我知道有问题，就是不会写啊，呵呵

shiyiming

data a;
input lung_cancer x1 x2 x3 x4 x5;
time=2-lung_cancer;*define time;
cards;
0 1 2 3 5 4
1 1 3 2 1 2
0 2 3 2 5 3
1 2 4 1 3 2
0 3 2 1 4 1
1 3 1 1 1 2
……/*共675条记录*/
；
proc phreg data=a;
model time*lung_cancer(0)=x1 x1 x2 x3 x4 x5/ties=discrete;*adjusting for x2 ...;
run;

shiyiming

sorry, only one x1 in last model

shiyiming

"且x1有3个等级，同时调整x2、x3、x3、x4等混杂因素变量（各有不同等级，有的有序，有的无序）。.."

Sorry, I did not pay attention to those conditions.
you know, phreg only can deal with contiunous variables; that is, if you do not regard those as continuous, you have to change those variables to dummies.for example, for x, you can write two new variables ,x1=(x=1) and x2=(x=2) and so on;
i believe that proc tphreg can deal with class directly, you can do some research on that.