R语言对数据进行非参数检验

逃离我推掉我的手 2022-05-24 06:55 674阅读 0赞

假设检验：参数检验运用样本的统计量来估计总体的参数，如用样本均值估计总体均值，用样本标准差估计总体标准差。

非参数检验则不考虑数据的具体值，而更多地运用了数据大小排序的信息，因此不可能以此估计总体的参数

### **1.原假设和备择假设** ###

原假设(null hypothesis)——原假设又称“ 0假设”，总是有符号 =， ≥ 或≤，表示为 H0。是研究者想收集证据予以反对的假设（生产实践中常对应正常情形，如均值与设计一致）；一般来说，原假设是一旦拒绝便要采取行动的假设。因此， 原假设总是“受到保护的假设” ，没有充分的证据是不能拒绝原假设的。例如，对一家信誉很好的工厂的产品进行检验，原假设一般是“ 产品合格”。备择假设(alternative hypothesis)——研究者想收集证据予以支持的假设， 一旦发生就要采取行动， 是与原假设对立的假设，也称“研究假设”，总是有符号 ≠， > 或 <，表示为 H1。

同时在实际应用中，我们有不同的需求，因此又有双侧检验和单侧检验的区分。

>  *  双侧检验——备择假设没有特定的方向性，并含有符号“=”的假设检验，称为双侧检验或双尾检验(two-tailed test)
>  *  单侧检验——备择假设具有特定的方向性，并含有符号“>”或“<”的假设检验，称为单侧检验或单尾检验(one-tailed test)。其中备择假设的方向为“<”，称为左侧检验，备择假设的方向为“>”，称为右侧检验。

原假设与备择假设形式：

![20180522220406215][]

双侧检验：

![format_png][]

左侧检验：

![format_png 1][]

### **2.非参数检验：** ###

1.正态W检验

例1. 已知15名学生体重如下，问是否服从正态分布

![format_png 2][]

解：

R语言代码：

w <- c(75.0, 64.0, 47.4, 66.9, 62.2, 62.2, 58.7, 63.5,  
66.6, 64.0, 57.0, 69.0, 56.9, 50.0, 72.0)  
shapiro.test(w)

P值>0.05，接受原假设，认为来自正态分布总体。

例2 ：H0:数据集服从正态分布

H1：数据集不服从正态分布

#是否服从正态分布
    #随机抽取3000个样本
    mysample <- error[sample(1:nrow(error), 3000, replace=FALSE),] 
    shapiro.test(mysample$B1)
    shapiro.test(mysample$Q1)

\#显示p值小于2.2e-16,因此认为不服从正太分布

![20180522203555288][]

2.Kolmogorov-Smirnov 检验

例3. 假定从分布函数未知的F（x）和G（x）的总体中分别抽出25个和20个观察值的随即样品，其数据由下表所示。现检验F（x）和G(x)是否相同。

![format_png 3][]

R语言代码：

X<-scan( )  
0.61 0.29 0.06 0.59 -1.73 -0.74 0.51 -0.56  
1.64 0.05 -0.06 0.64 -0.82 0.37 1.77  
2.36 1.31 1.05 -0.32 -0.40 1.06 -2.47  
0.39 1.09 -1.28  
Y<-scan( )  
2.20 1.66 1.38 0.20 0.36 0.00  
0.96 1.56 0.44 1.50 -0.30 0.66

2.31 3.29 -0.27 -0.37 0.38 0.70  
0.52 -0.71  
ks.test(X,Y)

P值>0.05,无法拒绝原假设，说明F（x）和G(x)分布函数相同。

例子：

> ks.test(error$B1,error$Q1)
    
    
    	Two-sample Kolmogorov-Smirnov test
    
    
    data:  error$B1 and error$Q1
    D = 0.11508, p-value < 2.2e-16
    alternative hypothesis: two-sided	Two-sample Kolmogorov-Smirnov test
    
    
    data:  error$B1 and error$Q1
    D = 0.11508, p-value < 2.2e-16
    alternative hypothesis: two-sided

P值<0.05,可以拒绝原假设，说明B1,Q1分布函数相同。

**3.列联表数据的检验**

例4. 为了研究吸烟是否与患肺癌相关，对63位肺癌患者及43名非肺癌患者（对照组）调查了其中的吸烟人数，得到2x2列联表，如下表所示

![format_png 4][]

解：

进行Pearson卡方检验

R语言代码：

x<-c(60, 3, 32, 11)  
dim(x)<- c(2,2)  
chisq.test(x,correct = F)

![format_png 5][]

P值<0.05，拒绝原假设，认为吸烟与患肺癌相关。

4.符号检验

例5. 用两种不同的饲料养猪，其增重情况如下表所示。试分析两种饲料养猪有无显著差异。

![format_png 6][]

R语言代码：

x<-scan()  
25 30 28 23 27 35 30 28 32 29 30 30 31 16  
  
