bakR (Q69251)

Several implementations of a novel Bayesian hierarchical statistical model of nucleotiderecoding RNA-seq experiments (NR-seq; TimeLapse-seq, SLAM-seq, TUC-seq, etc.)for analyzing and comparing NR-seq datasets (see 'Vock and Simon' (2023) <doi:10.1261/rna.079451.122>).NR-seq is a powerful extension of RNA-seq that provides information about the kineticsof RNA metabolism (e.g., RNA degradation rate constants), which is notably lackingin standard RNA-seq data. The statistical model makes maximal use of these high-throughputdatasets by sharing information across transcripts to significantly improveuncertainty quantification and increase statistical power. 'bakR' includes a maximallyefficient implementation of this model for conservative initial investigations of datasets. 'bakR'also provides more highly powered implementations using the probabilistic programming language'Stan' to sample from the full posterior distribution. 'bakR' performs multiple-testadjusted statistical inference with the output of these model implementations tohelp biologists separate signal from background. Methods to automatically visualize keyresults and detect batch effects are also provided.