[程式] [R] 用ks.test測試3 parameter weibull

[wertymin24 (werty)]

TITLE請著名軟體類別
R
[軟體程式類別]:
R

[程式問題]:
統計和軟體

[軟體熟悉度]:
高(1年以上)

[問題敘述]:
我已經用fitdistr(MASS)做出 3 parameter的weibull函數

想進一步的用ks.test測goodness of fit

但是ks.test好像只能測 2 parameter的? (google過沒有找到其他方法)

如果我先把自己的資料減掉fit出來的location可以再用一樣的shape和scale去測嗎?

p.s. 我是因為發現ks.test在測時 location跑掉會有不一樣的結果才這樣想的
     如下範例

[程式範例]:
#先做一筆資料
temp<-rweibull(1000,shape=1.75,scale=18)
#goodness of fit ,  pvalue=0.99 (當然..)
ks.test(temp,'pweibull',shape=1.75,scale=18)
#位移
temp<-temp+30
#pvalue<0.0001
ks.test(temp,'pweibull',shape=1.75,scale=18)



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

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 118.168.73.189
→ clickhere:y 不要給 'pweibull', 把相對應於x的cdf先求出來, 再餵
→ clickhere:給ks.test(x,y)就可以了. 例如, test location 28, 則
→ clickhere:y <- pweibull(temp - 28, 1.75, 18)
→ wertymin24:這樣子會不會其實不是測到3 parameter的goodness?
→ clickhere:眼殘了, y要給的是data不是cdf. sorry
→ clickhere:可能要改source code.
→ clickhere:x <- y(sort(x), ...) - (0:(n - 1))/n 這行要改成
→ clickhere:x <- pweibull(sort(x) - 28, ...) - (0:(n - 1))/n
→ wertymin24:樓上的意思是我要去改ks.test的ocde嗎?
→ wertymin24:我不懂(0:(n - 1))/n這是什麼?
→ clickhere:用ks.test為樣版,按定義覆寫它. emperical cdf of x.
→ clickhere:下面可能可行, 也比較容易. 寫個新函式叫 pweibull.new
→ clickhere:然後把 'pweibull.new' 當 y 丟給 ks.test 就行了.