Back to Multiple platform build/check report for BioC 3.8
ABCDEFGHIJKLMNOPQR[S]TUVWXYZ

CHECK report for SNPRelate on tokay1

This page was generated on 2019-04-13 11:24:16 -0400 (Sat, 13 Apr 2019).

Package 1476/1649HostnameOS / ArchINSTALLBUILDCHECKBUILD BIN
SNPRelate 1.16.0
Xiuwen Zheng
Snapshot Date: 2019-04-12 17:01:30 -0400 (Fri, 12 Apr 2019)
URL: https://git.bioconductor.org/packages/SNPRelate
Branch: RELEASE_3_8
Last Commit: 0e38e8d
Last Changed Date: 2018-10-30 11:41:55 -0400 (Tue, 30 Oct 2018)
malbec1 Linux (Ubuntu 16.04.6 LTS) / x86_64  OK  OK  NA 
tokay1 Windows Server 2012 R2 Standard / x64  OK  OK [ OK ] OK UNNEEDED, same version exists in internal repository
merida1 OS X 10.11.6 El Capitan / x86_64  OK  OK  OK  OK UNNEEDED, same version exists in internal repository

Summary

Package: SNPRelate
Version: 1.16.0
Command: C:\Users\biocbuild\bbs-3.8-bioc\R\bin\R.exe CMD check --force-multiarch --install=check:SNPRelate.install-out.txt --library=C:\Users\biocbuild\bbs-3.8-bioc\R\library --no-vignettes --timings SNPRelate_1.16.0.tar.gz
StartedAt: 2019-04-13 05:38:29 -0400 (Sat, 13 Apr 2019)
EndedAt: 2019-04-13 05:43:31 -0400 (Sat, 13 Apr 2019)
EllapsedTime: 302.3 seconds
RetCode: 0
Status:  OK  
CheckDir: SNPRelate.Rcheck
Warnings: 0

Command output

##############################################################################
##############################################################################
###
### Running command:
###
###   C:\Users\biocbuild\bbs-3.8-bioc\R\bin\R.exe CMD check --force-multiarch --install=check:SNPRelate.install-out.txt --library=C:\Users\biocbuild\bbs-3.8-bioc\R\library --no-vignettes --timings SNPRelate_1.16.0.tar.gz
###
##############################################################################
##############################################################################


* using log directory 'C:/Users/biocbuild/bbs-3.8-bioc/meat/SNPRelate.Rcheck'
* using R version 3.5.3 (2019-03-11)
* using platform: x86_64-w64-mingw32 (64-bit)
* using session charset: ISO8859-1
* using option '--no-vignettes'
* checking for file 'SNPRelate/DESCRIPTION' ... OK
* checking extension type ... Package
* this is package 'SNPRelate' version '1.16.0'
* checking package namespace information ... OK
* checking package dependencies ...Warning: unable to access index for repository https://CRAN.R-project.org/src/contrib:
  cannot open URL 'https://CRAN.R-project.org/src/contrib/PACKAGES'
 OK
* checking if this is a source package ... OK
* checking if there is a namespace ... OK
* checking for hidden files and directories ... OK
* checking for portable file names ... OK
* checking whether package 'SNPRelate' can be installed ... OK
* checking installed package size ... OK
* checking package directory ... OK
* checking 'build' directory ... OK
* checking DESCRIPTION meta-information ... OK
* checking top-level files ... OK
* checking for left-over files ... OK
* checking index information ... OK
* checking package subdirectories ... OK
* checking R files for non-ASCII characters ... OK
* checking R files for syntax errors ... OK
* loading checks for arch 'i386'
** checking whether the package can be loaded ... OK
** checking whether the package can be loaded with stated dependencies ... OK
** checking whether the package can be unloaded cleanly ... OK
** checking whether the namespace can be loaded with stated dependencies ... OK
** checking whether the namespace can be unloaded cleanly ... OK
* loading checks for arch 'x64'
** checking whether the package can be loaded ... OK
** checking whether the package can be loaded with stated dependencies ... OK
** checking whether the package can be unloaded cleanly ... OK
** checking whether the namespace can be loaded with stated dependencies ... OK
** checking whether the namespace can be unloaded cleanly ... OK
* checking dependencies in R code ... OK
* checking S3 generic/method consistency ... OK
* checking replacement functions ... OK
* checking foreign function calls ... OK
* checking R code for possible problems ... OK
* checking Rd files ... OK
* checking Rd metadata ... OK
* checking Rd cross-references ... OK
* checking for missing documentation entries ... OK
* checking for code/documentation mismatches ... OK
* checking Rd \usage sections ... OK
* checking Rd contents ... OK
* checking for unstated dependencies in examples ... OK
* checking contents of 'data' directory ... OK
* checking data for non-ASCII characters ... OK
* checking data for ASCII and uncompressed saves ... OK
* checking line endings in C/C++/Fortran sources/headers ... OK
* checking line endings in Makefiles ... OK
* checking compilation flags in Makevars ... OK
* checking for GNU extensions in Makefiles ... OK
* checking for portable use of $(BLAS_LIBS) and $(LAPACK_LIBS) ... OK
* checking compiled code ... NOTE
Note: information on .o files for i386 is not available
Note: information on .o files for x64 is not available
File 'C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/libs/i386/SNPRelate.dll':
  Found 'abort', possibly from 'abort' (C), 'runtime' (Fortran)
  Found 'exit', possibly from 'exit' (C), 'stop' (Fortran)
  Found 'printf', possibly from 'printf' (C)

Compiled code should not call entry points which might terminate R nor
write to stdout/stderr instead of to the console, nor use Fortran I/O
nor system RNGs. The detected symbols are linked into the code but
might come from libraries and not actually be called.

See 'Writing portable packages' in the 'Writing R Extensions' manual.
* checking files in 'vignettes' ... OK
* checking examples ...
** running examples for arch 'i386' ... OK
** running examples for arch 'x64' ... OK
* checking for unstated dependencies in 'tests' ... OK
* checking tests ...
** running tests for arch 'i386' ...
  Running 'runTests.R'
 OK
** running tests for arch 'x64' ...
  Running 'runTests.R'
 OK
* checking for unstated dependencies in vignettes ... OK
* checking package vignettes in 'inst/doc' ... OK
* checking running R code from vignettes ... SKIPPED
* checking re-building of vignette outputs ... SKIPPED
* checking PDF version of manual ... OK
* DONE

Status: 1 NOTE
See
  'C:/Users/biocbuild/bbs-3.8-bioc/meat/SNPRelate.Rcheck/00check.log'
for details.



Installation output

SNPRelate.Rcheck/00install.out

##############################################################################
##############################################################################
###
### Running command:
###
###   C:\cygwin\bin\curl.exe -O https://malbec1.bioconductor.org/BBS/3.8/bioc/src/contrib/SNPRelate_1.16.0.tar.gz && rm -rf SNPRelate.buildbin-libdir && mkdir SNPRelate.buildbin-libdir && C:\Users\biocbuild\bbs-3.8-bioc\R\bin\R.exe CMD INSTALL --merge-multiarch --build --library=SNPRelate.buildbin-libdir SNPRelate_1.16.0.tar.gz && C:\Users\biocbuild\bbs-3.8-bioc\R\bin\R.exe CMD INSTALL SNPRelate_1.16.0.zip && rm SNPRelate_1.16.0.tar.gz SNPRelate_1.16.0.zip
###
##############################################################################
##############################################################################


  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed

  0     0    0     0    0     0      0      0 --:--:-- --:--:-- --:--:--     0
100 1253k  100 1253k    0     0  25.3M      0 --:--:-- --:--:-- --:--:-- 27.8M

install for i386

* installing *source* package 'SNPRelate' ...
** libs
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c ConvToGDS.cpp -o ConvToGDS.o
C:/Rtools/mingw_32/bin/gcc  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O3 -Wall  -std=gnu99 -mtune=generic -c R_SNPRelate.c -o R_SNPRelate.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c SNPRelate.cpp -o SNPRelate.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c ThreadPool.cpp -o ThreadPool.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c dGenGWAS.cpp -o dGenGWAS.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c dVect.cpp -o dVect.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genBeta.cpp -o genBeta.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genEIGMIX.cpp -o genEIGMIX.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genFst.cpp -o genFst.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genHWE.cpp -o genHWE.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genIBD.cpp -o genIBD.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genIBS.cpp -o genIBS.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genKING.cpp -o genKING.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genLD.cpp -o genLD.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genPCA.cpp -o genPCA.o
C:/Rtools/mingw_32/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genSlideWin.cpp -o genSlideWin.o
C:/Rtools/mingw_32/bin/g++ -shared -s -static-libgcc -o SNPRelate.dll tmp.def ConvToGDS.o R_SNPRelate.o SNPRelate.o ThreadPool.o dGenGWAS.o dVect.o genBeta.o genEIGMIX.o genFst.o genHWE.o genIBD.o genIBS.o genKING.o genLD.o genPCA.o genSlideWin.o -LC:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/bin/i386 -lRlapack -LC:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/bin/i386 -lRblas -lgfortran -lm -lquadmath -LC:/extsoft/lib/i386 -LC:/extsoft/lib -LC:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/bin/i386 -lR
installing to C:/Users/biocbuild/bbs-3.8-bioc/meat/SNPRelate.buildbin-libdir/SNPRelate/libs/i386
** R
** data
** inst
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
  converting help for package 'SNPRelate'
    finding HTML links ... done
    SNPGDSFileClass-class                   html  
    SNPRelate-package                       html  
    hapmap_geno                             html  
    snpgdsAdmixPlot                         html  
    snpgdsAdmixProp                         html  
    snpgdsAlleleSwitch                      html  
    snpgdsApartSelection                    html  
    snpgdsBED2GDS                           html  
    finding level-2 HTML links ... done

    snpgdsClose                             html  
    snpgdsCombineGeno                       html  
    snpgdsCreateGeno                        html  
    snpgdsCreateGenoSet                     html  
    snpgdsCutTree                           html  
    snpgdsDiss                              html  
    snpgdsDrawTree                          html  
    snpgdsEIGMIX                            html  
    snpgdsErrMsg                            html  
    snpgdsExampleFileName                   html  
    snpgdsFst                               html  
    snpgdsGDS2BED                           html  
    snpgdsGDS2Eigen                         html  
    snpgdsGDS2PED                           html  
    snpgdsGEN2GDS                           html  
    snpgdsGRM                               html  
    snpgdsGetGeno                           html  
    snpgdsHCluster                          html  
    snpgdsHWE                               html  
    snpgdsIBDKING                           html  
    snpgdsIBDMLE                            html  
    snpgdsIBDMLELogLik                      html  
    snpgdsIBDMoM                            html  
    snpgdsIBDSelection                      html  
    snpgdsIBS                               html  
    snpgdsIBSNum                            html  
    snpgdsIndInb                            html  
    snpgdsIndInbCoef                        html  
    snpgdsIndivBeta                         html  
    snpgdsLDMat                             html  
    snpgdsLDpair                            html  
    snpgdsLDpruning                         html  
    snpgdsMergeGRM                          html  
    snpgdsOpen                              html  
    snpgdsOption                            html  
    snpgdsPCA                               html  
    snpgdsPCACorr                           html  
    snpgdsPCASNPLoading                     html  
    snpgdsPCASampLoading                    html  
    snpgdsPED2GDS                           html  
    snpgdsPairIBD                           html  
    snpgdsPairIBDMLELogLik                  html  
    snpgdsPairScore                         html  
    snpgdsSNPList                           html  
    snpgdsSNPListClass                      html  
    snpgdsSNPListIntersect                  html  
    snpgdsSNPRateFreq                       html  
    snpgdsSampMissRate                      html  
    snpgdsSelectSNP                         html  
    snpgdsSlidingWindow                     html  
    snpgdsSummary                           html  
    snpgdsTranspose                         html  
    snpgdsVCF2GDS                           html  
    snpgdsVCF2GDS_R                         html  
** building package indices
** installing vignettes
** testing if installed package can be loaded
In R CMD INSTALL

install for x64

* installing *source* package 'SNPRelate' ...
** libs
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c ConvToGDS.cpp -o ConvToGDS.o
C:/Rtools/mingw_64/bin/gcc  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -std=gnu99 -mtune=generic -c R_SNPRelate.c -o R_SNPRelate.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c SNPRelate.cpp -o SNPRelate.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c ThreadPool.cpp -o ThreadPool.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c dGenGWAS.cpp -o dGenGWAS.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c dVect.cpp -o dVect.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genBeta.cpp -o genBeta.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genEIGMIX.cpp -o genEIGMIX.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genFst.cpp -o genFst.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genHWE.cpp -o genHWE.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genIBD.cpp -o genIBD.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genIBS.cpp -o genIBS.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genKING.cpp -o genKING.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genLD.cpp -o genLD.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genPCA.cpp -o genPCA.o
C:/Rtools/mingw_64/bin/g++  -I"C:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/include" -DNDEBUG -DUSING_R -I. -I"C:/Users/biocbuild/bbs-3.8-bioc/R/library/gdsfmt/include"   -I"C:/extsoft/include"     -O2 -Wall  -mtune=generic -c genSlideWin.cpp -o genSlideWin.o
C:/Rtools/mingw_64/bin/g++ -shared -s -static-libgcc -o SNPRelate.dll tmp.def ConvToGDS.o R_SNPRelate.o SNPRelate.o ThreadPool.o dGenGWAS.o dVect.o genBeta.o genEIGMIX.o genFst.o genHWE.o genIBD.o genIBS.o genKING.o genLD.o genPCA.o genSlideWin.o -LC:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/bin/x64 -lRlapack -LC:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/bin/x64 -lRblas -lgfortran -lm -lquadmath -LC:/extsoft/lib/x64 -LC:/extsoft/lib -LC:/Users/BIOCBU˜1/BBS-3˜1.8-B/R/bin/x64 -lR
installing to C:/Users/biocbuild/bbs-3.8-bioc/meat/SNPRelate.buildbin-libdir/SNPRelate/libs/x64
** testing if installed package can be loaded
* MD5 sums
packaged installation of 'SNPRelate' as SNPRelate_1.16.0.zip
* DONE (SNPRelate)
In R CMD INSTALL
In R CMD INSTALL
* installing to library 'C:/Users/biocbuild/bbs-3.8-bioc/R/library'
package 'SNPRelate' successfully unpacked and MD5 sums checked
In R CMD INSTALL

Tests output

SNPRelate.Rcheck/tests_i386/runTests.Rout


