g******7 发帖数: 1433 | 1 假设我有500个row,有的variable是categorical的(male,female),想在其中做一
个proportion test,看看他们的proportion是不是significant different,我的sas
code该怎么写,是先把proportion算出来在run test吗,有没有比较简洁的code?
谢谢大家! |
k*******a 发帖数: 772 | 2 Proc freq
Chi square test |
g******7 发帖数: 1433 | 3 嗯,如果是binary的 Z test就可以了吧,再请教具体code怎么写。。,我搜了下网上
的,test 2个proportion的code没找到,倒有什么glm的在里面
【在 k*******a 的大作中提到】![](/moin_static193/solenoid/img/up.png) : Proc freq : Chi square test
|
k*******a 发帖数: 772 | 4 假设 你又一个categorical variable X 和 一个 categorical Gender
你要比较 X的每个level里面 Male 的proportioan是否 一样
proc freq;
table X*gender /chisq;
run; |
g******7 发帖数: 1433 | 5 不好意思我可能没表达清楚,等于我是有500个obs(全部都是病例),我想测试Male
gender是不是一个risk factor导致这个病, 我想的是把500个里面male的个数除以总
数500,跟female的test 看不是看significant different,求出P value,请问这样
SAS codes 该怎么写?
【在 k*******a 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 假设 你又一个categorical variable X 和 一个 categorical Gender : 你要比较 X的每个level里面 Male 的proportioan是否 一样 : proc freq; : table X*gender /chisq; : run;
|
k*******a 发帖数: 772 | 6 那你就是要测试 male的 p=0.5
proc ttest就可以了
【在 g******7 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 不好意思我可能没表达清楚,等于我是有500个obs(全部都是病例),我想测试Male : gender是不是一个risk factor导致这个病, 我想的是把500个里面male的个数除以总 : 数500,跟female的test 看不是看significant different,求出P value,请问这样 : SAS codes 该怎么写?
|
g******7 发帖数: 1433 | 7 嗯,是的,那这样SAS CODE怎么写比较简洁些,不需要自己算proportion吧
这个code应该不难,但我找了好多地方都没说。。
【在 k*******a 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 那你就是要测试 male的 p=0.5 : proc ttest就可以了
|
k*******a 发帖数: 772 | 8 proc freq;
table sex / testp=(0.5, 0.5);
run;
【在 g******7 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 嗯,是的,那这样SAS CODE怎么写比较简洁些,不需要自己算proportion吧 : 这个code应该不难,但我找了好多地方都没说。。
|
g******7 发帖数: 1433 | 9 谢谢! 用了你的以后run出来了,我在run之前的把sex中unknown和空白的obs drop掉
,只留male和female的,但这样出来写的是Chi-Square Test
for Specified Proportions,那这还是一个chisquare test? 那默认的是5% 的
significant level么?
【在 k*******a 的大作中提到】![](/moin_static193/solenoid/img/up.png) : proc freq; : table sex / testp=(0.5, 0.5); : run;
|
l****u 发帖数: 529 | 10 我怎么绝的你的设计有点问题,你这样求出的proportion是conditional proportion,
即p(x=gender/y=1)。
risk factor should compare p(y=1/x=gender)with p(y=0/x=gender)?
【在 g******7 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 不好意思我可能没表达清楚,等于我是有500个obs(全部都是病例),我想测试Male : gender是不是一个risk factor导致这个病, 我想的是把500个里面male的个数除以总 : 数500,跟female的test 看不是看significant different,求出P value,请问这样 : SAS codes 该怎么写?
|
|
|
g******7 发帖数: 1433 | 11 嗯,conditional proportion?我也不是很确定,但上面的code应该是male的
proportion跟0.5比,如果在significant level不同的话就reject h0,证明male 是一
个risk factor,我认为这样有道理。。
【在 l****u 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 我怎么绝的你的设计有点问题,你这样求出的proportion是conditional proportion, : 即p(x=gender/y=1)。 : risk factor should compare p(y=1/x=gender)with p(y=0/x=gender)?
|
l****u 发帖数: 529 | 12 如果你的数据全是病人,只能说明在这批病人数据中male和female比例均衡。不能说明
gender是这个病的risk factor
【在 g******7 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 嗯,conditional proportion?我也不是很确定,但上面的code应该是male的 : proportion跟0.5比,如果在significant level不同的话就reject h0,证明male 是一 : 个risk factor,我认为这样有道理。。
|
g******7 发帖数: 1433 | 13 有道理。。那该怎么test male gender 是个risk factor呢,还有类似于ethnicity这
种variable之类的?
【在 l****u 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 如果你的数据全是病人,只能说明在这批病人数据中male和female比例均衡。不能说明 : gender是这个病的risk factor
|
l****u 发帖数: 529 | 14 logistic regression is a very good model to define risk factor. of course
you can compare conditional proportion as p(y=1/gender) and p(y=0/gender).
In your case, it looks you only have the data of patients. If I were you, I
would change my objective to comparing the proportion of each variable
within the patients.
Thanks |
k*******a 发帖数: 772 | 15 你这个有点类似Case-Control study
你现在有Case, 你需要Control
然后比较Case里面Malde比例和Control里面Male的比例
【在 g******7 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 有道理。。那该怎么test male gender 是个risk factor呢,还有类似于ethnicity这 : 种variable之类的?
|