提问人:Steve 提问时间:6/20/2011 最后编辑:zx8754Steve 更新时间:10/6/2022 访问量:11038
将两个 DataFrame 逐个元素粘贴在一起
Paste together two dataframes element by element
问:
我需要逐个元素粘贴两个数据帧的内容,以便输入到另一个程序。我有一个均值数据框和一个均值标准误差数据框。
我尝试使用 R paste() 函数,但它似乎无法处理数据帧。当使用向量时,它似乎将第一个向量的所有元素连接成一个字符串,将第二个向量的所有元素连接成一个单独的字符串。相反,我需要将两个数据帧中的每个倒数元素连接在一起。
关于如何处理这个问题有什么建议吗?我包含了虚拟输入数据(datMean 和 datSE)和我想要的输出 (datNew)。我的真实数据框大小约为 10 行 x 150 列。
# means and SEM
datMean <- data.frame(a=rnorm(10, 3), b=rnorm(10, 3), d=rnorm(10, 3))
datSE <- data.frame(a=rnorm(10, 3)/100, b=rnorm(10, 3)/100, d=rnorm(10, 3)/100)
# what the output should look like
# i've chosen some arbitrary values here, and show only the first row.
datNew <- data.frame(a="2.889-2.926", b="1.342-1.389", d="2.569-2.576")
这个想法是让 datNew 中的每个元素都是一个由 'mean - se' 和 'mean + se' 组成的范围,用破折号 '-' 分隔。paste() 函数可以对一个元素执行此操作,如何在整个数据帧上执行此操作?
paste(datMean[1,1] - datSE[1,1], datMean[1,1] + datSE[1,1], sep="-")
编辑1:看着一些答案,我意识到我在问题中遗漏了一些重要的信息。原始数据框的每一行都有名称,我需要用这些名称重新构建最终数据框。例如:
rownames(datMean) <- LETTERS[1:10]
rownames(datSE) <- LETTERS[1:10]
我需要 datNew 最终再次拥有这 10 个行名。对于某些使用 melt() 的解决方案,这可能是有问题的。
答:
您可以一次对每一行执行此操作,但您要应用于两个 data.frame 之间的成对列。由于您每次都要执行特定的粘贴作业,因此请定义函数:
pfun <- function(x, y) paste(x - y, x + y, sep = "-")
然后使用以下函数构造新的 Data.Frame:
datNew <- data.frame(a = pfun(datMean$a, datSE$a), b = pfun(datMean$b, datSE$b), d = pfun(datMean$d, datSE$d))
会有更简单的方法来应用它,但也许这有助于你更好地理解。您可以将整个列传递到粘贴,但不能传递整个 data.frames。
使用循环匹配结果中的所有列,而无需单独指定它们。
首先创建一个列表来存储所有列,我们将转换为具有正确列名的 data.frame。
datNew <- vector("list", ncol(datMean))
命名确实假定列号、名称和顺序在两个输入 data.frame 之间完全匹配。
names(datNew) <- names(datMean)
for (i in 1:ncol(datMean)) {
datNew[[i]] <- pfun(datMean[[i]], datSE[[i]])
}
转换为 data.frame:
datNew <- as.data.frame(datNew)
评论
这是我如何理解你的问题。我使用 将均值和 SE 的数据从多列融化到一列。reshape2::melt
library(reshape2)
datMean <- melt(datMean)$value
datSE <- melt(datSE)$value
dat <- cbind(datMean, datSE)
apply(X = dat, MARGIN = 1, FUN = function(x) {
paste(x[1] - x[2], x[1] + x[2], sep = " - ")
})
结果
[1] "3.03886802467251 - 3.08551547263516"
[2] "3.01803172559258 - 3.05247871975711"
[3] "3.4609230722069 - 3.56097173966387"
[4] "1.35368243309618 - 1.45548512578821"
[5] "2.39936853846605 - 2.47570756724791"
[6] "3.21849170272184 - 3.29653660329785"
编辑
此解决方案尊重您的原始数据维度。我所做的是制作一个 3D 数组,并在保持三维 () 恒定的情况下一次处理每个单元格。[x,y, 1:2]
dat <- array(c(datMean, datSE), dim = c(10, 3, 2))
datNEW <- matrix(rep(NA, nrow(dat)*ncol(dat)), ncol = ncol(dat))
for (column in seq(ncol(dat))) {
cls <- rep(NA, nrow(dat))
for (rows in seq(nrow(dat))) {
tmp <- dat[rows, column, 1:2]
cls[rows] <- paste(tmp[1] - tmp[2], tmp[1] + tmp[2], sep = " - ")
}
datNEW[, column] <- cls
}
评论
下面是一种无需手动指定每列即可执行此操作的方法。首先,我们制作数据并使用包将它们放入数组中,四舍五入到 3,因为这样看起来更好:abind
datMean <- data.frame(a=rnorm(10, 3), b=rnorm(10, 3), d=rnorm(10, 3))
datSE <- data.frame(a=rnorm(10, 3)/100, b=rnorm(10, 3)/100, d=rnorm(10, 3)/100)
library(abind)
datArray <- round(abind(datMean,datSE,along=3),3)
然后我们可以将该函数应用于该数组的每个元素和列:paste
apply(datArray,1:2,function(x)paste(x[1]-x[2],"-",x[1]+x[2]))
a b d
[1,] "3.537 - 3.581" "3.358 - 3.436" "3.282 - 3.312"
[2,] "2.452 - 2.516" "1.372 - 1.44" "3.041 - 3.127"
[3,] "3.017 - 3.101" "3.14 - 3.228" "5.238 - 5.258"
[4,] "3.397 - 3.451" "2.783 - 2.839" "3.381 - 3.405"
[5,] "1.918 - 1.988" "2.978 - 3.02" "3.44 - 3.504"
[6,] "4.01 - 4.078" "3.014 - 3.068" "1.914 - 1.954"
[7,] "3.475 - 3.517" "2.117 - 2.159" "1.871 - 1.929"
[8,] "2.551 - 2.619" "3.907 - 3.975" "1.588 - 1.614"
[9,] "1.707 - 1.765" "2.63 - 2.678" "1.316 - 1.348"
[10,] "4.051 - 4.103" "3.532 - 3.628" "3.235 - 3.287"
评论
如果先转换为矩阵,则完全无需应用或循环即可完成。
MdatMean <- as.matrix(datMean)
MdatSE <- as.matrix(datSE)
matrix( paste(MdatMean - MdatSE, MdatMean + MdatSE, sep="-"),
nrow=nrow(MdatMean), dimnames=dimnames(MdatMean) )
您也可以考虑更好的格式。formatC
lo <- formatC(MdatMean - MdatSE, format="f", digits=3)
hi <- formatC(MdatMean + MdatSE, format="f", digits=3)
matrix( paste(lo, hi, sep="-"),
nrow=nrow(MdatMean), dimnames=dimnames(MdatMean) )
如果你想要最后,只需将最后一行括起来。data.frame
as.data.frame
评论
使用 mapply 粘贴和 cbind 保留行名:
x <- cbind(
datMean[, 0],
mapply(paste, round(datMean - datSE, 3), round(datMean + datSE, 3), sep = " - "))
x
# a b d
# A 3.268 - 3.321 5.226 - 5.308 2.3 - 2.358
# B 3.795 - 3.874 1.772 - 1.833 2.265 - 2.335
# C 1.305 - 1.346 1.238 - 1.291 2.812 - 2.874
# D 1.957 - 2.041 3.016 - 3.057 2.402 - 2.473
# E 4.73 - 4.786 2.909 - 2.963 2.245 - 2.297
# F 3.511 - 3.554 3.547 - 3.603 2.316 - 2.374
# G 3.601 - 3.689 3.073 - 3.144 3.145 - 3.215
# H 2.056 - 2.118 2.597 - 2.69 2.58 - 2.627
# I 1.802 - 1.835 2.794 - 2.895 2.452 - 2.5
# J 2.399 - 2.461 1.807 - 1.844 3.199 - 3.254
class(x)
# [1] "data.frame"
identical(rownames(datMean), rownames(x))
# [1] TRUE
上一个:在 R 中按两列对数据帧进行排序
评论