R version 3.5.3 (2019-03-11) -- "Great Truth"
Copyright (C) 2019 The R Foundation for Statistical Computing
Platform: i386-w64-mingw32/i386 (32-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> BiocGenerics:::testPackage("SNPRelate")
SNPRelate
Genetic Relationship Matrix (GRM, GCTA):
Excluding 8,088 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 1,000 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 282597
CPU capabilities:
Sat Apr 13 05:41:31 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:41:32 2019    Done.
Genetic Relationship Matrix (GRM, GCTA):
Excluding 7,088 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 2,000 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 559412
CPU capabilities:
Sat Apr 13 05:41:32 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:41:33 2019    Done.
Genetic Relationship Matrix (GRM, GCTA):
Excluding 5,288 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 3,800 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 1066957
CPU capabilities:
Sat Apr 13 05:41:33 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:41:34 2019    Done.
GRM merging:
    open 'tmp1.gds' (1,000 variants)
    open 'tmp2.gds' (2,000 variants)
    open 'tmp3.gds' (3,800 variants)
Weight: 0.147059, 0.294118, 0.558824
Output: tmp.gds

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Genetic Relationship Matrix (GRM, GCTA):
Excluding 2,288 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 6,800 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 1908966
CPU capabilities:
Sat Apr 13 05:41:35 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:41:36 2019    Done.
Genetic Relationship Matrix (GRM, IndivBeta):
Excluding 8,088 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 1,000 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 282597
CPU capabilities:
Sat Apr 13 05:41:36 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:41:37 2019    Done.
Genetic Relationship Matrix (GRM, IndivBeta):
Excluding 7,088 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 2,000 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 559412
CPU capabilities:
Sat Apr 13 05:41:37 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:41:39 2019    Done.
Genetic Relationship Matrix (GRM, IndivBeta):
Excluding 5,288 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 3,800 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 1066957
CPU capabilities:
Sat Apr 13 05:41:39 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:41:40 2019    Done.
GRM merging:
    open 'tmp1.gds' (1,000 variants)
    open 'tmp2.gds' (2,000 variants)
    open 'tmp3.gds' (3,800 variants)
Weight: 0.147059, 0.294118, 0.558824
Output: tmp.gds

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Writing ...

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Genetic Relationship Matrix (GRM, IndivBeta):
Excluding 2,288 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 6,800 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 1908966
CPU capabilities:
Sat Apr 13 05:41:41 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:41:42 2019    Done.
Linkage Disequilibrium (LD) estimation on genotypes:
Working space: 279 samples, 1000 SNPs
    using 1 (CPU) core.
    method: covariance
LD matrix:    the sum of all selected genotypes (0,1,2) = 283058
Linkage Disequilibrium (LD) estimation on genotypes:
Working space: 279 samples, 1000 SNPs
    using 1 (CPU) core.
    method: correlation
LD matrix:    the sum of all selected genotypes (0,1,2) = 283058
FUNCTION: SNPGDSFileClass
FUNCTION: SNPRelate-package
Start snpgdsBED2GDS ...
	BED file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.bed.gz" in the SNP-major mode (Sample X SNP)
	FAM file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.fam.gz", DONE.
	BIM file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.bim.gz", DONE.
Sat Apr 13 05:41:46 2019 	store sample id, snp id, position, and chromosome.
	start writing: 60 samples, 5000 SNPs ...
 	Sat Apr 13 05:41:46 2019	0%
 	Sat Apr 13 05:41:46 2019	100%
Sat Apr 13 05:41:46 2019 	Done.
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'HapMap.gds' (98.1K)
    # of fragments: 38
    save to 'HapMap.gds.tmp'
    rename 'HapMap.gds.tmp' (97.8K, reduced: 240B)
    # of fragments: 18
Principal Component Analysis (PCA) on genotypes:
Excluding 203 SNPs on non-autosomes
Excluding 28 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 60 samples, 4,769 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 124273
CPU capabilities:
Sat Apr 13 05:41:46 2019    (internal increment: 53996)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:41:47 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:41:47 2019    Done.
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:41:47 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:41:48 2019    Done.
Identity-By-State (IBS) analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
IBS:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:41:49 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:41:50 2019    Done.
Linkage Disequilibrium (LD) estimation on genotypes:
Working space: 279 samples, 200 SNPs
    using 1 (CPU) core.
    method: composite
LD matrix:    the sum of all selected genotypes (0,1,2) = 55417
FUNCTION: hapmap_geno
FUNCTION: snpgdsAdmixPlot
Eigen-analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Eigen-analysis:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities:
Sat Apr 13 05:41:51 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:41:51 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:41:51 2019    Done.
FUNCTION: snpgdsAdmixProp
Eigen-analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Eigen-analysis:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities:
Sat Apr 13 05:41:52 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:41:52 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:41:52 2019    Done.
FUNCTION: snpgdsAlleleSwitch
Strand-switching at 50 SNP locus/loci.
Unable to determine switching at 10 SNP locus/loci.
FUNCTION: snpgdsApartSelection
Sat Apr 13 05:41:53 2019	Chromosome 1, # of SNPs: 365
Sat Apr 13 05:41:53 2019	Chromosome 2, # of SNPs: 370
Sat Apr 13 05:41:53 2019	Chromosome 3, # of SNPs: 314
Sat Apr 13 05:41:53 2019	Chromosome 4, # of SNPs: 294
Sat Apr 13 05:41:53 2019	Chromosome 5, # of SNPs: 291
Sat Apr 13 05:41:53 2019	Chromosome 6, # of SNPs: 280
Sat Apr 13 05:41:53 2019	Chromosome 7, # of SNPs: 245
Sat Apr 13 05:41:53 2019	Chromosome 8, # of SNPs: 232
Sat Apr 13 05:41:53 2019	Chromosome 9, # of SNPs: 203
Sat Apr 13 05:41:53 2019	Chromosome 10, # of SNPs: 224
Sat Apr 13 05:41:53 2019	Chromosome 11, # of SNPs: 227
Sat Apr 13 05:41:53 2019	Chromosome 12, # of SNPs: 215
Sat Apr 13 05:41:53 2019	Chromosome 13, # of SNPs: 175
Sat Apr 13 05:41:53 2019	Chromosome 14, # of SNPs: 147
Sat Apr 13 05:41:53 2019	Chromosome 15, # of SNPs: 127
Sat Apr 13 05:41:53 2019	Chromosome 16, # of SNPs: 128
Sat Apr 13 05:41:53 2019	Chromosome 17, # of SNPs: 112
Sat Apr 13 05:41:53 2019	Chromosome 18, # of SNPs: 127
Sat Apr 13 05:41:53 2019	Chromosome 19, # of SNPs: 76
Sat Apr 13 05:41:53 2019	Chromosome 20, # of SNPs: 107
Sat Apr 13 05:41:53 2019	Chromosome 21, # of SNPs: 63
Sat Apr 13 05:41:53 2019	Chromosome 22, # of SNPs: 54
Sat Apr 13 05:41:53 2019	Chromosome 23, # of SNPs: 215
Total # of SNPs selected:4591
FUNCTION: snpgdsBED2GDS
Start snpgdsBED2GDS ...
	BED file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.bed.gz" in the SNP-major mode (Sample X SNP)
	FAM file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.fam.gz", DONE.
	BIM file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.bim.gz", DONE.
Sat Apr 13 05:41:53 2019 	store sample id, snp id, position, and chromosome.
	start writing: 60 samples, 5000 SNPs ...
 	Sat Apr 13 05:41:53 2019	0%
 	Sat Apr 13 05:41:53 2019	100%
Sat Apr 13 05:41:53 2019 	Done.
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'HapMap.gds' (98.1K)
    # of fragments: 38
    save to 'HapMap.gds.tmp'
    rename 'HapMap.gds.tmp' (97.8K, reduced: 240B)
    # of fragments: 18
FUNCTION: snpgdsClose
FUNCTION: snpgdsCombineGeno
Create a GDS genotype file:
The new dataset consists of 10 samples and 3000 SNPs
    write sample.id
    write snp.id
    write snp.rs.id
    write snp.position
    write snp.chromosome
    write snp.allele
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Create a GDS genotype file:
The new dataset consists of 20 samples and 3000 SNPs
    write sample.id
    write snp.id
    write snp.rs.id
    write snp.position
    write snp.chromosome
    write snp.allele
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Merge SNP GDS files:
    open 't1.gds' ...
        10 samples, 3000 SNPs
    open 't2.gds' ...
        20 samples, 3000 SNPs
Concatenating samples (mapping to the first GDS file) ...
    reference: 3000 SNPs (100.0%)
    file 2: 0 allele flips, 0 ambiguous locus/loci
        [no flip]: 3000
    create 'test.gds': 30 samples, 3000 SNPs
    FileFormat = SNP_ARRAY
    writing genotypes ...
Clean up the fragments of GDS file:
    open the file 'test.gds' (46.2K)
    # of fragments: 32
    save to 'test.gds.tmp'
    rename 'test.gds.tmp' (46.0K, reduced: 204B)
    # of fragments: 15
Done.
Create a GDS genotype file:
The new dataset consists of 279 samples and 100 SNPs
    write sample.id
    write snp.id
    write snp.rs.id
    write snp.position
    write snp.chromosome
    write snp.allele
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Create a GDS genotype file:
The new dataset consists of 279 samples and 200 SNPs
    write sample.id
    write snp.id
    write snp.rs.id
    write snp.position
    write snp.chromosome
    write snp.allele
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Merge SNP GDS files:
    open 't1.gds' ...
        279 samples, 100 SNPs
    open 't2.gds' ...
        279 samples, 200 SNPs
Concatenating SNPs ...
    create 'test.gds': 279 samples, 300 SNPs
    FileFormat = SNP_ARRAY
    writing genotypes ...
Clean up the fragments of GDS file:
    open the file 'test.gds' (19.1K)
    # of fragments: 32
    save to 'test.gds.tmp'
    rename 'test.gds.tmp' (18.9K, reduced: 204B)
    # of fragments: 15
Done.
FUNCTION: snpgdsCreateGeno
Principal Component Analysis (PCA) on genotypes:
Excluding 42 SNPs on non-autosomes
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 958 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 264760
CPU capabilities:
Sat Apr 13 05:41:53 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:41:54 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:41:54 2019    Done.
FUNCTION: snpgdsCreateGenoSet
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 75.42%, 540/716
Chromosome 2: 72.24%, 536/742
Chromosome 3: 74.71%, 455/609
Chromosome 4: 73.31%, 412/562
Chromosome 5: 77.03%, 436/566
Chromosome 6: 75.58%, 427/565
Chromosome 7: 75.42%, 356/472
Chromosome 8: 71.31%, 348/488
Chromosome 9: 77.88%, 324/416
Chromosome 10: 74.33%, 359/483
Chromosome 11: 77.40%, 346/447
Chromosome 12: 76.81%, 328/427
Chromosome 13: 75.58%, 260/344
Chromosome 14: 76.95%, 217/282
Chromosome 15: 76.34%, 200/262
Chromosome 16: 72.66%, 202/278
Chromosome 17: 74.40%, 154/207
Chromosome 18: 73.68%, 196/266
Chromosome 19: 85.00%, 102/120
Chromosome 20: 71.62%, 164/229
Chromosome 21: 76.98%, 97/126
Chromosome 22: 75.86%, 88/116
6,547 markers are selected in total.
Create a GDS genotype file:
The new dataset consists of 279 samples and 6547 SNPs
    write sample.id
    write snp.id
    write snp.rs.id
    write snp.position
    write snp.chromosome
    write snp.allele
SNP genotypes are stored in SNP-major mode (Sample X SNP).
FUNCTION: snpgdsCutTree
Individual dissimilarity analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Dissimilarity:    the sum of all selected genotypes (0,1,2) = 2446510
Dissimilarity:	Sat Apr 13 05:41:54 2019	0%
Dissimilarity:	Sat Apr 13 05:41:56 2019	100%
Determine groups by permutation (Z threshold: 15, outlier threshold: 5):
Create 3 groups.
Create 4 groups.
FUNCTION: snpgdsDiss
Individual dissimilarity analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Dissimilarity:    the sum of all selected genotypes (0,1,2) = 2446510
Dissimilarity:	Sat Apr 13 05:41:57 2019	0%
Dissimilarity:	Sat Apr 13 05:41:58 2019	100%
Determine groups by permutation (Z threshold: 15, outlier threshold: 5):
Create 3 groups.
FUNCTION: snpgdsDrawTree
Individual dissimilarity analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Dissimilarity:    the sum of all selected genotypes (0,1,2) = 2446510
Dissimilarity:	Sat Apr 13 05:41:59 2019	0%
Dissimilarity:	Sat Apr 13 05:42:01 2019	100%
Determine groups by permutation (Z threshold: 15, outlier threshold: 5):
Create 3 groups.
FUNCTION: snpgdsEIGMIX
Eigen-analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Eigen-analysis:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities:
Sat Apr 13 05:42:01 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:02 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:42:02 2019    Done.
FUNCTION: snpgdsErrMsg
FUNCTION: snpgdsExampleFileName
FUNCTION: snpgdsFst
Fst estimation on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
Method: Weir & Cockerham, 1984
# of Populations: 4
    CEU (92), HCB (47), JPT (47), YRI (93)
Fst estimation on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
Method: Weir & Hill, 2002
# of Populations: 4
    CEU (92), HCB (47), JPT (47), YRI (93)
FUNCTION: snpgdsGDS2BED
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: 0.95)
Converting from GDS to PLINK binary PED:
Working space: 279 samples, 8722 SNPs
Output a BIM file.
Output a BED file ...
		Sat Apr 13 05:42:03 2019	0%
		Sat Apr 13 05:42:03 2019	100%
Done.
FUNCTION: snpgdsGDS2Eigen
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: 0.95)
Converting from GDS to EIGENSOFT:
	save to *.snp: 8722 snps
	save to *.ind: 279 samples
	Output: 	Sat Apr 13 05:42:03 2019	0%
	Output: 	Sat Apr 13 05:42:03 2019	100%
Done.
FUNCTION: snpgdsGDS2PED
Converting from GDS to PLINK PED:
	Output a MAP file DONE.
	Output a PED file ...
		Output: 	Sat Apr 13 05:42:03 2019	0%
		Output: 	Sat Apr 13 05:42:04 2019	100%
