时间戳

首页 > 解决方案 > 计算 R 中多个数据帧的 z 分数

r - 计算 R 中多个数据帧的 z 分数

问题描述

我有十个行数和列数相等的数据框。它们看起来像这样：

df1 <- data.frame(geneID=c("AKT1","AKT2","AKT3","ALK",
                           "APC"),
                  CDKN2A=c(3490,9447,4368,908,204),
                  INPP4B=c(NA,9459,4395,1030,NA),
                  BCL2=c(NA,9480,4441,1209,NA),
                  IRS2=c(NA,NA,4639,1807,NA),
                  HRAS=c(3887,9600,4691,1936,1723))

df2 <- data.frame(geneID=c("AKT1","AKT2","AKT3","ALK",
                           "APC"),
                  CDKN2A=c(10892,17829,7156,1325,387),
                  INPP4B=c(NA,17840,7185,1474,NA),
                  BCL2=c(NA,17845,7196,1526,NA),
                  IRS2=c(NA,NA,12426,10244,NA),
                  HRAS=c(11152,17988,7545,2734,2423))

df3 <- data.frame(geneID=c("AKT1","AKT2","AKT3","ALK",
                           "APC"),
                  CDKN2A=c(11376,17103,8580,780,178),
                  INPP4B=c(NA,17318,9001,2829,NA),
                  BCL2=c(NA,17124,8621,1141,NA),
                  IRS2=c(NA,NA,8658,1397,NA),
                  HRAS=c(11454,17155,8683,1545,1345))

我想根据多个数据帧的均值和方差计算每个数据帧的 z 分数。z-score 应按如下方式计算：z-score=(x-mean(x))/sd(x)))。

我发现 plyr 的 ddply 函数可以完成这项工作，但解决方案是针对单个数据帧，而我有多个数据帧作为单独的文件，有 18214 行和 269 列。

我将不胜感激任何建议。

非常感谢您的帮助！奥尔哈

标签： rdataframestatistics

解决方案

这是一个选项，我们将数据集与bind_rows(from dplyr) 绑定在一起，然后按分组列分组并返回zscore转换后的数字列

library(dplyr)
bind_rows(df1, df2, df3, .id = 'grp') %>% 
     group_by(geneID) %>% 
     mutate(across(where(is.numeric), 
      ~(.- mean(., na.rm = TRUE))/sd(., na.rm = TRUE), .names = '{col}_zscore'))

注意：如果我们不需要新列，则删除该.names部分

如果我们需要在循环中执行此操作，而不绑定到单个 data.frame，可以循环list

library(purrr)
list(df1, df2, df3) %>% # // automatically => mget(ls('^df\\d+$'))
        map(~ .x %>%
          mutate(across(where(is.numeric), 
      ~(.- mean(., na.rm = TRUE))/sd(., na.rm = TRUE), .names = '{col}_zscore')))

推荐阅读