tests/benchmarking.R

In rcorty/mvGWAS: Uses the DGLM packge (double generalized linear model) to do a GWAS, assess significance of results, and plot results

# system.time(expr = y <- read.vcfR(file = 'data/chr22_JHS_2_10Mb.vcf.gz'), gcFirst = TRUE)
# dim(y)
# print(x = object.size(x = y), units = 'Gb')
#
# system.time(expr = y <- read.vcfR(file = 'data/chr22_JHS_3_10Mb.vcf.gz'), gcFirst = TRUE)
# dim(y)
# print(x = object.size(x = y), units = 'Gb')
#
# system.time(expr = y <- read.vcfR(file = 'data/chr22_JHS_4_10Mb.vcf.gz'), gcFirst = TRUE)
# dim(y)
# print(x = object.size(x = y), units = 'Gb')
#
# system.time(expr = y <- read.vcfR(file = 'data/chr22_JHS_5_10Mb.vcf.gz'), gcFirst = TRUE)
# dim(y)
# print(x = object.size(x = y), units = 'Gb')
#
# system.time(expr = y <- read.vcfR(file = 'data/chr22_JHS_6_10Mb.vcf.gz'), gcFirst = TRUE)
# dim(y)
# print(x = object.size(x = y), units = 'Gb')