******************************************************************* * Genome-wide Complex Trait Analysis (GCTA) * version 1.25.2 * (C) 2010-2013 Jian Yang, Hong Lee, Michael Goddard and Peter Visscher * The University of Queensland * MIT License ******************************************************************* Analysis started: Mon Aug 13 18:49:56 2018 Options: --mlma --bfile /home/mw489/Somo/AllChrom_TZ --keep /home/mw489/Somo/keeplist.txt --maf 0.05 --pheno /home/mw489/Somo/blups.txt --grm-bin /home/mw489/Somo/Kinship_TZ --out /home/mw489/Somo/HI_PYT_UKIRIGURU_TARI_81318 --thread-num 72 Note: the program will be running on 72 threads. Reading PLINK FAM file from [/home/mw489/Somo/AllChrom_TZ.fam]. 1344 individuals to be included from [/home/mw489/Somo/AllChrom_TZ.fam]. 84 individuals are kept from [/home/mw489/Somo/keeplist.txt]. Reading PLINK BIM file from [/home/mw489/Somo/AllChrom_TZ.bim]. 165816 SNPs to be included from [/home/mw489/Somo/AllChrom_TZ.bim]. Reading PLINK BED file from [/home/mw489/Somo/AllChrom_TZ.bed] in SNP-major format ... Genotype data for 84 individuals and 165816 SNPs to be included from [/home/mw489/Somo/AllChrom_TZ.bed]. Calculating allele frequencies ... Filtering SNPs with MAF > 0.05 ... After filtering SNPs with MAF > 0.05, there are 81917 SNPs (83899 SNPs with MAF < 0.05). Reading phenotypes from [/home/mw489/Somo/blups.txt]. Non-missing phenotypes of 84 individuals are included from [/home/mw489/Somo/blups.txt]. Reading IDs of the GRM from [/home/mw489/Somo/Kinship_TZ.grm.id]. 1344 IDs read from [/home/mw489/Somo/Kinship_TZ.grm.id]. Reading the GRM from [/home/mw489/Somo/Kinship_TZ.grm.bin]. Pairwise genetic relationships between 1344 individuals are included from [/home/mw489/Somo/Kinship_TZ.grm.bin]. 84 individuals are in common in these files. Performing MLM association analyses (including the candidate SNP) ... Performing REML analysis ... (Note: may take hours depending on sample size). 84 observations, 1 fixed effect(s), and 2 variance component(s)(including residual variance). Calculating prior values of variance components by EM-REML ... Updated prior values: 0.000648785 0.000608402 logL: 239.825 Running AI-REML algorithm ... Iter. logL V(G) V(e) 1 239.94 0.00076 0.00052 2 240.25 0.00105 0.00031 3 240.60 0.00111 0.00029 4 240.62 0.00112 0.00028 5 240.62 0.00112 0.00028 6 240.62 0.00112 0.00028 7 240.62 0.00112 0.00028 Log-likelihood ratio converged. Running association tests for 81917 SNPs ... Saving the results of the mixed linear model association analyses of 81917 SNPs to [/home/mw489/Somo/HI_PYT_UKIRIGURU_TARI_81318.mlma] ... Analysis finished: Mon Aug 13 18:49:59 2018 Computational time: 0:0:3