y<-scan()  
19 32 21 19 25 31 31 26 30 25 28 31 25 25  
  
binom.test(sum(x<y), length(x))

![format_png 7][]

sum(x < y)表示样品X小于样品Y的个数。计算出P值大于0.05，无法拒绝原假设，可以认为两种饲料养猪无显著差异。

例子6：H0：两组数据无显著的差别

H1：两组数据有明显的差别

binom.test(sum(error$B1<error$Q1), length(error$B1))

> binom.test(sum(error$B1<error$Q1), length(error$B1))
    
    	Exact binomial test
    
    data:  sum(error$B1 < error$Q1) and length(error$B1)
    number of successes = 34950, number of trials = 65535, p-value <2.2e-16
    alternative hypothesis: true probability of success is not equal to 0.5
    95 percent confidence interval:
     0.5294742 0.5371285
    sample estimates:
    probability of success 
                 0.5333028

P<0.05，所以拒绝原假设，接受两组数据有明显的差别

在符号检验法中，只计算符号的个数，而不考虑每个符号差所包含的绝对值的大小，因此常常使用弥补了这个缺点的wilcoxon符号秩检验。

###  ###

### 5.符号秩检验 ###

在R中，wilcox.test()函数可以用来做Wilcoxon秩和检验，也可以用于做Mann-Whitney U检验。当参数为单个样本，或者是两个样本相减，或者是两个参数，paired=F时，是Wilcoxon秩和检验。当paired = FALSE（独立样本）时，就是Mann-Whitney U检验，

在R语言中进行符号秩检验可以使用wilcox.test( )

wilcox.test(x, y = NULL,  
alternative = c("two.sided", "less", "greater"),  
mu = 0, paired = FALSE, exact = NULL, correct = TRUE,  
conf.int = FALSE, conf.level = 0.95, ...)

n为整数向量,表示试验的次数,当x为矩阵时,该值无效.  
p为向量,表示试验成功的概率,必须与x有相同的维数,且值在0至1之间,默认值为NULL.  
alternative为备择假设选项,取"two.sided"(默认值)表示双侧检验;取"less"表示备择假设为“<”的单侧检验，取"greater"表示备择假设为“>”的单侧检验.  
conf.level为0～1之间的数值(默认值为0.95),表示置信水平,它将用于计算比率p或比率差p1-p2的置信区间.  
correct为逻辑变量,表示是否对统计量作连续修正,默认值为TRUE.

其中x,y是观察数据构成的数据向量。alternative是备择假设，有单侧检验和双侧检验，mu待检参数，如中位数M0.paired是逻辑变量，说明变量x,y是否为成对数据。exact是逻辑变量，说明是否精确计算P值，当样本量较小时，此参数起作用，当样本两较大时，软件采用正态分布近似计算P值。correct是逻辑变量，说明是否对P值的计算采用连续性修正，相同秩次较多时，统计量要校正。conf.int是逻辑变量，说明是否给出相应的置信区间。

\------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

例7. 今测得10名非铅作业工人和7名铅作业工人的血铅值，如下表所示。试用Wilcoxon秩和检验分析两组工人血铅值有无差异。

![format_png 8][]

解：进行Wilcoxon秩和检验R语言同样可以使用wilcox.test( )

R语言代码：

x<-c(24, 26, 29, 34, 43, 58, 63, 72, 87, 101)

y<-c(82, 87, 97, 121, 164, 208, 213)

wilcox.test(x,y,alternative="less",exact=FALSE,correct=FALSE)

![format_png 9][]

P值小于0.05，拒绝原假设，即铅作业工人血铅值高于非作业工人。

例子2：H0：两组数据无显著的差别

H1：两组数据有明显的差别

> wilcox.test(error$B1,error$Q1,alternative="two.sided",exact=FALSE,correct=TRUE)
    
    	Wilcoxon rank sum test with continuity correction
    
    data:  error$B1 and error$Q1
    W = 2103600000, p-value = 1.606e-10
    alternative hypothesis: true location shift is not equal to 0

P<0.05,所以拒绝原假设，两组数据是有差别的

6.二元数据相关检验

例8 .某种矿石中两种有用成分A，B，取10个样品，每个样品中成分A的含量百分数x(%)，及B的含量百分数y(%)的数据下表所示，对两组数据进行相关性检验。

![format_png 10][]

解：进行相关性检验，在R语言中可以使用cor.test( )

cor.test(x, y,  
alternative = c("two.sided", "less", "greater"),  
method = c("pearson", "kendall", "spearman"),  
exact = NULL, conf.level = 0.95, ...)

\#其中x,y是数据长度相同的向量，alternative是备择假设，缺省值为"two.sided"，method是检验方法，缺省值是Pearson检验，conf.level是置信区间水平，缺省值为0.95

cor.test( )还有另一种使用格式

