#读取2.22_output_firststept.csv文件 getwd() setwd("/Analysis_B/20220326_JNILI/010_summary") data <- read.csv("output_summary12.csv", header = F) getwd() setwd("/home/Reference_database/new_taxdump") read.table("fullnamelineage.dmp", header = F, sep = '\t', fill = TRUE,stringsAsFactors=FALSE, quote = "") A <- read.table("fullnamelineage.dmp", header = F, sep = '\t', fill = TRUE,stringsAsFactors=FALSE, quote = "") [c(1,3,5)] B <- read.table("host.dmp", header = F, sep ='|') M1 <- merge(data,B,by='saccver_N',all.X=TRUE) M2 <- merge(M1,A,by='saccver_N',all.X=TRUE) write.csv(M2, file = "/Analysis_B/20220326_JNILI/010_summary/output_summary13.csv",row.names = T)这段代码什么意思
时间: 2023-11-02 13:04:56 浏览: 83
这段代码的主要作用是将两个数据文件进行合并,最后将合并后的数据写入一个新的csv文件中。具体实现过程如下:
1. 使用getwd函数获取当前的工作目录,然后使用setwd函数将工作目录更改为/Analysis_B/20220326_JNILI/010_summary/,并读取名为output_summary12.csv的csv文件到data数据框中。
2. 使用getwd函数获取当前的工作目录,然后使用setwd函数将工作目录更改为/home/Reference_database/new_taxdump/,并读取名为fullnamelineage.dmp的文本文件到A数据框中,同时只保留数据框中第1、3、5列的数据。
3. 使用read.table函数读取名为host.dmp的文本文件到B数据框中,文件中的列之间以"|"分隔。
4. 使用merge函数将data数据框和B数据框按“saccver_N”列进行合并,合并后的结果保存到数据框M1中,all.X=TRUE表示使用左连接方式,即只保留data数据框中的数据。
5. 使用merge函数将M1数据框和A数据框按“saccver_N”列进行合并,合并后的结果保存到数据框M2中,all.X=TRUE表示使用左连接方式,即只保留M1数据框中的数据。
6. 使用write.csv函数将合并后的数据框M2写入到/Analysis_B/20220326_JNILI/010_summary/output_summary13.csv文件中,row.names=T表示将行名写入文件中。
综上,这段代码主要的目的是将三个数据文件进行合并,并将合并后的数据写入一个新的csv文件中,方便后续的数据分析和处理。
阅读全文