DAEQTL mapping — daeqtl_mapping • daeqtlr

TODO

Usage

daeqtl_mapping(
  snp_pairs,
  zygosity,
  ae,
  fn = daeqtl_test,
  ...,
  .extra_cols = 2L,
  .n_chunks = future::nbrOfWorkers()
)

snp_pairs: A data frame with DAE SNP/candidate SNP pairs.
zygosity: A data frame with zygosity levels of samples. First column should be the SNP identity, and remaining columns should refer to biological samples. These should be factors, and their levels should be "hom" (homozygous) and "het" (heterozygous).
ae: A data frame with allelic expression (AE) ratios. Each row is for a locus. The locus identity is indicated in the first column. Remaining columns are samples.
fn: A function implementing the statistical association approach. This function needs to have two named arguments ae_hom and ae_het. Extra arguments will be read from ....
...: Extra arguments passed on to the call of the function passed in fn.
.extra_cols: The number of extra columns you are creating with fn.
.n_chunks: The number of chunks to divide snp_pairs for parallel processing. If you are going to run in parallel set this number to be equal to the total number of workers.

An updated version of the data frame snp_pairs. The update includes extra columns, typically pvalue and case.