样本量计算SAS程序大全.docx

估计总体均数时样本含量的计算1.1??SAS程序引用公式参照??颜红主编.医学统计学8年.人民卫生出版社,2561.2??例题，据文献报道，一个镇静药物中某有效成分含量为2.25，标准差为0.85，临床医师估计其有效成分的含量的均值在真值的正负0.10范围之内，问需要调查多大的样本？1.3??程序%let alpha=0.05; %let sigma=0.85; %let delta=0.1; data d1_1;n=(probit ((1-alpha/2))*sigma/delta)**2;n1=int(n)+1; proc print data=d1_1;run;2、估计总体率时样本含量的计算2.1??SAS程序引用公式参照??颜红主编.医学统计学8年.人民卫生出版社,2532.2??例题：对北京地区肠易激综合征的患病率进行调查，根据以往的文献资料，人群患病率约为15%，若将允许误差控制在3%，则样本量至少应为多少人？2.3??程序％let P＝0.15；％let delta＝0.03；％let alpha＝0.05；data d1_2;n=(probit((1-alpha/2))/delta)**2*p*(1-p);n1=int(n+1);proc print data=d1_2;run;3、成组设计定量资料样本含量计算（双侧）3.1??SAS程序引用公式参照??方积乾主编.医学统计学第6版.人民卫生出版社,2653.2??例题：在饮食中降低盐能否较低血压值的研究中，将受试者分为两个组别（低盐饮食组和高盐饮食组），预试验结果为两组血压值的标准差分别为12mmHg和10.3mmHg，假设两组样本比例为1:1，检验水准为a＝0.05（双侧），检验效能为1-B＝0.90，检测两组血压差为4mmHg，需要多大样本？3.3??程序%let alpha=0.05; %let beta=0.10; %let delta=12.33; %let std=25;%let q1=0.5;%let q2=0.5;data d1_3;n=(1/q1+1/q2)*((abs(probit((alpha/2)))+probit((1-beta)))*std/delta)**2;N=int(n)+1;n1=int(N*q1)+1;n2=int(N*q2)+1;proc print data=d1_3;run;说明：1）根据颜虹老师书中所述，无论单双侧t（B）均取单侧；2）此段程序中采用的Probit意为：标准正态分布左侧p分位数。probit（a）还是probit（B）均为负值，因此用probit(1-beta)以及abs(probit(alpha/2))，来分别表达t（B）和t（a）。3）此公式和程序同样适用于单组、配对或者交叉设计定量资料分析时的样本含量计算，std（标准差）为每对观察对象差值的标准差，得出的n为对子数即可。4）程序中q1和q2为样本比例，假如为1:1，那么q1＝q2=0.5。4、成组设计定量资料样本含量计算（单侧）4.1??SAS程序引用公式参照??方积乾主编.医学统计学第6版.人民卫生出版社,2654.2??例题：比较A、B两种降压药物降低收缩压的疗效，收缩压下降值相差5mmHg以上具有临床意义，假设两组药物标准差＝12mmHg，若使B组的样本量为A组的两倍，检验水准为a＝0.05，检验功效1－B＝0.90，问A、B各需要多少例？4.3??例题%let alpha=0.05; %let beta=0.10; %let delta=5; %let std=12;%let q1=0.333;%let q2=0.667;data d1_4;n=(1/q1+1/q2)*((abs(probit(alpha))+probit((1-beta)))*std/delta)**2;N=int(n)+1;n1=int(N*q1)+1;n2=int(N*q2)+1;proc print data=d1_4;run;说明：单侧检验和双侧检验的区别就在于用alpha替换alpha/2。5、配对二分类资料分析时的样本含量计算5.1??SAS程序引用公式参照??方积乾主编.医学统计学第6版.人民卫生出版社,2665.2??例题：拟比较甲、乙两种检测方法对腹泻婴幼儿乳糖不耐受的检出情况，初步估计甲法的阳性检出率为48％，乙法的阳性检出率为30%，两种方法一致的阳性检出率为25%。若检验水准a＝0.05（双侧），检验效能1-B＝0.10，问需要多少两本量？5.3??程序%let p1=0.48;%let p2=0.30;%let alpha=0.05;%let beta=0.10;%let pk=0.25;data d1_5;%let pc=(p1+p2-