SAS基本统计量及ttest检验

1、统计计算第一次上机练习学号：200930980106 姓名：何斌 年级专业：10级统计1班 指导老师：肖莉1、以下数据是变量X和Y的34次观测值：XYXYXYXY180200116100145165115120104100123110141135191205134135151180144160190220141125110130190190153145204235108110190210155160150170158130161145177185121125107115165195177205151135180240154150143160147155127135（1）分别画变量X和Y的等距频率直方图；（2）分别检验变量X和Y的观测数据是否可以认为来自正态总体；（3） 计算变量X和Y的描述统计量（样本均值、样本方差、中位数、众数、标准偏度 系数和标准峰度系数）。1解答如下：（1）SAS程序如下所示：data xy;input x y;cards;180 200 116 100 145 165 115 120104 100 123 110 141 135 191 205134 135 1

2、51 180 144 160 190 220141 125 110 130 190 190 153 145204 235 108 110 190 210 155 160150 170 158 130 161 145 177 185121 125 107 115 165 195 177 205151 135 180 240 154 150 143 160147 155 127 135Irun;proc print;run;proc gchart data=xy;vbar x y/type=percent;run;得到结果如下所示：图1变量X的等距频率直方图y MIDPOINT图2变量Y的等距频率直方图SAS程序如下所示：proc univariate data=xy normal; var x y;run;得到结果如下所示：正盍性检验检验yctfSF 值一-Shapiro-ffilkw0.960334Fi- IXJ.1500Cr =ijri eivon MisesW-Sq.052071Pr 啪-丸X).2500iLmler z 口口 一 l ：rl i ngA._Sq0.394098Ft

3、 A-Sq.为.500图3对变量X的正态性检验结果由图3的检验结果可知,P值005,故接受原假设，即认为变量X服从正态分布.Shap i i- 口-Ki lkwo. assessFr D0. 10E5Cr：iiTi er_Y0iL Mi iala.084862Fr ff-Sq0. 1786Ail de r e ori-ll:rl 1 np;A_S qa.504705Pr kSq . 19B5检验正态性检验轨计呈F值图4对变量Y的正态性检验结果由图4的检验结果可知,P值005,故接受原假设，即认为变量Y服从正态分布.SAS程序如下所示：proc univariate data=xy ;var x y;output out=XY_Tongjiliangmean=mean_X mean_Y var=Var_X Var_Y median=median_X median_Ymode=mode_X mode_Y skewness=S_X S_Y kurtosis=K_X K_Y ;run;proc print data=XY_Tongjiliang;run;得到结果如下所示：gyXrr.aorYE

4、_XKJfmaLrr.iiorYmodoji5*15LLO08150.2350.U7E.3610.45396771.59B151B. LIU -0.86635-0.74139150.5152.5190135图5变量X、Y的各统计量指标值2、已知一数据集SCORE中的数据如下：nosexmathphyeng9901F9893919902M8682859903M8794659904F8990789905F9698829906M928789(1) 产生一新数据集NEWSCORE，将该数据集存入d:example下,且转换为文本文件。(2) 计算MATH和ENG的均值，方差，偏度，峰度，并将这些结果存入数据集RESULT。(3) 产生一新变量，其值为MATH、PHY、ENG的平均值。2.解答如下：SAS程序如下所示： data score;input no sex$ math phy eng; cards;9901 F 98 93 919902M8682859903M8794659904F8990789905F9698829906M928789run;proc print;run;data n

5、ewscore; set score;file d:examplenewscore.txt; put no sex$ math phy eng;run;得到结果如下所示：newscoreLtxt文件编辑格式回琶(99 01F98939199 02M86S28599 03M87946599 04Fa?907899 05F9698829? 06M928789图6将数据存入记事本(2)SAS程序如下所示：proc means data=score;var math eng;output out=result mean=math_mean eng_mean var=math_var eng_var skew=math_skew eng_skew kurt=math_kurt eng_kurt;run;proc print data=result;run;得到结果如下所示：m at hm ；ELt hmatkTie：jiLerLg_Tiy!rL(Ti 江 th_varpng_V!irskewyrLg_zkewkurtyrig_kir t91. 333381. &66T23.866788.6B&T

6、0. 39223-1.23406-1.775461.5S04图7 math、eng的各统计量的值(3)SAS程序如下所示： data aa(drop=no sex); set score;average=(math+eng+phy)/3;run;proc print;run;得到结果如下所示：UtisTiathphyeiLE辽晳已 race93SI94.0000S2654.33339465S2.0000907885.66&T986292.0000ST89S9.3333图8各科成绩及其平均值3、某克山病区测得11名克山病病人与13名健康人的血磷值(mmol/L)如下表，问该地 区急性克山病病人与健康人的血磷值平均水平是否不同？血磷值(mmol/L)病人0.841.051.21.21.391.531.671.801.872.072.11健康人0.540.640.640.750.750.811.161.201.341.351.481.561.873.解答如下：(1)对数据进行正态性检验的SAS程序如下所示: data a;do type=healthy;do n=1 to 11; input

7、 xlzhi; output;end;end;cards;0.841.051.21.21.391.531.671.801.872.072.11run;data b;do type=sick; do n=1 to 13; input xlzhi; output;end;end; cards;0.540.640.640.750.750.811.161.201.341.351.481.561.87!run;data c;set a b;run;proc print;run;proc univariate data=c normal;class type;var xlzhi;run;对两种不同状况的人的数据进行正态性检验，得到结果如下所示:正态性检验检验轨计量-P值Shapir D-Wi lkW0.959147Pr D0.1500Cr；inier-vorL Wise5W-Sq 0. 027442Tr W-3q0.2500ArLiler e un-D :rt lmg1-Sq 0. 194403Pr A-5 ij0.2500图9健康人数据的正态性检验结果正态性检验统计量F值Sh4piro-Wlk0.926209Fr D0. 1471Ci-:aniei-vonW-Sq0.072543Pr0. 2422ArLilr z i：iTl_D :i_i-l l Rg0.429717PrXJ. 2500图10病人数据的正态性检验结果由图9、图10的检验结果可知,P值均大于005,故接受原假设，即认为两种情况下 的数据均服从正态分布.（2）两独立样本的均值的T检验，SAS程序如下所示：proc ttest data

《SAS基本统计量及ttest检验》由会员博****1分享，可在线阅读，更多相关《SAS基本统计量及ttest检验》请在金锄头文库上搜索。