非参数统计检验方法的应用

1、论文投稿领域：数理经济与计量经济学非参数统计检验方法的应用阮曙芬1 程娇翼1 张振中2（1. 中国地质大学数理学院，武汉 430074；2.中南大学数学科学与计算学院，长沙410075）摘要：本文对非参数统计中常用的三种假设检验方法进行了简单的介绍。运用Kruskal-Wallis检验方法对2002年前三季度的上海股市综合指数收益率数据进行了周末效应的检验，结果表明2002年上海股市综合指数收益率不具有周末效应。关键字：符号检验；Wilcoxon秩和检验；Kruskal-Wallis检验1引言非参数统计是统计分析的重要组成部分。非参数假设检验是在总体分布未知或者总体分布不满足参数统计对总体所做

2、的假定的时候，分析样本特点，寻找相应的非参数检验统计量。本文就是以此为出发点，介绍了非参数统计中假设检验常用的几个检验方法：符号检验、Wilcoxon秩和检验和Kruskal-Wallis检验，然后结合具体的问题和数据，在统计软件SAS中作相应的非参数检验。2非参数假设检验介绍2.1 配对样本的符号检验符号检验是根据正、负符号进行假设检验的方法。这种检验方法用于配对设计数值变量资料的假设检验，常常是差值不服从正态分布或者总体分布未知的情况下不能用检验的时候使用。其原理是对差值进行编制并冠以符号，然后对正负秩和进行比较检验。设随机变量相互独立同分布，分布为，在连续。假设检验问题： ：检验统计量可

3、取，即为中取正号的个数。在下，的分布是参数为和的二项分布。2.2 两独立样本的Wilcoxon秩和检验Wilcoxon秩和检验的理论背景如下：有两个总体，一个总体的样本为，相互独立同分布，分布为；另一个样本为，相互独立同分布，分布为，连续。问随机变量是否随机大于随机变量，即检验：，：，且有某些点不等号成立。将，共个随机变量一起排序，产生对应的秩。则Wilcoxon秩和检验统计量为：即在混合样本中的秩的和为Wilcoxon秩和检验统计量。2.3多样本的Kruskal-Wallis检验Kruskal-Wallis检验一般对多个总体的分布情况进行检验。其理论基础为：假设有种处理，对于第个检验体实行第

4、种处理产生的效果记为，其分布函数为。即：；：存在和，。 设观测值为。全体样本数为，的顺位记为。假定同顺位不会出现。考虑统计量：检验方法为：。 其中，。当较小时，可以查表得到的值；当该值较大时，近似服从自由度为的分布。因此，为自由度为的分布的右侧的分位数点。3 Kruskal-Wallis检验的应用股市的周末效应是指周一的收益率比其他交易日收益率低，且风险较大；周五的收益率比其他交易日高，且相对风险较小。下面分别对2002年的前三季度的上证综合指数进行周末效应的分析。本实证分析中，样本为2002年1月4日到2002年9月27日的上海股市综合指数（数据来源于）。指数收益率的计算公式为：，其中为第天

5、的指数，为第天的指数收益率。3.1收益率分布状况的分析首先计算收益率序列的方差，均值，偏度和峰度初步判断该序列是否服从正态分布。然后利用Kolmogorov-Smirnov等检验结果对收益率进行正态性检验。 SAS程序创建数据集：将excel数据导入SAS中，然后在分析家中利用数据计算得到:r0=p/lag1(p)和r=log(r0);data sasuser.chx1 sasuser.chx2 sasuser.chx3 sasuser.chx4 sasuser.chx5;set sasuser.ch01;select (w); when(1) output sasuser.chx1; whe

6、n(2) output sasuser.chx2; when(3) output sasuser.chx3; when(4) output sasuser.chx4; when(5) output sasuser.chx5;end;run; proc univariate data=sasuser.ch01; var r;run;SAS结果输出见表1汇总表1 上证指数收益率描述性统计分析星期周一周二周三周四周五全体数据均值-3.882 E-34.875 E-30.422 E-31.363 E-3-3.423 E-3-0.110 E-3t统计量-1.268（0.214）1.399（0.171）0

