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