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

看板Statistics作者 (werty)時間14年前 (2011/12/16 12:39), 編輯推噓0(0推 0噓 10→)
留言10則, 2人參與, 最新討論串2/2 (看更多)
※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 118.168.73.189
→
12/15 22:33,
y 不要給 'pweibull', 把相對應於x的cdf先求出來, 再餵
12/15 22:33
→
12/15 22:35,
給ks.test(x,y)就可以了. 例如, test location 28, 則
12/15 22:35
→
12/15 22:35,
y <- pweibull(temp - 28, 1.75, 18)
12/15 22:35
→
12/15 22:47,
這樣子會不會其實不是測到3 parameter的goodness?
12/15 22:47
→
12/15 22:52,
眼殘了, y要給的是data不是cdf. sorry
12/15 22:52
→
12/15 22:53,
可能要改source code.
12/15 22:53
→
12/15 22:54,
x <- y(sort(x), ...) - (0:(n - 1))/n 這行要改成
12/15 22:54
→
12/15 22:54,
x <- pweibull(sort(x) - 28, ...) - (0:(n - 1))/n
12/15 22:54
→
12/15 23:12,
樓上的意思是我要去改ks.test的ocde嗎?
12/15 23:12
→
12/15 23:12,
我不懂(0:(n - 1))/n這是什麼?
12/15 23:12
→
12/15 23:41,
用ks.test為樣版,按定義覆寫它. emperical cdf of x.
12/15 23:41
→
12/16 00:11,
下面可能可行, 也比較容易. 寫個新函式叫 pweibull.new
12/16 00:11
→
12/16 00:11,
然後把 'pweibull.new' 當 y 丟給 ks.test 就行了.
12/16 00:11
不好意思 可能是我沒描述清楚 也可能是我還沒看懂這部分的意思 我想知道的是，我可以直接平移我的weibull分佈讓他reduce為2 parameter的嗎? 或者我可以平移我的data讓我的data符合location=0,但相同scale和shape的weibull嗎? click大的意思看起來像是平移出一個新的weibull分佈是這樣子嗎? 謝謝 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.112.82.37
→
12/16 21:26, , 1F
這不就是location parameter的意思嘛!
12/16 21:26, 1F
→
12/16 21:27, , 2F
平移後law/distribution是相同的.
12/16 21:27, 2F
→
12/16 21:27, , 3F
www.weibull.com/hotwire/issue15/relbasics15.htm
12/16 21:27, 3F
→
12/16 21:29, , 4F
和Normal的mu是類似的意思.
12/16 21:29, 4F
→
12/16 21:59, , 5F
非常感謝你~"~ 我已經google一天了
12/16 21:59, 5F
→
12/17 03:11, , 6F
再想想,其實也沒必要寫新函式,直接從data減掉,再把新的
12/17 03:11, 6F
→
12/17 03:12, , 7F
data代進 x 就行了. 真的是腦殘了.
12/17 03:12, 7F
→
12/17 03:12, , 8F
sorry.
12/17 03:12, 8F
→
12/17 12:22, , 9F
不會~因為我其實也不是做統計的，我們系上統計很弱
12/17 12:22, 9F
→
12/17 12:22, , 10F
常常找不到人討論統計問題~有人願意跟我討論就很好了
12/17 12:22, 10F
‣ 更多 wertymin24 的文章
文章代碼(AID): #1EwikLCp (Statistics)
‣ 更多 wertymin24 的文章
文章代碼(AID): #1EwikLCp (Statistics)