Question

我有一个data.frame data：

Sector Var1 Var2 Var3
Abcd    1    1    1
Efgh    4    5    6
Ijkl    7    8    9

我想创建一个新的data.frame，它包含每个变量Sector的每个唯一配对之间的差异。所需结果的一个例子如下：

# result
Sector1 Sector2 Dif_Var1 Dif_Var2 Dif_Var3
Abcd    Efgh        3        4        5
Abcd    Ijkl        6        7        8
Efgh    Ijkl        3        3        3

我可以找到Sector的每个唯一对（代码如下），但我不确定这是否是继续进行的方式。下面这个最小的伪代码看起来是一个合适且可能成功的路线吗？

# pseudo code
result <- which(Sector1 == Sector) get Var1,2,3 values - which(Sector2 == Sector) get Var1,2,3 values

是否有一个r函数或包可以帮助达到预期的结果？

# unique pairings
# http://stackoverflow.com/q/23024059/1670053
df <- structure(list(Sector = structure(1:3, .Label = c("Abcd", "Efgh", 
"Ijkl"), class = "factor"), Var1 = c(1L, 4L, 7L), Var2 = c(1L, 
5L, 8L), Var3 = c(1L, 6L, 9L)), .Names = c("Sector", "Var1", 
"Var2", "Var3"), class = "data.frame", row.names = c(NA, -3L))
df.u <- expand.grid(df$Sector,df$Sector)
df.u2 <- as.data.frame(unique(t(apply(df.u, 1, function(x) sort(x)))))
data <- subset(df.u2, ! df.u2$V1 == df.u2$V2)

Answer 1

您可以使用您创建的V1和V2向量来定位df

df[data$V2, -1] - df[data$V1, -1]
#    Var1 Var2 Var3
#2      3    4    5
#3      6    7    8
#3.1    3    3    3

保留名称

cbind(data$V1, data$V2, df[data$V2, -1] - df[data$V1, -1])
#    data$V1 data$V2 Var1 Var2 Var3
#2      Abcd    Efgh    3    4    5
#3      Abcd    Ijkl    6    7    8
#3.1    Efgh    Ijkl    3    3    3

Answer 2

使用sqldf包：

library(sqldf)
sqldf("select A.Sector Sector1,
              B.Sector Sector2,
              B.Var1 - A.Var1 Var1, 
              B.Var2 - A.Var2 Var2, 
              B.Var3 - A.Var3 Var3
              from df A join df B
              on A.Sector < B.Sector")

，并提供：

  Sector1 Sector2 Var1 Var2 Var3
1    Abcd    Efgh    3    4    5
2    Abcd    Ijkl    6    7    8
3    Efgh    Ijkl    3    3    3

这也可以写成：

 nms <- names(df)[-1]
 sel <- toString( sprintf('B.%s - A.%s %s', nms, nms, nms) )
 fn$sqldf("select A.Sector Sector1, B.Sector Sector2, $sel 
           from df A join df B 
           on A.Sector < B.Sector")

已修订已修复并添加了变体。

计算data.frame变量的唯一配对之间的差异

2 个答案: