This function uses posterior samples to quantify uncertainty in the properties of fitted curves.

Usage

getStatsIntervals(posteriorStats, mass = 0.95, method = c("eti", "hdi"))

Arguments

posteriorStats

A data.table of statistics for posterior samples, as returned by getRhythmStats() or getDiffRhythmStats().

mass

Number between 0 and 1 indicating the probability mass for which to calculate the intervals.

method

String indicating the type of interval: 'eti' for equal-tailed using stats::quantile(), or 'hdi' for highest density using HDInterval::hdi().

Value

A data.table containing lower and upper bounds of various statistics for each feature or each feature-condition pair. For peak_trough_amp and rms_amp, a negative lower bound indicates a rhythm of the opposite phase.

getRhythmStats(), getDiffRhythmStats(), getExpectedMeasIntervals()

Examples

library('data.table')

y = GSE54650$y metadata = GSE54650$metadata