Applied Nonparametric Statistics-lec5
今天继续two-sample test
Ref: https://onlinecourses.science.psu.edu/stat464/print/book/export/html/6
- Mann-Whitney Test
前面说这个和Wilcoxon是identical的,只是统计量不同。现在我们来看一下它的统计量U。注意,现在检查的仍然是两个独立样本。
Treatment 1: x1, x2, ... , xm
Treatment 2: y1, y2, ... , yn
U = # of pairs of (Xi, Yj) for which Xi < Yj
H0 : the distributions are the same
H1 : the distributions are not the same
Table A4 contains lower tailed and upper tailed values for U under the null hypothesis. It can be shown that Uupper = mn - Ulower.
实际在R中的操作:
wilcox.test(new, trad, alternative="greater")
现在,我们来构造U的置信区间。好了,直接用R就好了啊:)
wilcox.test(ugrad, grad, conf.level=0.90, conf.int=T)
以上我们检查的都是两个样本的位置情况,下面要检查scale,也就是分布的形状,variability。
我们现在假设两个样本的均值相等,方差不同,现在想判断哪个方差更大。如果两个样本来自正态分布的总体,那么可以分别计算方差,然后
但是如果这两个不是来自正态分布的,就不可以了。此时,我们考虑非参数检验的方法。
- Siegel-Tukey Test (ST test):检查方差
- 把数据合在一起,从小到大排列。
- 最小的rank为1,最大的为2,次小的为3,次大的为4,蛇形排序。
- 做Wilcoxon rank-sum test. The smaller rank sums are associated with the treatment that has the largest variability.
如果把指定rank的方式改为:最小和最大为1,次小和次大为2,以此类推,这样就是Ansari-Bradley test。
R下,可以使用jmuOutlier包内的函数
siegel.test(x, y, alternative = c("two.sided", "less", "greater"), reverse = FALSE,all.perms = TRUE, num.sim = 20000)
Ansari-Bradley test:
ansari.test(x, y,
alternative = c("two.sided", "less", "greater"),
exact = NULL, conf.int = FALSE, conf.level = 0.95,
...)
事实上,检查方差一致性homogeneity of variance (homoscedasticity)还有其他方法,参考这篇文章:
http://www.cookbook-r.com/Statistical_analysis/Homogeneity_of_variance/
- Bartlett’s test:适合于数据是正态分布的;
- Levene’s test :在car包里,对于非正态分布的数据,比Bartlett's test更具健壮性;
- Fligner-Killeen test :非参数检验方法。
- Tests on Deviances
1. Obtain the deviances for the two treatments (devix and devjy) and compute RMD from the orginal data, denoted RMDobs.
2. Permute.样本数大的话,循环指定次数,否则循环所有。计算RMD。
- Kolmogorov-Smirnov Test
这个是检查general difference的,也就是说,会考虑location,scale,shape。如果两个样本分布的位置是否不同这一点未知,使用这个;如果已经
知道两个样本数据分布的位置是不同的,那么应该用wilcoxon。
- Calculate the observed test statistic, KSobs
^F1(W)是sample CDF
2. Find all the possible permutation of the data and calculate KS for each permutation.
3. The p-value is found by
R语言就简单了:
ks.test(a, b)
Applied Nonparametric Statistics-lec5的更多相关文章
- Applied Nonparametric Statistics-lec10
Ref:https://onlinecourses.science.psu.edu/stat464/print/book/export/html/14 估计CDF The Empirical CDF ...
- Applied Nonparametric Statistics-lec9
Ref:https://onlinecourses.science.psu.edu/stat464/print/book/export/html/12 前面我们考虑的情况是:response是连续的, ...
- Applied Nonparametric Statistics-lec8
Ref:https://onlinecourses.science.psu.edu/stat464/print/book/export/html/11 additive model value = t ...
- Applied Nonparametric Statistics-lec7
Ref: https://onlinecourses.science.psu.edu/stat464/print/book/export/html/9 经过前面的步骤,我们已经可以判断几个样本之间是否 ...
- Applied Nonparametric Statistics-lec6
Ref: https://onlinecourses.science.psu.edu/stat464/print/book/export/html/8 前面都是对一两个样本的检查,现在考虑k个样本的情 ...
- Applied Nonparametric Statistics-lec4
Ref: https://onlinecourses.science.psu.edu/stat464/print/book/export/html/5 Two sample test 直接使用R的t- ...
- Applied Nonparametric Statistics-lec3
Ref: https://onlinecourses.science.psu.edu/stat464/print/book/export/html/4 使用非参数方法的优势: 1. 对总体分布做的假设 ...
- Applied Nonparametric Statistics-lec2
Ref: https://onlinecourses.science.psu.edu/stat464/print/book/export/html/3 The Binomial Distributio ...
- Applied Nonparametric Statistics-lec1
参考网址: https://onlinecourses.science.psu.edu/stat464/node/2 Binomial Distribution Normal Distribution ...
随机推荐
- 如何安装使用windows自带的telnet服务
控制面板->程序和功能->打开或关闭Windows功能->Telnet 客户端 [ Telnet 服务器 ] 安装完成后重启cmd telnet ip port
- node.js安装Oracledb指导文档
https://community.oracle.com/docs/DOC-931127
- C语言-字符操作函数
1字符数组的初始化: 1.1 char string={'c','h','i','n','a'} 1.2char string={"china"}或者去掉{}即char strin ...
- 渣渣菜鸡的 ElasticSearch 源码解析 —— 启动流程(下)
关注我 转载请务必注明原创地址为:http://www.54tianzhisheng.cn/2018/08/12/es-code03/ 前提 上篇文章写完了 ES 流程启动的一部分,main 方法都入 ...
- wepy开发踩坑记录
与vue的不同 methods对象只存放tap等事件触发时的方法 events对象只存放$emit及$broadcast方法触发的事件 自定义方法及属性放在与methods平级的位置 props是动态 ...
- zuul的本地跳转
- 在IDEA中编辑struts国际化properties文件
在IDEA中编辑struts国际化properties文件 如果手工创建的web工程,struts的i18n属性文件,可以使用native2ascii工具转换(记得命令行的第二个文件名是要保存的文件名 ...
- SpringBoot 封装返回类以及session 添加获取
1.创建返回类Result public class Result<T>{ /*错误码*/ private Integer code; /*提示信息 */ private String m ...
- Redis数据库2
一.存储list 1.概述 list类型是按照插入顺序排序的字符串链表,可以在头部(left)和尾部(right)添加新的元素,如果不存在,则自动创建;如果list中所有元素被移除,则list ...
- 代码管理_Git中获取两个标签之间的变更代码行数
操作步骤: 1.先进入git的服务器,定位到 repositories 目录 2.再定位到具体的项目目录(登录bitbucket查看项目属性可以看到项目的存放目录),如下图: 3.浏览项目的标签,命令 ...