「統計学関連なんでもあり」の過去ログ---034

★ Bootstrap データセットのアウトプット ★

7140.　Bootstrap データセットのアウトプット　MM　2005/07/08 (金) 09:32
├7143.　Re: Bootstrap データセットのアウトプット　ひの　2005/07/08 (金) 16:29
└7141.　Re: Bootstrap データセットのアウトプット　青木繁伸　2005/07/08 (金) 10:07
　└7148.　Re^3: Bootstrap データセットのアウトプット　MM　2005/07/08 (金) 17:14
　　└7150.　Re^4: Bootstrap データセットのアウトプット　青木繁伸　2005/07/08 (金) 18:27
　　　└7152.　Re^5: Bootstrap データセットのアウトプット　青木繁伸　2005/07/08 (金) 18:35
　　　　└7153.　Re^6: Bootstrap データセットのアウトプット　mint　2005/07/08 (金) 20:30
　　　　　└7157.　Re^7: できました！　MM　2005/07/09 (土) 08:50

7140.　Bootstrap データセットのアウトプット　MM　　2005/07/08 (金) 09:32

集団内の独立したサンプルから，bootstrap（replacementありのresampling）を使って別のデータセットを作りたいのですがどうしたらよいか教えてください！サンプルは独立していて，drawの確立は同じだと仮定しています。簡単そうなのに，どうしたらいいのか分かりません。お願いします。

　　　　　[このページのトップへ]

7143.　Re: Bootstrap データセットのアウトプット　ひの　　2005/07/08 (金) 16:29

> 集団内の独立したサンプルから，bootstrap（replacementありのresampling）を使って別のデータセットを作りたいのですがどうしたらよいか教えてください！サンプルは独立していて，drawの確立は同じだと仮定しています。簡単そうなのに，どうしたらいいのか分かりません。お願いします。

いろいろ方法はありますが，できる人にお願いしてやってもらうというのが一番楽でしょうか。

　　　　　[このページのトップへ]

7141.　Re: Bootstrap データセットのアウトプット　青木繁伸　　2005/07/08 (金) 10:07

どういう環境で行おうとしているのですか？
C プログラムを書くんですか？VBA ですか。
R が使えますか。
SPSS とか SAS を使うんですか。

　　　　　[このページのトップへ]

7148.　Re^3: Bootstrap データセットのアウトプット　MM　　2005/07/08 (金) 17:14

> > どういう環境で行おうとしているのですか？
> > C プログラムを書くんですか？VBA ですか。
> > R が使えますか。
> > SPSS とか SAS を使うんですか。

お返事ありがとうございます。SASを使っています。Proc Multtestではアウトプットがテストをべースにベクターに変換したものなのですが，単純にデータをそのままリサンプリングしたいのです。Bootのマクロからもデータセットのみのアウトプットができますか？Rは使ったことありません。Mathematicaはありますが，まったく使いこなせていません。SASでもできるようならなにかクルーをください！Rなどが必要なら・・・誰かにたのんで書いてもらうしかないですね・・・。

　　　　　[このページのトップへ]

7150.　Re^4: Bootstrap データセットのアウトプット　青木繁伸　　2005/07/08 (金) 18:27

便利なやり方や，既にそのような機能が組み込み済みと言うこともあるかもしれませんが，どんな場合にもプログラムさえできればできる方法を。ただし，私はSASのプログラムは書けませんので，自然言語で。。

データ行列を読み取り data に記憶しておく。
サンプルサイズ（data の行数）を n とする。
以下を n 回，繰り返す。
　[0, 1) の一様乱数 x を発生させる。
　int(x*n)+1 行目の data を書き出す。
（int は整数部を採る関数）
以上

　　　　　[このページのトップへ]

7152.　Re^5: Bootstrap データセットのアウトプット　青木繁伸　　2005/07/08 (金) 18:35

ついつい R で書いてしまうが，ものすごく簡単

data = read.table("input.dat") # データを読んで
x = data[floor(runif(nrow(data))*nrow(data))+1,] # 抽出して
write.data(x, "output.dat") # 書き出す

　　　　　[このページのトップへ]

7153.　Re^6: Bootstrap データセットのアウトプット　mint　　2005/07/08 (金) 20:30

SASを知りませんが，以下のURLは，参考にならないでしょうか。

行動科学のデータ解析
http://www.littera.waseda.ac.jp/faculty/stok/
のSAS pgm

　　　　　[このページのトップへ]

7157.　Re^7: できました！　MM　　2005/07/09 (土) 08:50

教えていただいたことを参考に，SASでプログラムを書くことができました。青木さん，mintさん，どうもありがとうございました。

　　　　　[このページのトップへ]