小编use*_*129的帖子

重复重采样功能1000次?使用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
Run Code Online (Sandbox Code Playgroud)

从上面的代码我可以通过随机抽样数据得到1000 betas0&1.而且我想做1000次以获得不同的测试版.除了replicate()之外我该怎么做呢?

r lm statistics-bootstrap

10
推荐指数
1
解决办法
4万
查看次数

使用tapply作为子集数据组

我有一组汽车销售数据,我将我的数据分组到不同的组中,如下所示:

汽车品牌和销售年份.

toyota=subset(car, brand=="Toyota")
toyota.yr = cut(toyota$date, "year")
honda=subset(car, brand=="Honda")
honda.yr = cut(honda$date, "year")
Run Code Online (Sandbox Code Playgroud)

等等

所以现在我有6个小组的汽车品牌,然后我用tapply来获得每个品牌的销售均值:

tapply(toyota$price, toyota.yr, mean, na.rm=TRUE)
Run Code Online (Sandbox Code Playgroud)

我想对所有6个子组执行此操作,无论如何我可以同时执行此操作而不是键入tapply函数6次?

我感谢任何帮助,谢谢!

r

2
推荐指数
1
解决办法
3077
查看次数

标签 统计

r ×2

lm ×1

statistics-bootstrap ×1