GWAS曼哈顿图叠加代码

Question

求多个GWAS结果曼哈顿图合并叠加一起的R语言代码，非常感谢。

Answer 1

引自强大而免费的微信小程序：皆我百晓生。

---

以下是一个示例代码，用于将多个GWAS结果的曼哈顿图叠加在一起：

# 导入所需的R包
library(ggplot2)

# 创建绘图函数，用于绘制曼哈顿图
manhattan_plot <- function(data, title) {
  ggplot(data, aes(x = POSITION, y = -log10(PVALUE), color = GROUP)) +
    geom_point(size = 1) +
    scale_color_manual(values = c("blue", "red", "green")) +
    theme_bw() +
    theme(plot.title = element_text(hjust = 0.5),
          axis.title.x = element_blank(),
          axis.title.y = element_text(size=12),
          axis.text.x = element_blank(),
          axis.ticks.x = element_blank()) +
    ggtitle(title)
}

# 读取多个GWAS结果文件，并使用不同的组别进行标记
gwas_file_list <- c("gwas_result_1.txt", "gwas_result_2.txt", "gwas_result_3.txt")
group_list <- c("Group1", "Group2", "Group3")
gwas_data_list <- lapply(gwas_file_list, function(file){
  data <- read.table(file, header = TRUE, sep = "\t")
  data$GROUP <- file
  return(data)
})

# 合并所有GWAS结果的数据
combined_data <- do.call(rbind, gwas_data_list)

# 绘制合并后的曼哈顿图
manhattan_plot(combined_data, "Combined GWAS Manhattan Plot")

在上述代码中，你需要将gwas_file_list替换为你的GWAS结果文件名列表，并将group_list替换为与GWAS结果文件对应的组别列表。另外，gwas_result_1.txt等文件需要替换为你的实际GWAS结果文件名。