发简信
IP属地：陕西
#第一步计算admixtrue结果，写成admixtrue.sh。nohup bashadmixtrue.sh for i in {6..20}...
在vcf基因 型整理过程中，有不同平台不同样本的vcf一键合成，用perl脚本完成 #!/usr/bin/env perl use strict...
####根据个体ID,调取VCF文件 import vcf # 读取个体ID with open('171q.txt', 'r') as f: ...
import pysamdef extract_vcf_info(vcf_gz_filename, output_filename): vc...
import pysamdef add_ids_to_vcf(input_vcf, output_vcf): vcf = pysam.Var...
def vcf_to_012_without_pysam(vcf_filename, output_filename): data = []...
setwd("D:\\GWAS\\QTL") library(BGLR) library(plyr) library(data.table) w...
# 使用随机效应交互对遗传异质性建模 #上提到的问题，抽样之后准确率不稳定，有可能少数几次准确率很低，可能因为基因型频率出现问题（基因型异质性导...
library(BGLR) #载入BGLR包 nIter=2000 #### number of iteration设置迭代次数 burnIn=...