pairscan_null_kin: Generates a null distribution for the pairscan

Description Usage Arguments Value

View source: R/pairscan_null_kin.R

Description

This function generates a null distribution for the pairscan. For each permutation, it runs a single scan and selects markers in the same manner as for the true test.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
pairscan_null_kin(
  data_obj,
  geno_obj = NULL,
  kin_obj = NULL,
  scan_what = c("eigentraits", "raw_traits"),
  pairscan_null_size = NULL,
  max_pair_cor = NULL,
  min_per_geno = NULL,
  model_family = "gaussian",
  marker_selection_method = c("top_effects", "uniform", "effects_dist", "by_gene"),
  run_parallel = FALSE,
  n_cores = 4,
  verbose = FALSE
)

Arguments

data_obj

a Cape object

geno_obj

a genotype object

kin_obj

a kinship object

scan_what

A character string uniquely identifying whether eigentraits or raw traits should be scanned. Options are "eigentraits", "raw_traits"

pairscan_null_size

The total size of the null distribution. This is DIFFERENT than the number of permutations to run. Each permutation generates n choose 2 elements for the pairscan. So for example, a permutation that tests 100 pairs of markers will generate a null distribution of size 4950. This process is repeated until the total null size is reached. If the null size is set to 5000, two permutations of 100 markers would be done to get to a null distribution size of 5000.

max_pair_cor

A numeric value between 0 and 1 indicating the maximum Pearson correlation that two markers are allowed. If the correlation between a pair of markers exceeds this threshold, the pair is not tested. If this value is set to NULL, min_per_genotype must have a numeric value.

min_per_geno

The minimum number of individuals allowable per genotype. If for a given marker pair, one of the genotypes is underrepresented, the marker pair is not tested. If this value is NULL, max_pair_cor must have a numeric value.

model_family

Indicates the model family of the phenotypes. This can be either "gaussian" or "binomial".

marker_selection_method

options are "top_effects", "uniform", "effects_dist", "by_gene"

run_parallel

Whether to run the analysis on multiple CPUs.

n_cores

The number of CPUs to use if run_parallel is TRUE.

verbose

Whether to print progress to the screen. Defaults to FALSE.

Value

This function returns a list with two elements, one containing the results of the permutations, and the other containing the markers that were tested in the individual permutations.


cape documentation built on Feb. 10, 2021, 5:06 p.m.