01 2024 档案
EMMAX GWAS 分析后基因显著SNP位点
摘要:EMMAX 软件 library(tidyverse) # cat sig.R # GWAS output of emmax process_GWAS_file <- function(file_name, output_directory) { # 读取文件 data <- read.delim( 阅读全文
posted @ 2024-01-20 16:08 燕子南飞0415 阅读(194) 评论(0) 推荐(0)
合并CMplot曼哈顿和QQ-plot图
摘要:# cat merge.py import sys from PIL import Image def merge_images(image1_path, image2_path, output_path): # 打开两个图片 img1 = Image.open(image1_path) img2 阅读全文
posted @ 2024-01-20 13:24 燕子南飞0415 阅读(94) 评论(0) 推荐(0)
分割表型数据文件-(genotype,phenotype)
摘要:# cat split_file.py,分割表型文件 import csv import sys # 检查是否提供了输入文件名和trait文件名 if len(sys.argv) < 3: print("Usage: python script.py <input_file.txt> <trait_ 阅读全文
posted @ 2024-01-20 12:38 燕子南飞0415 阅读(44) 评论(0) 推荐(0)
CMplot 绘制曼哈顿图和QQ图
摘要:#!/usr/bin/env Rscript # cat ../manhattan_cmplot2.R # Rscript ../manhattan_cmplot2.R gwas_output.snp gwas_output.png library(tidyverse) args=commandAr 阅读全文
posted @ 2024-01-19 23:26 燕子南飞0415
emmax程序GWAS分析
摘要:EMMAX的主要特点: 线性混合模型:EMMAX使用线性混合模型(Linear Mixed Model, LMM),其中包括固定效应(如待检验的SNP)和随机效应(通常是亲缘关系矩阵或种群结构)来校正样本间的相关性。 高效计算:它通过使用近似方法来计算随机效应,使得分析大规模数据集(如人类或动植物基 阅读全文
posted @ 2024-01-19 23:16 燕子南飞0415