《医学统计学》计算题答案

暨南大学医学院医学统计学教研室1医学统计学计算分析题参考答案孙振球主编医学统计学第3版北京人民卫生出版社，2010第二章计量资料的统计描述计算分析题（P26）1根据某单位的体检资料，116名正常成年女子的血清甘油三酯测量结果如下，请据此资料（1）描述集中趋势应选择何指标并计算之。（2）描述离散趋势应选择何指标并计算之。（3）求该地正常成年女子血清甘油三酯的95参考值范围。（4）试估计该地正常成年女子血清甘油三脂在08MMOL/L以下者及15MMOL/L者各占正常成年女子总人数的百分比表21某单位116名正常成年女子的血清甘油三酯（MMOL/L）测量结果组段频数0610730890913101911251218131314915516171合计116（1）数据文件数值变量名组段，频数；用COMPUTE产生新变量“组中值”（也可直接输入组中值）。暨南大学医学院医学统计学教研室2（2）操作步骤ANALYZEDATAWEIGHTCASES；WEIGHTCASESBY频数。ANALYZEDESCRIPTIVESDESCRIPTIVES；将“组中值”选入VARIABLE框中；单击OK。（3）结果解释表21显示某单位116名正常成年女子的血清甘油三酯测量结果呈正态分布，故选择均数描述集中趋势，选择标准差描述离散趋势。某单位116名正常成年女子的血清甘油三酯测量结果的均数为116（MMOL/L），标准差为020（MMOL/L）；该地正常成年女子血清甘油三酯的95参考值范围是（077，155）MMOL/L。计算过程根据公式SX961，即116196020。该单位正常成年女子血清甘油三酯在08MMOL/L以下者估计占总人数的359，15MMOL/L以下者估计占总人数的9554。计算过程当X08时，U08116/02018，查附表1得曲线下面积为00359；当X15时，U15116/02017，查附表1得曲线下面积为10044609554。DESCRIPTIVESTATISTICS116651651156020442116组中值VALIDNLISTWISENMINIMUMMAXIMUMMEANSTDDEVIATION2某地微丝蚴血症患者42例治疗后7年，用间接荧光抗体试验测得抗体滴度如下。求平均抗体滴度。抗体滴度的倒数10204080160例数5121375（1）数据文件数值变量名抗体滴度，例数。（2）操作步骤ANALYZEDATAWEIGHTCASES；WEIGHTCASESBY例数。ANALYZEREPORTSCASESUMMARIZES；将“抗体滴度”选入VARIABLE框中，单击STATISTICS，将GEOMETRICMEAN从STATISTICS框中选入CELLSTATISTICS框。暨南大学医学院医学统计学教研室3（3）结果解释CASESUMMARIESGEOMETRICMEAN3683抗体滴度42例微丝幼血症患者治疗后7年平均抗体滴度为137。3测得某地300名正常人尿汞值，其频数表如下。试计算均数和中位数，何者的代表性较好并求正常人尿汞值95的参考值范围。表22某地300名正常人尿汞值（G/L）频数表组段例数049427858125016452022241628932936440544048352056260064068721合计300（1）数据文件数值变量名组中值，例数。组中值按2，6，1066，70输入。（2）操作步骤ANALYZEDATAWEIGHTCASES；WEIGHTCASESBY频数。表22显示该组正常人的尿汞值为明显的偏态分布，因而选用中位数描述平均水平，用百分位数法计算95参考值范围。由于尿汞值以过高为异常，故求单侧上限95参考值范围。ANALYZEDESCRIPTIVESFRIQUENCIES；将“组中值”选入VARIABLE框中；单击STATISTICS按钮，选择MEDIAN，选择PERCENTITILE计算P95，选择“VALUESAREGROUPEDMIDPOINT”，单击CONTINUE；单击OK。暨南大学医学院医学统计学教研室44结果解释该地正常人尿汞值的中位数为133G/L，95参考值范围为0368（G/L）。STATISTICS组中值30001333A3677BVALIDMISSINGNMEDIAN95PERCENTILESCALCULATEDFROMGROUPEDDATAAPERCENTILESARECALCULATEDFROMGROUPEDDATAB4（补充题）随机抽取某地110名正常成年男性，测得血清胆固醇MMOL/L结果如下。请据此资料（1）绘制频数分布图，判断该资料属于何种分布（2）描述集中趋势应选择何指标并计算之。（3）描述离散趋势应选择何指标并计算之。（4）求该地正常成年男性血清胆固醇的95参考值范围。表23某地110名正常成年男性血清胆固醇MMOL/L测量结果525536382645447575452546549603488429543543379354478447369460452437517523538429577460491481546593478608387426473294419499478507406382520452382608419340580387429445569554658546419439452458413463517452403655473476471450627387447481338437406442442455533426416322507476491497575359463582512369390598463380353364560578302465520530372575（1）数据文件数值变量X表示血清胆固醇。（2）操作步骤ANALYZEDESCRIPTIVESFRIQUENCIES（3）结果解释频数分布图显示某地110名正常成年男性血清胆固醇MMOL/L测量结果呈正态分布，故选择均数描述集中趋势，标准差描述离散趋势。样本均数为473（MMOL/L），标准差为080（MMOL/L）。该地正常成年男性血清胆固醇的95参考值范围是（077，155）MMOL/L。计算过程根据公式SX961，即473196080。暨南大学医学院医学统计学教研室5700600500400300200血清胆固醇MMOL/L20151050FREQUENCYMEAN47292STDDEV080073N110HISTOGRAMSTATISTICSX11004729280073VALIDMISSINGNMEANSTDDEVIATION5（补充题）随机抽取某地109名正常成年女性（1850岁），测得血清铁蛋白NG/ML结果如下。请据此资料（1）绘制频数分布图，判断该资料属于何种分布（2）描述集中趋势应选择何指标并计算之。（3）描述离散趋势应选择何指标并计算之。（4）求该地正常成年女性（1850岁）血清铁蛋白的95参考值范围。表24某地109名正常成年女性（1850岁）血清铁蛋白NG/ML测量结果441201651852302853304204605358805312516518523029033042047054092063126170185240300340420470560955751301721902403003444205006009858013217519024031036043050060010108514018020025031537544050068012009014718021525031538044050072014809015018422025532039044050077515201051531852302703204004405207801600113158185230280330400460530800暨南大学医学院医学统计学教研室6（1）数据文件数值变量X表示血清铁蛋白。（2）操作步骤ANALYZEDESCRIPTIVESFRIQUENCIES（3）结果解释频数分布图显示109名正常成年女性（1850岁）血清铁蛋白NG/ML测量结果呈明显偏态分布，故选择中位数描述集中趋势，四分位数间距描述离散趋势，用百分位数法计算95参考值范围。中位数为473（MMOL/L）；第25百分位数和第75百分位数分别为182NG/ML和470NG/ML，四分位数间距为288NG/ML；该地正常成年女性（1850岁）血清铁蛋白的95参考值范围是（605，1490）NG/ML。STATISTICSX10906050182003100047000149000VALIDMISSINGN25255075975PERCENTILES20001500100050000血清铁蛋白NG/ML302520151050FREQUENCY暨南大学医学院医学统计学教研室7第三章总体均数的估计与假设检验1某地随机抽样调查了部分健康成人的红细胞数和血红蛋白含量，结果如下表表31某年某地健康成年人的红细胞数和血红蛋白含量指标性别例数均数标准差标准值红细胞数1012/L）男360466058484女255418029433血红蛋白G/L男3601345711402女25511761021247请就上表资料（1）说明女性的红细胞数与血红蛋白的变异程度何者为大（2）计算男性两项指标的抽样误差。（3）试估计该地健康成年女性红细胞数的均数。（4）该地健康成年女性血红蛋白含量是否不同（5）该地男性两项血液指标是否均低于上表得标准值（若测定方法相同）参考答案如下（1）女性红细胞数和血红蛋白的变异系数分别为694和867，故血红蛋白比红细胞数变异度大。（2）男性两项指标的抽样误差分别为红细胞数的标准误为00311012/L；血红蛋白的标准误为0374G/L。（3）已知健康成人的红细胞数服从正态分布，而且例数较大为255，故采用正态近似法计算其总体均数的95可信区间。该地健康成年女性红细胞数总体均数的95可信区间为（414,422）1012/L。（4）经两样本T检验，可以推断血红蛋白含量的总体均数男性高于女性（T2423,P020）。4为了解某一新降血压药物的效果，将28名高血压患者随机等分到试验组和对照组，暨南大学医学院医学统计学教研室8试验组采用新降压药，对照组采用标准药物治疗，测得治疗前后舒张压（MMHG）的差值（前后）如下表。问新药和标准药的疗效是否不同表32两种药物治疗前后的舒张压（MMHG）之差新药1210784516321113481414标准药2910502108412345（1）数据文件数值变量GROUP表示分组（1新药；2标准药），BP表示舒张压差值。（2）SPSS操作步骤ANALYZECOMPAREMEANSINDEPENDENTSAMPLESTTEST（3）结果解释INDEPENDENTSAMPLESTEST1177353714260018786236539231364837142398300187862365390313668EQUALVARIANCESASSUMEDEQUALVARIANCESNOTASSUMEDBPFSIGLEVENESTESTFOREQUALITYOFVARIANCESTDFSIG2TAILEDMEANDIFFERENCESTDERRORDIFFERENCELOWERUPPER95CONFIDENCEINTERVALOFTHEDIFFERENCETTESTFOREQUALITYOFMEANS经两样本T检验，可以认为新药疗效优于标准药（T371,P0001），两药均数的差值为88（95CI39136）MMHG。5将钩端螺旋体病人的血清随机分为两组，分别用标准株作凝溶试验，测得稀释倍数如下表，问两组的平均效价是否不同表33钩端螺旋体病患者血清作凝溶试验测得的稀释倍数标准株（11人）1002004004004004008001600160016003200水生株（9人）100100100200200200200400400（1）数据文件数值变量GROUP表示分组（1标准株；2水生株），X表示稀释倍数。由于稀释倍数资料明显服从偏态分布，故将原始数据进行对数变换，然后进行两样本T检验。暨南大学医学院医学统计学教研室9（2）SPSS操作步骤ANALYZECOMPAREMEANSINDEPENDENTSAMPLESTTEST（3）结果解释INDEPENDENTSAMPLESTEST50630373149180065260416704175118769833451558500452604157251919786011EQUALVARIANCESASSUMEDEQUALVARIANCESNOTASSUMEDLGXFSIGLEVENESTESTFOREQUALITYOFVARIANCESTDFSIG2TAILEDMEANDIFFERENCESTDERRORDIFFERENCELOWERUPPER95CONFIDENCEINTERVALOFTHEDIFFERENCETTESTFOREQUALITYOFMEANS对原始数据取对数后，两总体方差仍然不等（F506，P0037），故采用两样本近似T检验。结果显示标准株的平均效价高于水生株（T335,P0004）。6某医师观察某新药治疗肺炎的疗效，将肺炎病人随机分为新药组和旧药组，得两组的退热天数如下表。表34新旧两药的退热天数分组例数平均退热天数退热天数的标准差新药353808旧药375209（1）两药平均退热天数之差的95CI（2）请对两药平均退热天数是否不同进行假设检验。（3）上述两种方法有何联系参考答案如下（1）两药平均退热天数之差总体均数的95可信区间为（10，18）天。（2）经两样本T检验，可以认为新药和旧药的平均退热天数不同（T696,P离子交换法B蒸馏法离子交换法CTESTSTATISTICSB474A635ZASYMPSIG2TAILED蒸馏法离子交换法BASEDONPOSITIVERANKSAWILCOXONSIGNEDRANKSTESTB暨南大学医学院医学统计学教研室292某实验室观察局部温热治疗小鼠移植肿瘤的疗效，以生存日数作为观察指标，试验结果如下表，问局部温热治疗小鼠移植肿瘤是否可延长小鼠生存日数实验组10121515161718202390以上对照组2345678910111213解生存时间一般服从偏态分布，故选用秩和检验。经两个独立样本比较的WILCOXON秩和检验，可认为小鼠生存日数实验组长于对照组，即温热治疗小鼠移植肿瘤可延长小鼠生存日数（Z363，P0000）。RANKS1017001700012692830022分组实验组对照组TOTAL生存日数NMEANRANKSUMOFRANKSTESTSTATISTICSB5000830003630000000AMANNWHITNEYUWILCOXONWZASYMPSIG2TAILEDEXACTSIG21TAILEDSIG生存日数NOTCORRECTEDFORTIESAGROUPINGVARIABLE分组B3下表资料是某药治疗两种不同病情的老年慢性支气管炎病人的疗效，问该药对两种病情的疗效是否不同某药对两种不同病情的老年慢性支气管炎病人的疗效疗效单纯型单纯型合并肺气肿控制6542显效186有效3023无效1311合计12682（1）数据文件数值变量病情（1单纯型；2单纯型合并肺气肿），疗效（3控制；2显效；1有效；0无效），频数。暨南大学医学院医学统计学教研室30（2）SPSS操作步骤ANALYZENONPARAMETRICTESTSINDEPENDENTSAMPLES解疗效分4个等级，故选用秩和检验。经两个独立样本比较的WILCOXON秩和检验，尚不能认为该药对两种病情的疗效不同（Z0543，P0587）。RANKS1261061813378508210192835750208病情单纯型单纯型合并肺气肿TOTAL疗效NMEANRANKSUMOFRANKSTESTSTATISTICSA49545008357500543587MANNWHITNEYUWILCOXONWZASYMPSIG2TAILED疗效GROUPINGVARIABLE病情A4据下表资料，问三种不同人群的血浆总皮质醇测定值有无差别如果有差别，作两两比较。三种不同人群的血浆总皮质醇测定值（102MMOL/L）正常人单纯性肥胖皮质醇增多症011017270052033281061055292069066359077086386086113408102138430108163430127204596192375662解本例为小样本资料，而且数据离散程度较大，如单纯肥胖组的标准差接近均值，故选用秩和检验。经多个样本比较的KRUSKALWALLISH检验，可认为三种不同人群的血浆总皮质醇测定值有差别（21813，P0000）。RANKS1096510117510251030人群正常人单纯性肥胖皮质醇增多症TOTAL皮质醇NMEANRANKTESTSTATISTICSA,B181302000CHISQUAREDFASYMPSIG皮质醇KRUSKALWALLISTESTAGROUPINGVARIABLE人群B3个独立样本两两比较的NEMENYI法检验用SPSS的SYNTAXEDITOR编程进行计算。程序第1行的变量名HC，R1，R2，R3，N1，N2，N3分别表示卡方值（1813），3个样本的平均秩（965，1175，2510）及样本含量（10,10,10），第3行为相应的数据。DATALISTFREE/HCR1R2R3N1N2N3BEGINDATA181396511752510101010ENDDATACOMPUTENN1N2N3暨南大学医学院医学统计学教研室31COMPUTEH12R1N12/N1R2N22/N2R3N32/N3/NN13N1COMPUTECH/HCCOMPUTEX12R1R22/NN1/121/N11/N2CCOMPUTEX13R1R32/NN1/121/N11/N3CCOMPUTEX23R2R32/NN1/121/N21/N3CCOMPUTEP121CDFCHISQX12,2COMPUTEP131CDFCHISQX13,2COMPUTEP231CDFCHISQX23,2FORMATSP12F83FORMATSP13F83FORMATSP23F83EXECUTE程序运行结果（下图）血浆总皮质醇测定值在正常人与单纯性肥胖者之间无差别（2028，P0867），在正常人与皮质醇增多症者之间有差别（21541，P0000），在单纯性肥胖者与皮质醇增多症者之间也有差别（21150，P0003）。附4个独立样本两两比较的NEMENYI法检验（第3行数据需要修改）DATALISTFREE/HCR1R2R3R4N1N2N3N4BEGINDATA1552435029102409222317151711ENDDATACOMPUTENN1N2N3N4COMPUTEH12R1N12/N1R2N22/N2R3N32/N3R4N42/N4/NN13N1COMPUTECH/HCCOMPUTEX12R1R22/NN1/121/N11/N2CCOMPUTEX13R1R32/NN1/121/N11/N3CCOMPUTEX14R1R42/NN1/121/N11/N4CCOMPUTEX23R2R32/NN1/121/N21/N3CCOMPUTEX24R2R42/NN1/121/N21/N4CCOMPUTEX34R3R42/NN1/121/N31/N4CCOMPUTEP121CDFCHISQX12,3COMPUTEP131CDFCHISQX13,3COMPUTEP141CDFCHISQX14,3暨南大学医学院医学统计学教研室32COMPUTEP231CDFCHISQX23,3COMPUTEP241CDFCHISQX24,3COMPUTEP341CDFCHISQX34,3FORMATSP12F83FORMATSP13F83FORMATSP14F83FORMATSP23F83FORMATSP24F83FORMATSP34F83EXECUTE5据下表资料，问三种产妇在产后1个月内的泌乳量有无差别三种产妇在产后1个月内的泌乳量乳量早产足月产过期产无3013210少3629214多3141434合计9783858解经多个样本比较的KRUSKALWALLISH检验，可认为三种产妇在产后1个月内的泌乳量有差别（21700，P0000）。RANKS9739521838505825853983993产妇种类早产足月产过期产TOTAL乳量NMEANRANKTESTSTATISTICSA,B169962000CHISQUAREDFASYMPSIG乳量KRUSKALWALLISTESTAGROUPINGVARIABLE产妇种类B610例食管癌病人在某种药物的保护下，做不同强度的放射线照射，观察血中淋巴细胞畸变百分数，结果如下表。问三者的淋巴细胞畸变百分数有无差别如果有差别，作两两比较。10例食管癌病人放射线照射后血中淋巴细胞畸变百分数病例号照射前照射6000照射90001100000210180120300679741200635102901606101701677105025081060259101009010407070暨南大学医学院医学统计学教研室33解本题为百分数资料，而且数据离散程度很大，3个组的标准差都接近均值，故选用秩和检验。经多个相关样本比较的FRIEDMANM检验，可认为三者的血中淋巴细胞畸变百分数有差别（2779，P0020）。RANKS130240230照射前照射6000照射9000MEANRANKTESTSTATISTICSA1077892020NCHISQUAREDFASYMPSIGFRIEDMANTESTA第九章双变量回归与相关1某地10名一年级女大学生的胸围（CM）与肺活量（L）数据如下表所示。试建立肺活量Y与胸围X的回归方程，并估计胸围为75CM时相应肺活量均数的95可信区间以及个体值的95可信区间。10名一年级女大学生的胸围（CM）与肺活量学生编号12345678910胸围X725839783884771817783748737794肺活量Y251311191328283286316191298328解建立的回归方程为XY0530431。但对回归系数进行T检验，尚不能认为胸围与肺活量之间有直线关系（T165，P0138）。故没有必要做进一步的直线回归分析分析。COEFFICIENTSA143025595595920530325041649138CONSTANT胸围XMODEL1BSTDERRORUNSTANDARDIZEDCOEFFICIENTSBETASTANDARDIZEDCOEFFICIENTSTSIGDEPENDENTVARIABLE肺活量YA2某医师测得10名正常成年男性的血浆清蛋白含量（G/L）及其血红蛋白含量（G/L）数据如下表所示。请对这两项指标作相关分析并计算相关系数的95可信区间。10名正常成年男性的血浆清蛋白含量（G/L）及其血红蛋白含量（G/L）编号12345678910血浆清蛋白含量355365385375365355345342346335血红蛋白含量1195120512751265120511851105109510851053解由原始数据及散点图的观察，可认为血浆清蛋白含量X及血红蛋白含量Y服从双变量正态分布，而且变量间呈直线趋势。经计算，PEARSON相关系数R0960（P0000），总体相关系数的95可信区间为（08352，09908）。可认为血浆清蛋白含量与血红蛋白含量之间有正的直线相关关系。暨南大学医学院医学统计学教研室34130012501200115011001050血红蛋白含量G/L390380370360350340330血浆清蛋白含量G/LCORRELATIONS1963000101096310001010PEARSONCORRELATIONSIG2TAILEDNPEARSONCORRELATIONSIG2TAILEDN血浆清蛋白含量血红蛋白含量血浆清蛋白含量血红蛋白含量CORRELATIONISSIGNIFICANTATTHE001LEVEL2TAILED可信区间的计算方法从计算机桌面选择程序附件计算器（科学型），输入相关系数0960，选择INV（反，逆）和HYP（双曲线函数），单击“TAN”，可得19459；输入19459，单击HYP和TAN可得0960；输入12051，单击HYP和TAN可得08352。3两名放射科医师对13张肺部X线片各自做出评定，评定方法是将X线片按病情严重程度给出等级，结果如下表所示。问他们的等级评定结果是否相关。两名放射科医师对13张肺部X线片的评定结果X线片编号12345678910111213甲医师乙医师两名放射科医师对13张肺部X线片的评定结果（用04分别表示5个等级）X线片编号12345678910111213甲医师2301023434032乙医师1322034344133解本资料为等级数据，故选用秩相关。经秩相关分析，可认为两名医师的评定结果存在正相关关系（RS0808，P0001）。CORRELATIONS1000808001131380810000011313CORRELATIONCOEFFICIENTSIG2TAILEDNCORRELATIONCOEFFICIENTSIG2TAILEDN甲医师乙医师SPEARMANSRHO甲医师乙医师CORRELATIONISSIGNIFICANTATTHE001LEVEL2TAILED暨南大学医学院医学统计学教研室35第十章统计表与统计图1对下列统计表提出意见，并绘制改进后的统计表。发病率职业性别例数例数工人男102329女77791合计10718168农民男8020250女86223合计9110110职员男18223126女163955合计19828141答该表没有标题，行的层次多，不应该有竖线，修改如下表1某年某地不同职业居民的某病发病率（）男女合计职业例数发病数发病率例数发病数发病率例数发病数发病率工人102329802025018223126农民7779186223163955职员107181689110110198281412请评价下表是否符合统计表的制表原则，并作改进。两组临床疗效观察花斑竹组（226例）维生素B12组（70例）治疗后治疗后消失消失项目治疗前例数平均天数改善无效治疗前例数平均天数改善无效食欲不振2532399451182866054900192851疲乏无力2472359511180665651911194341恶心呕吐1441379517174330301000937肝区疼痛171148865173919441276502441104肝脏肿大15914289323091433222688306864黄疸245232957168476595898323191麝絮211182863294122751295693724166GPT高246206837258432866538033038112答该表有竖线、标题不具体、层次多，疗效指标过多、例数显示不当。修改如下暨南大学医学院医学统计学教研室36表1花斑竹和维生素B12治疗某病的临床疗效观察花斑竹组维生素B12组观察指标治疗例数症状消失率（）症状消失平均天数治疗例数症状消失率（）症状消失平均天数食欲不振2539451182609001928疲乏无力2479511180569111943恶心呕吐144951717301000937肝区疼痛1718651739416502441肝脏肿大1598932309326883068黄疸2459571684599832319麝絮2118632941515693724GPT高24683725846680330383某地两年三种死因别死亡率资料如下表，请绘制合适的统计图描述该资料。某地两年三种死因别死亡率（1/10万）死因1952年1992年肺结核1632247心脏病725834恶性肿瘤5721563答该资料是比较两个年度三种死因的变化，适宜用复式直条图，绘制如下恶性肿瘤心脏病肺结核180160140120100806040200死亡率（1/10万）1992年1952年年份图1某地1952年与1992年三种死因的死亡率比较暨南大学医学院医学统计学教研室374某医生统计1992年某市机械工业生产性外伤例数见下表，请用适当的统计图描述。1992年某市机械工业生产性外伤分类的病例数外伤类型病例数百分比（）创伤3814057挫伤3053248眼外伤1181257烧伤92980其他43458合计9391000答该资料是描述生产性外伤的构成比，适宜用圆图，绘制圆图如下458其他98烧伤1257眼外伤3248挫伤4058创伤图21992年某市机械工业939例生产性外伤类型的百分比5某医生在冠心病药物的动物实验中得到以下结果，请用适当的统计图描述。家兔服药后血清胆固醇的变化服药后组别服药前2周4周6周8周10周12周实验组5368158898712989123251179410968对照组50881811081114645164551620114110答该资料是比较两组家兔的血清总胆固醇随时间变化的趋势，适宜用线图。该资料的不足之处是没有显示各组例数、血清胆固醇的单位和数据是否为均数，服药前后的血清胆固醇测量值相差15倍以上提示数据可能有误。因此，该资料不合适作为统计图的习题。暨南大学医学院医学统计学教研室386三种不同人群的血浆总皮质醇测定资料见下表，试选用合适类型的统计图描述以下资料。三种不同人群的血浆总皮质醇测定值（102MMOL/L）正常人单纯性肥胖皮质醇增多症011017270052033281061055292069066359077086386086113408102138430108163430127204596192375662答由于数据的离散程度很大，该资料合适用箱式图。皮质醇增多症单纯性肥胖正常人人群700600500400300200100000血浆总皮质醇测定值302010图6不同人群的血浆总皮质醇测定值（102MMOL/L）7补充习题中国卫生统计年鉴报告了20002010年我国甲乙类报告传染病中淋病和梅毒的发病率，请用适当的统计图描述。暨南大学医学院医学统计学教研室3920002010年我国甲乙类报告传染病淋病和梅毒发病率（1/10万）年淋病梅毒200018645082001148048020021328467200314094502004173471220051379967200612141280200711081588200899019492009902230720107912686答该资料合适用线图描述发病率随时间变化而变化的趋势。20102009200820072006200520042003200220012000年2500200015001000500发病率（1/10万）梅毒淋病图720002010年我国甲乙类报告传染病淋病和梅毒发病率8补充习题为研究药的抑癌作用，使一批小白鼠致癌后，按完全随机设计的方法随机分为四组，A、B、C和对照组，分别接受不同的处理。A、B、C三个实验组分别注射05、10和15ML30的注射液，对照组不用药。经一定时间以后，测定四组小白鼠的肿瘤重量（G），测量结果见下表。请分别用ERRORBAR和BOXPLOT描述该资料。暨南大学医学院医学统计学教研室40四组小白鼠的肿瘤重量（G）实验组对照组ABC36300433452318124224211344114525374036315637133270283206411921145026261345132321（1）ERRORBAR15ML10ML05ML对照组分组65432195CI肿瘤重量（G）图4某种药物对小白鼠抑癌作用实验结果（2）BOXPLOT（注对照组第7个测量值为离群值）15ML10ML05ML对照组分组7060504030201000肿瘤重量（G）7图5某种药物对小白鼠抑癌作用实验结果暨南大学医学院医学统计学教研室41第十一章多因素试验资料的方差分析1研究不同负荷（A因素）和不同航空供养装置（B因素）的呼吸阻力。A因素有三个水平，即在自行车功量计上静坐、250与600KG/MIN体力负荷；B因素有两个水平，即YX1与YX3供养系统。如有12名受试者，每人只能接受一种处理，应如何进行试验设计如果在一段时间内每人可重复接受所有的处理，又应如何设计试写出这两种条件下地随机化分组结果。解A因素有三个水平A1（静坐），A2（250KG/MIN体力负荷），A3（600KG/MIN体力负荷）；B因素有两个水平B1（YX1供氧系统），B2（YX3供氧系统）。G326，共有6个处理组T1（A1，B1），T2（A1，B2），T3（A2，B1），T4（A2，B2），T5（A3，B1），T6（A3，B2）。现有12名受试者，每人只能接受一种处理，用完全随机分组方法将12名受试者等分6组，每组2人。随机分组表如下（随机数字序号R为1，2分配T1，R为3，4分配T2，余类推）。受试者序号123456789101112随机数字221768658168952392358702随机数序号R327698124115101处理分组T2T1T4T3T5T4T6T2T6T3T5T1如果在一段时间内每个受试者可重复接受所有的处理，用区组随机化设计安排处理，即试验时间受试者序号1234561T2T4T3T6T5T12312各个试验时间的处理顺序随机确定（随机数序号R为对应的处理号）。如第一个受试者试验时间顺序123456随机数字193627594613随机数序号R243651处理分组T2T4T3T6T5T1第212个受试者的处理分配重复以上步骤。2下表为完全随机设计35析因试验结果，A因素为毫米波照射频率，B因素为照射时间，表内TK（K1,2，15）为各处理组的实验结果小计。试进行统计分析。暨南大学医学院医学统计学教研室42小鼠肝细胞DNA含量的合计（N5）（绝对单位AU）A因素（I3）B因素（J5）3604GHZ5005GHZ空白对照A因素小计（BI）即刻22035（T1）19380（T6）21820（T11）63235（B1）1天19155（T2）19140（T7）19875（T12）58170（B2）3天19705（T3）16630（T8）18825（T13）55160（B3）5天19120（T4）19810（T9）20615（T14）59545（B4）7天19204（T5）19755（T10）19090（T15）58085（B5）A因素小计（AI）99255（A1）94715（A2）100225（A3）294195（X）平方和X2117074解本题为完全随机设计35析因试验结果。没有给原始数据，不能用SPSS处理。结果略。3研究高锰酸盐处理后对甘蓝叶核黄素浓度测量结果的影响。试验分3天进行，每天安排1次22的析因处理，A因素为试样处理方式，B因素为试样重量，试验结果见下表。试进行析因设计资料的方差分析。甘蓝叶核黄素浓度测量结果（G/G）高锰酸盐处理不处理试验日期025G试样100G试样025G试样100G试样第1天272246395386第2天232242431395第3天248222452330解本题为有4个处理组的22析因设计试验。试验分3天进行，即每个试验重复3次。经析因设计资料的方差分析，可认为高锰酸盐处理（F10842，P0000）及其试样重量（F551，P0047）对甘蓝叶核黄素浓度测量结果有影响；还不能认为高锰酸盐和试样重量的交互作用对甘蓝叶核黄素浓度测量结果有影响（F197，P0198），即高锰酸盐处理对甘蓝叶核黄素浓度测量结果的作用与试样重量无关。TESTSOFBETWEENSUBJECTSEFFECTSDEPENDENTVARIABLE测量结果765529A32551763863400012358501112358501187108300071610817161081084190003640113640155110471302111302119711985284086605131768701281836911SOURCECORRECTEDMODELINTERCEPT高锰酸盐试样重量高锰酸盐试样重量ERRORTOTALCORRECTEDTOTALTYPEIIISUMOFSQUARESDFMEANSQUAREFSIGRSQUARED935ADJUSTEDRSQUARED911A暨南大学医学院医学统计学教研室43补充题将54只患肿瘤小鼠随机分为9组，进行33析因实验，观察微波照射时间和血卜啉（HPD）对小鼠肿瘤组织DNA含量的影响。微波照射时间有三个水平，照射15分钟（A1）、6小时（A2）、24小时（A3）；HPD的应用也有三个水平，即不用HPD（B1）、用HPD后微波照射（B2）、用微波照射后用HPD（B3）。实验结果见下表。试进行析因分析。小鼠肿瘤组织DNA含量（G/100）A1A2A3B1B2B3B1B2B3B1B2B34402022181412182302302901882753041881881451603302461022752463602303602302462187527526520216033013024621814526524634520226514524630410216027560218202145218202175经析因设计资料的方差分析，微波照射时间的主效应（F339，P0043）和血卜啉使用方法的主效应（F711，P0002）均有差别；尚不能认为微波照射时间与血卜啉使用方法之间有交互作用（F171，P0164），即微波照射时间对小鼠肿瘤组织DNA含量的作用与血卜啉使用方法无关；在9组小鼠中，对小鼠肿瘤组织DNA含量的影响最大的组是微波照射时间24小时与微波照射后用血卜啉，小鼠DNA含量的均值为1313（G/100），其它8组小鼠DNA含量的均值最小为1733（G/100），最大为2817（G/100）。TESTSOFBETWEENSUBJECTSEFFECTSDEPENDENTVARIABLEDNA含量113562333A81419529234790032702788167127027881676624030002764311121382155633870435800744422900372271080022791177846977944171016418361250045408027829999630005429717483353SOURCECORRECTEDMODELINTERCEPT微波照射时间血卜啉使用方法微波照射时间血卜啉使用方法ERRORTOTALCORRECTEDTOTALTYPEIIISUMOFSQUARESDFMEANSQUAREFSIGRSQUARED382ADJUSTEDRSQUARED272A暨南大学医学院医学统计学教研室44DESCRIPTIVESTATISTICSDEPENDENTVARIABLEDNA含量2816789759625633342386228831110946255618268718189833415562533381254617333449076205506432318252673957162463341960613117451646210066985818241396876718252005313218177788072718223727488054血卜啉使用方法不用HPD用HPD后微波照射用微波照射后用HPDTOTAL不用HPD用HPD后微波照射用微波照射后用HPDTOTAL不用HPD用HPD后微波照射用微波照射后用HPDTOTAL不用HPD用HPD后微波照射用微波照射后用HPDTOTAL微波照射时间15分钟6小时24小时TOTALMEANSTDDEVIATIONN4例115研究高频呼吸机A，B，C，D，E五个参数对通气量的影响，每个参数有高、低两个水平，试按正交设计安排试验。表1119高频呼吸机五个参数选择的正交设计与试验结果试验序号A1B2C4D8E15通气量升/分11111116262111221938311212236041122128435121122028612121348871221149108122224744921112183210211212485112121139451221222320813221114550142212250301522212552616222216664注表中（1），（2），（4），（8），（15）为L16（215）正交表的列号。暨南大学医学院医学统计学教研室45试对例115的表1119正交试验结果进行统计分析。分析A、B、C、D、E五个参数对通气量的影响，同时分析AB、AC、AE三个1阶交互作用。解首先对该资料的质量进行分析。该资料的研究对象应该是16个健康人或患者，但16次试验中每个试验没有重复。生物个体对相同一个试验的变异一般较大，应考虑重复每个试验。如将48名研究对象完全随机分为16组，每组3名受试者。假定该试验结果受研究对象个体差异的影响较小，可进行正交设计资料的方差分析。先用SPSS的MEAN过程分别计算出A、B、C、D、E五个参数两个水平的均值。可以看出A、B、C、D四个参数都是第2水平的均值高于第1水平，而参数E是第1水平的均值高于第2水平。这些观察到的差别可能由抽样误差或机会造成，需要进行假设检验。XAMEAN299213415500357356A12TOTALXXBMEAN252963461750357356B12TOTALXXCMEAN287213427500357356C12TOTALXXDMEAN334713380000357356D12TOTALXXEMEAN381388333325357356E12TOTALX经正交设计资料的方差分析，按005水准，呼吸机参数A、B、C对通气量的影响均有统计学意义（FA2473，P0002；FB7972，P0000；FC3599，P0001）；参数D、E对通气量的影响无统计学意义（FD375，P0094；FE422，P0079）；AB、AC、AE的交互作用也无统计学意义（FAB434，P0076；FAC003，P0865；FAE002，P0897）；A、B、C三个参数都是第2水平的均值高于第1水平，故应该取参数2。虽然参数D和E对通气量的影响无统计学意义，但样本提供的信息是参数D的2水平高于1水平，参数E的1水平高于2水平，故使通气量达到最大的参数选择为A、B、C、D均取2，E参数取1，这种参数水平组合的通气量为6664（L/MIN）,是16次试验中通气量最高的一种组合。TESTSOFBETWEENSUBJECTSEFFECTSDEPENDENTVARIABLEX3342257A84177821910100020432558