##############################################################################
##############################################################################
###
### Running command:
###
###   /home/biocbuild/R/R-4.3.1/bin/R CMD check --install=check:MungeSumstats.install-out.txt --library=/home/biocbuild/R/R-4.3.1/site-library --no-vignettes --timings MungeSumstats_1.10.1.tar.gz
###
##############################################################################
##############################################################################


* using log directory ‘/home/biocbuild/bbs-3.18-bioc/meat/MungeSumstats.Rcheck’
* using R version 4.3.1 (2023-06-16)
* using platform: aarch64-unknown-linux-gnu (64-bit)
* R was compiled by
    gcc (GCC) 10.3.1
    GNU Fortran (GCC) 10.3.1
* running under: openEuler 22.03 (LTS-SP1)
* using session charset: UTF-8
* using option ‘--no-vignettes’
* checking for file ‘MungeSumstats/DESCRIPTION’ ... OK
* checking extension type ... Package
* this is package ‘MungeSumstats’ version ‘1.10.1’
* package encoding: UTF-8
* checking package namespace information ... OK
* checking package dependencies ... 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 for sufficient/correct file permissions ... OK
* checking whether package ‘MungeSumstats’ 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
* 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 loading without being on the library search path ... OK
* checking startup messages can be suppressed ... 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 R/sysdata.rda ... OK
* checking files in ‘vignettes’ ... OK
* checking examples ... ERROR
Running examples in ‘MungeSumstats-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: format_sumstats
> ### Title: Check that summary statistics from GWAS are in a homogeneous
> ###   format
> ### Aliases: format_sumstats
> 
> ### ** Examples
> 
> # Pass path to Educational Attainment Okbay sumstat file to a temp directory
> 
> eduAttainOkbayPth <- system.file("extdata", "eduAttainOkbay.txt",
+     package = "MungeSumstats"
+ )
> 
> ## Call uses reference genome as default with more than 2GB of memory,
> ## which is more than what 32-bit Windows can handle so remove certain checks
> ## Using dbSNP = 144 for speed as it's smaller but you should use 155 unless
> ## you know what you are doing and need 144
> 
> is_32bit_windows <-
+     .Platform$OS.type == "windows" && .Platform$r_arch == "i386"
> if (!is_32bit_windows) {
+     reformatted <- format_sumstats(
+         path = eduAttainOkbayPth,
+         ref_genome = "GRCh37",
+         dbSNP = 144
+     )
+ } else {
+     reformatted <- format_sumstats(
+         path = eduAttainOkbayPth,
+         ref_genome = "GRCh37",
+         on_ref_genome = FALSE,
+         strand_ambig_filter = FALSE,
+         bi_allelic_filter = FALSE,
+         allele_flip_check = FALSE,
+         dbSNP=144
+     )
+ }


