應用在陣列第三維功能

我有一個大陣，這是一個可複製的例子：應用在陣列第三維功能

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

我需要一個快速的方式在我的大陣做到這一點。有什麼幫助嗎？謝謝。

來源

2016-09-01 Albert

我會'熔化'你的數組到data.table中並且處理它。但你也可以使用'apply'。 – Roland

適用方式：'數組（apply（-cube，3，frank，ties.method =「dense」），dim（cube））''但我不期望它快。我認爲羅蘭的想法可能會更好。 – Frank

功能apply就是爲此而設計的。第二個參數MARGIN指定應用該函數的維度。

res <- apply(-cube, 3, frank, ties.method = "dense") 
## reshape: 
array(res, dim=dim(cube))

來源

2016-09-01 15:25:16

謝謝@弗蘭克指出。我編輯了我的答案。 –

應用在陣列第三維功能

回答

相關問題