简单统计分析与过程.ppt

课本例6.7 data package; input sale1 sale2@@; dif=sale2-sale1; cards; 66 72 70 75 75 68 79 87 65 84 90 73 85 70 82 83 97 95 95 90 92 82 73 78 71 69 69 74 77 86 ; proc univariate normal; var dif; run; 第31页，共43页，2022年，5月20日，15点22分，星期五 关于简单统计分析与过程 第1页，共43页，2022年，5月20日，15点22分，星期五 一、假设检验与SAS过程 平均每户消费支出 500 600 700 800 900 1000 家庭数 8 15 30 25 13 9 例题6.1： 为了了解农村居民家庭消费水平是否有所提高，2008年，某市对其农村居民家庭进行了一次抽样调查，其中100户被抽样家庭的调查结果如下表： 表6.1 2008年某市农村居民家庭月均消费水平 若3年前该市农村居民家庭月均消费支出服从N（720，17580），假定2008年月均消费支出服从正态分布，问该市农村居民家庭月均消费支出是否有显著提高？（显著性水平0.05） 即在方差未知的情况下检验 第2页，共43页，2022年，5月20日，15点22分，星期五 统计量的计算值、临界值、显著性水平及检验概率之间的关系 假定 检验统计量Z服从正态分布 统计量的计算值：一次抽样观测值代入统计量Z后得到的数值Z0. 临界值：在给定的显著性水平下，由 第3页，共43页，2022年，5月20日，15点22分，星期五 统计量的计算值、临界值、显著性水平及检验概率之间的关系 检验概率： 由临界值和检验概率的计算公式，可知 因此，判断接受或拒绝H0只需看p大于还是小于 第4页，共43页，2022年，5月20日，15点22分，星期五 单样本和两样本下的假设检验 单样本的假设检验 （一）单样本的参数假设检验（正态分布总体） 总体均值的假设检验 检验统计量： 拒绝域： 第5页，共43页，2022年，5月20日，15点22分，星期五 检验统计量： 拒绝域： 第6页，共43页，2022年，5月20日，15点22分，星期五 对于总体均值的假设检验，可转化为均值是否为零的检验，可通过PROC MEANS过程实现，只需在选项中选择t,prt,和clm，alpha。 例6.1程序： data consume; input expend number @@; dif=expend-720; cards; 500 8 600 15 750 30 800 25 900 13 1000 9 ; proc means mean t prt; var dif; freq number; output out=meant t=tv; run; 第7页，共43页，2022年，5月20日，15点22分，星期五 由于检验变量dif=expend-720的t值=3.17， 概率pr|t|的值为0.0020，小于显著性水平0.05，故 在0.05的显著性水平下推断出dif的均值显著不为0，也即居民月均消费支出显著不等于720. 第8页，共43页，2022年，5月20日，15点22分，星期五 data a; set meant; k=_freq_-1; p=1-probt(tv,k); t1=tinv(0.95,k); proc print;run; p=1-probt(t,k) t1=tinv(0.95,k);计算t分布的0.95分位数 显然，tv的值t1且p值也0.05,故在0.05的显著性水平下拒绝原假设，也即接受居民月均消费支出显著大于720. 第9页，共43页，2022年，5月20日，15点22分，星期五 总体方差的假设检验 检验统计量： 拒绝域： 第10页，共43页，2022年，5月20日，15点22分，星期五 例6.2：检验例6.1中居民消费支出的方差是否有变化，即是否仍为17580。 proc means var; var expend; freq number; output out=test var=varex; run; data A(drop=_type_); set test; k=_freq_-1; chisq=k*varex/17580; p=1-probchi(chisq,k); ci1=cinv(0.025,k); ci2=cinv(0.975,k); proc print data=a noobs; run; 第11页，共43页，2022年，5月20日，15点22分，星期五 程序说明： ci1=cinv(0.025,k); ci2