我有一個 R 資料框,如下例所示。我希望計算觀察/行(所有組合)之間列值的差異。
my_df <- tibble(a=runif(5), b=runif(5), c=runif(5))
> my_df
# A tibble: 5 x 3
a b c
<dbl> <dbl> <dbl>
1 0.0513 0.267 0.846
2 0.614 0.683 0.937
3 0.230 0.700 0.0651
4 0.671 0.110 0.901
5 0.424 0.520 0.817
我已經嘗試了下面的代碼,它只給出了后續行之間的差異;我想要所有組合:row2 - row1; 第 3 行 - 第 1 行;行4 - 行1,行5 - 行1,行3 - 行2,行4 - 行2,依此類推...
另外,我寫的代碼對我來說似乎不是最好的(!),雖然它輸出了我想要的結果,但不是所有可能的組合!
my_diff <- as.data.frame(diff(as.matrix(my_df)))
> my_diff
a b c
1 0.5623574 0.41522579 0.09165630
2 -0.3837289 0.01755953 -0.87209740
3 0.4407068 -0.58982681 0.83540813
4 -0.2463205 0.40943495 -0.08358985
如果有人可以使用 R 解決我的問題,我很感激,如果可能的話,使用整潔的詩句選項。
謝謝。
uj5u.com熱心網友回復:
更新:一個整潔友好的解決方案:
library(tidyverse)
set.seed(1)
my_df <- tibble(a=runif(5), b=runif(5), c=runif(5))
給出:
# A tibble: 5 x 3
a b c
<dbl> <dbl> <dbl>
1 0.266 0.898 0.206
2 0.372 0.945 0.177
3 0.573 0.661 0.687
4 0.908 0.629 0.384
5 0.202 0.0618 0.770
從那里:
my_df %>%
mutate(ID = row_number()) %>%
slice(as.numeric(t(combn(1:nrow(.), 2)))) %>%
mutate(group = rep(1:(n()/2), 2)) %>%
group_by(group) %>%
summarize(comparison = paste0(ID[2], "-", ID[1]),
across(c(a, b, c), ~ .[2] - .[1])) %>%
select(-group)
這使:
# A tibble: 10 x 4
comparison a b c
<chr> <dbl> <dbl> <dbl>
1 2-1 0.107 0.0463 -0.0294
2 3-1 0.307 -0.238 0.481
3 4-1 0.643 -0.269 0.178
4 5-1 -0.0638 -0.837 0.564
5 3-2 0.201 -0.284 0.510
6 4-2 0.536 -0.316 0.208
7 5-2 -0.170 -0.883 0.593
8 4-3 0.335 -0.0317 -0.303
9 5-3 -0.371 -0.599 0.0828
10 5-4 -0.707 -0.567 0.386
uj5u.com熱心網友回復:
請讓我知道這是否是您的預期。
my_df <- tibble(a=runif(5), b=runif(5), c=runif(5))
# Generating the sequence to calculate the combinations
seq1 <- seq(1,nrow(my_df))
seq2 <- seq1
# Generating the Combinations
Combinations <- expand.grid(seq1, seq2)
# Removing the dupilicate Combinations
Combinations <- Combinations[which(Combinations$Var2 < Combinations$Var1),]
# Performing the subtraction
result <- my_df[Combinations$Var1,] - my_df[Combinations$Var2,]
根據評論更新:
result <- expand.grid(seq1,seq1)%>%
filter(Var1 > Var2)%>%
mutate(my_df[Var1,] - my_df[Var2,])
轉載請註明出處,本文鏈接:https://www.uj5u.com/yidong/334578.html
