24小时热门版块排行榜    

查看: 6392  |  回复: 11

漫天雕

金虫 (正式写手)

[求助] 数理统计高手请进:如何计算KS检验中的p-value 已有1人参与

如题,在假设检验中要求采用KS检验。按一些参考资料上所述,KS值为经验累计分布函数与假定累计分布函数之差的最大值,但我没有看到与该检验相关的p-value的定义。更不知道求解方法了?
不知道那位数理统计的高手能够指点指点?

[ 发自手机版 http://muchong.com/3g ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

宝贵青春,耗不起!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

2楼2013-12-13 15:18:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

漫天雕

金虫 (正式写手)

引用回帖:
2楼: Originally posted by feixiaolin at 2013-12-13 15:18:33
http://www.mathworks.cn/matlabcentral/newsreader/view_thread/249306

不是太明白,能不能简单表述一下?

[ 发自手机版 http://muchong.com/3g ]
宝贵青春,耗不起!
3楼2013-12-13 15:51:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
3楼: Originally posted by 漫天雕 at 2013-12-13 15:51:08
不是太明白,能不能简单表述一下?
...

公式
s = n*KSstatistic^2;
pValue = 2*exp(-(2.000071+.331/sqrt(n)+1.409/n)*s);
4楼2013-12-13 15:56:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

漫天雕

金虫 (正式写手)

引用回帖:
4楼: Originally posted by feixiaolin at 2013-12-13 15:56:26
公式
s = n*KSstatistic^2;
pValue = 2*exp(-(2.000071+.331/sqrt(n)+1.409/n)*s);...

这个计算公式与假设分布有关系吗?

[ 发自手机版 http://muchong.com/3g ]
宝贵青春,耗不起!
5楼2013-12-13 16:00:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
5楼: Originally posted by 漫天雕 at 2013-12-13 16:00:15
这个计算公式与假设分布有关系吗?
...

n样点数;KSstatistic=K-S值
6楼2013-12-13 16:02:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

漫天雕

金虫 (正式写手)

引用回帖:
6楼: Originally posted by feixiaolin at 2013-12-13 16:02:53
n样点数;KSstatistic=K-S值...

谢谢你,这个公式的计算我明白了,只是我有个疑问。一般来说,计算p-value需要知道统计量的分布,那KS检验的统计量服从什么分布?这种分布与检验中数据的原假设分布是否有关系,例如对于原假设为正态分布和Gamma分布,它们的统计量的分布是否一样?

[ 发自手机版 http://muchong.com/3g ]
宝贵青春,耗不起!
7楼2013-12-13 16:13:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
7楼: Originally posted by 漫天雕 at 2013-12-13 16:13:18
谢谢你,这个公式的计算我明白了,只是我有个疑问。一般来说,计算p-value需要知道统计量的分布,那KS检验的统计量服从什么分布?这种分布与检验中数据的原假设分布是否有关系,例如对于原假设为正态分布和Gamma分 ...

自己多: 谷歌学术,wiki百科,cnki.net,
8楼2013-12-13 16:49:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li8512420

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
R语言KS程序

X<-c(420, 500, 920, 1380, 1510, 1650, 1760, 2100, 2300, 2350)
ks.test(X, "pexp", 1/1500)   #pexp:指数分布

检验是否服从参数为 1/1500的指数分布


运行结果:
One-sample Kolmogorov-Smirnov test

data:  X
D = 0.3015, p-value = 0.2654
alternative hypothesis: two-sided
要不屈不挠,屡败屡战才能成为真正的强者——丘成桐
9楼2013-12-14 08:58:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li8512420

木虫 (正式写手)

【答案】应助回帖

MATLAB的统计工具箱里也有
[h,p,stat,cv]=kstest(x,cdf,alpha,tail)

h,决策值
p,P值
stat,测试统计量的值
cv拒绝域的临界值
x,样本
cdf,累积分布函数
alpha,显著性水平
tail,0双边检验,1右边检验,-1左边检验
要不屈不挠,屡败屡战才能成为真正的强者——丘成桐
10楼2013-12-14 09:05:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 漫天雕 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 16/800 2026-08-07 13:25 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +23 低垂的野花 2026-07-31 31/1550 2026-08-07 12:22 by speculiar
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 基金中了 +11 laoda193707 2026-08-06 11/550 2026-08-07 10:54 by TT和NN
[基金申请] 听说今天filecode变了 +21 布布和一二 2026-08-06 44/2200 2026-08-07 10:11 by 布布和一二
[基金申请] 固定端突然变了,今天 +4 archvillain 2026-08-06 8/400 2026-08-07 08:06 by 带你看大海
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了,论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见