作者lsshno1 (朝右邊鋒邁進)
看板Statistics
標題[程式] R 二元變數sampling
時間Mon Oct 19 22:01:07 2009
上面的問題感謝各位, 目前已經解決了
不過尚有一個問題需要處理
假設我現在的a為 seq(-10,10,length=1000) b為seq(0,40,length=1000)
我的posterior function是某個(a)%*%t(b)的function
所以我可以畫出 contour(posterior , a , b )的圖形
現在我希望做二元變數的sampling
老師提到的方式是先固定beta , 然後去算此beta下的alpha CDF , 接著
使用 -1
F (x) 服從 U(0,1)的方式抽樣,
最後把所有的alpha , beta 及 posterior的 density
不知道除了這樣固定一個變數 , 在抽另一個變數外
我還有什麼其他的處理辦法呢 感謝大家!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 60.250.86.244
→ clickhere:如果是貝氏,你的posterior已經是sample了. 10/20 01:42
→ clickhere:如果不是,inverse cdf是個gereral的觀念/出發點. 10/20 01:44