资源描述：
《随机抽样方法的SAS实现.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、中国卫生统计2007年2月第24卷第1期85随机抽样方法的SAS实现上海第二军医大学卫生统计学教研室(200433)王睿贺佳在实际调查过程中,我们往往采用抽样调查的方法factorsi=100of1000;outputout=b;run;获得有代表性的样本。目前流行病学调查中常用的随procprintdata=b;run;1机抽样方法有单纯随机抽样、系统抽样、分层抽样等。方法2:在实际应用中,面对大量的调查数据,如何使用统计软procplanseed=25070416;件快速有效地进行抽样,这是实际中经常碰到的问题,factorsa=100blo
2、ck=10;outputout=c;run;现介绍这些抽样方法如何在SAS82中实现。datac;setc;ifblock=1thensam=A;elsesam=B;no=n单纯随机抽样;run;单纯随机抽样是最简单、最基本的抽样方法。从datac;setc;ifsam=A;run;总体N个对象中,利用抽签或其他随机方法抽取nprocprintdata=c;varno;run;个,构成一个样本。一般用于数目不大的情况。在3UNIFORM过程SAS中有很多命令可用于单纯随机抽样。UNIFORM是一个产生均匀随机数的函数,可通11SURVEYSELE
3、CT过程过其产生的随机数进行抽样。基本命令语句如下:该过程提供了很多方法,用于基于概率的随机样dataa;set;name=uniform本的选择。基本命令语句如下:(seed);run;procsurveyselectdata=dataa;seta;ifname=<;run;method=srsn=procprintdata=a;run;out=seed=number;run;name是任意指定的变量名。samplerate即抽样31这
4、里使用的方法SRS即simplerandomsampling。N比例。SEED即种子数,为2-1的正整数。31是指抽样样本量。SEED指定一个2-1的正整举例:在1000个数据中随机抽取100个,即按数,用该值作为随机数发生器的初始数值,以便在下次01的比例抽样。运行该程序时得到的结果不变。dataa;doi=1to1000;output;end;run;/*产生1举例:在1000个数据中随机抽取100个。-1000共1000个数字*/dataa;doi=1to1000;output;end;run;/*产生1dataa;a=uniform(11111);run
5、;/*产生1000个-1000共1000个数字*/在0-1之间均匀分布的随机数*/procsurveyselectdata=amethod=srsn=100outdataa;seta;ifa=<01;run;=bseed=25070419;run;procprintdata=a;vari;run;procprintdata=b;run;因为UNIFORM语句产生的是均匀分布的随机数字,32PLAN过程所以使用(a=<01)不一定会正好得到100个数据,PLAN过程一般用于产生随机化设计方案,也可而是100左右个数据。以将其提供的随机数用以抽样。基本命令语句如
6、下:4RANUNI过程procplanseed=number;factorsname=mofn;RANUNI也是一个产生均匀随机数的函数,用法outputout=;run;与UNIFORM基本类似。基本命令语句如下:name是指这个设计中因子的名字,mofn规定从n个dataa;set;name=ranuni数值中选择m个,并且0;run;方法1:procprintdata=
7、a;run;procplanseed=25070419;举例:在1000个数据中随机抽取100个,即按01的比例抽样。通讯作者:贺佳(下转第93页)ChineseJournalofHealthStatistics,Feb2007,Vol.24,No.193实力。对该科加强管理,做好思想工作,调配人力,提高业务人员的思第二类科室妇科、按摩科、针灸科指标完成情况属于比较想素质,来促进该科完成好各项指标。好的科室。特别是中药使用率929%,是全院完成得最好的科第四类科室是儿科,儿科是我院完成各项业务指标最好的室,从全院来看这项指标完成得比较好难度
