All functions

count_by_maf()

count number of loci below and above a maf threshold

fix_allele_encoding()

Fixes allele encoding so that the minor allele homozygote is always encoded as 2

get_effects()

Get weights for markers using provided distribution function.

get_genotypes()

Get genotypes matrix (counts of reference allele) from a GenABEL gwaa-data, vcfR object or an object imported using seqminer

get_genotypes_gwaa_data()

Get genotypes matrix (counts of reference allele) from a GenABEL gwaa-data object

get_genotypes_seqminer_matrix()

Get genotypes matrix (counts of reference allele) from one of the seqminer read functions that return genotype matrix

get_genotypes_vcfR()

Get genotypes matrix (counts of reference allele) from a vcfR object

get_maf()

Get minor allele frequency vector

hush()

Silence too verbose code

impute_G()

Impute missing genotypes using minor allele frequency

plot_effects()

Plots effect sizes distribution dependent on minor allele frequency.

plot_mrk_density()

Plot density of markers in a region

read_region_vcf()

read specific region from VCF file to a gwasim type genotype matrix

read_regions_bed()

Read regions information from bed file.

recode_G()

Re-code genotypes according to the supplied genotype-phenotype map

rnorm_betas()

Effect sizes $$ from a normal distribution.

scan_CommonRare()

Run SKAT-O CommonRare model on all regions

sim_y_in_regions()

Simulate phenotypes in a number of regions

simulate_phenotype() simulate_phenotype()

Simulate phenotype given a set of common and rare variants.

tibble_to_gwaa()

Convert tibble with genotypes along with a vector of phenotypes to GenABEL internal raw gwaa.data format.

tibble_to_raw_genotypes()

Convert tibble with genotypes to GenABEL internal raw genotype format and writ it into a file.

validate_maf()

Validate minor allele frequency vector