******::NOTE::******
 - Formatted results will be saved to `tempdir()` by default.
 - This means all formatted summary stats will be deleted upon ending the R session.
 - To keep formatted summary stats, change `save_path`  ( e.g. `save_path=file.path('./formatted',basename(path))` ),   or make sure to copy files elsewhere after processing  ( e.g. `file.copy(save_path, './formatted/' )`.
 ******************** 

Formatted summary statistics will be saved to ==>  /home/biocbuild/tmp/RtmpsNIurT/file1a967a112d71d8.tsv.gz
Warning: replacing previous import ‘utils::findMatches’ by ‘S4Vectors::findMatches’ when loading ‘SNPlocs.Hsapiens.dbSNP144.GRCh37’
Importing tabular file: /home/biocbuild/R/R-4.3.1/site-library/MungeSumstats/extdata/eduAttainOkbay.txt
Checking for empty columns.
Infer Effect Column
First line of summary statistics file: 
MarkerName	CHR	POS	A1	A2	EAF	Beta	SE	Pval	
Allele columns are ambiguous, attempting to infer direction
Standardising column headers.
First line of summary statistics file: 
MarkerName	CHR	POS	A1	A2	EAF	Beta	SE	Pval	
Loading SNPlocs data.
Loading reference genome data.
Preprocessing RSIDs.
Validating RSIDs of 93 SNPs using BSgenome::snpsById...
* checking for unstated dependencies in ‘tests’ ... OK
* checking tests ...
  Running ‘testthat.R’/home/biocbuild/R/R-4.3.1/bin/BATCH: line 60: 1745512 Killed                  ${R_HOME}/bin/R -f ${in} ${opts} ${R_BATCH_OPTIONS} > ${out} 2>&1

 ERROR
Running the tests in ‘tests/testthat.R’ failed.
Last 13 lines of output:
     - 93 unique variants
     - 70 genome-wide significant variants (P<5e-8)
     - 20 chromosomes
  Checking for multi-GWAS.
  Checking for multiple RSIDs on one row.
  Checking SNP RSIDs.
  Checking for merged allele column.
  Checking A1 is uppercase
  Checking A2 is uppercase
  Checking for incorrect base-pair positions
  Checking for correct direction of A1 (reference) and A2 (alternative allele).
  Loading SNPlocs data.
  Loading reference genome data.
  Preprocessing RSIDs.
  Validating RSIDs of 93 SNPs using BSgenome::snpsById...
* 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: 2 ERRORs
See
  ‘/home/biocbuild/bbs-3.18-bioc/meat/MungeSumstats.Rcheck/00check.log’
for details.

##############################################################################
##############################################################################
###
### Running command:
###
###   /home/biocbuild/R/R-4.3.1/bin/R CMD INSTALL MungeSumstats
###
##############################################################################
##############################################################################


* installing to library ‘/home/biocbuild/R/R-4.3.1/site-library’
* installing *source* package ‘MungeSumstats’ ...
** using staged installation
** R
** data
** inst
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
** building package indices
** installing vignettes
** testing if installed package can be loaded from temporary location
** testing if installed package can be loaded from final location
** testing if installed package keeps a record of temporary installation path
* DONE (MungeSumstats)

R version 4.3.1 (2023-06-16) -- "Beagle Scouts"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: aarch64-unknown-linux-gnu (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.

> library(testthat)
> library(MungeSumstats)
> 
> test_check("MungeSumstats")
Collecting metadata from Open GWAS.
Filtering metadata by substring criteria.
Filtering metadata by sample/case/control/SNP size criteria.
Excluding sample/case/control size with NAs.
Found 3 GWAS datasets matching search criteria across:
   - 3 trait(s)
   - 1 population(s)
   - 2 category(ies)
   - 2 subcategory(ies)
   - 2 publication(s)
   - 2 consortia(ium)
   - 1 genome build(s)
Collecting metadata from Open GWAS.
Filtering metadata by substring criteria.
Found 49 GWAS datasets matching search criteria across:
   - 44 trait(s)
   - 4 population(s)
   - 2 category(ies)
   - 2 subcategory(ies)
   - 9 publication(s)
   - 5 consortia(ium)
   - 1 genome build(s)
Downloading VCF ==> /home/biocbuild/tmp/RtmpkeCfF1/ieu-a-298.vcf.gz
Downloading with download.file.
trying URL 'https://gwas.mrcieu.ac.uk/files/ieu-a-298/ieu-a-298.vcf.gz'
Content type 'application/gzip' length 234480 bytes (228 KB)
==================================================
downloaded 228 KB

Downloading VCF index ==> https://gwas.mrcieu.ac.uk/files/ieu-a-298/ieu-a-298.vcf.gz.tbi
Downloading with download.file.
trying URL 'https://gwas.mrcieu.ac.uk/files/ieu-a-298/ieu-a-298.vcf.gz.tbi'
Content type 'application/gzip' length 37803 bytes (36 KB)
==================================================
downloaded 36 KB

Processing 1 datasets from Open GWAS.

========== Processing dataset : a-fake-id ==========

Downloading VCF ==> /home/biocbuild/tmp/RtmpkeCfF1/a-fake-id.vcf.gz
Downloading with download.file.
trying URL 'https://gwas.mrcieu.ac.uk/files/a-fake-id/a-fake-id.vcf.gz'
Processing 1 datasets from Open GWAS.

========== Processing dataset : ieu-a-298 ==========

Using previously downloaded VCF.
Formatted summary statistics will be saved to ==>  /home/biocbuild/tmp/RtmpkeCfF1/ieu-a-298/ieu-a-298.tsv.gz


******::NOTE::******
 - Formatted results will be saved to `tempdir()` by default.
 - This means all formatted summary stats will be deleted upon ending the R session.
 - To keep formatted summary stats, change `save_path`  ( e.g. `save_path=file.path('./formatted',basename(path))` ),   or make sure to copy files elsewhere after processing  ( e.g. `file.copy(save_path, './formatted/' )`.
 ******************** 

Formatted summary statistics will be saved to ==>  /home/biocbuild/tmp/RtmpkeCfF1/file1aa268672a946b.tsv.gz
Reading header.
Tabular format detected.
Importing tabular file: /home/biocbuild/tmp/RtmpkeCfF1/file1aa26863d5a5b4
Checking for empty columns.
Infer Effect Column
First line of summary statistics file: 
MarkerName	CHR	POS	A0	A1	EAF	Beta	SE	Pval	
Standardising column headers.
First line of summary statistics file: 
MarkerName	CHR	POS	A0	A1	EAF	Beta	SE	Pval	
Summary statistics report:
   - 93 rows
   - 93 unique variants
   - 70 genome-wide significant variants (P<5e-8)
   - 20 chromosomes
Checking for multi-GWAS.
Checking for multiple RSIDs on one row.
Checking SNP RSIDs.
Checking for merged allele column.
Checking A1 is uppercase
Checking A2 is uppercase
Checking for incorrect base-pair positions
Checking for missing data.
Checking for duplicate columns.
Checking for duplicated rows.
INFO column not available. Skipping INFO score filtering step.
Filtering SNPs, ensuring SE>0.
Ensuring all SNPs have N<5 std dev above mean.
47 SNPs (50.5%) have FRQ values > 0.5. Conventionally the FRQ column is intended to show the minor/effect allele frequency.
The FRQ column was mapped from one of the following from the inputted  summary statistics file:
FRQ, EAF, FREQUENCY, FRQ_U, F_U, MAF, FREQ, FREQ_TESTED_ALLELE, FRQ_TESTED_ALLELE, FREQ_EFFECT_ALLELE, FRQ_EFFECT_ALLELE, EFFECT_ALLELE_FREQUENCY, EFFECT_ALLELE_FREQ, EFFECT_ALLELE_FRQ, A2FREQ, A2FRQ, ALLELE_FREQUENCY, ALLELE_FREQ, ALLELE_FRQ, AF, MINOR_AF, EFFECT_AF, A2_AF, EFF_AF, ALT_AF, ALTERNATIVE_AF, INC_AF, A_2_AF, TESTED_AF, ALLELEFREQ, ALT_FREQ, EAF_HRC, EFFECTALLELEFREQ, FREQ.B, FREQ_EUROPEAN_1000GENOMES, FREQ_HAPMAP, FREQ_TESTED_ALLELE_IN_HRS, FRQ_U_113154, FRQ_U_31358, FRQ_U_344901, FRQ_U_43456, POOLED_ALT_AF, AF_ALT, AF.ALT, AF-ALT, ALT.AF, ALT-AF, A2.AF, A2-AF, AF.EFF, AF_EFF, ALL_AF
As frq_is_maf=TRUE, the FRQ column will not be renamed. If the FRQ values were intended to represent major allele frequency,
set frq_is_maf=FALSE to rename the column as MAJOR_ALLELE_FRQ and differentiate it from minor/effect allele frequency.
Sorting coordinates with 'data.table'.
Writing in tabular format ==> /home/biocbuild/tmp/RtmpkeCfF1/file1aa268672a946b.tsv.gz
Summary statistics report:
   - 93 rows (100% of original 93 rows)
   - 93 unique variants
   - 70 genome-wide significant variants (P<5e-8)
   - 20 chromosomes
Done munging in 0.079 minutes.
Successfully finished preparing sumstats file, preview:
Reading header.
          SNP CHR       BP A1 A2     FRQ   BETA    SE         P
1:   rs301800   1  8490603  T  C 0.17910  0.019 0.003 1.794e-08
2: rs11210860   1 43982527  A  G 0.36940  0.017 0.003 2.359e-10
3: rs34305371   1 72733610  A  G 0.08769  0.035 0.005 3.762e-14
4:  rs2568955   1 72762169  T  C 0.23690 -0.017 0.003 1.797e-08
Returning path to saved data.


******::NOTE::******
 - Formatted results will be saved to `tempdir()` by default.
 - This means all formatted summary stats will be deleted upon ending the R session.
 - To keep formatted summary stats, change `save_path`  ( e.g. `save_path=file.path('./formatted',basename(path))` ),   or make sure to copy files elsewhere after processing  ( e.g. `file.copy(save_path, './formatted/' )`.
 ******************** 

Formatted summary statistics will be saved to ==>  /home/biocbuild/tmp/RtmpkeCfF1/file1aa2685aebfa0f.tsv.gz
Reading header.
Tabular format detected.
Importing tabular file: /home/biocbuild/tmp/RtmpkeCfF1/file1aa26863d5a5b4
Checking for empty columns.
Infer Effect Column
First line of summary statistics file: 
MarkerName	CHR	POS	A1	A2	EAF	Beta	SE	Pval	
Allele columns are ambiguous, attempting to infer direction
Can't infer allele columns from sumstats
Standardising column headers.
First line of summary statistics file: 
MarkerName	CHR	POS	A1	A2	EAF	Beta	SE	Pval	
Summary statistics report:
   - 93 rows
   - 93 unique variants
   - 70 genome-wide significant variants (P<5e-8)
   - 20 chromosomes
Checking for multi-GWAS.
Checking for multiple RSIDs on one row.
Checking SNP RSIDs.
Checking for merged allele column.
Checking A1 is uppercase
Checking A2 is uppercase
Checking for incorrect base-pair positions
Checking for missing data.
Checking for duplicate columns.
Checking for duplicated rows.
INFO column not available. Skipping INFO score filtering step.
Filtering SNPs, ensuring SE>0.
Ensuring all SNPs have N<5 std dev above mean.
47 SNPs (50.5%) have FRQ values > 0.5. Conventionally the FRQ column is intended to show the minor/effect allele frequency.
The FRQ column was mapped from one of the following from the inputted  summary statistics file:
FRQ, EAF, FREQUENCY, FRQ_U, F_U, MAF, FREQ, FREQ_TESTED_ALLELE, FRQ_TESTED_ALLELE, FREQ_EFFECT_ALLELE, FRQ_EFFECT_ALLELE, EFFECT_ALLELE_FREQUENCY, EFFECT_ALLELE_FREQ, EFFECT_ALLELE_FRQ, A2FREQ, A2FRQ, ALLELE_FREQUENCY, ALLELE_FREQ, ALLELE_FRQ, AF, MINOR_AF, EFFECT_AF, A2_AF, EFF_AF, ALT_AF, ALTERNATIVE_AF, INC_AF, A_2_AF, TESTED_AF, ALLELEFREQ, ALT_FREQ, EAF_HRC, EFFECTALLELEFREQ, FREQ.B, FREQ_EUROPEAN_1000GENOMES, FREQ_HAPMAP, FREQ_TESTED_ALLELE_IN_HRS, FRQ_U_113154, FRQ_U_31358, FRQ_U_344901, FRQ_U_43456, POOLED_ALT_AF, AF_ALT, AF.ALT, AF-ALT, ALT.AF, ALT-AF, A2.AF, A2-AF, AF.EFF, AF_EFF, ALL_AF
As frq_is_maf=TRUE, the FRQ column will not be renamed. If the FRQ values were intended to represent major allele frequency,
set frq_is_maf=FALSE to rename the column as MAJOR_ALLELE_FRQ and differentiate it from minor/effect allele frequency.
Sorting coordinates with 'data.table'.
Writing in tabular format ==> /home/biocbuild/tmp/RtmpkeCfF1/file1aa2685aebfa0f.tsv.gz
Summary statistics report:
   - 93 rows (100% of original 93 rows)
   - 93 unique variants
   - 70 genome-wide significant variants (P<5e-8)
   - 20 chromosomes
Done munging in 0.073 minutes.
Successfully finished preparing sumstats file, preview:
Reading header.
          SNP CHR       BP A1 A2     FRQ   BETA    SE         P
1:   rs301800   1  8490603  T  C 0.17910  0.019 0.003 1.794e-08
2: rs11210860   1 43982527  A  G 0.36940  0.017 0.003 2.359e-10
3: rs34305371   1 72733610  A  G 0.08769  0.035 0.005 3.762e-14
4:  rs2568955   1 72762169  T  C 0.23690 -0.017 0.003 1.797e-08
Returning path to saved data.


******::NOTE::******
 - Formatted results will be saved to `tempdir()` by default.
 - This means all formatted summary stats will be deleted upon ending the R session.
 - To keep formatted summary stats, change `save_path`  ( e.g. `save_path=file.path('./formatted',basename(path))` ),   or make sure to copy files elsewhere after processing  ( e.g. `file.copy(save_path, './formatted/' )`.
 ******************** 

Formatted summary statistics will be saved to ==>  /home/biocbuild/tmp/RtmpkeCfF1/file1aa268fcd72ea.tsv.gz
Reading header.
Tabular format detected.
Importing tabular file: /home/biocbuild/tmp/RtmpkeCfF1/file1aa2686bafb8cf
Checking for empty columns.
Infer Effect Column
First line of summary statistics file: 
MarkerName	CHR	POS	A2	A1	EAF	Beta	SE	Pval	
Allele columns are ambiguous, attempting to infer direction
Can't infer allele columns from sumstats
Standardising column headers.
First line of summary statistics file: 
MarkerName	CHR	POS	A2	A1	EAF	Beta	SE	Pval	
Summary statistics report:
   - 93 rows
   - 93 unique variants
   - 70 genome-wide significant variants (P<5e-8)
   - 20 chromosomes
Checking for multi-GWAS.
Checking for multiple RSIDs on one row.
Checking SNP RSIDs.
Checking for merged allele column.
Checking A1 is uppercase
Checking A2 is uppercase
Checking for incorrect base-pair positions
Checking for correct direction of A1 (reference) and A2 (alternative allele).
Loading SNPlocs data.
Loading reference genome data.
Preprocessing RSIDs.
Validating RSIDs of 93 SNPs using BSgenome::snpsById...