cor.test(formula, data, subset, na.action, ...) \#其中formula是公式，形如'~u+v' , 'u', 'v' 必须是具有相同长度的数值向量，data是数据框，subset是可选择向量，表示观察值的子集。

假设此例中两组数据均来自正态分布，使用pearson相关性检验，

R语言代码：

ore<-data.frame(  
x=c(67, 54, 72, 64, 39, 22, 58, 43, 46, 34),  
y=c(24, 15, 23, 19, 16, 11, 20, 16, 17, 13)  
)  
cor.test(ore$x,ore$y)

![format_png 11][]

可见P值<0.05,拒绝原假设，认为X与Y相关。

--------------------

### 项目推荐： ###

**[2000多G的计算机各行业电子资源分享（持续更新）][2000_G]**

**[2020年微信小程序全栈项目之喵喵交友【附课件和源码】][2020]**

**[Spring Boot开发小而美的个人博客【附课件和源码】][Spring Boot]**

**[Java微服务实战296集大型视频-谷粒商城【附代码和课件】][Java_296_-]**

**[Java开发微服务畅购商城实战【全357集大项目】-附代码和课件][Java_357_-]**

**[最全最详细数据结构与算法视频-【附课件和源码】][-]**

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNTg3NTc1_size_16_color_FFFFFF_t_70_pic_center]

\------------------------- end -----------------------

--------------------

[20180522220406215]: /images/20220524/c9b8e382af0a4d3bbc2d605c9aba02df.png
[format_png]: /images/20220524/76f9d580ea524592b1fb644c38f4bef0.png
[format_png 1]: /images/20220524/2eec3c74c37c419787be64d75c162721.png
[format_png 2]: /images/20220524/25e9fff949784be4af4bcb05c1d61eee.png
[20180522203555288]: /images/20220524/db4035e041b14ebaa99a0c9e04faae83.png
[format_png 3]: /images/20220524/5644c25737ea45588c3868e7ae50d5db.png
[format_png 4]: /images/20220524/7b9b1c583f4b4c5a8e47b84f84fc9060.png
[format_png 5]: /images/20220524/01ffb19877ba4fa6b6d32e340c8c1d1c.png
[format_png 6]: /images/20220524/176b08efd1d346d080045e289bcb9dc7.png
[format_png 7]: /images/20220524/4bd271ee76f4485bbbae3c3395d346af.png
[format_png 8]: /images/20220524/65178deccfbd4249b51f66194f01ff97.png
[format_png 9]: /images/20220524/df4e8d7eca904b0382eb4c25b43e3e4d.png
[format_png 10]: /images/20220524/181feb5f5e7d41559b0ec9e67d792b05.png
[format_png 11]: /images/20220524/7e437bbfa6634ef5bb871e29bf92bc4d.png
[2000_G]: https://mp.weixin.qq.com/s/sP4JgGWkCzpgwKr9sAV2_Q
[2020]: https://mp.weixin.qq.com/s?__biz=MzIyNTI3NDQ4NQ==&mid=2247487704&idx=1&sn=5f4b2127c4d49fd07ae072a0721424a2&chksm=e8036fc2df74e6d489c4aa9b06f917ef7cee6027f13e150fca53cf79d5d188d4ccc1af49e098&scene=21#wechat_redirect
[Spring Boot]: https://mp.weixin.qq.com/s?__biz=MzIyNTI3NDQ4NQ==&mid=2247487798&idx=2&sn=ac0293b996521b872a9dba5fbb3e65e6&chksm=e8036e2cdf74e73aba104a9a994a5b2e31483e8dcbe0f1d9936f6d5173b887e1560f59d2819c&scene=21#wechat_redirect
[Java_296_-]: https://mp.weixin.qq.com/s?__biz=MzIyNTI3NDQ4NQ==&mid=2247487674&idx=1&sn=7aff0bdf2bb727303f3d3618995aef21&chksm=e8036fa0df74e6b6d872c7e6ece179c524ed463a4a6b74c96875475c9a3d5ddb903427dd993b&scene=21#wechat_redirect
[Java_357_-]: https://mp.weixin.qq.com/s?__biz=MzIyNTI3NDQ4NQ==&mid=2247486376&idx=1&sn=d1fef270c463ea8ac663f6fbfedd70a0&chksm=e80374b2df74fda4d3bafba878a106a19e18c5fcda266008f4f37975847a21bc612ffcd5ff39&scene=21#wechat_redirect
[-]: https://mp.weixin.qq.com/s?__biz=MzIyNTI3NDQ4NQ==&mid=2247487750&idx=1&sn=747bccbb5f5ea6b58915198de40da777&chksm=e8036e1cdf74e70ae97a5e8e265b49d7236d904d291203309159d07ba1724033062c0e370843&scene=21#wechat_redirect
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNTg3NTc1_size_16_color_FFFFFF_t_70_pic_center]: /images/20220524/5d9c0c5fa2f541f0b0d344ceccf2e5fe.png