当前位置: 首页 > 知识库问答 >
问题:

通过切片行名-R将2列数据帧转换为矩阵

酆鸿哲
2023-03-14

我在R中有两个数据帧,我想把它们转换成矩阵。这是我手上的示例数据

一个data.framedat有2列,一列有名称,第二列有一些分数。

> head(dat, n=20)
                                                     V1           V2
1                                      4star_Active_TSS  10.99561503
2                               4star_Bivalent_Enhancer   0.42425920
3                             4star_Bivalent_Poised_TSS   0.31163730
4                                       4star_Enhancers 139.64713405
5                             4star_Flanking_Active_TSS  23.11961717
6                       4star_Flanking_Bivalent_TSS_Enh   0.17654506
7                                 4star_Genic_enhancers  44.42394542
8                                 4star_Heterochromatin -34.61099049
9                                   4star_Quiescent_Low -28.54240987
10                             4star_Repressed_PolyComb  -0.42096698
11                           4star_Strong_transcription  12.70895605
12                        4star_Transcr_at_gene_5_and_3   0.00000000
13                        4star_Weak_Repressed_PolyComb   0.08902141
14                             4star_Weak_transcription  19.28859369
15                          4star_ZNF_genes_and_repeats  -1.58340662
16          A549_EtOH_0.02pct_Lung_Carcinoma_Active_TSS  14.01552989
17   A549_EtOH_0.02pct_Lung_Carcinoma_Bivalent_Enhancer   1.00466761
18 A549_EtOH_0.02pct_Lung_Carcinoma_Bivalent_Poised_TSS   1.20607773
19           A549_EtOH_0.02pct_Lung_Carcinoma_Enhancers  63.36004048
20 A549_EtOH_0.02pct_Lung_Carcinoma_Flanking_Active_TSS  48.27400816

其他数据帧为一列,包含上述数据帧的row.names的半名称。

> states=read.delim("states.txt", header=FALSE)
> states
                          V1
1                 Active_TSS
2          Bivalent_Enhancer
3        Bivalent_Poised_TSS
4                  Enhancers
5        Flanking_Active_TSS
6  Flanking_Bivalent_TSS_Enh
7            Genic_enhancers
8            Heterochromatin
9              Quiescent_Low
10        Repressed_PolyComb
11      Strong_transcription
12   Transcr_at_gene_5_and_3
13   Weak_Repressed_PolyComb
14        Weak_transcription
15     ZNF_genes_and_repeats

我想要结果矩阵/数据。框架类似于下面的内容。我要切掉第一个数据。框架的第一列包含第二个数据。框架的列,并创建一个如下所示的矩阵。

>dd_matrix
                                 Active_TSS Bivalent_Enhancer
4star                              10.99562         0.4242592
A549_EtOH_0.02pct_Lung_Carcinoma   14.01553         1.0046676
                                 Bivalent_Poised_TSS Enhancers
4star                                      0.3116373 139.64713
A549_EtOH_0.02pct_Lung_Carcinoma           1.2060777  63.36004
                                 Flanking_Active_TSS Flanking_Bivalent_TSS_Enh
4star                                       23.11962                 0.1765451
A549_EtOH_0.02pct_Lung_Carcinoma            48.27401                 1.1449923
                                 Genic_enhancers Heterochromatin Quiescent_Low
4star                                  44.423945      -34.610990     -28.54241
A549_EtOH_0.02pct_Lung_Carcinoma        5.976754       -1.274768     -31.68228
                                 Repressed_PolyComb Strong_transcription
4star                                     -0.420967            12.708956
A549_EtOH_0.02pct_Lung_Carcinoma          -0.331186             3.375022
                                 Transcr_at_gene_5_and_3
4star                                           0.000000
A549_EtOH_0.02pct_Lung_Carcinoma                1.501412
                                 Weak_Repressed_PolyComb Weak_transcription
4star                                         0.08902141           19.28859
A549_EtOH_0.02pct_Lung_Carcinoma             -0.05151471           11.19855
                                 ZNF_genes_and_repeats
4star                                        -1.583407
A549_EtOH_0.02pct_Lung_Carcinoma              0.000000

任何帮助都将不胜感激。任何Runix解决方案都可以工作。

非常感谢。

共有2个答案

闻人宝
2023-03-14

棘手的事情似乎是,状态和前缀之间没有分隔符,有些状态是其他状态的一部分。不管怎样,这就是我想到的:

library(tidyverse)
dat  %>%
  rowwise() %>%
  mutate(postfix = max(states$V1[str_detect(V1, states$V1)])) %>%
  mutate(prefix = str_replace(V1, str_c("_", postfix), "")) %>%
  melt(id.vars = c("postfix", "prefix"), measure.vars = "V2") %>%
  dcast(prefix ~ postfix) 

在第一个mutate中,识别最长的匹配状态以形成后缀。在第二个mutate中,删除该后缀以形成前缀。

连曜灿
2023-03-14

下面是一个使用dplyr/tidyr的解决方案:

require(tidyverse);
df %>%
    separate(V1, into = c("what","states"), "_", extra = "merge") %>%
    spread(states, V2) %>%
    column_to_rownames("what");
