群体遗传进化PCA分析

群体遗传进化PCA分析

#方法1:

## plink分析PCA

plink --vcf  $workdir/00.filter/clean.sorted.vcf.gz --pca 10 --out  plink_pca   \

    --allow-extra-chr --set-missing-var-ids @:#    --vcf-half-call missing


#绘图


pca_plink_plot.r -i plink_pca.eigenvec -f $GROUP -g group --name plink_pca





#方法2:

## smartpca分析PCA


vcf2smartpca.py --vcf  $workdir/00.filter/clean.sorted.vcf.gz -q 20 -g $GROUP -p passing_snps_q20

smartpca -p par.passing_snps_q20 > passing_snps_q20_logfile.txt 

#绘图展示


smartpca_plot.r -i passing_snps_q20_loci.evec -f $GROUP -g group --name smartpca

smartpca_plot.py -i passing_snps_q20_loci.evec -v passing_snps_q20_loci.eval -p smartpca3d




  • 发表于 2021-09-01 13:47
  • 阅读 ( 2901 )
  • 分类:遗传进化

0 条评论

请先 登录 后评论
omicsgene
omicsgene

生物信息

702 篇文章

作家榜 »

  1. omicsgene 702 文章
  2. 安生水 351 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. xun 82 文章
  6. rzx 78 文章
  7. 红橙子 78 文章
  8. CORNERSTONE 72 文章