いずれかまたは特定の折り目がある行をフィルタリングする

Nov 28 2020

私はまだRを学習しています。このデータセットがあり、5つの列があり、最初の列はtracking_idで、次の4つの列には4つのグループの値があります。

列( "CD44low_rep")と比較して8倍高いまたは4倍低い最後の3つの列( "CD44hi_CD69low_rep"、 "CD44hi_CD69hi_CD103low_rep"、 "CD44hi_CD69hi_CD103hi_rep")を比較した後、行をフィルター処理したいと思います。

これを達成する方法は?

回答

1 akrun Nov 28 2020 at 05:10

'CD44low_rep'列に8と4を掛けてから、とを使用して対象の列>=<=それぞれ比較し、TRUE値の行方向の合計をで取得し、rowSums3に等しいかどうか(つまり、比較される列の数)を確認します。&両方の比較から単一の論理ベクトルを返し、それを使用して行をサブセット化します

nm1 <- c("CD44hi_CD69low_rep",  "CD44hi_CD69hi_CD103low_rep", 
         "CD44hi_CD69hi_CD103hi_rep")
i1 <- (rowSums(df1[nm1]  >= (df1$CD44low_rep * 8)) == 3) & (rowSums(df1[nm1] <= (df1$CD44low_rep * 4)) == 3)

df1[i1,]
# A tibble: 798 x 5
#   tracking_id   CD44low_rep CD44hi_CD69low_rep CD44hi_CD69hi_CD103low_rep CD44hi_CD69hi_CD103hi_rep
#   <chr>               <dbl>              <dbl>                      <dbl>                     <dbl>
# 1 1600014C23Rik           0                  0                          0                         0
# 2 1600019K03Rik           0                  0                          0                         0
# 3 1700006E09Rik           0                  0                          0                         0
# 4 1700010M22Rik           0                  0                          0                         0
# 5 1700011A15Rik           0                  0                          0                         0
# 6 1700016P04Rik           0                  0                          0                         0
# 7 1700018G05Rik           0                  0                          0                         0
# 8 1700019A02Rik           0                  0                          0                         0
# 9 1700024B18Rik           0                  0                          0                         0
#10 1700024G13Rik           0                  0                          0                         0
# … with 788 more rows

またはdplyr、を使用して、対象の列をループすることで同じ式を使用acrossします(デフォルトでは、all列のチェックを行います)

library(dplyr)
df1 %>%
     filter(across(contains('hi'), ~ (. >= (CD44low_rep * 8)) & 
                (. <= (CD44low_rep * 4))))

-出力

# A tibble: 798 x 5
#   tracking_id   CD44low_rep CD44hi_CD69low_rep CD44hi_CD69hi_CD103low_rep CD44hi_CD69hi_CD103hi_rep
#   <chr>               <dbl>              <dbl>                      <dbl>                     <dbl>
# 1 1600014C23Rik           0                  0                          0                         0
# 2 1600019K03Rik           0                  0                          0                         0
# 3 1700006E09Rik           0                  0                          0                         0
# 4 1700010M22Rik           0                  0                          0                         0
# 5 1700011A15Rik           0                  0                          0                         0
# 6 1700016P04Rik           0                  0                          0                         0
# 7 1700018G05Rik           0                  0                          0                         0
# 8 1700019A02Rik           0                  0                          0                         0
# 9 1700024B18Rik           0                  0                          0                         0
#10 1700024G13Rik           0                  0                          0                         0
# … with 788 more rows