Post by K R Patil on Jan 13, 2014 18:06:43 GMT
Hi,
I am running gcta 1.24 on PLINK imputed data. I want to calculated GRMs (in gz format) for each of the chromosomes but it runs into an error. The details are at the end of the post. I also tried to run windows version which also runs into an error.
I will really appreciate any help.
Best regards,
Kaustubh
Details on gcta64 run on ubuntu 64 bit with 3G of RAM (on vmware host windows 7)
(run_gcta_chr.sh is a simple shell script iterating over chromosomes 1..22)
--
*******************************************************************
* Genome-wide Complex Trait Analysis (GCTA)
* version 1.24
* (C) 2010-2013 Jian Yang, Hong Lee, Michael Goddard and Peter Visscher
* The University of Queensland
* MIT License
*******************************************************************
Analysis started: Mon Jan 13 10:01:32 2014
Options:
--bfile PLINK/Imputed/chr1-FIN
--chr 1
--make-grm
--out chr1-FIN-grm
Note: This is a multi-thread program. You could specify the number of threads by the --thread-num option to speed up the computation if there are multiple processors in your machine.
Reading PLINK FAM file from [PLINK/Imputed/chr1-FIN.fam].
5654 individuals to be included from [PLINK/Imputed/chr1-FIN.fam].
Reading PLINK BIM file from [PLINK/Imputed/chr1-FIN.bim].
238819 SNPs to be included from [PLINK/Imputed/chr1-FIN.bim].
238819 SNPs on chromosome 1 are included in the analysis.
Reading PLINK BED file from [PLINK/Imputed/chr1-FIN.bed] in SNP-major format ...
Genotype data for 5654 individuals and 238819 SNPs to be included from [PLINK/Imputed/chr1-FIN.bed].
terminate called after throwing an instance of 'std::bad_alloc'
what(): std::bad_alloc
./run_gcta_chr.sh: line 2: 5245 Aborted
I am running gcta 1.24 on PLINK imputed data. I want to calculated GRMs (in gz format) for each of the chromosomes but it runs into an error. The details are at the end of the post. I also tried to run windows version which also runs into an error.
I will really appreciate any help.
Best regards,
Kaustubh
Details on gcta64 run on ubuntu 64 bit with 3G of RAM (on vmware host windows 7)
(run_gcta_chr.sh is a simple shell script iterating over chromosomes 1..22)
--
*******************************************************************
* Genome-wide Complex Trait Analysis (GCTA)
* version 1.24
* (C) 2010-2013 Jian Yang, Hong Lee, Michael Goddard and Peter Visscher
* The University of Queensland
* MIT License
*******************************************************************
Analysis started: Mon Jan 13 10:01:32 2014
Options:
--bfile PLINK/Imputed/chr1-FIN
--chr 1
--make-grm
--out chr1-FIN-grm
Note: This is a multi-thread program. You could specify the number of threads by the --thread-num option to speed up the computation if there are multiple processors in your machine.
Reading PLINK FAM file from [PLINK/Imputed/chr1-FIN.fam].
5654 individuals to be included from [PLINK/Imputed/chr1-FIN.fam].
Reading PLINK BIM file from [PLINK/Imputed/chr1-FIN.bim].
238819 SNPs to be included from [PLINK/Imputed/chr1-FIN.bim].
238819 SNPs on chromosome 1 are included in the analysis.
Reading PLINK BED file from [PLINK/Imputed/chr1-FIN.bed] in SNP-major format ...
Genotype data for 5654 individuals and 238819 SNPs to be included from [PLINK/Imputed/chr1-FIN.bed].
terminate called after throwing an instance of 'std::bad_alloc'
what(): std::bad_alloc
./run_gcta_chr.sh: line 2: 5245 Aborted