add_DRD: Add deletion rate difference

In SchwartzLab/brainSTORM: brainSTORM. A package to process STORM-seq based data

Add deletion rate difference

Description

Calculated the deletion rate to nucleotide reads for each position and subtracts it for the rate in the control group.

Usage

add_DRD(
  STORM,
  group_A,
  group_B,
  newColName = "auto",
  onNucs = c("A", "C", "G", "T"),
  minReadCounts = 20
)

Arguments

STORM	list. STORM object as output by storm_STORM
group_A	character. Name of group to be compared as found in STORM$META$group
group_B	character. Name of group to be compared as found in STORM$META$group
newColName	character. Name of calculated metric to be stored in STORM$RES, assigned by default based on group_A and group_B
onNucs	character. Nucleotide(s) in which the metric will be calculated
minReadCounts	integer. Minimal number of nucleotide reads (NOT coverage) in position for metric to be calculated

SchwartzLab/brainSTORM documentation built on May 14, 2022, 5:14 p.m.