文件怎么处理不重要,需要关注的如何将文件读入和写出
library(tidyverse)
# 读取数据文件
res <- read.table("test.txt", header = TRUE)
name <- read.csv("name.csv", header = FALSE)
head(name)
# 设置文件路径和文件模式
folder_path <- "D:/1.1BIO/Test/"
file_pattern <- "*motif"
# 获取满足文件模式的文件列表
temp <- list.files(path = folder_path, pattern = file_pattern)
temp
# 循环处理每个文件
for (i in 1:length(temp)) {
# 打印进度
progress_percent <- round(i / length(temp) * 100, 2)
print(paste0("Progress:", progress_percent, "%"))
# 读取当前文件
mo <- read.table(paste0(folder_path, temp[i]), header = TRUE)
colnames(mo) <- c("s", "e")
# 合并数据框
df <- merge(mo, res, by.x = "s", by.y = "Gr")
colnames(df)[3:4] <- c("s_Gh", "s_At")
df <- merge(df, res, by.x = "e", by.y = "Gr")
colnames(df)[5:6] <- c("e_Gh", "e_At")
fr <- df[, c(4, 6)]
# 将结果写入文件
write.table(fr, paste0(temp[i], "_At"), sep = "\t", row.names = FALSE, quote = FALSE)
# 合并 name 数据框以获取名称信息
fr <- merge(fr, name, by.x = "s_At", by.y = "V1")
colnames(fr)[3] <- "s_name"
head(fr)
fr <- merge(fr, name, by.x = "e_At", by.y = "V1")
colnames(fr)[4] <- "e_name"
head(fr)
fr <- fr[, c(3, 4)]
# 写入带有名称信息的结果文件
write.table(fr, paste0(temp[i], "_At_name"), sep = "\t", row.names = FALSE, quote = FALSE)
}