卫生统计学（教学讲义）第五章参数估计和假设检验的SAS实现.docVIP

第五章　参数估计和假设检验的SAS实现 例5－1 随机抽取某地25名正常成年男子，测得146 139 153 138 137 125 142 134 133 122 137 128 140 137 139 128 131 158 138 151 147 144 151 117 118 该样本的均数为，标准差为，求该地正常成年男子总体均数95%可信区间。SAS分析程序 data ch5_1; 定义数据集名 n=25; 变量n为样本量 mean=137.32; 变量mean为样本均数 std=10.63; 变量std为样本标准差 t=tinv(0.975,n-1); 给出自由度为n-1的t分布的0.975的分位点 in=t*std/sqrt(n); lclm=mean-in; lclm为下可信限 uclm=mean+in; uclm为上可信限 run; proc print; var lclm uclm; SAS软件输出结果 Obs lclm uclm 1 132.932 141.708 SAS软件输出结果的解释 SAS计算的正常成年男子总体均数95%可信区间132.932 ~ 141.708 (g/L)，而132.932和141.708分别为下可信限和上可信限。本例自由度(=25-1=24，经查t0.05/2,24=2.064，则： 即该地正常成年男子总体均数的95％可信区间为：~141.7(g/L)。本例～为可信区间，而和分别为下可信限和上可信限。例5－2 某市200年120名7岁男童的身高=123.62(cm)，标准差s=4.75(cm)，计算该市7岁男童总体均数90%的可信区间。SAS分析程序 data ch5_2; 定义数据集名 n=120; 变量n为样本量 mean=123.62; 变量mean为样本均数 std=4.75; 变量std为样本标准差 t=probit(0.95); 给出标准正态0.95的分位点 in=t*std/sqrt(n); lclm=mean-in; lclm为下可信限 uclm=mean+in; uclm为上可信限 run; proc print; var lclm uclm; SAS软件输出结果 Obs lclm uclm 1 122.907 124.333 SAS软件输出结果的解释 SAS计算的该市7岁男童平均身高的90％可信区间为：122.907～124.333 (cm)，而122.907和124.333分别为下可信限和上可信限。 教材中的说明 因n=120100，故可以用标准正态分布t分布，0.10=1.645， 即该市7岁男童平均身高的90％可信区间为：122.91～124.33(cm)，可认为该市7岁男童平均身高在122.91～124.33(cm)之间。 例5－ 为研究铅暴露对儿童智商(IQ)的影响，某研究调查了78名铅暴露(其血铅水平≥40(g/100ml)的6岁儿童，测得其平均IQ为88.02，标准差为12.21；同时选择了78名铅非暴露的6岁儿童作为对照，测得其平均IQ为92.89，标准差为13.34。试估计铅暴露的儿童智商IQ的平均水平与铅非暴露儿童相差多少，并估计两个人群IQ的总体均数之差的95%可信区间。 SAS分析程序 data ch5_3; 定义数据集名 n1=78; 铅暴露儿童的样本量 n2=78; 铅非暴露儿童的样本量 m1=92.89; 铅暴露的儿童平均IQ m2=88.02; 铅非暴露的儿童平均IQ s1=13.34;