FC2ブログ

 日々の研究生活での出来事,考えたことなどのメモ。あるいはお知らせ。

  お知らせのナビゲーター   トップページ > スポンサー広告> 統計 > 多変量行列のブートストラップ  

スポンサーサイト

-- - --/-- [--] - --:--

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

多変量行列のブートストラップ

2012 - 04/21 [Sat] - 08:48


 1変量データをブートストラップリサンプリングする方法はよく見かけるのですが,多変量データをリサンプリングする方法はあまり説明されていません。なので,ここではそれをするためのRのコードをメモしておくことにします。基本は1変量の時と変わりません。ただ行全体でブートストラップをすればよいだけです。


x<- c(1:40)
(dat <- matrix(x, 10,4))
d.boot <- matrix(0,10,4)

n<-1
for(i in sample(10,replace=T)){
d.boot[n,] <- dat[i,]
n <- n+1
}
d.boot


 上記を実行すると,以下のようになります。(ただし,乱数による結果なので,実行するたびに並び順は異なります。)

リサンプリング前

[,1] [,2] [,3] [,4]
[1,] 1 11 21 31
[2,] 2 12 22 32
[3,] 3 13 23 33
[4,] 4 14 24 34
[5,] 5 15 25 35
[6,] 6 16 26 36
[7,] 7 17 27 37
[8,] 8 18 28 38
[9,] 9 19 29 39
[10,] 10 20 30 40


リサンプリング後

[,1] [,2] [,3] [,4]
[1,] 4 14 24 34
[2,] 3 13 23 33
[3,] 8 18 28 38
[4,] 10 20 30 40
[5,] 3 13 23 33
[6,] 2 12 22 32
[7,] 8 18 28 38
[8,] 10 20 30 40
[9,] 5 15 25 35
[10,] 8 18 28 38
スポンサーサイト

承認待ちコメント

このコメントは管理者の承認待ちです

コメントの投稿





管理者にだけ表示を許可する

トラックバック

http://tarotwork.blog40.fc2.com/tb.php/102-9141db34

 | HOME | 

プロフィール

管理者

 専門は心理統計学で構造方程式モデリングを中心に探索的分析法の研究をしています。

 ・著書:5冊
 ・翻訳:2冊
 ・査読論文:4本
 ・紀要論文:2本
 ・報告書等:2本
 ・国際発表:5回
 ・国内発表:12回

最近の記事+コメント

カテゴリ

全記事一覧表示

過去の記事

カウンタ

現在の閲覧者数

カレンダー

08 | 2018/09 | 10
- - - - - - 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 - - - - - -

ブログ内検索

リンク

上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。