Fig. 3

Example output summary from ReTrOS-Switch applied to protein data. Data is luciferase tagged protein for the LHY clock gene in wild-type Arabidopsis plants (ROBuST Project) [15]. Panel a) shows the raw input data (blue circles, extremes of replicates shown in shaded region and median shown with line), the fitted protein expression (black line) and the estimated switch events (μ±1.96σ as vertical red/green lines). Panel b) shows the back-calculated mRNA expression profiles (black line) with the estimated switch times. Panel c) shows the baseline-removed estimated switch time probability density (black line) and the estimated switch times Gaussian mixture model (blue line, with shaded red/green region indicating μ±1.96σ). Panel d) shows the accepted switch time samples, along with the switch events (red/green). The burn-in period of the chain is shown by the dashed black line. Panel e) shows the accepted samples from the precision (σ ²), δ _m and δ _p parameter chains in blue and a histogram of the accepted model size. Summaries of the parameter chains are shown in black (median and lower/upper quartiles) and red (μ and 1.96 σ). If mRNA expression data is used, specific output panels are removed accordingly. The model in Eq. (9) is used here. The default switch parameter values in the parameter setting GUI are used, with the prior mRNA and protein degradation rates set by choosing luc reporter