作者pcshdan (我的人生就像冒險一樣)
看板Statistics
標題[問題] 關於[R]的問題
時間Sun Jul 20 14:51:07 2008
我是R的新手
所以內容敘述不清楚請見諒
我想請問關於bootstrap resampling的問題
今天我有一組data 裡面有1000筆資料
資料的尺度是intreval
請問我要如何抽出新的100組data (抽後放回)
讓每一組data裡面同樣有1000筆資料的呢
又一組裡面如果只要有500筆資料又要怎麼寫呢
我看R的PDF檔查了很久都沒查到所以才來問問題
我有看到指令boot()和boot.ci()
但裡面似乎都要寫一個函數
我不懂到底是要如何使用?
試了很久都不行
因為我的最終目的是要形成一個新的矩陣
也就是那100組data我最後會寫一個函數使每一組data都產生一個係數
然後做此係數的信賴區間
請問這樣是不是使用loop會比較好
我是第一次學統計軟體
問題敘述有點長
麻煩大家了
謝謝^^
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 219.91.78.147
※ 編輯: pcshdan 來自: 219.91.78.147 (07/20 14:55)
※ 編輯: pcshdan 來自: 219.91.78.147 (07/20 14:58)
1F：推 mrk:囧，有這個函數喔，我直接用runif取。Orz 07/20 18:02
2F：推 tew:用google搜尋 bootstrap R 前面幾個網頁就有範例 07/20 18:49
3F：推 tew:我常這樣搜尋 學會自己要的東西 07/20 18:50
4F：推 waltervic:我也在學R= = 07/20 23:01
5F：→ clickhere:bootstrap要resampling 1000筆, 500是不對的. 07/21 02:23
6F：→ clickhere:把你寫好的產生係數的函數丟給boot(),boot.ci()就行了. 07/21 02:24
7F：→ clickhere:sample的問題, boot()會處理. 07/21 02:25
8F：→ clickhere:用sample(..., replace=T)可以直接rexampling. 07/21 02:26
9F：→ pcshdan:感恩喔~謝謝各位的回答^^但我把產生係數的函數丟進boot() 07/21 18:08
10F：→ pcshdan:它一直說找不到引數"m" 不知道是哪步出問題@@ 07/21 18:09
11F：→ pcshdan:似乎在寫boot()裡面的函數跟平常的寫法不一樣 07/21 23:13