Get SNP group counts

Tally locus-specific SNP group sizes.

get_SNPgroup_counts(merged_DT, grouping_vars = "Locus", verbose = TRUE)

Arguments

merged_DT

Merged fine-mapping results generated by echolocatoR::finemap_loci.

grouping_vars

Column names in merged_DT to group counts by.

verbose

Print messages.

Value

Named list of report tables.

See also

Other summarise: get_CS_bins(), get_CS_counts(), results_report()

Examples

merged_DT <- echodata::get_Nalls2019_merged()
snp_groups <- get_SNPgroup_counts(merged_DT = merged_DT)
#> All loci (75) : 
#>            Total.SNPs          nom.sig.GWAS              sig.GWAS 
#>               4948.16                924.68                 82.36 
#>                    CS             Consensus          topConsensus 
#>                  7.88                  2.69                  1.47 
#> topConsensus.leadGWAS 
#>                  0.41 
#> Loci with at least one Consensus SNP (69) : 
#>            Total.SNPs          nom.sig.GWAS              sig.GWAS 
#>               5019.07                911.41                 84.28 
#>                    CS             Consensus          topConsensus 
#>                  7.77                  2.93                  1.59 
#> topConsensus.leadGWAS 
#>                  0.45