创建一个二维数组

Question

我想创建一个100,000乘100,000的二维数组，其中包含Cp值。所以我得到Cp [1，]，Cp [2，]，Cp [3，]等等 然后找到Cp [1，]的所有均值并将其存储为meanCp [1]，Cp [2，]的平均值并存储为meanCp [2]。 有没有比这里更好更快的方式呢？

K=0.5
x=numeric()
stdCp=numeric()
meanCp=numeric()
u0=0
Cp <- matrix(0,100000,100000)
for (j in 1:100000) {
  x=rnorm(100000,0,1)
  for (i in 1:100000) {    
    if (j==1){
     Cp[j,i]= max(0,x[i]-(u0+K))
    } else {
       Cp[j,i]= max(0,x[i]-(u0+K)+Cp[j-1,i])
    }
  }
 meanCp[j]=mean(Cp[j,])
 stdCp[j]=sd(Cp[j,])
}
meanCp

所以例如meanCp [1]给出的值介于0.19和0.20之间

Answer 1

您可以参考 colMeans 功能。同一帮助html中的其他一些功能也很有用。