楼主: ihsihs
4613 2

[问答] ks.test参数设置疑问 [推广有奖]

  • 1关注
  • 1粉丝

博士生

21%

还不是VIP/贵宾

-

威望
0
论坛币
612 个
通用积分
7.4619
学术水平
6 点
热心指数
6 点
信用等级
1 点
经验
2918 点
帖子
127
精华
0
在线时间
306 小时
注册时间
2008-6-18
最后登录
2024-4-17

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
我用一直分布生成了一组数据,我想看看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吧?
二维码

扫码加我 拉你入群

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

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

关键词:参数设置 test Est Exponential expo

离2200还有多远
沙发
ihsihs 发表于 2017-10-17 05:31:45 |只看作者 |坛友微信交流群
那里是没有指定指数分布的参数lamada

使用道具

藤椅
ihsihs 发表于 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)

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-24 16:27