统计计算课程设计.doc

统计计算 课程设计报告学院 专业 姓名 学号 评语：分数 题型一： 1、某医生随机抽取正常人和脑病病人各11例，测定尿中类固醇排出量（mg/dl），结果如表1。该医生根据此资料算得正常人尿中类固醇排出量的均数=4.266mg/dl，标准差S1=0.985mg/dl；脑病病人尿中类固醇排出量的均数=5.732mg/dl，标准差S2=1.626mg/dl，配对t检验结果，t = 3.098，P 0.05，所以我们可以认为正常人和脑病病人的尿中类固醇排出量方差齐性。正态性检验的结果如下：p=1p=2因为是小样本，所以观察Shapiro-Wilk检验结果：说明在0.05水平下，正常人的尿中类固醇排出量的正态性统计量W的P值为：0.22590.05，脑病病人的尿中类固醇排出量的正态性统计量W的P值为：0.54040.05，所以正常人和脑病病人的尿中类固醇排出量都通过正态性检验，即正常人和脑病病人的尿中类固醇排出量的分布都服从正态分布。所以选择采用独立组两样本T检验，输出结果如下：结果表明：正常人和脑病病人的尿中类固醇排出量数据是方差齐性的（Equal），所以观察相应Pooled方法的检验统计量的P值=0.01880.05，结论为正常人和脑病病人的尿中类固醇排出量有显著差异。2、某研究者为研究核黄素缺乏对尿中氨基氮的影响，将60只Wistar大白鼠随机分为核黄素缺乏、限食量、不限食量三组不同饲料组。每组20只大白鼠。一周后测尿中氨基氮的三天排出量，结果如表2。该研究者对上述资料采用了两样本均数t检验进行两两比较，得出结论：三组之间均数差异均有统计学意义（P0.05）。检验进行两两比较，得出结论：三组之间均数差异均有统计学意义（P0.05，所以我们可以认为三组不同饲料组的大白鼠进食一周后尿中氨基氮的三天排出量（mg）的方差是齐性的。正态性检验的结果如下：p=1p=2p=3因为是小样本，所以观察Shapiro-Wilk检验结果：说明在0.05水平下，核黄素缺乏饲料组的大白鼠进食一周后尿中氨基氮的三天排出量（mg）的正态性统计量W的P值为：0.25260.05；限食量饲料组的大白鼠进食一周后尿中氨基氮的三天排出量（mg）的正态性统计量W的P值为：0.00010.05；不限食量饲料组的大白鼠进食一周后尿中氨基氮的三天排出量（mg）的正态性统计量W的P值为：0.02560.05。所以除核黄素缺乏饲料组外，限食量饲料组和不限食量饲料组的大白鼠进食一周后尿中氨基氮的三天排出量（mg）都不通过正态性检验，即核黄素缺乏饲料组的大白鼠进食一周后尿中氨基氮的三天排出量（mg）是服从正态分布的，而限食量饲料组和不限食量饲料组的大白鼠进食一周后尿中氨基氮的三天排出量（mg）的分布不是正态分布。所以选择采用非参数方差分析，输出结果如下：三组不同饲料组的大白鼠进食一周后尿中氨基氮的三天排出量（mg）的均值分别为：29.500，17.950和44.050，可以初步看出三组不同饲料组的大白鼠进食一周后尿中氨基氮的三天排出量（mg）是存在差异的。用Bonferroni法进行多重比较，输出结果说明：三组不同饲料组的大白鼠进食一周后尿中氨基氮的三天排出量（mg）之间有显著差异。3、某医师用改良的Seldingers插管技术对8例经临床和病理证实的恶性滋养细胞肿瘤进行选择性盆腔动脉插管灌注化疗。测定治疗前后血中的HCG含量如表3。该医师考虑到数据相差较大，采用对数变换后进行两样本均数比较的t检验，得t =2.460，P0.05，差异有统计学意义，故认为治疗前后血中HCG的含量有差别。 表3 灌注化疗前后HCG含量测定结果（pmol/L）病例号 灌注前（X1）灌注后（X2）lgX1 lgX2 1 1280000 210000 6.1072 5.3222 2 75500 3300 4.8779 3.5185 3 12450 2210 4.0952 3.3444 4 1500000 9.3 6.1761 0.9685 5 10000 2500 4.0000 3.3979 6 9700 1203 3.9868 3.0803 7 15588 4825 4.1928 3.6835 8 4223 914 3.6256 2.9609 【问题3】 1、这是什么资料? 2、该实验属于何种设计方案？ 3、该医师统计方法是否正确？为什么？请用SAS程序和SAS结果解释原因。答：1、这是成对组（配对）数据的资料。2、该实验属于配对比较设计方案。3、该医师统计方法不正确；因为两独立组样本均数t检验要求每组观察值是来自正态分布的总体的，而经检验证明，上述的两组数据并不都服从正态分布，所以不可以采用两独立组样本均数t检验，而应考虑采用非参数检验；而且上述数据的内部有明显地进行了“以前”和“以后”情况的对比，所以不应该采用简单的两样本均数比较的t检验，而应考虑采用成对组（配对）检验，所以综上所述：上述数据的统计分析应考虑采用成对组Wilcoxon符号秩检验。SAS源程序：data tk13;input x y;diff=x-y;cards;6.10725.32224.87793.51854.09523.34446.17610.968543.39793.98683.08034.19283.68353.62562.9609;proc print;run;proc univariate data=tk13 normal;var diff;run;SAS结果及分析：（1）因为两组数据相差较大，所以首先对数据进行对数变换，来保证两组数据的方差齐性，然后再对对数变换后的数据进行相关的检验分析。原数据的方差齐性检验：输出结果中，统计量的p值=0.0000110150.05，所以我们可以认为为对数化后的灌注化疗前后HCG含量测定结果（pmol/L）的方差是齐性的。（2）录入数据：变量解释：x为灌注前HCG含量测定结果的对数化结果（pmol/L）；y为灌注后HCG含量测定结果的对数化结果（pmol/L）；diff为灌注前后HCG含量测定结果的对数化结果的差值。因为是小样本，所以观察Shapiro-Wilk检验结果：正态性检验的结果来看，P值0.0001 0.05，所以在0.05显著性水平下拒绝灌注化疗前后HCG含量测定结果的差值服从正态分布，故可采用非参数的符号检验和符号秩和检验。而从符号检验统计量M=4, P=0.0078c1 then do; w1=0; c1=t+fs+w1; end; else do; w1=c1-t; c1=t+fs+w1; n=n+1; end; ss=ranexp(seed)*0.9; if c1c2 then do; w2=0; c2=c1+w2+ss; end; if c1c1 then do; w2=c2-c1; c2=c1+w2+ss; n=n+1; end; if c1c2 & tc1否是否是否计算该顾客在两服务台的总逗留时间：wz=c2-t计算该顾客离开第二个服务台的时间：c2=w2+t+ss计算在第二个服务台前的等待时间：w2=0队列中的等待人数：n=n+1计算在第一个服务台前的等待时间：w1=0c1c2计算并输出服务台服务系统的主要指标，如顾客在两服务台的总逗留时间，排队中的顾客人数；并计算顾客在两服务台的平均总逗留时间和排队中的顾客平均数模拟时钟t600（或1000）计算在第一个服务台前的等待时间：w1=0是计算该顾客离开第一个服务台的时间：c1=w1+t+fs队列中的等待人数：n=n+1计算在第一个服务台前的等待时间：w1=c1-t否tc1推进模拟时钟：置t=t+ct服务台服务模拟：按相应统计规律产生顾客到服