发简信
IP属地：河南
df=read.csv("19en.csv", header=T, row.names=1)head(df)data.pca <- prcomp(df[,1:19])summ...
ggplot(data = df, mapping = aes(x=year, y=value,linetype=type, colour=type, shape=type,...
下载参考基因组wget https://ftp.ncbi.nih.gov/genomes/all/GCA/000/001/405/GCA_000001405.28_GRCh3...
1 SNP名称替换 替换固定列的内容 awk -F "\t" '{ OFS="\t" }{1"chr"2=4"'snp'"}1' goat232-1.map > goat23...
cd /mnt/d/peng/climate/paper-climate/climate-1/NJtreeplink --file 72breeds --indep-pai...
E:\ > cd e:E:E:\ > cd plink-1E:\plink-1>plink –file test Map 更新Plink --goat --file data...
library(LEA) 结果会生成test.geno文件的数据. output = ped2lfmm("72breeds-ro.ped") 使用LEA进行structure...
for line in cat chr_griddochr={line##*,}./SweeD -name B1.chr{i}.vcf -grid ${grid} -mins...
pophelper包的安装 pophelper包在Github中，所以需要用到devtools包来进行安装 输入文件 sfiles <- list.files("struct...
plink 做PCA --vcf vcf 或者vcf.gz --recode 输出格式 --out 输入前缀 --const-fid 添加群体信息 --allow-extr...
