全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
5173 2
2017-10-17
我用一直分布生成了一组数据,我想看看ks检验的效力
这是我的程序:
x=rnorm(100)
c=abs(x)
z=fitdistr(c,"exponential")

ks.test(c,"pexp")
ks.test(c,"pexp",z$e)

我发现这两次检验的结果不一样,我想知道的是在ks.test(c,"pexp")没有指定指数,那么它用的是什么呢,及然后和后面的结果不一样,就应该不是MLE吧?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2017-10-17 05:31:45
那里是没有指定指数分布的参数lamada
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-10-17 06:03:37
通过以下程序的试探,我想我明白了,没有给定参数的,应该默认的是标准指数,比如指数分布,默认lemada=1,正态分布默认两个参数为0,1

y=rexp(100,2)
z=fitdistr(y,"exponential")
a=mean(y)
b=1/a
ks.test(y,"pexp")
ks.test(y,"pexp",b)
ks.test(y,"pexp",z$e)
ks.test(y,"pexp",2)

y=rexp(1000,2)
z=fitdistr(y,"exponential")
a=mean(y)
b=1/a
ks.test(y,"pexp")
ks.test(y,"pexp",b)
ks.test(y,"pexp",z$e)
ks.test(y,"pexp",2)

y=rexp(1000,1)
z=fitdistr(y,"exponential")
a=mean(y)
b=1/a
ks.test(y,"pexp")
ks.test(y,"pexp",b)
ks.test(y,"pexp",z$e)
ks.test(y,"pexp",1)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群