2014-05-09 33 views
0

我試圖在河上運行,從線性迴歸的引導我到目前爲止的代碼是引導線性迴歸

hprice<-lm(dat[,1]~dat[,3]+dat[,4]+dat[,5]+dat[,6]) 
print (hprice) 
pricefunc<-function(data,ind) lm(data[ind,1]~data[ind,3]+data[ind,4]+data[ind,5]+data[ind,6]) 
hpboot<-boot(dat,pricefunc, 1000) 

這似乎並沒有工作。

我真的不明白統計的論點,我會說這是我出錯的地方。

感謝

+0

到底是什麼「不工作「 意思?您沒有提供足夠的代碼供其他人運行(即沒有樣本數據)。看看[這個問題](http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example)如何製作一個可重複的例子。 – MrFlick

+0

你是對的,我沒有提供足夠的代碼。我最大的問題是我無法確定boot()中的'statistic'參數。 我得到了一個答案,使它的工作,現在我只需要解決如何。謝謝 – Jamzy

回答

2

如果需要的係數估計你必須添加$coeflm功能

pricefunc<-function(data,ind) lm(data[ind,1]~data[ind,3]+data[ind,4]+data[ind,5]+data[ind,6])$coef 

然後你可以運行:

boot(dat,pricefunc, 1000) 
+0

好。我試過了,我認爲它工作。謝謝。 $ coef做什麼? – Jamzy

+0

命令'lm'給出了關於線性模型的信息列表,但'boot'函數只是想對模型的係數進行推理,因此您必須使用'$ coef'來提取關於線性模型的信息係數估計。 –