楼主: shmilyhbf
1434 0

[问答] 请高手帮忙解决一下一下几个问题,跪求了。。 [推广有奖]

  • 0关注
  • 0粉丝

初中生

38%

还不是VIP/贵宾

-

威望
0
论坛币
1 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
226 点
帖子
6
精华
0
在线时间
14 小时
注册时间
2016-5-29
最后登录
2017-5-9

楼主
shmilyhbf 发表于 2016-6-11 22:14:00 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

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

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币

1本文是利用PROC MIXed程序来实现多个研究之间生存率的合并分析。主要数据操作过程如下:

1、第一步、资料的录入,原始获得的数据见图1(见附件)所示,


我通过自学,将上述资料录入到SAS中。代码如下:

data survival; /* 注:输入Rfmr(f=1-4;m=1-3;r=1,2)(两因素+重复)*/

/*顺序行标f→列标m→重复r*/

do i=1 to 10;/*注:定义试验编号,共10个试验*/

    do j=1 to 2 by 1; /*注:定义治疗组,1代表对照组(即control组),2代表实验组(即exp组)步长为1时by 语句可省略*/

        dok=1 to 7; /*注:定义时间点,步长为1时by 语句可省略*/

                inputsr se @@; /* 注: @@可以使数据连续输入*/

                output;

        end;

     end;

end;

cards;

0.74 0.04 0.580.05 0.52 0.05 0.50 0.05 0.47 0.05 . . 0.43 0.05 0.80 0.04 0.74 0.04 0.66 0.050.62 0.05 0.52 0.05 . . 0.48 0.05

. . 0.5 0.05 .. 0.25 0.04 . . . . 0.16 0.04 . . 0.67 0.05 0.45 0.05 . . . . 0.37 0.05

0.77 0.04 . .0.52 0.05 . . 0.29 0.05 . . . . 0.80 0.04 . . 0.53 0.05 . . 0.40 0.05 . . . .

. . 0.51 0.05. . . . . . . . 0.09 0.03 . . 0.68 0.05 . . . . . . . . 0.31 0.05

. . 0.63 0.05. . 0.45 0.05 . . . . . .  . . . . . . .. . . . . . .

. . 0.57 0.05. . . . . . . . 0.23 0.04 . . 0.78 0.04 . . 0.52 0.05 . . . . 0.42 0.05

. . . . . . .. . . . . . . 0.88 0.04 0.58 0.05 0.43 0.05 0.28 0.05 0.19 0.04 . . . .

. . 0.69 0.05. . . . . . . . 0.23 0.04 . . 0.78 0.04 . . . . . . . . 0.47 0.05

. . 0.68 0.05. . . . . . . . . . . . . . . . . . . . . . 0.42 0.05

. . . . . . .. . . 0.19 0.04 . . . . . . . . . . . . 0.47 0.05 . .

; /*可以看到输入数据时与表格形式一致*/

run;

proc print; /*在结果窗口输出*/

run;

第二步是作图,作图要求如下

1、首先按实验组和对照组分开画图,实验组为图片1,对照组为图片2;

2、然后所有实验组中,同一个试验中不同时点的生存率以时间点为横坐标,生存率为纵坐标,作图。

这步的代码我不懂,请帮忙编程。

第三步:利用上述资料进行meta分析。这个分析我已经通过作者获得SAS程序,可是我看不懂代码,请帮忙看一下代码,并帮忙解释代码的意义。SAS见附件survival_macro.sas。

第四步:由于我没有使用survival_macro.sas进行运算,所以我不知道论文中的Finger2和finger3是不是由survival_macro.sas获得的,如果不是,那么这个图要通过什么代码实现。

另外一组数据录入如下:

data survival; /* 注:输入Rfmr(f=1-4;m=1-3;r=1,2)(两因素+重复)*/

/*顺序行标f→列标m→重复r*/

do i=1 to 14;/*注:定义试验编号*/

    do j=1 to 2 by 1; /*注:定义治疗组,步长为1时by 语句可省略*/

        dok=1 to 5; /*注:定义时间点,步长为1时by 语句可省略*/

                inputsr se @@; /* 注: @@可以使数据连续输入*/

                output;

        end;

     end;

end;

cards;

49 12 46 1242 12 40 12 40 12               54 8 25 823 7 23 7 23 7

55 10 50 1036  9 . . . .                    40 8 23 7 23 7 23 7  .  .

54 10 47 1340 13 40 13 . .                 54 9 42 828 8 28 8 .  .

70 23 70 2370 23 70 23  . .                48 17 48 17 17 13  . .  . .

54 4 46 5 426  . . . .                      40 5 21 4 16 4 16 4 . .

54  2 43 3 40 3 39 3 .  .                    50 4 32 4 24 4 18 4 . .

59 8 49 9 479 47 9 47 9                    . . . . .. . . . .

61 8 53 8 538 53 8 53 8                    . . . . .. . . . .

. . . . . .. . . .                         60 9 48 932 9 32 9 32 9

. . . . . .. . . .                        44 5  26 4 17 5 16 4 .  .

. . . . . .. . . .                        50 3 33 326 3 22 3 19 3

. . . . . .. . . .                        62 3 38 329 3 24 3 22 3

. . . . . .. . . .                        50 10 24 8 16 7 12 6 . .

. . . . . .. . . .                        76 7 53 853 8 50 8 50 8

; /*可以看到输入数据时与表格形式一致*/

run;

proc print; /*在结果窗口输出*/

run;




二维码

扫码加我 拉你入群

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

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

关键词:proc mixed Survival control inputs output survival control 程序 资料

survival_macro.zip
下载链接: https://bbs.pinggu.org/a-2048493.html

1.39 KB

本附件包括:

  • survival_macro.sas

另一组模拟数据的代码.doc

202 KB

Meta-analysis of summary survival curve data.pdf

195.47 KB

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

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-29 06:03