#      Active_TSS Bivalent_Enhancer Bivalent_Poised_TSS Enhancers
#4star   10.99562         0.4242592           0.3116373  139.6471
#A549          NA                NA                  NA        NA
#      EtOH_0.02pct_Lung_Carcinoma_Active_TSS
#4star                                     NA
#A549                                14.01553
#      EtOH_0.02pct_Lung_Carcinoma_Bivalent_Enhancer
#4star                                            NA
#A549                                       1.004668
#      EtOH_0.02pct_Lung_Carcinoma_Bivalent_Poised_TSS
#4star                                              NA
#A549                                         1.206078
#      EtOH_0.02pct_Lung_Carcinoma_Enhancers
#4star                                    NA
#A549                               63.36004
#      EtOH_0.02pct_Lung_Carcinoma_Flanking_Active_TSS Flanking_Active_TSS
#4star                                              NA            23.11962
#A549                                         48.27401                  NA
#      Flanking_Bivalent_TSS_Enh Genic_enhancers Heterochromatin Quiescent_Low
#4star                 0.1765451        44.42395       -34.61099     -28.54241
#A549                         NA              NA              NA            NA
#      Repressed_PolyComb Strong_transcription Transcr_at_gene_5_and_3
#4star          -0.420967             12.70896                       0
#A549                  NA                   NA                      NA
#      Weak_Repressed_PolyComb Weak_transcription ZNF_genes_and_repeats
#4star              0.08902141           19.28859             -1.583407
#A549                       NA                 NA                    NA

说明:通过在第一个“quot上拆分,将V1分为两列;然后将states用作键,将V2用作值,将what列转换为行名称,将扩展为宽格式。

df <- read.table(text =
    "                                                    V1           V2
1                                      4star_Active_TSS  10.99561503
2                               4star_Bivalent_Enhancer   0.42425920
3                             4star_Bivalent_Poised_TSS   0.31163730
4                                       4star_Enhancers 139.64713405
5                             4star_Flanking_Active_TSS  23.11961717
6                       4star_Flanking_Bivalent_TSS_Enh   0.17654506
7                                 4star_Genic_enhancers  44.42394542
8                                 4star_Heterochromatin -34.61099049
9                                   4star_Quiescent_Low -28.54240987
10                             4star_Repressed_PolyComb  -0.42096698
11                           4star_Strong_transcription  12.70895605
12                        4star_Transcr_at_gene_5_and_3   0.00000000
13                        4star_Weak_Repressed_PolyComb   0.08902141
14                             4star_Weak_transcription  19.28859369
15                          4star_ZNF_genes_and_repeats  -1.58340662
16          A549_EtOH_0.02pct_Lung_Carcinoma_Active_TSS  14.01552989
17   A549_EtOH_0.02pct_Lung_Carcinoma_Bivalent_Enhancer   1.00466761
18 A549_EtOH_0.02pct_Lung_Carcinoma_Bivalent_Poised_TSS   1.20607773
19           A549_EtOH_0.02pct_Lung_Carcinoma_Enhancers  63.36004048
20 A549_EtOH_0.02pct_Lung_Carcinoma_Flanking_Active_TSS  48.27400816", header = T)

 类似资料:
  • 问题内容: 我有一个要转换为json格式的数据框: 我的数据帧称为res1: 当我做: 我得到这个: 我需要这个json输出像这样,有什么想法吗? 问题答案: 怎么样 通过使用,我们实际上将大的data.frame分解为每一行的单独的data.frame。通过从结果列表中删除名称,该函数将结果包装在数组中,而不是命名对象中。

  • 我有一个尺寸为3695 X 20的数据框。第一列包含字母数字标识符,其他19列都是数字。因此,rownames(df)提供数字1-3695,colnames(df)提供列的名称。df[,1]提供字母数字标识符。 我希望将数据框转换为矩阵,并使用现有数据框的第1列作为新矩阵的行名称,并将数据框的列名保留为矩阵的列名。 我还想自动化这个过程,用于类似但不同维度的数据帧。因此,如果解决这个问题需要知道行

  • 本文向大家介绍如何通过R中的行将数据帧值转换为向量?,包括了如何通过R中的行将数据帧值转换为向量?的使用技巧和注意事项,需要的朋友参考一下 数据可以任何形式提供给我们,但有可能不适合用于分析。有时数据记录在数据框中,但我们可能需要将其作为向量。在这种情况下,我们必须更改向量中数据帧的值。这可以通过在将数据帧与t换位后将它们读取为矢量来读取数据帧值来完成。 示例 请看以下数据帧- 让我们看另一个例子

  • 我有一个数据框,行名称中有单元格ID,列x中有分类。 我还有一个矩阵,列中有细胞ID,行中有基因。 我想粘贴来自数据框的矩阵分类的名称,请参见下面的示例。

  • 本文向大家介绍如何在R中将数据帧转换为data.table?,包括了如何在R中将数据帧转换为data.table?的使用技巧和注意事项,需要的朋友参考一下 由于对data.table的操作有时比数据帧快,因此我们可能希望将数据帧转换为data.table对象。数据框和data.table之间的主要区别在于,数据框可在基础R中使用,但是要使用data.table,我们必须安装软件包data.tabl

  • 我有这样的数据框 输出: . 输出: 我想分组的数据'ChannelPartnerID',因为我想合并数据的唯一值的唯一值在'ChannelPartnerID'和转换'品牌'列值为列包含销售价格为每一列 我希望它像这样输出: