我长期以来一直reshape2::melt
以一种相当非标准的方式使用它:我正在运行数字实验并得到一个矩阵。然后我融化它并产生一些图像。
受 和 之间的相似性的启发,我现在正试图在类矩阵的对象上实现相同的输出。到目前为止没有运气:reshape2
tidyr
library(reshape2)
library(tidyr)
set.seed(42)
mat <- matrix(runif(6), 3)
mat2 <- mat
colnames(mat2) <- letters[1:2]
rownames(mat2) <- letters[3:5]
melt(mat)
melt(mat2)
gather(mat) # fails
gather(mat2) # fails
请注意,如果它们存在,它melt
很聪明并保留。dimnames
我已经了解了它是如何工作的,因此我可以将以下函数添加到方法调度中:
gather.matrix <- function(mat) {
if (is.null(dimnames(mat))) {
grid <- expand.grid(seq.int(nrow(mat)), seq.int(ncol(mat)))
} else {
grid <- expand.grid(dimnames(mat))
}
cbind(grid, value = as.vector(mat))
}
all.equal(melt(mat),
gather.matrix(mat))
#[1] TRUE
all.equal(melt(mat2),
gather.matrix(mat2))
#[1] TRUE
但问题是,我可以强制gather
采取与我的情况相同的方式melt
吗?是否有任何参数组合可以在mat
和上产生所需的输出mat2
?