如何写1：2 配对卡方sas程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

如题，如何写1：2的配对卡方分析sas程序，有无简单的代码？
模拟数据如下（不知道什么数据格式比较好，先给个整理过的吧，非原始数据）

a 为病例（1表示暴露），b为对照（0为无暴露，1表示两个对照有一个暴露，2表示两个暴露），x为对子数；

data n;
input a b x;
cards;
0 0 5
0 1 6
0 2 8
1 0 6
1 1 8
1 2 9
;
run;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：sas程序 如何写 cards Input 模拟数据 程序 如何

按照CMH的思路.京剧

http://blog.sina.com.cn/s/blog_a3a926360101kqtc.html

jingju11 发表于 2014-6-25 05:18
按照CMH的思路.京剧

http://blog.sina.com.cn/s/blog_a3a926360101kqtc.html

京大神，我没太看懂你的数据库建库程序，我想要的是1：2配对的，你的结果感觉像1：1的，另外CMH可以做1：2的么？

可能写错了。你自己的数据自己应该更了解。也或许那个总结的数据表根本就是不对的。CMH可能根本不能做1：2的情形。那么什么样的匹配可以呢？也就是说，你的想法是什么呢？
京剧

jingju11 发表于 2014-6-25 08:50
可能写错了。你自己的数据自己应该更了解。也或许那个总结的数据表根本就是不对的。CMH可能根本不能做1：2的 ...

京大神， 这是我朋友的一个数据，他做的是1：2病例对照研究，请我给他做卡方分析。病例是高血压，对照是正常，暴露因素是喝酒300g；现在就是想了解病例和对照的 暴露率差异是否不同。 模拟原始数据如下：
data n;
input case control1 contro2;
cards;
0 0 1
0 0 1
0 1 1
1 1 1
0 0 0
0 1 0
1 1 1
1 0 1
1 1 0
0 0 0
0 0 0
0 0 1
0 0 1
1 1 1
1 0 1
1 0 1
;
run;
*1表示暴露;

if data set has a format above, use the code like

1. input Case1 ctrl1 ctrl2;
   
2. id ++1;
   
3. case =1; exp =1; n =0+Case1; output;
   
4. case =1; exp =0; n =1-Case1; output;
   
5. case =0; exp =1; n =0+ctrl1; output;
   
6. case =0; exp =0; n =1-ctrl1; output;
   
7. case =0; exp =1; n =0+ctrl2; output;
   
8. case =0; exp =0; n =1-ctrl2; output;

复制代码

to get the input data.

The way I am doing should be consistent with the method form Breslow & Day.

JingJu

jingju11 发表于 2014-6-25 11:38
if data set has a format above, use the code like

the results are basically
OR from MLE = 0.5679. from CMH = 0.5714. The chisquare without correction is t =0.45, with correction of 0.5 = 0.125 with df =1. that is, p-value > 0.5.
JingJu

Emails flood in. 他们都询问我结果是怎么来的。其实，我是用公式手算出来的。 现在特别忙。有空在 sas里把公式写出来。OR_MLE 需要解出高次方程。
京剧

http://blog.sina.com.cn/s/blog_a3a926360101kqtc.html

jingju11 发表于 2014-6-26 06:38
Emails flood in. 他们都询问我结果是怎么来的。其实，我是用公式手算出来的。 现在特别忙。有空在 sas里把 ...

I have updated the calculaiton bases on Breslow & Day. The formual is little bit tricky. I do some comments on it.

JingJu

http://blog.sina.com.cn/s/blog_a3a926360101kqtc.html

配对卡方可以使用网页版SPSSAU进行操作与分析。