全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
4934 2
2010-05-07
我在sas上实现广义估计方程时得到结果和书上给的结果不一致。请各位帮忙看看。例子如下:
        在研究某化学物质的致畸试验中,将32只孕鼠分为对照组和试验组,对照组用正常的饲料喂养,
实验组的饲料添加了一定剂量的可疑物。待孕鼠分娩后,观察某畸形仔鼠的发生数。
i是组别变量,j是孕鼠变量,x是畸形仔鼠数,y是孕鼠所生仔鼠数。数据如下:怎么在sas中实现GEE的分析,
请各位给出程序,谢了。
                  i       j      x     y
111313
121212
1399
1488
1588
161213
171112
18910
19910
11089
1111113
11245
11357
114710
115710
11699
2171212
2181111
2191010
22099
2211011
222910
223910
22489
22589
22645
22779
22847
229510
23036
231310
23207
二维码

扫码加我 拉你入群

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

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

全部回复
2020-7-17 15:11:49
我是觉得看您的响应变量要是使用畸形数目的话,可以用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;






二维码

扫码加我 拉你入群

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

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

2021-6-1 15:03:33
subject 和 within是设错了吗
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群