我在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
任何帮助都将不胜感激。任何R
或unix
解决方案都可以工作。
非常感谢。
棘手的事情似乎是,状态和前缀之间没有分隔符,有些状态是其他状态的一部分。不管怎样,这就是我想到的:
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
中,删除该后缀
以形成前缀。
下面是一个使用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'和转换'品牌'列值为列包含销售价格为每一列 我希望它像这样输出: