[英]Repeat the re-sampling function for 1000 times ? Using lapply?
請我出去! 我感謝任何幫助! 謝謝!
重復進行1000次重復采樣我遇到了麻煩。 我嘗試使用replicate()來做到這一點,但它不起作用。 有沒有其他方法可以做到這一點? 任何人都可以告訴我,如果這可能通過使用lapply完成? 以下是我的代碼:
#sampling 1000 betas0 & 1 (coefficients) from the data
get.beta=function(data,indices){
data=data[indices,] #let boot to select sample
lm.out=lm(y ~ x,data=data)
return(lm.out$coefficients)
}
n=nrow(data)
get.beta(data,1:n)
bootcoe=boot(data,get.beta,R=1000) #generate 1000 random samples
head(bootcoe$t) #look at the betas
從上面的代碼我可以通過隨機抽樣數據得到1000 betas0&1。 我想做1000次以獲得不同的測試版。 除了replicate()之外我該怎么做呢?
這更像是一個擴展注釋,我證明replicate
應該有效。 這是CLT的一個例子。 只需更換花括號之間的線條即可。
x <- replicate(1000, {
mm <- runif(10)
mean(mm)
})
hist(x)
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.