Код ниже для отбора проб с заменой.
random_nm<- function(n,m,v){
x <- sample(seq(length(v)),n*m,replace=TRUE)
rowMeans(matrix(v[x],nrow=n,ncol=m))
}
dat <- rnorm(100)
random_nm(5,20,dat)
[1] 0.05149423 -0.16239828 0.45485354 -0.10758954 -0.13597127
Если вы хотите сэмплы без замены, необходимо изменить первую строку функции random_rn.