来自多个不同大小的重采样的组平均值

Question

我需要通过重新采样不同的样本大小（例如：2、3、5）来计算每组（年份）的变量（FRic）的平均值。这只是我的桌子的一个例子：

   FRic <-c(6.644638923, 5.490853966, 4.948740617, 7.134710176, 5.167193087, 6.964663197, 4.157164504, 5.338576065, 1.463129916, 0.143467626, 4.962991887, 7.439127662, 7.127384236, 0.068020357, 2.733224992, 1.089417593)

   Year <- c(2002, 2002, 2002, 2002, 2002, 2003, 2003, 2003, 2003, 2003, 2003, 2004, 2004, 2004, 2004, 2004)

   Veg<- c(FRic, Year)

我每年都使用下面的脚本，这很烦人，因为我的表有9年700行。我想改进脚本，以获取每个样本量的重采样事件（N = 999）每年平均值的最终表。有什么更好的脚本可以使用吗？

    V2002<- Veg[Veg$Year == "2002",]
    class=c(V2002$FRic) 
    n = 3 
    N =999 
    statV2002 = numeric(N)
    for (i in 1:N){
    #bootstrap sample counterparts to observed samples are denoted with "B"
    classB= sample (class, n, replace=T)
    statV2002[i] = mean(classB)
    }

谢谢您的帮助。

最好

拍子