R:SNP数据筛选.R
# 清空工作环境,删除所有对象
rm(list = ls())
library(data.table)
# 加载data.table包
library(data.table)
# 设置工作目录
setwd("C:\\Users\\Administrator\\Desktop") # 设置为你的工作目录路径
# 读取大文件
file_path <- "AMP_lines456_maf0.05_het2N_geno0.6.hmp.txt" # 替换为你的文件路径
data <- fread(file_path, header = TRUE, sep = "\t") # 使用fread高效读取大文件
# 筛选并在指定范围上下50KB
start_range <- 136606701 - 50000 # 计算范围下限
end_range <- 136613008 + 50000 # 计算范围上限
# 根据条件筛选chr6数据
filtered_data <- data[chrom == "6" & pos >= start_range & pos <= end_range]
# 将筛选结果写入新的文件,确保NA值被正确写出,不添加引号
output_path <- "filtered_data.txt" # 输出文件路径
fwrite(filtered_data, output_path, sep = "\t", na = "NA", quote = FALSE) # 写入文件
# 显示筛选后的结果
print(filtered_data) # 打印筛选后的数据
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· DeepSeek “源神”启动!「GitHub 热点速览」
· 我与微信审核的“相爱相杀”看个人小程序副业
· 微软正式发布.NET 10 Preview 1:开启下一代开发框架新篇章
· C# 集成 DeepSeek 模型实现 AI 私有化(本地部署与 API 调用教程)
· spring官宣接入deepseek,真的太香了~
2023-06-19 PSO算法