在sas上实现广义估计方程时遇到的问题。

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我在sas上实现广义估计方程时得到结果和书上给的结果不一致。请各位帮忙看看。例子如下：
        在研究某化学物质的致畸试验中，将32只孕鼠分为对照组和试验组，对照组用正常的饲料喂养，
实验组的饲料添加了一定剂量的可疑物。待孕鼠分娩后，观察某畸形仔鼠的发生数。
i是组别变量，j是孕鼠变量，x是畸形仔鼠数，y是孕鼠所生仔鼠数。数据如下：怎么在sas中实现GEE的分析，
请各位给出程序，谢了。

i	j	x	y
1	1	13	13
1	2	12	12
1	3	9	9
1	4	8	8
1	5	8	8
1	6	12	13
1	7	11	12
1	8	9	10
1	9	9	10
1	10	8	9
1	11	11	13
1	12	4	5
1	13	5	7
1	14	7	10
1	15	7	10
1	16	9	9
2	17	12	12
2	18	11	11
2	19	10	10
2	20	9	9
2	21	10	11
2	22	9	10
2	23	9	10
2	24	8	9
2	25	8	9
2	26	4	5
2	27	7	9
2	28	4	7
2	29	5	10
2	30	3	6
2	31	3	10
2	32	0	7

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：广义估计方程 请各位帮忙 化学物质 对照组 GEE 程序 化学 饲料

我是觉得看您的响应变量要是使用畸形数目的话，可以用poisson分布的gee,  要是畸形率的话，可以采用线性gee，也请大家积极讨论
大体上程序（gee/genmod 似乎都可以，还请大神指点）
proc gee data=aa;
model y=x j; /*j为个体变量，根据需求设置*/
repeated subject=i/within=j /type=ind/exch /ar(1)/un  corrw;/*根据QIC指标确定作业相关矩阵 ，越小越好
run;

subject 和 within是设错了吗