SAS基本统计量和ttest检验.docx

《统计计算》第一次上机练习学号：200930980106 姓名：何斌 年级专业：10级统计1班 指导老师：肖莉1、以下数据是变量X和Y的34次观测值：XYXYXYXY180200116100145165115120104100123110141135191205134135151180144160190220141125110130190190153145204235108110190210155160150170158130161145177185121125107115165195177205151135180240154150143160147155127135分别画变量X和Y的等距频率直方图；分别检验变量X和Y的观测数据是否可以认为来自正态总体；计算变量X和Y的描述统计量(样本均值、样本方差、中位数、众数、标准偏度系数和标准峰度系数)。1.解答如下:(1)SAS程序如下所示:data xy;input x y@@;cards;180200116100145165115120104100123110141135191205134135151180144160190220141125110130190190153145204235108110190210155160150170158130161145177185121125107115165195177205151135180240154150143160147155127135;run;procprint;run;procgchartdata=xy;vbar x y/type=percent;run;得到结果如下所示: 图1变量X的等距频率直方图图2 变量Y的等距频率直方图(2)SAS程序如下所示:procunivariatedata=xy normal;var x y;run;得到结果如下所示:图3 对变量X的正态性检验结果由图3的检验结果可知,P值0.05,故接受原假设,即认为变量X服从正态分布.图4 对变量Y的正态性检验结果由图4的检验结果可知,P值0.05,故接受原假设,即认为变量Y服从正态分布.(3) SAS程序如下所示:procunivariatedata=xy ;var x y;outputout=XY_Tongjiliang mean=mean_X mean_Y var=Var_X Var_Y median=median_X median_Y mode=mode_X mode_Y skewness=S_X S_Y kurtosis=K_X K_Y ;run;procprintdata=XY_Tongjiliang;run;得到结果如下所示:图5 变量X、Y的各统计量指标值2、已知一数据集SCORE中的数据如下：nosexmathphyeng9901F9893919902M8682859903M8794659904F8990789905F9698829906M928789产生一新数据集NEWSCORE，将该数据集存入d:\example下，且转换为文本文件。计算MATH和ENG的均值，方差，偏度，峰度，并将这些结果存入数据集RESULT。产生一新变量，其值为MATH、PHY、ENG的平均值。2.解答如下：(1)SAS程序如下所示：data score;input no sex$ math phy eng;cards;9901F9893919902M8682859903M8794659904F8990789905F9698829906M928789;run;procprint;run;data newscore;set score;filed:\example\newscore.txt;put no sex$ math phy eng;run;得到结果如下所示：图6 将数据存入记事本（2）SAS程序如下所示：procmeansdata=score;var math eng;outputout=result mean=math_mean eng_mean var=math_var eng_varskew=math_skew eng_skew kurt=math_kurt eng_kurt;run;procprintdata=result;run;得到结果如下所示：图7 math、eng的各统计量的值（3）SAS程序如下所示：data aa(drop=no sex);set score;average=(math+eng+phy)/3;run;procprint;run;得到结果如下所示：图8 各科成绩及其平均值3、某克山病区测得11名克山病病人与13名健康人的血磷值(mmol/L)如下表，问该地区急性克山病病人与健康人的血磷值平均水平是否不同？血磷值(