plink Error: Failed to extract eigenvector(s) from GRM.

Error: Failed to extract eigenvector(s) from GRM.

001、

复制代码
root@DESKTOP-1N42TVH:/home/test5# ls
test.map  test.ped
root@DESKTOP-1N42TVH:/home/test5# cat test.map
1       rs11260595      0       1028961
1       rs6671356       0       1029889
1       rs6604968       0       1031563
1       rs4970405       0       1038818
1       rs12726255      0       1039813
1       rs4970409       0       1043315
1       rs7548798       0       1050037
1       rs17160824      0       1050471
1       rs11807848      0       1051029
1       rs12757754      0       1051125
root@DESKTOP-1N42TVH:/home/test5# cat test.ped
1444 NA12739 0 0 1 -9 C C C T G A G A G A A G C C A G C T G A
1344 NA10850 0 k 2 -9 A C C T G A A A G A G G T C G G C T A A
1328 NA06984 0 0 1 2 C C C T G A G A G A A G C C A G C T G A
1463 NA12877 0 0 1 -9 C C T T A A A A A A G G C C G G T T A A
1418 NA12275 0 0 2 1 C C T T A A A A A A G G C C G G T T A A
13291 NA06986 0 0 1 1 C C C T G A G A G A A G C C A G C T G A
root@DESKTOP-1N42TVH:/home/test5# plink --file test --pca 3
PLINK v1.90b6.26 64-bit (2 Apr 2022)           www.cog-genomics.org/plink/1.9/
(C) 2005-2022 Shaun Purcell, Christopher Chang   GNU General Public License v3
Logging to plink.log.
Options in effect:
  --file test
  --pca 3

16007 MB RAM detected; reserving 8003 MB for main workspace.
.ped scan complete (for binary autoconversion).
Performing single-pass .bed write (10 variants, 6 people).
--file: plink-temporary.bed + plink-temporary.bim + plink-temporary.fam
written.
10 variants loaded from .bim file.
6 people (4 males, 2 females) loaded from .fam.
3 phenotype values loaded from .fam.
Using up to 11 threads (change this with --threads).
Before main variant filters, 5 founders and 1 nonfounder present.
Calculating allele frequencies... done.
Total genotyping rate is exactly 1.
10 variants and 6 people pass filters and QC.
Among remaining phenotypes, 1 is a case and 2 are controls.  (3 phenotypes are
missing.)
Relationship matrix calculation complete.
[extracting eigenvalues and eigenvectors]
Error: Failed to extract eigenvector(s) from GRM.
复制代码

 

 

002、

复制代码
root@DESKTOP-1N42TVH:/home/test5# ls
test.map  test.ped
root@DESKTOP-1N42TVH:/home/test5# cat test.map
1       rs11260595      0       1028961
1       rs6671356       0       1029889
1       rs6604968       0       1031563
1       rs4970405       0       1038818
1       rs12726255      0       1039813
1       rs4970409       0       1043315
1       rs7548798       0       1050037
1       rs17160824      0       1050471
1       rs11807848      0       1051029
1       rs12757754      0       1051125
root@DESKTOP-1N42TVH:/home/test5# cat test.ped
1444 NA12739 0 0 1 -9 C C C T G A G A G A A G C C A G C T G A
1344 NA10850 0 0 2 -9 A C C T G A A A G A G G T C G G C T A A
1328 NA06984 0 0 1 2 C C C T G A G A G A A G C C A G C T G A
1463 NA12877 0 0 1 -9 C C T T A A A A A A G G C C G G T T A A
1418 NA12275 0 0 2 1 C C T T A A A A A A G G C C G G T T A A
13291 NA06986 0 0 1 1 C C C T G A G A G A A G C C A G C T G A
root@DESKTOP-1N42TVH:/home/test5# plink --file test --pca 3
PLINK v1.90b6.26 64-bit (2 Apr 2022)           www.cog-genomics.org/plink/1.9/
(C) 2005-2022 Shaun Purcell, Christopher Chang   GNU General Public License v3
Logging to plink.log.
Options in effect:
  --file test
  --pca 3

16007 MB RAM detected; reserving 8003 MB for main workspace.
.ped scan complete (for binary autoconversion).
Performing single-pass .bed write (10 variants, 6 people).
--file: plink-temporary.bed + plink-temporary.bim + plink-temporary.fam
written.
10 variants loaded from .bim file.
6 people (4 males, 2 females) loaded from .fam.
3 phenotype values loaded from .fam.
Using up to 11 threads (change this with --threads).
Before main variant filters, 6 founders and 0 nonfounders present.
Calculating allele frequencies... done.
Total genotyping rate is exactly 1.
10 variants and 6 people pass filters and QC.
Among remaining phenotypes, 1 is a case and 2 are controls.  (3 phenotypes are
missing.)
Relationship matrix calculation complete.
--pca: Results saved to plink.eigenval and plink.eigenvec .
复制代码

 

posted @   小鲨鱼2018  阅读(579)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 单元测试从入门到精通
· 上周热点回顾(3.3-3.9)
· winform 绘制太阳,地球,月球 运作规律
历史上的今天:
2021-07-12 c语言为什么设计这么多种数据类型
2021-07-12 c语言中不同进制数的表示
2021-07-12 计算机的字长是指什么
点击右上角即可分享
微信分享提示