R中验证正太分布检验及求取其方差和均值

àì夳堔傛蜴生んèń 2022-08-12 12:00 117阅读 0赞

**什么是正太分布检验？**  
判断一样本所代表的背景总体与理论正态分布是否没有显著差异的检验。  
  
**方法一****概率密度曲线比较法**  
看样本与正太分布概率密度曲线的拟合程度，R代码如下：

**\[plain\]** [view plain][] [copy][view plain] [![在CODE上查看代码片][CODE]][CODE_CODE] [![派生到我的代码片][ico_fork.svg]][ico_fork.svg 1]

1.  norm\_expression <- function(x) (1/sqrt(2\*pi))\*exp(-0.5\*x^2) 
2.  \#curve(norm\_expression, -4, 4, col="red") \#标准正太分布概率密度曲线 
3.  
4.  \#画样本概率密度图 
5.  s <- rnorm(100) \#产生样本 
6.  d <- density(s) 
7.  plot(d, col="green", ylim=c(0, 0.5)) 
8.  
9.  \#添加正太分布概率密度图 
10. s2 <- seq(from=-4, to=4, length.out=100) 
11. lines(s2, norm\_expression(s2), col="red")

画图结果如下：

![20140412174333234][]

**方法二****正太Q-Q图法**

使用Q-Q图来判断数据是否服从正太分布，R代码如下：

**\[plain\]** [view plain][] [copy][view plain] [![在CODE上查看代码片][CODE]][CODE_CODE] [![派生到我的代码片][ico_fork.svg]][ico_fork.svg 1]

1.  s <- rnorm(100) \#产生样本 
2.  qqnorm(s) 
3.  qqline(s)

画图结果如下，可见数据分布集中在对角线上，可以认为总体服从正太分布：

![20140412174406718][]

**方法三经验法则**

约68.3%数值分布在距离平均值有1个标准差之内的范围，约95.4%数值分布在距离平均值有2个标准差之内的范围，以及约99.7%数值分布在距离平均值有3个标准差之内的范围。称为“68-95-99.7法则”或“经验法则”。

使用R的验证代码如下：

**\[plain\]** [view plain][] [copy][view plain] [![在CODE上查看代码片][CODE]][CODE_CODE] [![派生到我的代码片][ico_fork.svg]][ico_fork.svg 1]

1.  s <- rnorm(10000) \#产生样本 
2.  sum(abs(s - mean(s)) < sd(s)) / length(s) 
3.  sum(abs(s - mean(s)) < 2\*sd(s)) / length(s)

程序结果：

> s <- rnorm(10000) \#产生样本  
> sum(abs(s - mean(s)) < sd(s)) / length(s)  
\[1\] 0.6871  
> sum(abs(s - mean(s)) < 2\*sd(s)) / length(s)  
\[1\] 0.9538

**方法四****统计检验方法**

使用样本偏度和样本峰度来估计总体偏度和峰度，在正太分布的假定下，样本偏度和峰度均服从均值为零、方差分别为6/T和24/T的正太分布，可以分别检验偏度和峰度，也可以将两个统计量结合起来生成一个服从自由度为2的卡方分布的统计量，再进行检验【参见《金融时间序列分析》第三版P8~P9】。理论方面的东西略...

可以使用夏皮罗-威尔克（Shapiro-Wilk）检验，代码如下：

**\[plain\]** [view plain][] [copy][view plain] [![在CODE上查看代码片][CODE]][CODE_CODE] [![派生到我的代码片][ico_fork.svg]][ico_fork.svg 1]

1.  s <- rnorm(1000) \#产生样本 
2.  shapiro.test(s)

检验结果:

> shapiro.test(s)  
  
  
Shapiro-Wilk normality test  
  
  
data: s  
W = 0.9987, p-value = 0.6716

shapiro.test函数输出一个p值，照惯例，p<0.05说明总体不太可能是正太分布，否则不能提供这么个证据，也就是说这个检验比较保守，倾向于错误的过分证明正态性。

R的nortest包能提供其他正太检验方法【参见《R语言经典实例》P221】。

**方法五经验法**

见[智库百科][Link 1]，使用R实现一遍感觉有误，仅作参考。

用R软件来实现求均值 顺序统计量 中位数 方差标准差 峰度系数 偏度系数

avg <- mean(test)
    mid <- median(test)
    var <- var(test)
    std <- sd(test)
    
    library(e1071)
    kurtosis <- kurtosis(test)
    skewness <- skewness(test)

[view plain]: http://blog.csdn.net/yucan1001/article/details/23539639#
[CODE]: https://code.csdn.net/assets/CODE_ico.png
[CODE_CODE]: https://code.csdn.net/snippets/288555
[ico_fork.svg]: https://code.csdn.net/assets/ico_fork.svg
[ico_fork.svg 1]: https://code.csdn.net/snippets/288555/fork
[20140412174333234]: /images/20220810/416c5261cb0e4a01bdcaed5227a03e65.png
[20140412174406718]: /images/20220810/e852c4bfb67f4710b879637224bf56f6.png
[Link 1]: http://wiki.mbalib.com/wiki/%E6%AD%A3%E6%80%81%E5%88%86%E5%B8%83%E6%A3%80%E9%AA%8C