全部版块 我的主页
论坛 提问 悬赏 求职 新闻 读书 功能一区 悬赏大厅 文献求助专区
892 1
2016-01-30
悬赏 10 个论坛币 未解决
不知道为什么我老是得不出我想要的模拟程序 悬赏跪求大神帮助
/*我所要建立瓮模型模型规则简介:*/
/*本翁模型规定,在试验开始时,翁中放两类球,红球1个,黑球1个。*/
/*模拟开始:*/
/*(1)由试验组和对照组中发生病情进展的先后顺序而生成的组别数据集sorrt*/
/*(2)刚开始在翁中随机抽取一个球观察颜色后放回翁中,*/
/*     如果第一例发生进展的病人来自对照组,且第一次摸出的是红球,则该病人转向试验组治疗,并且给翁中增加10个红球;*/
/*           如果第一例发生进展的病人来自对照组,且第一次摸出的是黑球,则该病人继续在对照组中治疗,同样给翁中增加10个红球;*/
/*     如果第一例发生进展的病人来自试验组,且第一次摸出的是红球,则该病人继续在试验组治疗,并且给翁中增加1个黑球;*/
/*           如果第一例发生进展的病人来自试验组,且第一次摸出的是黑球,则该病人继续在试验组治疗,同样给翁中增加1个黑球;*/
/*(3)由上一步可知翁中球的个数和比例发生了变化,接着进行第二次抽球,抽出球后观察其颜色后放回,此时*/
/*     如果第二例发生进展的病人来自对照组,且第二次摸出的是红球,则该病人转向试验组治疗,并且给翁中增加10个红球;*/
/*           如果第二例发生进展的病人来自对照组,且第二次摸出的是黑球,则该病人继续在对照组中治疗,同样给翁中增加10个红球;*/
/*     如果第二例发生进展的病人来自试验组,且第二次摸出的是红球,则该病人继续在试验组治疗,并且给翁中增加1个黑球;*/
/*           如果第二例发生进展的病人来自试验组,且第二次摸出的是黑球,则该病人继续在试验组治疗,同样给翁中增加1个黑球;*/
/*(4)依次类推*/

data sortt;/*建立一组在治疗过程中发生进展的病人的分组情况的数据集*/
input num group$@@;
cards;
1        对照组
2        试验组
3        对照组
4        对照组
5        试验组
6        对照组
7        对照组
8        对照组
9        对照组
10        对照组
11        试验组
12        对照组
13        对照组
14        试验组
15        对照组
16        对照组
17        对照组
18        对照组
19        对照组
20        对照组
;
run;
/*options nonotes nosource;*/
%macro ballurn (nred,nblack);/*建立一个瓮模型,每次仅仅随机抽一个球,nred为翁中红球个数,nblack为翁中黑球个数*/
data red;
%do x1=1 %to &nred;
RorB='红球';
output;
%end;
data black;
%do x1=1 %to &nblack;
RorB='黑球';
output;
%end;
run;
data conjoin;
set red black;
run;
proc surveyselect  data=conjoin out=dong  method=SRS n=1 seed=2052014476;/*用proc surveyselect随机抽取一个球*/
/*id   RorB;*/
run;
/*说明:method用来指定随机抽样方法的,其中SRS是指不放回简单随机抽样(Simple Random Samping);*/
/*urs是指放回简单随机抽样(Unrestricted Random Sampling);sys是指系统抽样(Systematic Sampling)。*/
/*seed用来指定
二维码

扫码加我 拉你入群

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

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

全部回复
2016-1-30 00:48:17
不好意思  这个帖子发错地方了 应该放到SAS版块的
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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