發(fā)簡信
IP屬地：陜西
#第一步計算admixtrue結(jié)果，寫成admixtrue.sh疗疟。nohup bashadmixtrue.sh for i in {6..20}...
在vcf基因 型整理過程中如失，有不同平臺不同樣本的vcf一鍵合成，用perl腳本完成 #!/usr/bin/env perl use strict...
####根據(jù)個體ID,調(diào)取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...
# 使用隨機(jī)效應(yīng)交互對遺傳異質(zhì)性建模 #上提到的問題限嫌，抽樣之后準(zhǔn)確率不穩(wěn)定靴庆，有可能少數(shù)幾次準(zhǔn)確率很低，可能因為基因型頻率出現(xiàn)問題（基因型異質(zhì)性導(dǎo)...
library(BGLR) #載入BGLR包 nIter=2000 #### number of iteration設(shè)置迭代次數(shù) burnIn=...