In bailey-lab/miplicorn: A Framework for Molecular Inversion Probe and Amplicon Analysis

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(miplicorn)

File Sizes

In the world of genomic sequencing, files are often several gigabytes large containing millions of data points. Reading in such files to local machines, such as your laptop, can take an excruciating amount of time.

While there are programs that can handle large amounts of data, an easy and simple solution is to process your data in chunks. For instance, instead of looking at ten chromosomes simultaneously, it may be simpler to focus on two or three at a time.

Filters

The entire read_tbl_*() family of functions provide the ability to filter data so that data may load and run faster. This works by filtering even before objects are loaded into R. Data can be filtered using any of the information present in the metadata, and you may even filter on multiple conditions.

cov_file <- miplicorn_example("coverage_AA_table.csv")

read_tbl_coverage(cov_file)

read_tbl_coverage(cov_file, gene == "atp6")

read_tbl_coverage(cov_file, gene == "atp6", targeted == "Yes")

bailey-lab/miplicorn documentation built on March 19, 2023, 7:40 p.m.