3
我有一個大陣,這是一個可複製的例子:應用在陣列第三維功能
cube <- array(c(1:10,5:15,-5:+5,1:18), c(4, 4, 3))
cube
, , 1
[,1] [,2] [,3] [,4]
[1,] 1 5 9 7
[2,] 2 6 10 8
[3,] 3 7 5 9
[4,] 4 8 6 10
, , 2
[,1] [,2] [,3] [,4]
[1,] 11 15 -2 2
[2,] 12 -5 -1 3
[3,] 13 -4 0 4
[4,] 14 -3 1 5
, , 3
[,1] [,2] [,3] [,4]
[1,] 1 5 9 13
[2,] 2 6 10 14
[3,] 3 7 11 15
[4,] 4 8 12 16
我想申請這個功能對於我的每個陣列(第三維)的「地板」:
m1 <- cube[,,1]
m2 <- cube[,,2]
m3 <- cube[,,3]
library(data.table)
m1[] <- frank(-m1, ties.method = "dense")
m2[] <- frank(-m2, ties.method = "dense")
m3[] <- frank(-m3, ties.method = "dense")
然後,因爲可複製的例子非常小。我可以做很容易通過結合矩陣的最後工作:
z <- array(c(m1 , m2, m3) , dim = c(4 , 4 , 3))
z
, , 1
[,1] [,2] [,3] [,4]
[1,] 10 6 2 4
[2,] 9 5 1 3
[3,] 8 4 6 2
[4,] 7 3 5 1
, , 2
[,1] [,2] [,3] [,4]
[1,] 5 1 13 9
[2,] 4 16 12 8
[3,] 3 15 11 7
[4,] 2 14 10 6
, , 3
[,1] [,2] [,3] [,4]
[1,] 16 12 8 4
[2,] 15 11 7 3
[3,] 14 10 6 2
[4,] 13 9 5 1
我需要一個快速的方式在我的大陣做到這一點。有什麼幫助嗎?謝謝。
我會'熔化'你的數組到data.table中並且處理它。但你也可以使用'apply'。 – Roland
適用方式:'數組(apply(-cube,3,frank,ties.method =「dense」),dim(cube))''但我不期望它快。我認爲羅蘭的想法可能會更好。 – Frank