use*_*129 10 r lm statistics-bootstrap
请我出去!我感谢任何帮助!谢谢!
重复进行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
Run Code Online (Sandbox Code Playgroud)
从上面的代码我可以通过随机抽样数据得到1000 betas0&1.而且我想做1000次以获得不同的测试版.除了replicate()之外我该怎么做呢?
Rom*_*rik 24
这更像是一个扩展的评论,我证明它replicate应该工作.这是CLT的一个例子.只需更换花括号之间的线条即可.
x <- replicate(1000, {
mm <- runif(10)
mean(mm)
})
hist(x)
Run Code Online (Sandbox Code Playgroud)
