******************************************************************* * 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:47:55 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/MCBSDS_CET_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]. 327 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 327 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 87801 SNPs (78015 SNPs with MAF < 0.05). Reading phenotypes from [/home/mw489/Somo/blups.txt]. Non-missing phenotypes of 327 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]. 327 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). 327 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.028012 0.0321194 logL: 306.959 Running AI-REML algorithm ... Iter. logL V(G) V(e) 1 309.80 0.02160 0.03585 2 311.64 0.01787 0.03839 3 312.47 0.01560 0.04012 4 312.85 0.01111 0.04384 5 313.18 0.01136 0.04389 6 313.18 0.01137 0.04388 7 313.18 0.01137 0.04388 Log-likelihood ratio converged. Running association tests for 87801 SNPs ... Saving the results of the mixed linear model association analyses of 87801 SNPs to [/home/mw489/Somo/MCBSDS_CET_UKIRIGURU_TARI_81318.mlma] ... Analysis finished: Mon Aug 13 18:48:02 2018 Computational time: 0:0:7