摘要： ###1.交互状态 #!/bin/bash #count genes in region,need three para #author lee if [ $# -ne 3 ] then echo "Usage chr start end" else num=$(awk -v chr=$1 -v s 阅读全文

摘要： if(!require("qqman"))install.packages("qqman") if(!require("CMplot"))install.packages("CMplot") library(qqman) #library(CMplot) library(tidyverse) set 阅读全文

摘要： #!/bin/bash grep -v ^# $i|awk '{sum[$1+=1]END{for(i in sum)print i"\t" sum[i]}' >chr_info #按大小排序 sort -n chr_info #某列求和 awk '{sum +=$2;END {print sum} 阅读全文

摘要： .: read.table("Fig_chr1_wild_africa.bin")->data1read.table("Fig_chr2_wild_africa.bin")->data2read.table("Fig_chr3_wild_africa.bin")->data3read.table(" 阅读全文

摘要： #author lee #date 2021.8.26 library(tidyverse) #1创建数据 df <- data.frame( "grammer" = c("Python","C","Java","GO",NA,"SQL","PHP","Python"), "score" = c(1 阅读全文

摘要： #!/bin/bash #vlookup #author lee #date 2021.8.26 if [ $# -ne 2 ] then echo "File number supplied is not right" echo "Input File1 that have your info" 阅读全文

摘要： #！/bin/bash plink --bfile ../cold/cold.filter --pca 2 --out 2pc https://zhuanlan.zhihu.com/p/357759550 向mrMLM包中添加主成分，需要加两行header 阅读全文

摘要： #PBS -N admixture #PBS -l nodes=1:ppn=12 #PBS -q middle cd $PBS_O_WORKDIR for K in {1..15}; do ~/biosoft/dist/admixture_linux-1.3.0/admixture -j12 --c 阅读全文

摘要： #!/bin/bash #use id to keep vcf and hmp #author leee #time 2021.8.23 echo "Input ID_file as args1 and Phenofile as args2!" >>${1}.log echo "Your ID_fi 阅读全文

摘要： 注意表型首行使用标签<Trait> 脚本mrmlm.r #!/public/home/caisl/biosoft/miniconda3/envs/r4.1/bin/R library('mrMLM') mrMLM(fileGen='184.hmp.txt',filePhe='pheno',Genfo 阅读全文