[問題] boot函數和optim函數搭配使用的問題

看板R_Language作者phil5566 (5566)時間7年前 (2017/01/10 00:14)推噓2(2推 0噓 30→)

留言32則, 5人參與討論串1/2 (看更多)

[問題類型]: 程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來) [軟體熟悉度]: 新手 [問題敘述]: 就是關於在R上用bootstrap這個統計方法時, 我會用boot這個function,而這個boot(), 括弧裡我只要給 1.T(假設T是data所形成的矩陣資料形式), 2.statistic(是自己定義的一個函數,用來計算原始dataT,bootstrapp後的資料,所要計算出統計量的數學式子) 3.R (要bootstrap的次數) 假設上述的第2點statistic的部分我是自己定義一個函數叫"theta.boot", 這個函數運算過程會使用到optim()這個function, 但optim()運算出的結果會有"$par"和"$convergence"的這兩部分我希望convergence=0的情形時,此時$par的解才是我想要的簡單說就是我希望算出的解"要收斂" ...................問題來了........................, 每bootstrap一次,當下那一次bootstrap過後的data是固定的,代入theta.boot運算, 由於給optim()設的初始值也是固定的,所以算出最後的解也是固定的, 然而這個解可能是收斂也可能不收斂,又我設定bootstrap1000次, 勢必這1000個解中有收斂也有沒收斂的,我希望要有1000個收斂解,要怎麼做到? 在程式範例中附上code 求各位大大高手先進解答~感激不盡~ [程式範例]: http://codepad.org/ongXCSxq [環境敘述]: win7, R x64 3.3.1, "boot"package [關鍵字]: boot,optim -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 36.224.8.250 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1483978448.A.5A4.html

→

andrew43

01/10 01:29, , 1^F

01/10 01:29, 1^F

→

andrew43

01/10 01:30, , 2^F

01/10 01:30, 2^F

→

andrew43

01/10 01:31, , 3^F

01/10 01:31, 3^F