reedssorenson/RNAdecay: Maximum Likelihood Decay Modeling of RNA Degradation Data
Version 0.99.20

RNA degradation is monitored through measurement of RNA abundance after inhibiting RNA synthesis. This package has functions and example scripts to facilitate (1) data normalization, (2) data modeling using constant decay rate or time-dependent decay rate models, (3) the evaluation of treatment or genotype effects, and (4) plotting of the the data and models. Data Normalization: functions and scripts make easy the normalization to the initial (T0) RNA abundance, as well as a method to correct for artificial inflation of Reads per Million (RPM) abundance in global assesements as the total size of the RNA pool deacreases. Modeling: Normalized data is then modeled using maximum likelihood to fit parameters. For making treatment or genotype comparisons (up to four), the modeling step models all possible treatement effects on each gene by repeating the modeling with constraints on the model parameters (i.e., the decay rate of treatments A and B are modeled once with them being equal and again allowing them to both vary independently). Model Selection: The AICc value is calculated for each model, and the model with the lowest AICc is chosen. Modeling results of selected models are then compiled into a single data frame. Graphical Plotting: a function is provided to easily visualize the data and the selected model using ggplot2 package functions.

Getting started

Package details

Bioconductor views DifferentialExpression GeneExpression GeneRegulation Normalization RNASeq Regression Software TimeCourse Transcription Transcriptomics WorkflowStep
Package repositoryView on GitHub
Installation Install the latest version of this package by entering the following in R:
reedssorenson/RNAdecay documentation built on April 21, 2018, 10:29 p.m.