の任意の数に分割マトリックスは、私はほぼ同じ大きさの、このマトリックス内にほぼ正方形の地域のためのいくつかの統計を計算したい私は値R:ブロック
set.seed(1)
A <- matrix(runif(25),ncol=5)
のマトリックスを有していると言います。出力のこれらの種類のいずれかが行うだろう:
N1 <- matrix(c(rep(c("A","A","B","B","B"),2),rep(c("C","C","D","D","D"),3)),ncol=5)
N2 <- matrix(c(rep(c("A","A","A","B","B"),3),rep(c("C","C","D","D","D"),2)),ncol=5)
N1
[,1] [,2] [,3] [,4] [,5]
[1,] "A" "A" "C" "C" "C"
[2,] "A" "A" "C" "C" "C"
[3,] "B" "B" "D" "D" "D"
[4,] "B" "B" "D" "D" "D"
[5,] "B" "B" "D" "D" "D"
N2
[,1] [,2] [,3] [,4] [,5]
[1,] "A" "A" "A" "C" "C"
[2,] "A" "A" "A" "C" "C"
[3,] "A" "A" "A" "D" "D"
[4,] "B" "B" "B" "D" "D"
[5,] "B" "B" "B" "D" "D"
他の近似を、私はいつも行列を回転させることができるので、もOKです。それから私はこのように、tapply()
を使用して統計を計算するためにこれらの近所行列を使用することができます。
tapply(A,N1,mean)
A B C D
0.6201744 0.5057402 0.4574495 0.5594227
私が欲しい私にN1
などブロック状の地域の任意の数の任意の次元の行列を作ることができる機能ですN2
。私は、そのような関数が、希望のブロック数が偶数でない状況をどのように扱うのか把握しようとするのは苦労している。 N1
とN2
は4つの地域がありますが、私はこのようないくつかの出力何かのために5をしたかったと言う:
N3 <- matrix(c("A","A","B","B","B","A","A","C","C","C","D","D","C","C","C",
"D","D","E","E","E","D","D","E","E","E"),ncol=5)
[,1] [,2] [,3] [,4] [,5]
[1,] "A" "A" "D" "D" "D"
[2,] "A" "A" "D" "D" "D"
[3,] "B" "C" "C" "E" "E"
[4,] "B" "C" "C" "E" "E"
[5,] "B" "C" "C" "E" "E"
は、誰もが、既存の機能を知っていスプリットのこの種を行う、または1つを作成する方法上の任意のアイデアを持つことができます?ありがとうございました!
[[編集]] 私の最後の関数、考慮にヴィンセントのアドバイス取っ:youreeの中央付近に
DecideBLocks <- function(A,nhoods){
nc <- ncol(A)
nr <- nrow(A)
nhood_side <- floor(sqrt((nc*nr)/nhoods))
Neighborhoods <- matrix(paste(ceiling(col(A)/nhood_side), ceiling(row(A)/nhood_side), sep="-"), nc=ncol(A))
nhoods.out <- length(unique(c(Neighborhoods)))
if (nhoods.out != nhoods){
cat(nhoods.out,"neighborhoods created.\nThese were on average",nhood_side,"by",nhood_side,"cells\nit's a different number than that stated the function tries to round things to square neighborhoods\n")
}
return(Neighborhoods)
}
A <- matrix(rnorm(120),12)
B <- DecideBLocks(A,13)
ヴィンセントに感謝します。私は、私が必要とするものを、 –