Generate plots for COG analysis

tli14 · Aug 3, 2022 · 6c23171 · 6c23171
1 parent d33061c
commit 6c23171
Showing 1 changed file with 6 additions and 6 deletions.
diff --git a/R_code/Fig6.COG_analysis.R → R_code/Fig5.COG_analysis.R b/R_code/Fig6.COG_analysis.R → R_code/Fig5.COG_analysis.R
@@ -6,12 +6,12 @@ library(reshape2)
 ## Here, the COG analysis of  E. coli fragmentation groups are compared (based on the core gene representatives generated by using Roary with 100% CG threshold).
 
 ## Load COG analysis summary for E. coli datasets.
-EC_rad1_ori <- read.table(file = "R_code/Fig6.COG_analysis/EC_ori.txt", header = FALSE)
-EC_rad1_50cut <- read.table(file = "R_code/Fig6.COG_analysis/EC_50cut.txt", header = FALSE)
-EC_rad1_100cut <- read.table(file = "R_code/Fig6.COG_analysis/EC_100cut.txt", header = FALSE)
-EC_rad1_200cut <- read.table(file = "R_code/Fig6.COG_analysis/EC_200cut.txt", header = FALSE)
-EC_rad1_300cut <- read.table(file = "R_code/Fig6.COG_analysis/EC_300cut.txt", header = FALSE)
-EC_rad1_400cut <- read.table(file = "R_code/Fig6.COG_analysis/EC_400cut.txt", header = FALSE)
+EC_rad1_ori <- read.table(file = "R_code/Fig5.COG_analysis/EC_ori.txt", header = FALSE)
+EC_rad1_50cut <- read.table(file = "R_code/Fig5.COG_analysis/EC_50cut.txt", header = FALSE)
+EC_rad1_100cut <- read.table(file = "R_code/Fig5.COG_analysis/EC_100cut.txt", header = FALSE)
+EC_rad1_200cut <- read.table(file = "R_code/Fig5.COG_analysis/EC_200cut.txt", header = FALSE)
+EC_rad1_300cut <- read.table(file = "R_code/Fig5.COG_analysis/EC_300cut.txt", header = FALSE)
+EC_rad1_400cut <- read.table(file = "R_code/Fig5.COG_analysis/EC_400cut.txt", header = FALSE)
 
 ## Reformat and combined the data.
 group <- rep("ori", length(EC_rad1_ori$V1))