老师好,我在做群里遗传进化算pca时,得到的文件都是0。但是没有报错

attachments-2024-06-X23cNj5666791dcd8e850.png这是输入的vcf文件,200万个SNP

nohup文件信息:PLINK v1.90b4.6 64-bit (15 Aug 2017)           www.cog-genomics.org/plink/1.9/

(C) 2005-2017 Shaun Purcell, Christopher Chang   GNU General Public License v3

Logging to plink_pca.log.

Options in effect:

  --allow-extra-chr

  --out plink_pca

  --pca 10

  --set-missing-var-ids @:#

  --vcf /work/00.filter/merged_all.vcf

  --vcf-half-call missing


30310 MB RAM detected; reserving 15155 MB for main workspace.


--vcf: 1k variants complete.

--vcf: 2k variants complete.

--vcf: 3k variants complete.

--vcf: 4k variants complete.

--vcf: 5k variants complete.

--vcf: 6k variants complete.

--vcf: 7k variants complete.

--vcf: 8k variants complete.

--vcf: 9k variants complete.

--vcf: 10k variants complete.

--vcf: 11k variants complete.

--vcf: 12k variants complete.

...

2044980 markers complete.

2045040 markers complete.

2045100 markers complete.

2045160 markers complete.

2045220 markers complete.

2045265 markers complete.

Relationship matrix calculation complete.

[extracting eigenvalues and eigenvectors]

--pca: Results saved to plink_pca.eigenval and plink_pca.eigenvec .

报错:
PLINK v1.90b4.6 64-bit (15 Aug 2017)           www.cog-genomics.org/plink/1.9/
(C) 2005-2017 Shaun Purcell, Christopher Chang   GNU General Public License v3
Logging to plink_pca.log.
Options in effect:
  --allow-extra-chr
  --out plink_pca
  --pca 10
  --set-missing-var-ids @:#
  --vcf /work/00.filter/merged_all.vcf.gz
  --vcf-half-call missing

30310 MB RAM detected; reserving 15155 MB for main workspace.
--vcf: plink_pca-temporary.bed + plink_pca-temporary.bim +
plink_pca-temporary.fam written.
2045265 variants loaded from .bim file.
2045265 missing IDs set.
684 people (0 males, 0 females, 684 ambiguous) loaded from .fam.
Ambiguous sex IDs written to plink_pca.nosex .
Using up to 23 threads (change this with --threads).
Before main variant filters, 684 founders and 0 nonfounders present.
Calculating allele frequencies... done.
Total genotyping rate is 0.391328.
2045265 variants and 684 people pass filters and QC.
Note: No phenotypes present.
Relationship matrix calculation complete.
--pca: Results saved to plink_pca.eigenval and plink_pca.eigenvec .
请先 登录 后评论

2 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

那是不是有包错信息,你可以截图看看的额,不然不知道什么原因;

可能是你的输入文件有问题吧,你吧你的输入文件都打开看看的;

请先 登录 后评论
马杰宇

attachments-2024-06-Nn8lmnGD66791ce29aa32.png

请先 登录 后评论