7、.180（0.858）0.423（0.675）-1.776（0.085）-0.084（0.933）自由度3434343433169方差0.319 E-30.413 E-30.186 E-30.352 E-30.123 E-30.283 E-3偏度-1.0692.2922.6960.992-0.6301.144峰度3.6667.85012.4333.6460.3897.084由上表可知，上证指数收益率序列的偏度和峰度分别为1.144和 7.084，而正态分布的偏度和峰度分别为0和3，所以我们可以初步断定指数收益率序列为非正态分布。为了进一步证实这一论断，我们对收益率序列进行Kolmogorov-

8、Smirnov检验。SAS自动输出包括Kolmogorov-Smirnov检验统计量在内的四种检验正态分布的检验统计量。SAS程序proc univariate data=sasuser.chx1 normal; var r;histogram r; probplot r; run;SAS结果输出见表1汇总，图1和图2。表2 上证指数收益率的正态性检验Kolmogorov-Smirnov检验对应的p值Shapino-Wilk检验对应p值Cramer-von Mises检验对应的p值Anderson-Darling检验对应的p值自由度周一0.04640.01060.02720.021134周二0

9、.01000.00010.00500.005034周三0.01000.00010.00500.005034周四0.09100.01260.01640.017234周五0.15000.10100.15280.096233全体数据0.01000.00010.00500.0050169图1 上证综合指数收益率分布的直方图图2 上证综合指数收益率分布的概率图包括Kolmogorov-Smirnov检验统计量在内的四种检验正态分布的检验统计量均表明上海综合指数收益率序列不服从正态分布，图1和图2也说明了这一点。所以要采用非参数方法进行以后的周末效应的检验。3.2周末效应存在性的Kruskal-Walli

10、s检验我们利用Kruskal-Wallis检验2002年前三季度上证综合指数收益率的周末效应的存在性。SAS程序proc npar1way wilcoxon data=sasuser.ch01; class w; var r;run;SAS结果输出-The NPAR1WAY ProcedureWilcoxon Scores (Rank Sums) for Variable rClassified by Variable ww N Sum of Scores Expected Under H0 Std Dev Under H0 Mean Score5 33 2576.0 2805.0 252.1

11、50749 78.0606061 34 2610.0 2890.0 255.000000 76.7647062 34 3206.0 2890.0 255.000000 94.2941183 34 2996.0 2890.0 255.000000 88.1176474 34 2977.0 2890.0 255.000000 87.558824Kruskal-Wallis TestChi-Square 3.0846DF 4Pr > Chi-Square 0.5438- KW检验得3.086，4，0.5348，所以不能拒绝,即周一到周五得上证综合指数收益率得分布，所以我们认为在2002年的前三

12、季度中，上海市股市综合指数收益率不存在周末效应。参考文献：1 Damodar N. Gujarati. Basic Econometrics.北京：中国人民大学出版社，2005. p791-p800.2 George E. P. Box, Gwilym M. Jenkins, Gregory C. Reinsel. Time Series Analysis Forecasting And Control. 3 何书元. 应用时间序列分析. 北京：北京大学出版社, 2003. p218-p226.4 张卓. SAS软件的应用. 统计与信息论坛(2005),Vol.20, No.4. p104-p

13、106.5 樊欣，邵谦谦.SAS 8.X 经济统计.北京：北京希望电子出版社，2003. p28-p60.6 岳朝龙，黄永兴，严钟. SAS系统与经济统计分析. 合肥：中国科学技术大学出版社, 2004. p469-p4877 李彦萍.发达与非发达地区收入与消费非参数统计分析.山西农业大学学报(2005)，Vol.4， No.4. p334-p339.8 刘彤.利用非参数方法对上海股市周末效应的研究.数理统计与管理(2003)，Vol.22，No.1. p69-p71.Application of Nonparametric statistical MethodRUAN Shu-fen，CHE

14、NG Jiao-yi，ZHANG Zhen-zhong(School of Mathematics and Physics, China University of Geosciences, Wuhan 430074)Abstract: In this paper, we simply introduce  three common hypothesis tests. Using the Kruskal- -Wallis test, we do week effect test about the shanghai synthetic  index of which we take the first three quarte