gpt4 book ai didi

r - 3D阵列->应用-> 3D阵列

转载 作者:行者123 更新时间:2023-12-02 22:53:47 25 4
gpt4 key购买 nike

似乎apply仅靠一个边距操作就不会重新组装3D阵列。考虑:

 arr <- array(
runif(2*4*3),
dim=c(2, 4, 3),
dimnames=list(a=paste0("a", 1:2), b=paste0("b", 1:4), c=paste0("c", 1:3))
)
# , , c = c1
#
# b
# a b1 b2 b3 b4
# a1 0.7321399 0.8851802 0.2469866 0.9307044
# a2 0.5896138 0.6183046 0.7732842 0.6652637
#
# , , c = c2
# b
# a b1 b2 b3 b4
# a1 0.5894680 0.7839048 0.3854357 0.56555024
# a2 0.6158995 0.6530224 0.8401427 0.04044974
#
# , , c = c3
# b
# a b1 b2 b3 b4
# a1 0.3500653 0.7052743 0.42487635 0.5689287
# a2 0.4097346 0.4527939 0.07192528 0.8638655


现在,制作一个4 x 4矩阵,以随机排列每个 arr[, , i]中的列,并使用 apply矩阵乘以 a*b中的每个 arr子矩阵以重新排列其列。重要的是,每个 apply迭代的结果都是一个矩阵

cols.shuf.mx <- matrix(c(0,1,0,0,1,0,0,0,0,0,0,1,0,0,1,0), ncol=4)
apply(arr, 3, `%*%`, cols.shuf.mx)
# c
# c1 c2 c3
# [1,] 0.8851802 0.78390483 0.70527431
# [2,] 0.6183046 0.65302236 0.45279387
# [3,] 0.7321399 0.58946800 0.35006532
# [4,] 0.5896138 0.61589947 0.40973463
# [5,] 0.9307044 0.56555024 0.56892870
# [6,] 0.6652637 0.04044974 0.86386552
# [7,] 0.2469866 0.38543569 0.42487635
# [8,] 0.7732842 0.84014275 0.07192528


而我希望结果是:

# , , c = c1
#
# a 1 2 3 4
# a1 0.8851802 0.7321399 0.9307044 0.2469866
# a2 0.6183046 0.5896138 0.6652637 0.7732842
#
# , , c = c2
#
# a 1 2 3 4
# a1 0.7839048 0.5894680 0.56555024 0.3854357
# a2 0.6530224 0.6158995 0.04044974 0.8401427
#
# , , c = c3
#
# a 1 2 3 4
# a1 0.7052743 0.3500653 0.5689287 0.42487635
# a2 0.4527939 0.4097346 0.8638655 0.07192528


我可以通过 plyr::aaply获得预期的结果:

aperm(aaply(arr, 3, `%*%`, cols.shuf.mx), c(2, 3, 1))


但我想知道是否有一种简单的基本方法来实现此结果(即我是否在这里缺少明显的东西以获得期望的结果)。

我意识到这里发生的是记录在案的( If each call to FUN returns a vector of length n, then apply returns an array of dimension c(n, dim(X)[MARGIN]) if n > 1),但是对我来说仍然很奇怪,如果函数返回具有尺寸的对象,则基本上将它们忽略。

最佳答案

这是一个不太理想的解决方案,需要预先知道函数结果矩阵的维数:

vapply(
1:dim(arr)[3],
function(x, y) arr[,,x] %*% y,
FUN.VALUE=arr[,,1],
y=cols.shuf.mx
)

关于r - 3D阵列->应用-> 3D阵列,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/21196818/

25 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com