FUNCTION: snpgdsGEN2GDS
running snpgdsGEN2GDS ...
FUNCTION: snpgdsGRM
Genetic Relationship Matrix (GRM, GCTA):
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities:
Sat Apr 13 05:42:04 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:05 2019    Done.
Genetic Relationship Matrix (GRM, GCTA):
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities:
Sat Apr 13 05:42:05 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:06 2019    Done.
FUNCTION: snpgdsGetGeno
Genotype matrix: 1000 SNPs X 279 samples
Genotype matrix: 279 samples X 1000 SNPs
FUNCTION: snpgdsHCluster
Individual dissimilarity analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Dissimilarity:    the sum of all selected genotypes (0,1,2) = 2446510
Dissimilarity:	Sat Apr 13 05:42:06 2019	0%
Dissimilarity:	Sat Apr 13 05:42:08 2019	100%
Determine groups by permutation (Z threshold: 15, outlier threshold: 5):
Create 3 groups.
FUNCTION: snpgdsHWE
Keeping 716 SNPs according to chromosome 1
Excluding 160 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
FUNCTION: snpgdsIBDKING
IBD analysis (KING method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
No family is specified, and all individuals are treated as singletons.
Relationship inference in the presence of population stratification.
KING IBD:    the sum of all selected genotypes (0,1,2) = 702139
CPU capabilities:
Sat Apr 13 05:42:09 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:09 2019    Done.
IBD analysis (KING method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
No family is specified, and all individuals are treated as singletons.
Relationship inference in the presence of population stratification.
KING IBD:    the sum of all selected genotypes (0,1,2) = 702139
CPU capabilities:
Sat Apr 13 05:42:09 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:10 2019    Done.
IBD analysis (KING method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
# of families: 20, and within- and between-family relationship are estimated differently.
Relationship inference in the presence of population stratification.
KING IBD:    the sum of all selected genotypes (0,1,2) = 702139
CPU capabilities:
Sat Apr 13 05:42:11 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:11 2019    Done.
IBD analysis (KING method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
Relationship inference in a homogeneous population.
KING IBD:    the sum of all selected genotypes (0,1,2) = 702139
Sat Apr 13 05:42:11 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:12 2019    Done.
IBD analysis (KING method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
Relationship inference in a homogeneous population.
KING IBD:    the sum of all selected genotypes (0,1,2) = 702139
Sat Apr 13 05:42:12 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:12 2019    Done.
FUNCTION: snpgdsIBDMLE
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1,581 SNPs (monomorphic: TRUE, MAF: 0.05, missing rate: 0.05)
Working space: 30 samples, 7,142 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 54.75%, 392/716
Chromosome 2: 54.45%, 404/742
Chromosome 3: 55.67%, 339/609
Chromosome 4: 56.94%, 320/562
Chromosome 5: 56.71%, 321/566
Chromosome 6: 52.74%, 298/565
Chromosome 7: 55.72%, 263/472
Chromosome 8: 50.41%, 246/488
Chromosome 9: 54.33%, 226/416
Chromosome 10: 50.10%, 242/483
Chromosome 11: 54.81%, 245/447
Chromosome 12: 55.27%, 236/427
Chromosome 13: 53.49%, 184/344
Chromosome 14: 55.67%, 157/282
Chromosome 15: 54.58%, 143/262
Chromosome 16: 54.68%, 152/278
Chromosome 17: 55.56%, 115/207
Chromosome 18: 55.64%, 148/266
Chromosome 19: 66.67%, 80/120
Chromosome 20: 52.84%, 121/229
Chromosome 21: 50.79%, 64/126
Chromosome 22: 50.86%, 59/116
4,755 markers are selected in total.
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 30 samples, 250 SNPs
    using 1 (CPU) core
MLE IBD:    the sum of all selected genotypes (0,1,2) = 7222
MLE IBD:	Sat Apr 13 05:42:12 2019	0%
MLE IBD:	Sat Apr 13 05:42:14 2019	100%
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.482, sd: 0.283
MLE IBD:    the sum of all selected genotypes (0,1,2) = 6018
MLE IBD:	Sat Apr 13 05:42:14 2019	0%
MLE IBD:	Sat Apr 13 05:42:14 2019	100%
FUNCTION: snpgdsIBDMLELogLik
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1,581 SNPs (monomorphic: TRUE, MAF: 0.05, missing rate: 0.05)
Working space: 30 samples, 7,142 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 54.75%, 392/716
Chromosome 2: 54.45%, 404/742
Chromosome 3: 55.67%, 339/609
Chromosome 4: 56.94%, 320/562
Chromosome 5: 56.71%, 321/566
Chromosome 6: 52.74%, 298/565
Chromosome 7: 55.72%, 263/472
Chromosome 8: 50.41%, 246/488
Chromosome 9: 54.33%, 226/416
Chromosome 10: 50.10%, 242/483
Chromosome 11: 54.81%, 245/447
Chromosome 12: 55.27%, 236/427
Chromosome 13: 53.49%, 184/344
Chromosome 14: 55.67%, 157/282
Chromosome 15: 54.58%, 143/262
Chromosome 16: 54.68%, 152/278
Chromosome 17: 55.56%, 115/207
Chromosome 18: 55.64%, 148/266
Chromosome 19: 66.67%, 80/120
Chromosome 20: 52.84%, 121/229
Chromosome 21: 50.79%, 64/126
Chromosome 22: 50.86%, 59/116
4,755 markers are selected in total.
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 30 samples, 250 SNPs
    using 1 (CPU) core
MLE IBD:    the sum of all selected genotypes (0,1,2) = 7222
MLE IBD:	Sat Apr 13 05:42:14 2019	0%
MLE IBD:	Sat Apr 13 05:42:15 2019	100%
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.482, sd: 0.283
MLE IBD:    the sum of all selected genotypes (0,1,2) = 6018
MLE IBD:	Sat Apr 13 05:42:15 2019	0%
MLE IBD:	Sat Apr 13 05:42:15 2019	100%
FUNCTION: snpgdsIBDMoM
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 702139
Sat Apr 13 05:42:15 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:16 2019    Done.
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 563 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 93 samples, 8,160 SNPs
    using 1 (CPU) core
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 755648
Sat Apr 13 05:42:16 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:16 2019    Done.
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 563 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 93 samples, 8,160 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.500, sd: 0.315
*** A correction factor based on allele count is not used, since the allele frequencies are specified.
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 755648
Sat Apr 13 05:42:16 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:16 2019    Done.
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 563 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 8,160 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.500, sd: 0.315
*** A correction factor based on allele count is not used, since the allele frequencies are specified.
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 203285
Sat Apr 13 05:42:16 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:16 2019    Done.
FUNCTION: snpgdsIBDSelection
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 563 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 93 samples, 8,160 SNPs
    using 1 (CPU) core
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 755648
Sat Apr 13 05:42:16 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:16 2019    Done.
FUNCTION: snpgdsIBS
Identity-By-State (IBS) analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
IBS:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:42:16 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 2s
Sat Apr 13 05:42:18 2019    Done.
FUNCTION: snpgdsIBSNum
Identity-By-State (IBS) analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
IBS:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:42:18 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:19 2019    Done.
FUNCTION: snpgdsIndInb
Estimating individual inbreeding coefficients:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
FUNCTION: snpgdsIndInbCoef
FUNCTION: snpgdsIndivBeta
Individual Inbreeding and Relatedness (beta estimator):
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Individual Beta:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities:
Sat Apr 13 05:42:20 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:21 2019    Done.
FUNCTION: snpgdsLDMat
Linkage Disequilibrium (LD) estimation on genotypes:
Working space: 279 samples, 203 SNPs
    using 1 (CPU) core.
    method: composite
LD matrix:    the sum of all selected genotypes (0,1,2) = 56582
Linkage Disequilibrium (LD) estimation on genotypes:
Working space: 279 samples, 203 SNPs
    using 1 (CPU) core.
    sliding window size: 203 
    method: composite
LD matrix:    the sum of all selected genotypes (0,1,2) = 56582
FUNCTION: snpgdsLDpair
FUNCTION: snpgdsLDpruning
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 75.42%, 540/716
Chromosome 2: 72.24%, 536/742
Chromosome 3: 74.71%, 455/609
Chromosome 4: 73.31%, 412/562
Chromosome 5: 77.03%, 436/566
Chromosome 6: 75.58%, 427/565
Chromosome 7: 75.42%, 356/472
Chromosome 8: 71.31%, 348/488
Chromosome 9: 77.88%, 324/416
Chromosome 10: 74.33%, 359/483
Chromosome 11: 77.40%, 346/447
Chromosome 12: 76.81%, 328/427
Chromosome 13: 75.58%, 260/344
Chromosome 14: 76.95%, 217/282
Chromosome 15: 76.34%, 200/262
Chromosome 16: 72.66%, 202/278
Chromosome 17: 74.40%, 154/207
Chromosome 18: 73.68%, 196/266
Chromosome 19: 85.00%, 102/120
Chromosome 20: 71.62%, 164/229
Chromosome 21: 76.98%, 97/126
Chromosome 22: 75.86%, 88/116
6,547 markers are selected in total.
FUNCTION: snpgdsMergeGRM
Genetic Relationship Matrix (GRM, GCTA):
Excluding 2,288 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 6,800 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 1908966
CPU capabilities:
Sat Apr 13 05:42:22 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:22 2019    Done.
Genetic Relationship Matrix (GRM, GCTA):
Excluding 5,688 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 3,400 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 951558
CPU capabilities:
Sat Apr 13 05:42:22 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:23 2019    Done.
Genetic Relationship Matrix (GRM, GCTA):
Excluding 5,688 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 3,400 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 957408
CPU capabilities:
Sat Apr 13 05:42:23 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:24 2019    Done.
GRM merging:
    open 'tmp1.gds' (3,400 variants)
    open 'tmp2.gds' (3,400 variants)
Weight: 0.5, 0.5
Output: tmp.gds

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
GRM merging:
    open 'tmp1.gds' (3,400 variants)
    open 'tmp2.gds' (3,400 variants)
Weight: 0.5, 0.5

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
FUNCTION: snpgdsOpen
FUNCTION: snpgdsOption
FUNCTION: snpgdsPCA
Principal Component Analysis (PCA) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities:
Sat Apr 13 05:42:25 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:26 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:42:26 2019    Done.
FUNCTION: snpgdsPCACorr
Principal Component Analysis (PCA) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities:
Sat Apr 13 05:42:26 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:26 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:42:27 2019    Done.
SNP Correlation:
Working space: 279 samples, 9088 SNPs
    using 1 (CPU) core
    using the top 4 eigenvectors
Correlation:    the sum of all selected genotypes (0,1,2) = 2553065
Sat Apr 13 05:42:27 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:27 2019    Done.
SNP Correlation:
Working space: 279 samples, 9088 SNPs
    using 1 (CPU) core
    using the top 4 eigenvectors
Creating 'test.gds' ...
Correlation:    the sum of all selected genotypes (0,1,2) = 2553065
Sat Apr 13 05:42:27 2019

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:27 2019    Done.
FUNCTION: snpgdsPCASNPLoading
Principal Component Analysis (PCA) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities:
Sat Apr 13 05:42:27 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:27 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:42:27 2019    Done.
SNP loading:
Working space: 279 samples, 8722 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
SNP Loading:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:42:27 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:27 2019    Done.
FUNCTION: snpgdsPCASampLoading
Principal Component Analysis (PCA) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities:
Sat Apr 13 05:42:28 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:28 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:42:28 2019    Done.
SNP loading:
Working space: 279 samples, 8722 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
SNP Loading:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:42:28 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:28 2019    Done.
Sample loading:
Working space: 100 samples, 8722 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
Sample Loading:    the sum of all selected genotypes (0,1,2) = 878146
Sat Apr 13 05:42:28 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:28 2019    Done.
FUNCTION: snpgdsPED2GDS
Converting from GDS to PLINK PED:
	Output a MAP file DONE.
	Output a PED file ...
		Output: 	Sat Apr 13 05:42:28 2019	0%
		Output: 	Sat Apr 13 05:42:29 2019	100%
PLINK PED/MAP to GDS Format:
Import 9088 variants from 'tmp.map'
Chromosome:
  1  10  11  12  13  14  15  16  17  18  19   2  20  21  22   3   4   5   6   7 
716 483 447 427 344 282 262 278 207 266 120 742 229 126 116 609 562 566 565 472 
  8   9   X 
488 416 365 
Reading 'tmp.ped'
Output: 'test.gds'
Import 279 samples
Transpose the genotypic matrix ...
Done.
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test.gds' (1.3M)
    # of fragments: 50
    save to 'test.gds.tmp'
    rename 'test.gds.tmp' (711.4K, reduced: 618.7K)
    # of fragments: 26
FUNCTION: snpgdsPairIBD
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1,646 SNPs (monomorphic: TRUE, MAF: 0.05, missing rate: 0.05)
Working space: 93 samples, 7,077 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 62.29%, 446/716
Chromosome 2: 62.40%, 463/742
Chromosome 3: 60.76%, 370/609
Chromosome 4: 64.59%, 363/562
Chromosome 5: 62.54%, 354/566
Chromosome 6: 60.00%, 339/565
Chromosome 7: 62.71%, 296/472
Chromosome 8: 58.61%, 286/488
Chromosome 9: 62.98%, 262/416
Chromosome 10: 60.66%, 293/483
Chromosome 11: 62.86%, 281/447
Chromosome 12: 63.00%, 269/427
Chromosome 13: 63.08%, 217/344
Chromosome 14: 63.48%, 179/282
Chromosome 15: 63.36%, 166/262
Chromosome 16: 61.87%, 172/278
Chromosome 17: 65.70%, 136/207
Chromosome 18: 59.02%, 157/266
Chromosome 19: 68.33%, 82/120
Chromosome 20: 67.25%, 154/229
Chromosome 21: 61.11%, 77/126
Chromosome 22: 57.76%, 67/116
5,429 markers are selected in total.
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.506, sd: 0.300
MLE IBD:    the sum of all selected genotypes (0,1,2) = 6339
MLE IBD:	Sat Apr 13 05:42:30 2019	0%
MLE IBD:	Sat Apr 13 05:42:31 2019	100%
IBD analysis (PLINK method of moment) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.506, sd: 0.300
*** A correction factor based on allele count is not used, since the allele frequencies are specified.
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 6339
Sat Apr 13 05:42:31 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:31 2019    Done.
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.506, sd: 0.300
MLE IBD:    the sum of all selected genotypes (0,1,2) = 6339
MLE IBD:	Sat Apr 13 05:42:31 2019	0%
MLE IBD:	Sat Apr 13 05:42:31 2019	100%
Genotype matrix: 250 SNPs X 25 samples
[1] -364.4667
[1] -377.9771
[1] -390.3299
[1] -380.5342
[1] -378.2092
[1] -383.641
[1] -364.5406
[1] -375.196
[1] -383.5752
[1] -385.7771
[1] -368.6237
[1] -390.7766
[1] -367.5082
[1] -383.595
[1] -385.7602
[1] -378.5134
[1] -388.3292
[1] -376.4419
[1] -378.3509
[1] -386.8318
[1] -385.3803
[1] -382.988
[1] -369.9328
[1] -383.3912
FUNCTION: snpgdsPairIBDMLELogLik
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1,646 SNPs (monomorphic: TRUE, MAF: 0.05, missing rate: 0.05)
Working space: 93 samples, 7,077 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 62.29%, 446/716
Chromosome 2: 62.40%, 463/742
Chromosome 3: 60.76%, 370/609
Chromosome 4: 64.59%, 363/562
Chromosome 5: 62.54%, 354/566
Chromosome 6: 60.00%, 339/565
Chromosome 7: 62.71%, 296/472
Chromosome 8: 58.61%, 286/488
Chromosome 9: 62.98%, 262/416
Chromosome 10: 60.66%, 293/483
Chromosome 11: 62.86%, 281/447
Chromosome 12: 63.00%, 269/427
Chromosome 13: 63.08%, 217/344
Chromosome 14: 63.48%, 179/282
Chromosome 15: 63.36%, 166/262
Chromosome 16: 61.87%, 172/278
Chromosome 17: 65.70%, 136/207
Chromosome 18: 59.02%, 157/266
Chromosome 19: 68.33%, 82/120
Chromosome 20: 67.25%, 154/229
Chromosome 21: 61.11%, 77/126
Chromosome 22: 57.76%, 67/116
5,429 markers are selected in total.
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.506, sd: 0.300
MLE IBD:    the sum of all selected genotypes (0,1,2) = 6339
MLE IBD:	Sat Apr 13 05:42:32 2019	0%
MLE IBD:	Sat Apr 13 05:42:33 2019	100%
IBD analysis (PLINK method of moment) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.506, sd: 0.300
*** A correction factor based on allele count is not used, since the allele frequencies are specified.
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 6339
Sat Apr 13 05:42:33 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:33 2019    Done.
Genotype matrix: 250 SNPs X 25 samples
[1] -364.4667
[1] -377.9771
[1] -390.3299
[1] -380.5342
[1] -378.2092
[1] -383.641
[1] -364.5406
[1] -375.196
[1] -383.5752
[1] -385.7771
[1] -368.6237
[1] -390.7766
[1] -367.5082
[1] -383.595
[1] -385.7602
[1] -378.5134
[1] -388.3292
[1] -376.4419
[1] -378.3509
[1] -386.8318
[1] -385.3803
[1] -382.988
[1] -369.9328
[1] -383.3912
FUNCTION: snpgdsPairScore
Excluding 365 SNPs on non-autosomes
Working space: 120 samples, 8723 SNPs
Method: IBS
Genotype Score:    the sum of all selected genotypes (0,1,2) = 1050236
Working space: 120 samples, 8723 SNPs
Method: IBS
Genotype Score:    the sum of all selected genotypes (0,1,2) = 1050236
Working space: 120 samples, 8723 SNPs
Method: IBS
Genotype Score:    the sum of all selected genotypes (0,1,2) = 1050236
Working space: 120 samples, 8723 SNPs
Method: IBS
Genotype Score:    the sum of all selected genotypes (0,1,2) = 1050236
Working space: 120 samples, 8723 SNPs
Method: IBS
Output: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\tmp.gds
Genotype Score:    the sum of all selected genotypes (0,1,2) = 1050236
FUNCTION: snpgdsSNPList
FUNCTION: snpgdsSNPListClass
FUNCTION: snpgdsSNPListIntersect
FUNCTION: snpgdsSNPRateFreq
FUNCTION: snpgdsSampMissRate
FUNCTION: snpgdsSelectSNP
Excluding 365 SNPs on non-autosomes
Excluding 1,221 SNPs (monomorphic: TRUE, MAF: 0.05, missing rate: 0.95)
FUNCTION: snpgdsSlidingWindow
Sliding Window Analysis:
Excluding 8 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 9,080 SNPs
    using 1 (CPU) core
    window size: 500000, shift: 100000 (basepair)
Chromosome Set: 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23
Sat Apr 13 05:42:33 2019, Chromosome 1 (716 SNPs), 2448 windows
Sat Apr 13 05:42:33 2019, Chromosome 2 (742 SNPs), 2416 windows
Sat Apr 13 05:42:34 2019, Chromosome 3 (609 SNPs), 1985 windows
Sat Apr 13 05:42:34 2019, Chromosome 4 (562 SNPs), 1894 windows
Sat Apr 13 05:42:34 2019, Chromosome 5 (566 SNPs), 1797 windows
Sat Apr 13 05:42:34 2019, Chromosome 6 (565 SNPs), 1694 windows
Sat Apr 13 05:42:34 2019, Chromosome 7 (472 SNPs), 1573 windows
Sat Apr 13 05:42:34 2019, Chromosome 8 (488 SNPs), 1445 windows
Sat Apr 13 05:42:34 2019, Chromosome 9 (416 SNPs), 1393 windows
Sat Apr 13 05:42:34 2019, Chromosome 10 (483 SNPs), 1343 windows
Sat Apr 13 05:42:34 2019, Chromosome 11 (447 SNPs), 1338 windows
Sat Apr 13 05:42:34 2019, Chromosome 12 (427 SNPs), 1316 windows
Sat Apr 13 05:42:34 2019, Chromosome 13 (344 SNPs), 948 windows
Sat Apr 13 05:42:34 2019, Chromosome 14 (281 SNPs), 847 windows
Sat Apr 13 05:42:34 2019, Chromosome 15 (262 SNPs), 774 windows
Sat Apr 13 05:42:34 2019, Chromosome 16 (278 SNPs), 873 windows
Sat Apr 13 05:42:34 2019, Chromosome 17 (207 SNPs), 773 windows
Sat Apr 13 05:42:34 2019, Chromosome 18 (266 SNPs), 753 windows
Sat Apr 13 05:42:34 2019, Chromosome 19 (120 SNPs), 627 windows
Sat Apr 13 05:42:34 2019, Chromosome 20 (229 SNPs), 602 windows
Sat Apr 13 05:42:34 2019, Chromosome 21 (126 SNPs), 311 windows
Sat Apr 13 05:42:34 2019, Chromosome 22 (116 SNPs), 312 windows
Sat Apr 13 05:42:34 2019, Chromosome 23 (358 SNPs), 1507 windows
Sat Apr 13 05:42:34 2019 	Done.
FUNCTION: snpgdsSummary
The file name: C:\Users\biocbuild\bbs-3.8-bioc\R\library\SNPRelate\extdata\hapmap_geno.gds 
The total number of samples: 279 
The total number of SNPs: 9088 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
FUNCTION: snpgdsTranspose
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test.gds 
The total number of samples: 279 
The total number of SNPs: 9088 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
SNP genotypes: 279 samples, 9088 SNPs
Genotype matrix is being transposed ...
Clean up the fragments of GDS file:
    open the file 'test.gds' (1.3M)
    # of fragments: 28
    save to 'test.gds.tmp'
    rename 'test.gds.tmp' (709.6K, reduced: 619.1K)
    # of fragments: 26
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test.gds 
The total number of samples: 279 
The total number of SNPs: 9088 
SNP genotypes are stored in individual-major mode (SNP X Sample).
FUNCTION: snpgdsVCF2GDS
##fileformat=VCFv4.1
##fileDate=20090805
##source=myImputationProgramV3.1
##reference=file:///seq/references/1000GenomesPilot-NCBI36.fasta
##contig=<ID=20,length=62435964,assembly=B36,md5=f126cdf8a6e0c7f379d618ff66beb2da,species="Homo sapiens",taxonomy=x>
##phasing=partial
##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data">
##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth">
##INFO=<ID=AF,Number=A,Type=Float,Description="Allele Frequency">
##INFO=<ID=AA,Number=1,Type=String,Description="Ancestral Allele">
##INFO=<ID=DB,Number=0,Type=Flag,Description="dbSNP membership, build 129">
##INFO=<ID=H2,Number=0,Type=Flag,Description="HapMap2 membership">
##FILTER=<ID=q10,Description="Quality below 10">
##FILTER=<ID=s50,Description="Less than 50% of samples have data">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality">
##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">
##FORMAT=<ID=HQ,Number=2,Type=Integer,Description="Haplotype Quality">
#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	NA00001	NA00002	NA00003
20	14370	rs6054257	G	A	29	PASS	NS=3;DP=14;AF=0.5;DB;H2	GT:GQ:DP:HQ	0|0:48:1:51,51	1|0:48:8:51,51	1/1:43:5:.,.
20	17330	.	T	A	3	q10	NS=3;DP=11;AF=0.017	GT:GQ:DP:HQ	0|0:49:3:58,50	0|1:3:5:65,3	0/0:41:3
20	1110696	rs6040355	A	G,T	67	PASS	NS=2;DP=10;AF=0.333,0.667;AA=T;DB	GT:GQ:DP:HQ	1|2:21:6:23,27	2|1:2:0:18,2	2/2:35:4
20	1230237	.	T	.	47	PASS	NS=3;DP=13;AA=T	GT:GQ:DP:HQ	0|0:54:7:56,60	0|0:48:4:51,51	0/0:61:2
20	1234567	microsat1	GTC	G,GTCT	50	PASS	NS=3;DP=9;AA=G	GT:GQ:DP	0/1:35:4	0/2:17:2	1/1:40:3
VCF Format ==> SNP GDS Format
Method: exacting biallelic SNPs
Number of samples: 3
Parsing "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf" ...
	import 2 variants.
+ genotype   { Bit2 3x2, 2B } *
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test1.gds' (2.9K)
    # of fragments: 46
    save to 'test1.gds.tmp'
    rename 'test1.gds.tmp' (2.6K, reduced: 312B)
    # of fragments: 20
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test1.gds 
The total number of samples: 3 
The total number of SNPs: 2 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
VCF Format ==> SNP GDS Format
Method: exacting biallelic SNPs
Number of samples: 3
Parsing "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf" ...
	import 2 variants.
+ genotype   { Bit2 3x2, 2B } *
SNP genotypes: 3 samples, 2 SNPs
Genotype matrix is being transposed ...
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test2.gds' (3.0K)
    # of fragments: 48
    save to 'test2.gds.tmp'
    rename 'test2.gds.tmp' (2.6K, reduced: 417B)
    # of fragments: 20
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test2.gds 
The total number of samples: 3 
The total number of SNPs: 2 
SNP genotypes are stored in individual-major mode (SNP X Sample).
VCF Format ==> SNP GDS Format
Method: dosage (0,1,2) of reference allele for all variant sites
Number of samples: 3
Parsing "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf" ...
	import 5 variants.
+ genotype   { Bit2 3x5, 4B } *
SNP genotypes: 3 samples, 5 SNPs
Genotype matrix is being transposed ...
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test3.gds' (3.1K)
    # of fragments: 48
    save to 'test3.gds.tmp'
    rename 'test3.gds.tmp' (2.7K, reduced: 419B)
    # of fragments: 20
Some of 'snp.allele' are not standard (e.g., A/G,T).
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test3.gds 
The total number of samples: 3 
The total number of SNPs: 5 
SNP genotypes are stored in individual-major mode (SNP X Sample).
The number of valid samples: 3 
The number of biallelic unique SNPs: 2 
VCF Format ==> SNP GDS Format
Method: dosage (0,1,2) of reference allele for all variant sites
Number of samples: 3
Parsing "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf" ...
	import 5 variants.
+ genotype   { Bit2 3x5, 4B } *
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test4.gds' (3.0K)
    # of fragments: 46
    save to 'test4.gds.tmp'
    rename 'test4.gds.tmp' (2.7K, reduced: 312B)
    # of fragments: 20
Some of 'snp.allele' are not standard (e.g., A/G,T).
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test4.gds 
The total number of samples: 3 
The total number of SNPs: 5 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
The number of valid samples: 3 
The number of biallelic unique SNPs: 2 
VCF Format ==> SNP GDS Format
Method: dosage (0,1,2) of reference allele for all variant sites
Number of samples: 3
Parsing "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf" ...
	import 5 variants.
+ genotype   { Bit2 3x5, 4B } *
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test5.gds' (3.0K)
    # of fragments: 46
    save to 'test5.gds.tmp'
    rename 'test5.gds.tmp' (2.7K, reduced: 312B)
    # of fragments: 20
Some of 'snp.allele' are not standard (e.g., T/A,G).
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test5.gds 
The total number of samples: 3 
The total number of SNPs: 5 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
The number of valid samples: 3 
The number of biallelic unique SNPs: 2 
FUNCTION: snpgdsVCF2GDS_R
##fileformat=VCFv4.1
##fileDate=20090805
##source=myImputationProgramV3.1
##reference=file:///seq/references/1000GenomesPilot-NCBI36.fasta
##contig=<ID=20,length=62435964,assembly=B36,md5=f126cdf8a6e0c7f379d618ff66beb2da,species="Homo sapiens",taxonomy=x>
##phasing=partial
##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data">
##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth">
##INFO=<ID=AF,Number=A,Type=Float,Description="Allele Frequency">
##INFO=<ID=AA,Number=1,Type=String,Description="Ancestral Allele">
##INFO=<ID=DB,Number=0,Type=Flag,Description="dbSNP membership, build 129">
##INFO=<ID=H2,Number=0,Type=Flag,Description="HapMap2 membership">
##FILTER=<ID=q10,Description="Quality below 10">
##FILTER=<ID=s50,Description="Less than 50% of samples have data">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality">
##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">
##FORMAT=<ID=HQ,Number=2,Type=Integer,Description="Haplotype Quality">
#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	NA00001	NA00002	NA00003
20	14370	rs6054257	G	A	29	PASS	NS=3;DP=14;AF=0.5;DB;H2	GT:GQ:DP:HQ	0|0:48:1:51,51	1|0:48:8:51,51	1/1:43:5:.,.
20	17330	.	T	A	3	q10	NS=3;DP=11;AF=0.017	GT:GQ:DP:HQ	0|0:49:3:58,50	0|1:3:5:65,3	0/0:41:3
20	1110696	rs6040355	A	G,T	67	PASS	NS=2;DP=10;AF=0.333,0.667;AA=T;DB	GT:GQ:DP:HQ	1|2:21:6:23,27	2|1:2:0:18,2	2/2:35:4
20	1230237	.	T	.	47	PASS	NS=3;DP=13;AA=T	GT:GQ:DP:HQ	0|0:54:7:56,60	0|0:48:4:51,51	0/0:61:2
20	1234567	microsat1	GTC	G,GTCT	50	PASS	NS=3;DP=9;AA=G	GT:GQ:DP	0/1:35:4	0/2:17:2	1/1:40:3
Start snpgdsVCF2GDS ...
	Extracting bi-allelic and polymorhpic SNPs.
	Scanning ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
	content: 5 rows x 12 columns
Sat Apr 13 05:42:34 2019 	store sample id, snp id, position, and chromosome.
	start writing: 3 samples, 2 SNPs ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
[1] 1
Sat Apr 13 05:42:35 2019 	Done.
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test1.gds 
The total number of samples: 3 
The total number of SNPs: 2 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Start snpgdsVCF2GDS ...
	Extracting bi-allelic and polymorhpic SNPs.
	Scanning ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
	content: 5 rows x 12 columns
Sat Apr 13 05:42:35 2019 	store sample id, snp id, position, and chromosome.
	start writing: 3 samples, 2 SNPs ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
[1] 1
Sat Apr 13 05:42:35 2019 	Done.
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test2.gds 
The total number of samples: 3 
The total number of SNPs: 2 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Start snpgdsVCF2GDS ...
	Storing dosage of the reference allele for all variant sites, including bi-allelic SNPs, multi-allelic SNPs, indels and structural variants.
	Scanning ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
	content: 5 rows x 12 columns
Sat Apr 13 05:42:35 2019 	store sample id, snp id, position, and chromosome.
	start writing: 3 samples, 5 SNPs ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
Sat Apr 13 05:42:35 2019 	Done.
Some of 'snp.allele' are not standard (e.g., A/G,T).
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test3.gds 
The total number of samples: 3 
The total number of SNPs: 5 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
The number of valid samples: 3 
The number of biallelic unique SNPs: 2 
Start snpgdsVCF2GDS ...
	Storing dosage of the reference allele for all variant sites, including bi-allelic SNPs, multi-allelic SNPs, indels and structural variants.
	Scanning ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
	content: 5 rows x 12 columns
Sat Apr 13 05:42:35 2019 	store sample id, snp id, position, and chromosome.
	start writing: 3 samples, 5 SNPs ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
Sat Apr 13 05:42:35 2019 	Done.
Some of 'snp.allele' are not standard (e.g., A/G,T).
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_i386\test4.gds 
The total number of samples: 3 
The total number of SNPs: 5 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
The number of valid samples: 3 
The number of biallelic unique SNPs: 2 
SNP Correlation:
Working space: 90 samples, 9088 SNPs
    using 1 (CPU) core
    using the top 2 eigenvectors
Correlation:    the sum of all selected genotypes (0,1,2) = 824424
Sat Apr 13 05:42:37 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:37 2019    Done.
SNP Correlation:
Working space: 90 samples, 9088 SNPs
    using 1 (CPU) core
    using the top 2 eigenvectors
Creating 'test.gds' ...
Correlation:    the sum of all selected genotypes (0,1,2) = 824424
Sat Apr 13 05:42:37 2019

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:37 2019    Done.
SNP loading:
Working space: 90 samples, 8695 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
SNP Loading:    the sum of all selected genotypes (0,1,2) = 787449
Sat Apr 13 05:42:37 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:37 2019    Done.
Sample loading:
Working space: 100 samples, 8695 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
Sample Loading:    the sum of all selected genotypes (0,1,2) = 875255
Sat Apr 13 05:42:37 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:37 2019    Done.
SNP Correlation:
Working space: 90 samples, 9088 SNPs
    using 2 (CPU) cores
    using the top 2 eigenvectors
Correlation:    the sum of all selected genotypes (0,1,2) = 824424
Sat Apr 13 05:42:37 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:37 2019    Done.
SNP Correlation:
Working space: 90 samples, 9088 SNPs
    using 2 (CPU) cores
    using the top 2 eigenvectors
Creating 'test.gds' ...
Correlation:    the sum of all selected genotypes (0,1,2) = 824424
Sat Apr 13 05:42:37 2019

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:37 2019    Done.
SNP loading:
Working space: 90 samples, 8695 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
SNP Loading:    the sum of all selected genotypes (0,1,2) = 787449
Sat Apr 13 05:42:37 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:37 2019    Done.
Sample loading:
Working space: 100 samples, 8695 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
Sample Loading:    the sum of all selected genotypes (0,1,2) = 875255
Sat Apr 13 05:42:37 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:37 2019    Done.


RUNIT TEST PROTOCOL -- Sat Apr 13 05:42:38 2019 
*********************************************** 
Number of test functions: 13 
Number of errors: 0 
Number of failures: 0 

 
1 Test Suite : 
SNPRelate RUnit Tests - 13 test functions, 0 errors, 0 failures
Number of test functions: 13 
Number of errors: 0 
Number of failures: 0 
> 
> proc.time()
   user  system elapsed 
  65.21    2.10   67.64 

SNPRelate.Rcheck/tests_x64/runTests.Rout


R version 3.5.3 (2019-03-11) -- "Great Truth"
Copyright (C) 2019 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64 (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> BiocGenerics:::testPackage("SNPRelate")
SNPRelate -- supported by Streaming SIMD Extensions 2 (SSE2)
Genetic Relationship Matrix (GRM, GCTA):
Excluding 8,088 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 1,000 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 282597
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:39 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:40 2019    Done.
Genetic Relationship Matrix (GRM, GCTA):
Excluding 7,088 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 2,000 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 559412
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:40 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:41 2019    Done.
Genetic Relationship Matrix (GRM, GCTA):
Excluding 5,288 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 3,800 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 1066957
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:41 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:42 2019    Done.
GRM merging:
    open 'tmp1.gds' (1,000 variants)
    open 'tmp2.gds' (2,000 variants)
    open 'tmp3.gds' (3,800 variants)
Weight: 0.147059, 0.294118, 0.558824
Output: tmp.gds

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Genetic Relationship Matrix (GRM, GCTA):
Excluding 2,288 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 6,800 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 1908966
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:42 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:43 2019    Done.
Genetic Relationship Matrix (GRM, IndivBeta):
Excluding 8,088 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 1,000 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 282597
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:43 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:43 2019    Done.
Genetic Relationship Matrix (GRM, IndivBeta):
Excluding 7,088 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 2,000 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 559412
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:43 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:44 2019    Done.
Genetic Relationship Matrix (GRM, IndivBeta):
Excluding 5,288 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 3,800 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 1066957
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:44 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:45 2019    Done.
GRM merging:
    open 'tmp1.gds' (1,000 variants)
    open 'tmp2.gds' (2,000 variants)
    open 'tmp3.gds' (3,800 variants)
Weight: 0.147059, 0.294118, 0.558824
Output: tmp.gds

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Writing ...

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Genetic Relationship Matrix (GRM, IndivBeta):
Excluding 2,288 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 6,800 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 1908966
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:45 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:46 2019    Done.
Linkage Disequilibrium (LD) estimation on genotypes:
Working space: 279 samples, 1000 SNPs
    using 1 (CPU) core.
    method: covariance
LD matrix:    the sum of all selected genotypes (0,1,2) = 283058
Linkage Disequilibrium (LD) estimation on genotypes:
Working space: 279 samples, 1000 SNPs
    using 1 (CPU) core.
    method: correlation
LD matrix:    the sum of all selected genotypes (0,1,2) = 283058
FUNCTION: SNPGDSFileClass
FUNCTION: SNPRelate-package
Start snpgdsBED2GDS ...
	BED file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.bed.gz" in the SNP-major mode (Sample X SNP)
	FAM file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.fam.gz", DONE.
	BIM file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.bim.gz", DONE.
Sat Apr 13 05:42:50 2019 	store sample id, snp id, position, and chromosome.
	start writing: 60 samples, 5000 SNPs ...
 	Sat Apr 13 05:42:50 2019	0%
 	Sat Apr 13 05:42:50 2019	100%
Sat Apr 13 05:42:50 2019 	Done.
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'HapMap.gds' (98.1K)
    # of fragments: 38
    save to 'HapMap.gds.tmp'
    rename 'HapMap.gds.tmp' (97.8K, reduced: 240B)
    # of fragments: 18
Principal Component Analysis (PCA) on genotypes:
Excluding 203 SNPs on non-autosomes
Excluding 28 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 60 samples, 4,769 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 124273
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:50 2019    (internal increment: 53996)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:50 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:42:50 2019    Done.
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:42:50 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:42:51 2019    Done.
Identity-By-State (IBS) analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
IBS:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:42:51 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:51 2019    Done.
Linkage Disequilibrium (LD) estimation on genotypes:
Working space: 279 samples, 200 SNPs
    using 1 (CPU) core.
    method: composite
LD matrix:    the sum of all selected genotypes (0,1,2) = 55417
FUNCTION: hapmap_geno
FUNCTION: snpgdsAdmixPlot
Eigen-analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Eigen-analysis:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:52 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:52 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:42:53 2019    Done.
FUNCTION: snpgdsAdmixProp
Eigen-analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Eigen-analysis:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:53 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:53 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:42:53 2019    Done.
FUNCTION: snpgdsAlleleSwitch
Strand-switching at 50 SNP locus/loci.
Unable to determine switching at 10 SNP locus/loci.
FUNCTION: snpgdsApartSelection
Sat Apr 13 05:42:53 2019	Chromosome 1, # of SNPs: 365
Sat Apr 13 05:42:53 2019	Chromosome 2, # of SNPs: 370
Sat Apr 13 05:42:53 2019	Chromosome 3, # of SNPs: 314
Sat Apr 13 05:42:53 2019	Chromosome 4, # of SNPs: 294
Sat Apr 13 05:42:53 2019	Chromosome 5, # of SNPs: 291
Sat Apr 13 05:42:53 2019	Chromosome 6, # of SNPs: 280
Sat Apr 13 05:42:53 2019	Chromosome 7, # of SNPs: 245
Sat Apr 13 05:42:53 2019	Chromosome 8, # of SNPs: 232
Sat Apr 13 05:42:53 2019	Chromosome 9, # of SNPs: 203
Sat Apr 13 05:42:53 2019	Chromosome 10, # of SNPs: 224
Sat Apr 13 05:42:53 2019	Chromosome 11, # of SNPs: 227
Sat Apr 13 05:42:53 2019	Chromosome 12, # of SNPs: 215
Sat Apr 13 05:42:53 2019	Chromosome 13, # of SNPs: 175
Sat Apr 13 05:42:53 2019	Chromosome 14, # of SNPs: 147
Sat Apr 13 05:42:53 2019	Chromosome 15, # of SNPs: 127
Sat Apr 13 05:42:53 2019	Chromosome 16, # of SNPs: 128
Sat Apr 13 05:42:53 2019	Chromosome 17, # of SNPs: 112
Sat Apr 13 05:42:53 2019	Chromosome 18, # of SNPs: 127
Sat Apr 13 05:42:53 2019	Chromosome 19, # of SNPs: 76
Sat Apr 13 05:42:53 2019	Chromosome 20, # of SNPs: 107
Sat Apr 13 05:42:53 2019	Chromosome 21, # of SNPs: 63
Sat Apr 13 05:42:53 2019	Chromosome 22, # of SNPs: 54
Sat Apr 13 05:42:53 2019	Chromosome 23, # of SNPs: 215
Total # of SNPs selected:4591
FUNCTION: snpgdsBED2GDS
Start snpgdsBED2GDS ...
	BED file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.bed.gz" in the SNP-major mode (Sample X SNP)
	FAM file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.fam.gz", DONE.
	BIM file: "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/plinkhapmap.bim.gz", DONE.
Sat Apr 13 05:42:53 2019 	store sample id, snp id, position, and chromosome.
	start writing: 60 samples, 5000 SNPs ...
 	Sat Apr 13 05:42:53 2019	0%
 	Sat Apr 13 05:42:53 2019	100%
Sat Apr 13 05:42:53 2019 	Done.
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'HapMap.gds' (98.1K)
    # of fragments: 38
    save to 'HapMap.gds.tmp'
    rename 'HapMap.gds.tmp' (97.8K, reduced: 240B)
    # of fragments: 18
FUNCTION: snpgdsClose
FUNCTION: snpgdsCombineGeno
Create a GDS genotype file:
The new dataset consists of 10 samples and 3000 SNPs
    write sample.id
    write snp.id
    write snp.rs.id
    write snp.position
    write snp.chromosome
    write snp.allele
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Create a GDS genotype file:
The new dataset consists of 20 samples and 3000 SNPs
    write sample.id
    write snp.id
    write snp.rs.id
    write snp.position
    write snp.chromosome
    write snp.allele
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Merge SNP GDS files:
    open 't1.gds' ...
        10 samples, 3000 SNPs
    open 't2.gds' ...
        20 samples, 3000 SNPs
Concatenating samples (mapping to the first GDS file) ...
    reference: 3000 SNPs (100.0%)
    file 2: 0 allele flips, 0 ambiguous locus/loci
        [no flip]: 3000
    create 'test.gds': 30 samples, 3000 SNPs
    FileFormat = SNP_ARRAY
    writing genotypes ...
Clean up the fragments of GDS file:
    open the file 'test.gds' (46.2K)
    # of fragments: 32
    save to 'test.gds.tmp'
    rename 'test.gds.tmp' (46.0K, reduced: 204B)
    # of fragments: 15
Done.
Create a GDS genotype file:
The new dataset consists of 279 samples and 100 SNPs
    write sample.id
    write snp.id
    write snp.rs.id
    write snp.position
    write snp.chromosome
    write snp.allele
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Create a GDS genotype file:
The new dataset consists of 279 samples and 200 SNPs
    write sample.id
    write snp.id
    write snp.rs.id
    write snp.position
    write snp.chromosome
    write snp.allele
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Merge SNP GDS files:
    open 't1.gds' ...
        279 samples, 100 SNPs
    open 't2.gds' ...
        279 samples, 200 SNPs
Concatenating SNPs ...
    create 'test.gds': 279 samples, 300 SNPs
    FileFormat = SNP_ARRAY
    writing genotypes ...
Clean up the fragments of GDS file:
    open the file 'test.gds' (19.1K)
    # of fragments: 32
    save to 'test.gds.tmp'
    rename 'test.gds.tmp' (18.9K, reduced: 204B)
    # of fragments: 15
Done.
FUNCTION: snpgdsCreateGeno
Principal Component Analysis (PCA) on genotypes:
Excluding 42 SNPs on non-autosomes
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 958 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 264760
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:42:54 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:42:54 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:42:54 2019    Done.
FUNCTION: snpgdsCreateGenoSet
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 75.42%, 540/716
Chromosome 2: 72.24%, 536/742
Chromosome 3: 74.71%, 455/609
Chromosome 4: 73.31%, 412/562
Chromosome 5: 77.03%, 436/566
Chromosome 6: 75.58%, 427/565
Chromosome 7: 75.42%, 356/472
Chromosome 8: 71.31%, 348/488
Chromosome 9: 77.88%, 324/416
Chromosome 10: 74.33%, 359/483
Chromosome 11: 77.40%, 346/447
Chromosome 12: 76.81%, 328/427
Chromosome 13: 75.58%, 260/344
Chromosome 14: 76.95%, 217/282
Chromosome 15: 76.34%, 200/262
Chromosome 16: 72.66%, 202/278
Chromosome 17: 74.40%, 154/207
Chromosome 18: 73.68%, 196/266
Chromosome 19: 85.00%, 102/120
Chromosome 20: 71.62%, 164/229
Chromosome 21: 76.98%, 97/126
Chromosome 22: 75.86%, 88/116
6,547 markers are selected in total.
Create a GDS genotype file:
The new dataset consists of 279 samples and 6547 SNPs
    write sample.id
    write snp.id
    write snp.rs.id
    write snp.position
    write snp.chromosome
    write snp.allele
SNP genotypes are stored in SNP-major mode (Sample X SNP).
FUNCTION: snpgdsCutTree
Individual dissimilarity analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Dissimilarity:    the sum of all selected genotypes (0,1,2) = 2446510
Dissimilarity:	Sat Apr 13 05:42:55 2019	0%
Dissimilarity:	Sat Apr 13 05:42:56 2019	100%
Determine groups by permutation (Z threshold: 15, outlier threshold: 5):
Create 3 groups.
Create 4 groups.
FUNCTION: snpgdsDiss
Individual dissimilarity analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Dissimilarity:    the sum of all selected genotypes (0,1,2) = 2446510
Dissimilarity:	Sat Apr 13 05:42:57 2019	0%
Dissimilarity:	Sat Apr 13 05:42:58 2019	100%
Determine groups by permutation (Z threshold: 15, outlier threshold: 5):
Create 3 groups.
FUNCTION: snpgdsDrawTree
Individual dissimilarity analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Dissimilarity:    the sum of all selected genotypes (0,1,2) = 2446510
Dissimilarity:	Sat Apr 13 05:42:59 2019	0%
Dissimilarity:	Sat Apr 13 05:43:00 2019	100%
Determine groups by permutation (Z threshold: 15, outlier threshold: 5):
Create 3 groups.
FUNCTION: snpgdsEIGMIX
Eigen-analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Eigen-analysis:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:01 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:01 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:43:01 2019    Done.
FUNCTION: snpgdsErrMsg
FUNCTION: snpgdsExampleFileName
FUNCTION: snpgdsFst
Fst estimation on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
Method: Weir & Cockerham, 1984
# of Populations: 4
    CEU (92), HCB (47), JPT (47), YRI (93)
Fst estimation on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
Method: Weir & Hill, 2002
# of Populations: 4
    CEU (92), HCB (47), JPT (47), YRI (93)
FUNCTION: snpgdsGDS2BED
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: 0.95)
Converting from GDS to PLINK binary PED:
Working space: 279 samples, 8722 SNPs
Output a BIM file.
Output a BED file ...
		Sat Apr 13 05:43:02 2019	0%
		Sat Apr 13 05:43:02 2019	100%
Done.
FUNCTION: snpgdsGDS2Eigen
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: 0.95)
Converting from GDS to EIGENSOFT:
	save to *.snp: 8722 snps
	save to *.ind: 279 samples
	Output: 	Sat Apr 13 05:43:02 2019	0%
	Output: 	Sat Apr 13 05:43:02 2019	100%
Done.
FUNCTION: snpgdsGDS2PED
Converting from GDS to PLINK PED:
	Output a MAP file DONE.
	Output a PED file ...
		Output: 	Sat Apr 13 05:43:02 2019	0%
		Output: 	Sat Apr 13 05:43:02 2019	100%
FUNCTION: snpgdsGEN2GDS
running snpgdsGEN2GDS ...
FUNCTION: snpgdsGRM
Genetic Relationship Matrix (GRM, GCTA):
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:02 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:43:03 2019    Done.
Genetic Relationship Matrix (GRM, GCTA):
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:03 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:04 2019    Done.
FUNCTION: snpgdsGetGeno
Genotype matrix: 1000 SNPs X 279 samples
Genotype matrix: 279 samples X 1000 SNPs
FUNCTION: snpgdsHCluster
Individual dissimilarity analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Dissimilarity:    the sum of all selected genotypes (0,1,2) = 2446510
Dissimilarity:	Sat Apr 13 05:43:04 2019	0%
Dissimilarity:	Sat Apr 13 05:43:05 2019	100%
Determine groups by permutation (Z threshold: 15, outlier threshold: 5):
Create 3 groups.
FUNCTION: snpgdsHWE
Keeping 716 SNPs according to chromosome 1
Excluding 160 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
FUNCTION: snpgdsIBDKING
IBD analysis (KING method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
No family is specified, and all individuals are treated as singletons.
Relationship inference in the presence of population stratification.
KING IBD:    the sum of all selected genotypes (0,1,2) = 702139
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:06 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:06 2019    Done.
IBD analysis (KING method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
No family is specified, and all individuals are treated as singletons.
Relationship inference in the presence of population stratification.
KING IBD:    the sum of all selected genotypes (0,1,2) = 702139
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:06 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:06 2019    Done.
IBD analysis (KING method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
# of families: 20, and within- and between-family relationship are estimated differently.
Relationship inference in the presence of population stratification.
KING IBD:    the sum of all selected genotypes (0,1,2) = 702139
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:07 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:07 2019    Done.
IBD analysis (KING method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
Relationship inference in a homogeneous population.
KING IBD:    the sum of all selected genotypes (0,1,2) = 702139
Sat Apr 13 05:43:07 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:07 2019    Done.
IBD analysis (KING method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
Relationship inference in a homogeneous population.
KING IBD:    the sum of all selected genotypes (0,1,2) = 702139
Sat Apr 13 05:43:07 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:43:08 2019    Done.
FUNCTION: snpgdsIBDMLE
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1,581 SNPs (monomorphic: TRUE, MAF: 0.05, missing rate: 0.05)
Working space: 30 samples, 7,142 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 54.75%, 392/716
Chromosome 2: 54.45%, 404/742
Chromosome 3: 55.67%, 339/609
Chromosome 4: 56.94%, 320/562
Chromosome 5: 56.71%, 321/566
Chromosome 6: 52.74%, 298/565
Chromosome 7: 55.72%, 263/472
Chromosome 8: 50.41%, 246/488
Chromosome 9: 54.33%, 226/416
Chromosome 10: 50.10%, 242/483
Chromosome 11: 54.81%, 245/447
Chromosome 12: 55.04%, 235/427
Chromosome 13: 53.49%, 184/344
Chromosome 14: 55.67%, 157/282
Chromosome 15: 54.58%, 143/262
Chromosome 16: 54.68%, 152/278
Chromosome 17: 55.56%, 115/207
Chromosome 18: 55.64%, 148/266
Chromosome 19: 66.67%, 80/120
Chromosome 20: 52.84%, 121/229
Chromosome 21: 50.79%, 64/126
Chromosome 22: 50.86%, 59/116
4,754 markers are selected in total.
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 30 samples, 250 SNPs
    using 1 (CPU) core
MLE IBD:    the sum of all selected genotypes (0,1,2) = 7195
MLE IBD:	Sat Apr 13 05:43:08 2019	0%
MLE IBD:	Sat Apr 13 05:43:09 2019	100%
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.480, sd: 0.281
MLE IBD:    the sum of all selected genotypes (0,1,2) = 5991
MLE IBD:	Sat Apr 13 05:43:09 2019	0%
MLE IBD:	Sat Apr 13 05:43:09 2019	100%
FUNCTION: snpgdsIBDMLELogLik
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1,581 SNPs (monomorphic: TRUE, MAF: 0.05, missing rate: 0.05)
Working space: 30 samples, 7,142 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 54.75%, 392/716
Chromosome 2: 54.45%, 404/742
Chromosome 3: 55.67%, 339/609
Chromosome 4: 56.94%, 320/562
Chromosome 5: 56.71%, 321/566
Chromosome 6: 52.74%, 298/565
Chromosome 7: 55.72%, 263/472
Chromosome 8: 50.41%, 246/488
Chromosome 9: 54.33%, 226/416
Chromosome 10: 50.10%, 242/483
Chromosome 11: 54.81%, 245/447
Chromosome 12: 55.04%, 235/427
Chromosome 13: 53.49%, 184/344
Chromosome 14: 55.67%, 157/282
Chromosome 15: 54.58%, 143/262
Chromosome 16: 54.68%, 152/278
Chromosome 17: 55.56%, 115/207
Chromosome 18: 55.64%, 148/266
Chromosome 19: 66.67%, 80/120
Chromosome 20: 52.84%, 121/229
Chromosome 21: 50.79%, 64/126
Chromosome 22: 50.86%, 59/116
4,754 markers are selected in total.
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 30 samples, 250 SNPs
    using 1 (CPU) core
MLE IBD:    the sum of all selected genotypes (0,1,2) = 7195
MLE IBD:	Sat Apr 13 05:43:09 2019	0%
MLE IBD:	Sat Apr 13 05:43:10 2019	100%
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.480, sd: 0.281
MLE IBD:    the sum of all selected genotypes (0,1,2) = 5991
MLE IBD:	Sat Apr 13 05:43:10 2019	0%
MLE IBD:	Sat Apr 13 05:43:10 2019	100%
FUNCTION: snpgdsIBDMoM
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1,217 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 92 samples, 7,506 SNPs
    using 1 (CPU) core
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 702139
Sat Apr 13 05:43:11 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:11 2019    Done.
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 563 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 93 samples, 8,160 SNPs
    using 1 (CPU) core
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 755648
Sat Apr 13 05:43:11 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:11 2019    Done.
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 563 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 93 samples, 8,160 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.500, sd: 0.315
*** A correction factor based on allele count is not used, since the allele frequencies are specified.
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 755648
Sat Apr 13 05:43:11 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:11 2019    Done.
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 563 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 8,160 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.500, sd: 0.315
*** A correction factor based on allele count is not used, since the allele frequencies are specified.
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 203285
Sat Apr 13 05:43:11 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:11 2019    Done.
FUNCTION: snpgdsIBDSelection
IBD analysis (PLINK method of moment) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 563 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 93 samples, 8,160 SNPs
    using 1 (CPU) core
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 755648
Sat Apr 13 05:43:11 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:11 2019    Done.
FUNCTION: snpgdsIBS
Identity-By-State (IBS) analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
IBS:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:43:11 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:11 2019    Done.
FUNCTION: snpgdsIBSNum
Identity-By-State (IBS) analysis on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
IBS:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:43:11 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:43:12 2019    Done.
FUNCTION: snpgdsIndInb
Estimating individual inbreeding coefficients:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
FUNCTION: snpgdsIndInbCoef
FUNCTION: snpgdsIndivBeta
Individual Inbreeding and Relatedness (beta estimator):
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
Individual Beta:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:12 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:12 2019    Done.
FUNCTION: snpgdsLDMat
Linkage Disequilibrium (LD) estimation on genotypes:
Working space: 279 samples, 203 SNPs
    using 1 (CPU) core.
    method: composite
LD matrix:    the sum of all selected genotypes (0,1,2) = 56582
Linkage Disequilibrium (LD) estimation on genotypes:
Working space: 279 samples, 203 SNPs
    using 1 (CPU) core.
    sliding window size: 203 
    method: composite
LD matrix:    the sum of all selected genotypes (0,1,2) = 56582
FUNCTION: snpgdsLDpair
FUNCTION: snpgdsLDpruning
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 75.42%, 540/716
Chromosome 2: 72.24%, 536/742
Chromosome 3: 74.71%, 455/609
Chromosome 4: 73.31%, 412/562
Chromosome 5: 77.03%, 436/566
Chromosome 6: 75.58%, 427/565
Chromosome 7: 75.42%, 356/472
Chromosome 8: 71.31%, 348/488
Chromosome 9: 77.88%, 324/416
Chromosome 10: 74.33%, 359/483
Chromosome 11: 77.40%, 346/447
Chromosome 12: 76.81%, 328/427
Chromosome 13: 75.58%, 260/344
Chromosome 14: 76.95%, 217/282
Chromosome 15: 76.34%, 200/262
Chromosome 16: 72.66%, 202/278
Chromosome 17: 74.40%, 154/207
Chromosome 18: 73.68%, 196/266
Chromosome 19: 85.00%, 102/120
Chromosome 20: 71.62%, 164/229
Chromosome 21: 76.98%, 97/126
Chromosome 22: 75.86%, 88/116
6,547 markers are selected in total.
FUNCTION: snpgdsMergeGRM
Genetic Relationship Matrix (GRM, GCTA):
Excluding 2,288 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 6,800 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 1908966
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:13 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:13 2019    Done.
Genetic Relationship Matrix (GRM, GCTA):
Excluding 5,688 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 3,400 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 951558
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:13 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:43:14 2019    Done.
Genetic Relationship Matrix (GRM, GCTA):
Excluding 5,688 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 3,400 SNPs
    using 1 (CPU) core
GRM Calculation:    the sum of all selected genotypes (0,1,2) = 957408
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:14 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Saving to the GDS file:

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:14 2019    Done.
GRM merging:
    open 'tmp1.gds' (3,400 variants)
    open 'tmp2.gds' (3,400 variants)
Weight: 0.5, 0.5
Output: tmp.gds

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
GRM merging:
    open 'tmp1.gds' (3,400 variants)
    open 'tmp2.gds' (3,400 variants)
Weight: 0.5, 0.5

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
FUNCTION: snpgdsOpen
FUNCTION: snpgdsOption
FUNCTION: snpgdsPCA
Principal Component Analysis (PCA) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:15 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:15 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:43:15 2019    Done.
FUNCTION: snpgdsPCACorr
Principal Component Analysis (PCA) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:15 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:43:16 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:43:16 2019    Done.
SNP Correlation:
Working space: 279 samples, 9088 SNPs
    using 1 (CPU) core
    using the top 4 eigenvectors
Correlation:    the sum of all selected genotypes (0,1,2) = 2553065
Sat Apr 13 05:43:16 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:16 2019    Done.
SNP Correlation:
Working space: 279 samples, 9088 SNPs
    using 1 (CPU) core
    using the top 4 eigenvectors
Creating 'test.gds' ...
Correlation:    the sum of all selected genotypes (0,1,2) = 2553065
Sat Apr 13 05:43:16 2019

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:16 2019    Done.
FUNCTION: snpgdsPCASNPLoading
Principal Component Analysis (PCA) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:16 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:16 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:43:16 2019    Done.
SNP loading:
Working space: 279 samples, 8722 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
SNP Loading:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:43:16 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 1s
Sat Apr 13 05:43:17 2019    Done.
FUNCTION: snpgdsPCASampLoading
Principal Component Analysis (PCA) on genotypes:
Excluding 365 SNPs on non-autosomes
Excluding 1 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 8,722 SNPs
    using 1 (CPU) core
PCA:    the sum of all selected genotypes (0,1,2) = 2446510
CPU capabilities: Double-Precision SSE2
Sat Apr 13 05:43:17 2019    (internal increment: 11612)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:17 2019    Begin (eigenvalues and eigenvectors)
Sat Apr 13 05:43:17 2019    Done.
SNP loading:
Working space: 279 samples, 8722 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
SNP Loading:    the sum of all selected genotypes (0,1,2) = 2446510
Sat Apr 13 05:43:17 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:17 2019    Done.
Sample loading:
Working space: 100 samples, 8722 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
Sample Loading:    the sum of all selected genotypes (0,1,2) = 878146
Sat Apr 13 05:43:17 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:17 2019    Done.
FUNCTION: snpgdsPED2GDS
Converting from GDS to PLINK PED:
	Output a MAP file DONE.
	Output a PED file ...
		Output: 	Sat Apr 13 05:43:17 2019	0%
		Output: 	Sat Apr 13 05:43:17 2019	100%
PLINK PED/MAP to GDS Format:
Import 9088 variants from 'tmp.map'
Chromosome:
  1  10  11  12  13  14  15  16  17  18  19   2  20  21  22   3   4   5   6   7 
716 483 447 427 344 282 262 278 207 266 120 742 229 126 116 609 562 566 565 472 
  8   9   X 
488 416 365 
Reading 'tmp.ped'
Output: 'test.gds'
Import 279 samples
Transpose the genotypic matrix ...
Done.
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test.gds' (1.3M)
    # of fragments: 50
    save to 'test.gds.tmp'
    rename 'test.gds.tmp' (711.4K, reduced: 618.7K)
    # of fragments: 26
FUNCTION: snpgdsPairIBD
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1,646 SNPs (monomorphic: TRUE, MAF: 0.05, missing rate: 0.05)
Working space: 93 samples, 7,077 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 62.29%, 446/716
Chromosome 2: 62.40%, 463/742
Chromosome 3: 60.76%, 370/609
Chromosome 4: 64.59%, 363/562
Chromosome 5: 62.54%, 354/566
Chromosome 6: 60.00%, 339/565
Chromosome 7: 62.71%, 296/472
Chromosome 8: 58.61%, 286/488
Chromosome 9: 62.98%, 262/416
Chromosome 10: 60.66%, 293/483
Chromosome 11: 62.86%, 281/447
Chromosome 12: 63.00%, 269/427
Chromosome 13: 63.08%, 217/344
Chromosome 14: 63.48%, 179/282
Chromosome 15: 63.36%, 166/262
Chromosome 16: 61.87%, 172/278
Chromosome 17: 65.70%, 136/207
Chromosome 18: 59.02%, 157/266
Chromosome 19: 68.33%, 82/120
Chromosome 20: 67.25%, 154/229
Chromosome 21: 61.11%, 77/126
Chromosome 22: 57.76%, 67/116
5,429 markers are selected in total.
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.506, sd: 0.300
MLE IBD:    the sum of all selected genotypes (0,1,2) = 6339
MLE IBD:	Sat Apr 13 05:43:18 2019	0%
MLE IBD:	Sat Apr 13 05:43:19 2019	100%
IBD analysis (PLINK method of moment) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.506, sd: 0.300
*** A correction factor based on allele count is not used, since the allele frequencies are specified.
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 6339
Sat Apr 13 05:43:19 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:19 2019    Done.
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.506, sd: 0.300
MLE IBD:    the sum of all selected genotypes (0,1,2) = 6339
MLE IBD:	Sat Apr 13 05:43:19 2019	0%
MLE IBD:	Sat Apr 13 05:43:20 2019	100%
Genotype matrix: 250 SNPs X 25 samples
[1] -364.4667
[1] -377.9771
[1] -390.3299
[1] -380.5342
[1] -378.2092
[1] -383.641
[1] -364.5406
[1] -375.196
[1] -383.5752
[1] -385.7771
[1] -368.6237
[1] -390.7766
[1] -367.5082
[1] -383.595
[1] -385.7602
[1] -378.5134
[1] -388.3292
[1] -376.4419
[1] -378.3509
[1] -386.8318
[1] -385.3803
[1] -382.988
[1] -369.9328
[1] -383.3912
FUNCTION: snpgdsPairIBDMLELogLik
SNP pruning based on LD:
Excluding 365 SNPs on non-autosomes
Excluding 1,646 SNPs (monomorphic: TRUE, MAF: 0.05, missing rate: 0.05)
Working space: 93 samples, 7,077 SNPs
    using 1 (CPU) core
    sliding window: 500,000 basepairs, Inf SNPs
    |LD| threshold: 0.2
    method: composite
Chromosome 1: 62.29%, 446/716
Chromosome 2: 62.40%, 463/742
Chromosome 3: 60.76%, 370/609
Chromosome 4: 64.59%, 363/562
Chromosome 5: 62.54%, 354/566
Chromosome 6: 60.00%, 339/565
Chromosome 7: 62.71%, 296/472
Chromosome 8: 58.61%, 286/488
Chromosome 9: 62.98%, 262/416
Chromosome 10: 60.66%, 293/483
Chromosome 11: 62.86%, 281/447
Chromosome 12: 63.00%, 269/427
Chromosome 13: 63.08%, 217/344
Chromosome 14: 63.48%, 179/282
Chromosome 15: 63.36%, 166/262
Chromosome 16: 61.87%, 172/278
Chromosome 17: 65.70%, 136/207
Chromosome 18: 59.02%, 157/266
Chromosome 19: 68.33%, 82/120
Chromosome 20: 67.25%, 154/229
Chromosome 21: 61.11%, 77/126
Chromosome 22: 57.76%, 67/116
5,429 markers are selected in total.
Identity-By-Descent analysis (MLE) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.506, sd: 0.300
MLE IBD:    the sum of all selected genotypes (0,1,2) = 6339
MLE IBD:	Sat Apr 13 05:43:20 2019	0%
MLE IBD:	Sat Apr 13 05:43:21 2019	100%
IBD analysis (PLINK method of moment) on genotypes:
Excluding 8,838 SNPs (non-autosomes or non-selection)
Excluding 0 SNP (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 25 samples, 250 SNPs
    using 1 (CPU) core
Specifying allele frequencies, mean: 0.506, sd: 0.300
*** A correction factor based on allele count is not used, since the allele frequencies are specified.
PLINK IBD:    the sum of all selected genotypes (0,1,2) = 6339
Sat Apr 13 05:43:21 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:21 2019    Done.
Genotype matrix: 250 SNPs X 25 samples
[1] -364.4667
[1] -377.9771
[1] -390.3299
[1] -380.5342
[1] -378.2092
[1] -383.641
[1] -364.5406
[1] -375.196
[1] -383.5752
[1] -385.7771
[1] -368.6237
[1] -390.7766
[1] -367.5082
[1] -383.595
[1] -385.7602
[1] -378.5134
[1] -388.3292
[1] -376.4419
[1] -378.3509
[1] -386.8318
[1] -385.3803
[1] -382.988
[1] -369.9328
[1] -383.3912
FUNCTION: snpgdsPairScore
Excluding 365 SNPs on non-autosomes
Working space: 120 samples, 8723 SNPs
Method: IBS
Genotype Score:    the sum of all selected genotypes (0,1,2) = 1050236
Working space: 120 samples, 8723 SNPs
Method: IBS
Genotype Score:    the sum of all selected genotypes (0,1,2) = 1050236
Working space: 120 samples, 8723 SNPs
Method: IBS
Genotype Score:    the sum of all selected genotypes (0,1,2) = 1050236
Working space: 120 samples, 8723 SNPs
Method: IBS
Genotype Score:    the sum of all selected genotypes (0,1,2) = 1050236
Working space: 120 samples, 8723 SNPs
Method: IBS
Output: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\tmp.gds
Genotype Score:    the sum of all selected genotypes (0,1,2) = 1050236
FUNCTION: snpgdsSNPList
FUNCTION: snpgdsSNPListClass
FUNCTION: snpgdsSNPListIntersect
FUNCTION: snpgdsSNPRateFreq
FUNCTION: snpgdsSampMissRate
FUNCTION: snpgdsSelectSNP
Excluding 365 SNPs on non-autosomes
Excluding 1,221 SNPs (monomorphic: TRUE, MAF: 0.05, missing rate: 0.95)
FUNCTION: snpgdsSlidingWindow
Sliding Window Analysis:
Excluding 8 SNPs (monomorphic: TRUE, MAF: NaN, missing rate: NaN)
Working space: 279 samples, 9,080 SNPs
    using 1 (CPU) core
    window size: 500000, shift: 100000 (basepair)
Chromosome Set: 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23
Sat Apr 13 05:43:22 2019, Chromosome 1 (716 SNPs), 2448 windows
Sat Apr 13 05:43:22 2019, Chromosome 2 (742 SNPs), 2416 windows
Sat Apr 13 05:43:22 2019, Chromosome 3 (609 SNPs), 1985 windows
Sat Apr 13 05:43:22 2019, Chromosome 4 (562 SNPs), 1894 windows
Sat Apr 13 05:43:22 2019, Chromosome 5 (566 SNPs), 1797 windows
Sat Apr 13 05:43:22 2019, Chromosome 6 (565 SNPs), 1694 windows
Sat Apr 13 05:43:22 2019, Chromosome 7 (472 SNPs), 1573 windows
Sat Apr 13 05:43:22 2019, Chromosome 8 (488 SNPs), 1445 windows
Sat Apr 13 05:43:22 2019, Chromosome 9 (416 SNPs), 1393 windows
Sat Apr 13 05:43:22 2019, Chromosome 10 (483 SNPs), 1343 windows
Sat Apr 13 05:43:22 2019, Chromosome 11 (447 SNPs), 1338 windows
Sat Apr 13 05:43:22 2019, Chromosome 12 (427 SNPs), 1316 windows
Sat Apr 13 05:43:22 2019, Chromosome 13 (344 SNPs), 948 windows
Sat Apr 13 05:43:22 2019, Chromosome 14 (281 SNPs), 847 windows
Sat Apr 13 05:43:22 2019, Chromosome 15 (262 SNPs), 774 windows
Sat Apr 13 05:43:22 2019, Chromosome 16 (278 SNPs), 873 windows
Sat Apr 13 05:43:22 2019, Chromosome 17 (207 SNPs), 773 windows
Sat Apr 13 05:43:22 2019, Chromosome 18 (266 SNPs), 753 windows
Sat Apr 13 05:43:22 2019, Chromosome 19 (120 SNPs), 627 windows
Sat Apr 13 05:43:22 2019, Chromosome 20 (229 SNPs), 602 windows
Sat Apr 13 05:43:22 2019, Chromosome 21 (126 SNPs), 311 windows
Sat Apr 13 05:43:22 2019, Chromosome 22 (116 SNPs), 312 windows
Sat Apr 13 05:43:22 2019, Chromosome 23 (358 SNPs), 1507 windows
Sat Apr 13 05:43:22 2019 	Done.
FUNCTION: snpgdsSummary
The file name: C:\Users\biocbuild\bbs-3.8-bioc\R\library\SNPRelate\extdata\hapmap_geno.gds 
The total number of samples: 279 
The total number of SNPs: 9088 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
FUNCTION: snpgdsTranspose
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test.gds 
The total number of samples: 279 
The total number of SNPs: 9088 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
SNP genotypes: 279 samples, 9088 SNPs
Genotype matrix is being transposed ...
Clean up the fragments of GDS file:
    open the file 'test.gds' (1.3M)
    # of fragments: 28
    save to 'test.gds.tmp'
    rename 'test.gds.tmp' (709.6K, reduced: 619.1K)
    # of fragments: 26
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test.gds 
The total number of samples: 279 
The total number of SNPs: 9088 
SNP genotypes are stored in individual-major mode (SNP X Sample).
FUNCTION: snpgdsVCF2GDS
##fileformat=VCFv4.1
##fileDate=20090805
##source=myImputationProgramV3.1
##reference=file:///seq/references/1000GenomesPilot-NCBI36.fasta
##contig=<ID=20,length=62435964,assembly=B36,md5=f126cdf8a6e0c7f379d618ff66beb2da,species="Homo sapiens",taxonomy=x>
##phasing=partial
##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data">
##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth">
##INFO=<ID=AF,Number=A,Type=Float,Description="Allele Frequency">
##INFO=<ID=AA,Number=1,Type=String,Description="Ancestral Allele">
##INFO=<ID=DB,Number=0,Type=Flag,Description="dbSNP membership, build 129">
##INFO=<ID=H2,Number=0,Type=Flag,Description="HapMap2 membership">
##FILTER=<ID=q10,Description="Quality below 10">
##FILTER=<ID=s50,Description="Less than 50% of samples have data">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality">
##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">
##FORMAT=<ID=HQ,Number=2,Type=Integer,Description="Haplotype Quality">
#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	NA00001	NA00002	NA00003
20	14370	rs6054257	G	A	29	PASS	NS=3;DP=14;AF=0.5;DB;H2	GT:GQ:DP:HQ	0|0:48:1:51,51	1|0:48:8:51,51	1/1:43:5:.,.
20	17330	.	T	A	3	q10	NS=3;DP=11;AF=0.017	GT:GQ:DP:HQ	0|0:49:3:58,50	0|1:3:5:65,3	0/0:41:3
20	1110696	rs6040355	A	G,T	67	PASS	NS=2;DP=10;AF=0.333,0.667;AA=T;DB	GT:GQ:DP:HQ	1|2:21:6:23,27	2|1:2:0:18,2	2/2:35:4
20	1230237	.	T	.	47	PASS	NS=3;DP=13;AA=T	GT:GQ:DP:HQ	0|0:54:7:56,60	0|0:48:4:51,51	0/0:61:2
20	1234567	microsat1	GTC	G,GTCT	50	PASS	NS=3;DP=9;AA=G	GT:GQ:DP	0/1:35:4	0/2:17:2	1/1:40:3
VCF Format ==> SNP GDS Format
Method: exacting biallelic SNPs
Number of samples: 3
Parsing "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf" ...
	import 2 variants.
+ genotype   { Bit2 3x2, 2B } *
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test1.gds' (2.9K)
    # of fragments: 46
    save to 'test1.gds.tmp'
    rename 'test1.gds.tmp' (2.6K, reduced: 312B)
    # of fragments: 20
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test1.gds 
The total number of samples: 3 
The total number of SNPs: 2 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
VCF Format ==> SNP GDS Format
Method: exacting biallelic SNPs
Number of samples: 3
Parsing "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf" ...
	import 2 variants.
+ genotype   { Bit2 3x2, 2B } *
SNP genotypes: 3 samples, 2 SNPs
Genotype matrix is being transposed ...
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test2.gds' (3.0K)
    # of fragments: 48
    save to 'test2.gds.tmp'
    rename 'test2.gds.tmp' (2.6K, reduced: 417B)
    # of fragments: 20
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test2.gds 
The total number of samples: 3 
The total number of SNPs: 2 
SNP genotypes are stored in individual-major mode (SNP X Sample).
VCF Format ==> SNP GDS Format
Method: dosage (0,1,2) of reference allele for all variant sites
Number of samples: 3
Parsing "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf" ...
	import 5 variants.
+ genotype   { Bit2 3x5, 4B } *
SNP genotypes: 3 samples, 5 SNPs
Genotype matrix is being transposed ...
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test3.gds' (3.1K)
    # of fragments: 48
    save to 'test3.gds.tmp'
    rename 'test3.gds.tmp' (2.7K, reduced: 419B)
    # of fragments: 20
Some of 'snp.allele' are not standard (e.g., A/G,T).
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test3.gds 
The total number of samples: 3 
The total number of SNPs: 5 
SNP genotypes are stored in individual-major mode (SNP X Sample).
The number of valid samples: 3 
The number of biallelic unique SNPs: 2 
VCF Format ==> SNP GDS Format
Method: dosage (0,1,2) of reference allele for all variant sites
Number of samples: 3
Parsing "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf" ...
	import 5 variants.
+ genotype   { Bit2 3x5, 4B } *
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test4.gds' (3.0K)
    # of fragments: 46
    save to 'test4.gds.tmp'
    rename 'test4.gds.tmp' (2.7K, reduced: 312B)
    # of fragments: 20
Some of 'snp.allele' are not standard (e.g., A/G,T).
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test4.gds 
The total number of samples: 3 
The total number of SNPs: 5 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
The number of valid samples: 3 
The number of biallelic unique SNPs: 2 
VCF Format ==> SNP GDS Format
Method: dosage (0,1,2) of reference allele for all variant sites
Number of samples: 3
Parsing "C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf" ...
	import 5 variants.
+ genotype   { Bit2 3x5, 4B } *
Optimize the access efficiency ...
Clean up the fragments of GDS file:
    open the file 'test5.gds' (3.0K)
    # of fragments: 46
    save to 'test5.gds.tmp'
    rename 'test5.gds.tmp' (2.7K, reduced: 312B)
    # of fragments: 20
Some of 'snp.allele' are not standard (e.g., T/A,G).
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test5.gds 
The total number of samples: 3 
The total number of SNPs: 5 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
The number of valid samples: 3 
The number of biallelic unique SNPs: 2 
FUNCTION: snpgdsVCF2GDS_R
##fileformat=VCFv4.1
##fileDate=20090805
##source=myImputationProgramV3.1
##reference=file:///seq/references/1000GenomesPilot-NCBI36.fasta
##contig=<ID=20,length=62435964,assembly=B36,md5=f126cdf8a6e0c7f379d618ff66beb2da,species="Homo sapiens",taxonomy=x>
##phasing=partial
##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data">
##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth">
##INFO=<ID=AF,Number=A,Type=Float,Description="Allele Frequency">
##INFO=<ID=AA,Number=1,Type=String,Description="Ancestral Allele">
##INFO=<ID=DB,Number=0,Type=Flag,Description="dbSNP membership, build 129">
##INFO=<ID=H2,Number=0,Type=Flag,Description="HapMap2 membership">
##FILTER=<ID=q10,Description="Quality below 10">
##FILTER=<ID=s50,Description="Less than 50% of samples have data">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality">
##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">
##FORMAT=<ID=HQ,Number=2,Type=Integer,Description="Haplotype Quality">
#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	NA00001	NA00002	NA00003
20	14370	rs6054257	G	A	29	PASS	NS=3;DP=14;AF=0.5;DB;H2	GT:GQ:DP:HQ	0|0:48:1:51,51	1|0:48:8:51,51	1/1:43:5:.,.
20	17330	.	T	A	3	q10	NS=3;DP=11;AF=0.017	GT:GQ:DP:HQ	0|0:49:3:58,50	0|1:3:5:65,3	0/0:41:3
20	1110696	rs6040355	A	G,T	67	PASS	NS=2;DP=10;AF=0.333,0.667;AA=T;DB	GT:GQ:DP:HQ	1|2:21:6:23,27	2|1:2:0:18,2	2/2:35:4
20	1230237	.	T	.	47	PASS	NS=3;DP=13;AA=T	GT:GQ:DP:HQ	0|0:54:7:56,60	0|0:48:4:51,51	0/0:61:2
20	1234567	microsat1	GTC	G,GTCT	50	PASS	NS=3;DP=9;AA=G	GT:GQ:DP	0/1:35:4	0/2:17:2	1/1:40:3
Start snpgdsVCF2GDS ...
	Extracting bi-allelic and polymorhpic SNPs.
	Scanning ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
	content: 5 rows x 12 columns
Sat Apr 13 05:43:23 2019 	store sample id, snp id, position, and chromosome.
	start writing: 3 samples, 2 SNPs ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
[1] 1
Sat Apr 13 05:43:23 2019 	Done.
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test1.gds 
The total number of samples: 3 
The total number of SNPs: 2 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Start snpgdsVCF2GDS ...
	Extracting bi-allelic and polymorhpic SNPs.
	Scanning ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
	content: 5 rows x 12 columns
Sat Apr 13 05:43:23 2019 	store sample id, snp id, position, and chromosome.
	start writing: 3 samples, 2 SNPs ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
[1] 1
Sat Apr 13 05:43:23 2019 	Done.
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test2.gds 
The total number of samples: 3 
The total number of SNPs: 2 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
Start snpgdsVCF2GDS ...
	Storing dosage of the reference allele for all variant sites, including bi-allelic SNPs, multi-allelic SNPs, indels and structural variants.
	Scanning ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
	content: 5 rows x 12 columns
Sat Apr 13 05:43:23 2019 	store sample id, snp id, position, and chromosome.
	start writing: 3 samples, 5 SNPs ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
Sat Apr 13 05:43:23 2019 	Done.
Some of 'snp.allele' are not standard (e.g., A/G,T).
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test3.gds 
The total number of samples: 3 
The total number of SNPs: 5 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
The number of valid samples: 3 
The number of biallelic unique SNPs: 2 
Start snpgdsVCF2GDS ...
	Storing dosage of the reference allele for all variant sites, including bi-allelic SNPs, multi-allelic SNPs, indels and structural variants.
	Scanning ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
	content: 5 rows x 12 columns
Sat Apr 13 05:43:23 2019 	store sample id, snp id, position, and chromosome.
	start writing: 3 samples, 5 SNPs ...
	file: C:/Users/biocbuild/bbs-3.8-bioc/R/library/SNPRelate/extdata/sequence.vcf
Sat Apr 13 05:43:23 2019 	Done.
Some of 'snp.allele' are not standard (e.g., A/G,T).
The file name: C:\Users\biocbuild\bbs-3.8-bioc\meat\SNPRelate.Rcheck\tests_x64\test4.gds 
The total number of samples: 3 
The total number of SNPs: 5 
SNP genotypes are stored in SNP-major mode (Sample X SNP).
The number of valid samples: 3 
The number of biallelic unique SNPs: 2 
SNP Correlation:
Working space: 90 samples, 9088 SNPs
    using 1 (CPU) core
    using the top 2 eigenvectors
Correlation:    the sum of all selected genotypes (0,1,2) = 824424
Sat Apr 13 05:43:24 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:24 2019    Done.
SNP Correlation:
Working space: 90 samples, 9088 SNPs
    using 1 (CPU) core
    using the top 2 eigenvectors
Creating 'test.gds' ...
Correlation:    the sum of all selected genotypes (0,1,2) = 824424
Sat Apr 13 05:43:24 2019

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:24 2019    Done.
SNP loading:
Working space: 90 samples, 8695 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
SNP Loading:    the sum of all selected genotypes (0,1,2) = 787449
Sat Apr 13 05:43:24 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:24 2019    Done.
Sample loading:
Working space: 100 samples, 8695 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
Sample Loading:    the sum of all selected genotypes (0,1,2) = 875255
Sat Apr 13 05:43:24 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:24 2019    Done.
SNP Correlation:
Working space: 90 samples, 9088 SNPs
    using 2 (CPU) cores
    using the top 2 eigenvectors
Correlation:    the sum of all selected genotypes (0,1,2) = 824424
Sat Apr 13 05:43:24 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:24 2019    Done.
SNP Correlation:
Working space: 90 samples, 9088 SNPs
    using 2 (CPU) cores
    using the top 2 eigenvectors
Creating 'test.gds' ...
Correlation:    the sum of all selected genotypes (0,1,2) = 824424
Sat Apr 13 05:43:24 2019

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:24 2019    Done.
SNP loading:
Working space: 90 samples, 8695 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
SNP Loading:    the sum of all selected genotypes (0,1,2) = 787449
Sat Apr 13 05:43:24 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:24 2019    Done.
Sample loading:
Working space: 100 samples, 8695 SNPs
    using 1 (CPU) core
    using the top 8 eigenvectors
Sample Loading:    the sum of all selected genotypes (0,1,2) = 875255
Sat Apr 13 05:43:24 2019    (internal increment: 65536)

[..................................................]  0%, ETC: ---        
[==================================================] 100%, completed in 0s
Sat Apr 13 05:43:24 2019    Done.


RUNIT TEST PROTOCOL -- Sat Apr 13 05:43:24 2019 
*********************************************** 
Number of test functions: 13 
Number of errors: 0 
Number of failures: 0 

 
1 Test Suite : 
SNPRelate RUnit Tests - 13 test functions, 0 errors, 0 failures
Number of test functions: 13 
Number of errors: 0 
Number of failures: 0 
> 
> proc.time()
   user  system elapsed 
  44.21    2.17   46.17 

Example timings

SNPRelate.Rcheck/examples_i386/SNPRelate-Ex.timings

nameusersystemelapsed
SNPGDSFileClass-class0.030.000.03
SNPRelate-package3.970.174.14
snpgdsAdmixPlot0.960.020.98
snpgdsAdmixProp0.860.000.86
snpgdsAlleleSwitch0.110.020.13
snpgdsApartSelection0.080.010.09
snpgdsBED2GDS0.100.020.11
snpgdsClose0.010.000.02
snpgdsCombineGeno0.240.030.26
snpgdsCreateGeno0.870.000.88
snpgdsCreateGenoSet0.230.010.25
snpgdsCutTree2.660.082.73
snpgdsDiss202
snpgdsDrawTree2.020.002.02
snpgdsEIGMIX0.850.020.86
snpgdsErrMsg000
snpgdsExampleFileName000
snpgdsFst0.040.000.05
snpgdsGDS2BED0.100.020.10
snpgdsGDS2Eigen0.230.000.24
snpgdsGDS2PED0.520.060.58
snpgdsGEN2GDS000
snpgdsGRM2.110.062.17
snpgdsGetGeno0.070.020.09
snpgdsHCluster2.720.012.74
snpgdsHWE0.030.000.03
snpgdsIBDKING2.660.082.73
snpgdsIBDMLE1.280.051.33
snpgdsIBDMLELogLik1.350.001.34
snpgdsIBDMoM0.680.010.70
snpgdsIBDSelection0.190.020.21
snpgdsIBS1.450.011.47
snpgdsIBSNum1.80.01.8
snpgdsIndInb0.040.000.04
snpgdsIndInbCoef0.020.000.02
snpgdsIndivBeta1.280.021.30
snpgdsLDMat0.390.030.42
snpgdsLDpair000
snpgdsLDpruning0.060.000.06
snpgdsMergeGRM3.250.143.45
snpgdsOpen0.020.000.02
snpgdsOption000
snpgdsPCA0.860.030.89
snpgdsPCACorr0.920.050.97
snpgdsPCASNPLoading0.780.010.80
snpgdsPCASampLoading0.830.020.84
snpgdsPED2GDS1.640.121.77
snpgdsPairIBD1.840.001.85
snpgdsPairIBDMLELogLik1.270.001.26
snpgdsPairScore0.190.050.24
snpgdsSNPList0.010.000.01
snpgdsSNPListIntersect0.100.000.09
snpgdsSNPRateFreq0.030.000.04
snpgdsSampMissRate0.010.000.01
snpgdsSelectSNP0.020.000.02
snpgdsSlidingWindow0.670.040.72
snpgdsSummary0.050.000.04
snpgdsTranspose0.150.000.16
snpgdsVCF2GDS0.160.140.30
snpgdsVCF2GDS_R0.050.070.11

SNPRelate.Rcheck/examples_x64/SNPRelate-Ex.timings

nameusersystemelapsed
SNPGDSFileClass-class0.030.000.04
SNPRelate-package1.690.091.78
snpgdsAdmixPlot0.420.000.42
snpgdsAdmixProp0.400.030.44
snpgdsAlleleSwitch0.110.000.11
snpgdsApartSelection0.070.020.08
snpgdsBED2GDS0.060.050.11
snpgdsClose0.010.000.01
snpgdsCombineGeno0.190.030.36
snpgdsCreateGeno0.380.030.41
snpgdsCreateGenoSet0.180.000.19
snpgdsCutTree2.280.092.37
snpgdsDiss1.80.01.8
snpgdsDrawTree1.780.011.80
snpgdsEIGMIX0.390.020.41
snpgdsErrMsg000
snpgdsExampleFileName000
snpgdsFst0.050.000.04
snpgdsGDS2BED0.080.000.11
snpgdsGDS2Eigen0.190.000.22
snpgdsGDS2PED0.360.090.46
snpgdsGEN2GDS000
snpgdsGRM1.060.101.15
snpgdsGetGeno0.060.010.08
snpgdsHCluster2.020.022.03
snpgdsHWE0.010.000.02
snpgdsIBDKING1.600.041.64
snpgdsIBDMLE1.320.001.33
snpgdsIBDMLELogLik1.330.001.32
snpgdsIBDMoM0.240.000.24
snpgdsIBDSelection0.070.000.08
snpgdsIBS0.360.000.36
snpgdsIBSNum0.530.000.53
snpgdsIndInb0.020.000.02
snpgdsIndInbCoef0.010.000.01
snpgdsIndivBeta0.30.00.3
snpgdsLDMat0.340.020.36
snpgdsLDpair0.020.000.02
snpgdsLDpruning0.040.000.04
snpgdsMergeGRM1.960.142.10
snpgdsOpen0.010.000.01
snpgdsOption000
snpgdsPCA0.580.020.60
snpgdsPCACorr0.610.090.70
snpgdsPCASNPLoading0.590.030.62
snpgdsPCASampLoading0.570.000.57
snpgdsPED2GDS1.260.051.31
snpgdsPairIBD1.780.001.78
snpgdsPairIBDMLELogLik1.190.011.20
snpgdsPairScore0.160.080.24
snpgdsSNPList0.010.000.01
snpgdsSNPListIntersect0.10.00.1
snpgdsSNPRateFreq0.010.000.01
snpgdsSampMissRate0.020.000.02
snpgdsSelectSNP000
snpgdsSlidingWindow0.720.000.72
snpgdsSummary0.040.000.04
snpgdsTranspose0.130.000.13
snpgdsVCF2GDS0.080.110.19
snpgdsVCF2GDS_R0.060.050.17