Loading required namespace: GenomicFiles Using local VCF. File already tabix-indexed. Finding empty VCF columns based on first 10,000 rows. 1 sample detected: EBI-a-GCST005060 Constructing ScanVcfParam object. VCF contains: 23,315,353 variant(s) x 1 sample(s) Reading VCF file: multi-threaded (4 threads) Renaming ID as SNP. VCF file has -log10 P-values; these will be converted to unadjusted p-values in the 'P' column. No INFO (SI) column detected. Standardising column headers. First line of summary statistics file: SNP chr BP end REF ALT FILTER ReverseComplementedAlleles ES LP SE P Summary statistics report: - 22,718,565 rows - 22,717,995 unique variants - 132 genome-wide significant variants (P<5e-8) - 22 chromosomes Checking for multi-GWAS. Checking for multiple RSIDs on one row. Inferring genome build. Loading SNPlocs data. Loading reference genome data. Preprocessing RSIDs. Validating RSIDs of 10,000 SNPs using BSgenome::snpsById... Loading required package: BiocGenerics Attaching package: ‘BiocGenerics’ The following objects are masked from ‘package:stats’: IQR, mad, sd, var, xtabs The following objects are masked from ‘package:base’: anyDuplicated, aperm, append, as.data.frame, basename, cbind, colnames, dirname, do.call, duplicated, eval, evalq, Filter, Find, get, grep, grepl, intersect, is.unsorted, lapply, Map, mapply, match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank, rbind, Reduce, rownames, sapply, setdiff, sort, table, tapply, union, unique, unsplit, which.max, which.min Loading required package: S4Vectors Loading required package: stats4 Attaching package: ‘S4Vectors’ The following objects are masked from ‘package:base’: expand.grid, I, unname BSgenome::snpsById done in 30 seconds. Loading SNPlocs data. Loading reference genome data. Preprocessing RSIDs. Validating RSIDs of 10,000 SNPs using BSgenome::snpsById... BSgenome::snpsById done in 41 seconds. Inferred genome build: GRCH37 Checking SNP RSIDs. 4 SNP IDs are not correctly formatted. These will be corrected from the reference genome. Found Indels. These won't be checked against the reference genome as it does not contain Indels. WARNING If your sumstat doesn't contain Indels, set the indel param to FALSE & rerun MungeSumstats::format_sumstats() Loading SNPlocs data. Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/logs/snp_not_found_from_chr_bp.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. Checking for merged allele column. Checking A1 is uppercase Checking A2 is uppercase Checking for incorrect base-pair positions Ensuring all SNPs are on the reference genome. Loading SNPlocs data. Loading reference genome data. Preprocessing RSIDs. Validating RSIDs of 22,717,991 SNPs using BSgenome::snpsById... BSgenome::snpsById done in 361 seconds. Found 1,596,689 Indels. These won't be checked against the reference genome as it does not contain Indels. WARNING If your sumstat doesn't contain Indels, set the indel param to FALSE & rerun MungeSumstats::format_sumstats() Checking for correct direction of A1 (reference) and A2 (alternative allele). There are 285 SNPs where neither A1 nor A2 match the reference genome. These will be removed. Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/logs/alleles_dont_match_ref_gen.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. There are 171 SNPs where A1 doesn't match the reference genome. These will be flipped with their effect columns. Reordering so first three column headers are SNP, CHR and BP in this order. Reordering so the fourth and fifth columns are A1 and A2. Checking for missing data. Checking for duplicate columns. Checking for duplicate SNPs from SNP ID. 457 RSIDs are duplicated in the sumstats file. These duplicates will be removed Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/logs/dup_snp_id.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. Checking for SNPs with duplicated base-pair positions. 2 base-pair positions are duplicated in the sumstats file. These duplicates will be removed. Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/logs/dup_base_pair_position.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. INFO column not available. Skipping INFO score filtering step. Filtering SNPs, ensuring SE>0. Ensuring all SNPs have N<5 std dev above mean. Checking for bi-allelic SNPs. 499,519 SNPs are non-biallelic. These will be removed. Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/logs/snp_bi_allelic.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. Computing Z-score from P using formula: `sign(BETA)*sqrt(stats::qchisq(P,1,lower=FALSE)` Warning: When method is an integer, must be >0. Sorting coordinates with 'data.table'. Loading required namespace: GenomicFiles Using local VCF. File already tabix-indexed. Finding empty VCF columns based on first 10,000 rows. 1 sample detected: EBI-a-GCST005060 Constructing ScanVcfParam object. VCF contains: 23,315,353 variant(s) x 1 sample(s) Reading VCF file: multi-threaded (4 threads) Renaming ID as SNP. VCF file has -log10 P-values; these will be converted to unadjusted p-values in the 'P' column. No INFO (SI) column detected. Standardising column headers. First line of summary statistics file: SNP chr BP end REF ALT FILTER ReverseComplementedAlleles ES LP SE P Summary statistics report: - 22,718,565 rows - 22,717,995 unique variants - 132 genome-wide significant variants (P<5e-8) - 22 chromosomes Checking for multi-GWAS. Checking for multiple RSIDs on one row. Inferring genome build. Loading SNPlocs data. Loading reference genome data. Preprocessing RSIDs. Validating RSIDs of 10,000 SNPs using BSgenome::snpsById... Loading required package: BiocGenerics Attaching package: ‘BiocGenerics’ The following objects are masked from ‘package:stats’: IQR, mad, sd, var, xtabs The following objects are masked from ‘package:base’: anyDuplicated, aperm, append, as.data.frame, basename, cbind, colnames, dirname, do.call, duplicated, eval, evalq, Filter, Find, get, grep, grepl, intersect, is.unsorted, lapply, Map, mapply, match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank, rbind, Reduce, rownames, sapply, setdiff, sort, table, tapply, union, unique, unsplit, which.max, which.min Loading required package: S4Vectors Loading required package: stats4 Attaching package: ‘S4Vectors’ The following objects are masked from ‘package:base’: expand.grid, I, unname BSgenome::snpsById done in 65 seconds. Loading SNPlocs data. Loading reference genome data. Preprocessing RSIDs. Validating RSIDs of 10,000 SNPs using BSgenome::snpsById... BSgenome::snpsById done in 83 seconds. Inferred genome build: GRCH37 Checking SNP RSIDs. 4 SNP IDs are not correctly formatted. These will be corrected from the reference genome. Found Indels. These won't be checked against the reference genome as it does not contain Indels. WARNING If your sumstat doesn't contain Indels, set the indel param to FALSE & rerun MungeSumstats::format_sumstats() Loading SNPlocs data. Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/logs/snp_not_found_from_chr_bp.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. Checking for merged allele column. Checking A1 is uppercase Checking A2 is uppercase Checking for incorrect base-pair positions Coercing BP column to numeric. Ensuring all SNPs are on the reference genome. Loading SNPlocs data. Loading reference genome data. Preprocessing RSIDs. Validating RSIDs of 22,717,991 SNPs using BSgenome::snpsById... BSgenome::snpsById done in 318 seconds. Found 1,596,689 Indels. These won't be checked against the reference genome as it does not contain Indels. WARNING If your sumstat doesn't contain Indels, set the indel param to FALSE & rerun MungeSumstats::format_sumstats() Checking for correct direction of A1 (reference) and A2 (alternative allele). There are 285 SNPs where neither A1 nor A2 match the reference genome. These will be removed. Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/logs/alleles_dont_match_ref_gen.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. There are 171 SNPs where A1 doesn't match the reference genome. These will be flipped with their effect columns. Reordering so first three column headers are SNP, CHR and BP in this order. Reordering so the fourth and fifth columns are A1 and A2. Checking for missing data. Checking for duplicate columns. Checking for duplicate SNPs from SNP ID. 457 RSIDs are duplicated in the sumstats file. These duplicates will be removed Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/logs/dup_snp_id.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. Checking for SNPs with duplicated base-pair positions. 2 base-pair positions are duplicated in the sumstats file. These duplicates will be removed. Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/logs/dup_base_pair_position.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. INFO column not available. Skipping INFO score filtering step. Filtering SNPs, ensuring SE>0. Ensuring all SNPs have N<5 std dev above mean. Checking for bi-allelic SNPs. 499,519 SNPs are non-biallelic. These will be removed. Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/logs/snp_bi_allelic.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. Computing Z-score from P using formula: `sign(BETA)*sqrt(stats::qchisq(P,1,lower=FALSE)` Sorting coordinates with 'data.table'. Sorting coordinates with 'data.table'. Writing in tabular format ==> /rds/general/project/neurogenomics-lab/ephemeral/MAGMA_Files_Public/data/GWAS_munged/ebi-a-GCST005060/ebi-a-GCST005060.tsv Writing uncompressed instead of gzipped to enable tabix indexing. Converting full summary stats file to tabix format for fast querying... Reading header. Ensuring file is bgzipped. Tabix-indexing file. Removing temporary .tsv file. Summary statistics report: - 20,525,481 rows (90.3% of original 22,718,565 rows) - 20,525,481 unique variants - 129 genome-wide significant variants (P<5e-8) - 22 chromosomes Done munging in 46.03 minutes. Successfully finished preparing sumstats file, preview: Reading header. SNP CHR BP A1 A2 END FILTER REVERSECOMPLEMENTEDALLELES BETA 1: rs190291950 1 52144 T A 52144 PASS FALSE 0.95 2: rs185273034 1 85622 A T 85622 PASS FALSE 1.04 3: rs140735660 1 87755 G A 87755 PASS FALSE -0.22 4: rs183898652 1 99719 C T 99719 PASS FALSE 0.38 5: rs185940535 1 258589 C A 258589 PASS FALSE -0.36 LP SE P Z 1: 0.646970 0.79 0.22543949 1.2121904 2: 0.810775 0.73 0.15460552 1.4234508 3: 0.126603 0.68 0.74713142 -0.3224241 4: 1.298940 0.19 0.05024120 1.9579047 5: 1.229130 0.19 0.05900244 -1.8881751