jiangbeilu 发表于 2018-4-7 19:06 
您好,我现在要对一个名为ddt的数据表进行处理,运行您编写的程序返回错误,请问哪里出错了呢
##构建一个矩阵 30*36
set.seed(201811)
origin <- matrix(sample(ddt[ , ],36*30,replace=T))
##编写1个每行取8个数,取100次的函数,返回这100次的平均值
sample8_100th <- function(x){
mat <- matrix(NA,nrow=8,ncol=100)
for(i in 1:100){
mat[,i] <- sample(x,8)
}
return(mean(mat))
}
##用apply对矩阵行进行计算,对每一行执行上述函数的操作,然后求个平均值,就可以得到30个值的平均值,也就是半径
## 圆的面积只与半径有关,所以只需要知道半径就好
mean(apply(origin,1,sample8_100th))