計算矩陣中每個單元的迭代後方差

Question

如何計算循環後矩陣中每個單元的方差？所以，現在如果我做

m = matrix(0,10,10) 
    n = 100 
    v = 1 

    rad2 <- function(matrix, repeats, v) { 
    idx <- sample(length(matrix), repeats, replace = TRUE) # indices 
    flip <- sample(c(-1, 1), repeats, replace = TRUE) # subtract or add 
    newVal <- aggregate(v * flip ~ idx, FUN = sum) # calculate new values for indices 
    matrix[newVal[[1]]] <- matrix[newVal[[1]]] + newVal[[2]] # add new values 
    variance = M2/(n-1) 
    return(matrix) 
} 

：

我有這個迄今爲止

rad2(m, n, v) 
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] 
[1,] 0 0 0 1 0 -2 1 -2 1  0 
[2,] 1 1 0 0 0 0 -1 2 0 -2 
[3,] -1 1 -1 1 0 0 0 1 -1 -1 
[4,] 0 0 1 0 1 0 1 1 1 -1 
[5,] 0 1 1 -2 0 0 1 0 -1 -1 
[6,] -2 -3 0 1 1 0 0 1 0 -2 
[7,] 0 0 0 1 0 2 -1 0 -1  1 
[8,] 2 0 -1 0 -1 -1 -1 0 -1  0 
[9,] 0 0 1 1 -1 1 1 0 0  1 
[10,] 0 -3 1 0 -2 0 0 -2 -1  0 

我想經過100次迭代這個功能來計算每個單元的方差。輸出可以在一個表中或一個向量中。最後應該有100個值。我怎樣才能做到這一點？

編輯：

如果我這樣做，而不是：

n=10 
for (i in 1:n) { 
tmp <- rad(m) 
m <- tmp 
outv <- unlist(sapply(m, function(x) var(m))) 
finalv <- outv 
} 

我得到finalv輸出。但是，如何在矩陣中的每個循環之後替換每個單元格的方差值，而不是一遍又一遍地重寫？

Answer 1

我相信這應該做你想做的。只是一些簡單的功能。

newdata=matrix(NA,100,100) 
output=NULL 
for (i in 1:100){newdata[i,]=as.numeric(rad2(m, n, v))} 
for (n in 1:100){output[n]=mean(newdata[,n])} 
matrix(output,10,10) 

一種更普遍的方式

newdata=matrix(NA,length(m),length(m)) 
output=NULL 
for (i in 1:length(m)){newdata[i,]=as.numeric(rad2(m, n, v))} 
for (b in 1:length(m)){output[b]=mean(newdata[,b])} 
matrix(output,dim(m)[1],dim(m)[2])