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

No.03051　Rのbootstrap　　【エミリ】　2007/03/23(Fri) 11:44

Rのbootについてお伺いしたいです。
目的：回帰分析のサンプルの数が不足ですので，リサンプリングで回帰分析をしたいです。
Rのbootのドキュメントの例に基づいて，以下のソースを書きました。

nuke <- nuclear[,c(1,2,5,7,8,10,11)]
nuke.lm <- lm(log(cost)~date+log(cap)+ne+ ct+log(cum.n)+pt, data=nuke)
nuke.lm.diag=ls.diag(nuke.lm)
nuke.lm.diag.res=nuke.lm.diag$std.res
resid=nuke.lm.diag.res
fit=fitted(nuke.lm)

myfun=function(dat,inds,i.pred){
	assign(".inds",inds,envir=.GlobalEnv)
	lm.b=lm(log(cost)[.inds]~date+log(cap)+ne+ ct+log(cum.n)+pt, data=dat)
	remove(".inds",envir=.GlobalEnv)
	c(coef(lm.b))
}
myboot=boot(nuke,myfun,R=999,m=1)

しかし，結果の偏差はすごく大きいです。myfunの書き方は正しいでしょうか？
間違ったら，正しいソースを教えていただけないでしょうか？

以上，よろしくお願いいたします。

No.03053　Re: Rのbootstrap　　【青木繁伸】　2007/03/23(Fri) 18:41

自分で書いたプログラムが正しいかどうか知る一番の方法は，答えの分かっている問題を解いてみることです。