楼主: kongqunyu
2010 0

SAS help 中 catmod 例28.8程序中疑问 [推广有奖]

  • 10关注
  • 2粉丝

副教授

37%

还不是VIP/贵宾

-

威望
0
论坛币
2977 个
通用积分
0.4469
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
14735 点
帖子
183
精华
0
在线时间
1274 小时
注册时间
2008-7-16
最后登录
2019-1-7

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
疑问在于以下程序中,类似矩阵的东西的意思是什么? 如何构建? 我有一份资料和这个和类似,只是时间上的重复是两次 ,题目中这个是三次,想套用程序,该如何改动?
请高手帮忙~~  我从事卫生统计,要学习的东西很多~~

SAS help 中 catmod process 例子28.8 中:

title 'Growth Curve Analysis';
   data growth2;
      input Diagnosis $ Treatment $ week1 $ week2 $ week4 $ count @@;
      datalines;
   mild std n n n 16    severe std n n n  2
   mild std n n a 13    severe std n n a  2
   
   mild new n n n 31    severe new n n n  7
   mild new n n a  0    severe new n n a  2
   ;
\\Diagnosis (mild or severe)
Treatment (std or new).
The times used in the design matrix (0, 1, 2) correspond to the logarithms (base 2) of the actual times (1, 2, 4). \\

proc catmod data=growth2 order=data;
      title2 'Reduced Logistic Model';
      weight count;
      population Diagnosis Treatment;
      response logit;
      model week1*week2*week4=(1 0 0 0,  /* mild, std */
                                                     1 0 1 0,
                                                     1 0 2 0,
   
                                                     1 0 0 0,  /* mild, new */
                                                     1 0 0 1,
                                                     1 0 0 2,
   
                                                     0 1 0 0,  /* severe, std */
                                                     0 1 1 0,
                                                     0 1 2 0,
   
                                                    0 1 0 0,  /* severe, new */
                                                    0 1 0 1,
                                                    0 1 0 2)
             (1='Mild diagnosis, week 1',
              2='Severe diagnosis, week 1',
              3='Time effect for std trt',
              4='Time effect for new trt')
              / freq design;
      contrast 'Diagnosis effect, week 1' all_parms 1 -1 0 0;
      contrast 'Equal time effects' all_parms 0 0 1 -1;
   quit;

二维码

扫码加我 拉你入群

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

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

关键词:CATMOD Help ATM CAT elp process title count 程序 如何

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-6-17 09:49