Re: [問題] R 兩個問題

看板Statistics作者 (ideal)時間17年前 (2008/07/08 11:09), 編輯推噓2(207)
留言9則, 2人參與, 最新討論串2/2 (看更多)
1. If you can calulate the inverse function of the cumulative function easily, then randomly sample U1,...,Un from Uniform distribution and compute values of the inverse function at U1,...,Un. Or try gibb's sampling. 2. I think there should be a function for arctan in R. ※ 引述《laba1014 (嘿喲~加油!!!)》之銘言: : 不好意思有兩個R的問題要請教 : : 1. 請問若已知某分配的 p.d.f 型態 , : 譬如 f(x) =1/a(√π)^(-1) * gamma(m)/gamma(m-1) *(1+ x^2/a^2)^(-m) : 要怎麼寫程式才能從該分配中進行抽樣本的動作 ? : 2. 反函數要怎麼寫 ? : 譬如 已知x 但想知道 tan^(-1)(x) 應該怎麼做 ? : 希望版上前輩能不吝指教 <(_ _)> -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 208.125.17.166

07/08 16:39, , 1F
感謝 <(_ _)> 我去找找 gibb's sampling要怎麼用
07/08 16:39, 1F

07/09 08:13, , 2F
這跟Gibb無關吧?!
07/09 08:13, 2F

07/09 08:21, , 3F
如果m是正整數,那這是cauchy的變型,在R用rcauchy就行了.
07/09 08:21, 3F

07/09 08:23, , 4F
上面是錯的, m=2才是cauchy
07/09 08:23, 4F

07/09 08:24, , 5F
sorry, m=1才是
07/09 08:24, 5F

07/09 08:29, , 6F
所以你這是t分佈的變型. m應該是a的函數吧?
07/09 08:29, 6F

07/09 08:30, , 7F
在R中,用 rt(), 再套用rejection method就行了.
07/09 08:30, 7F

07/09 17:25, , 8F
C大眼光真是有夠利 我是想從pearson type 4 的分配取
07/09 17:25, 8F

07/09 17:26, , 9F
r.v 這分配的確包含上述所提的分配只是參數多不少
07/09 17:26, 9F
文章代碼(AID): #18SjfaiW (Statistics)
討論串 (同標題文章)
本文引述了以下文章的的內容:
完整討論串 (本文為第 2 之 2 篇):
文章代碼(AID): #18SjfaiW (Statistics)