Provides statistical tests and algorithms for the detection of change points in time series and point processes - particularly for changes in the mean in time series and for changes in the rate and in the variance in point processes. References - Michael Messer, Marietta Kirchner, Julia Schiemann, Jochen Roeper, Ralph Neininger and Gaby Schneider (2014), A multiple filter test for the detection of rate changes in renewal processes with varying variance <doi:10.1214/14-AOAS782>. Stefan Albert, Michael Messer, Julia Schiemann, Jochen Roeper, Gaby Schneider (2017), Multi-scale detection of variance changes in renewal processes in the presence of rate change points <doi:10.1111/jtsa.12254>. Michael Messer, Kaue M. Costa, Jochen Roeper and Gaby Schneider (2017), Multi-scale detection of rate changes in spike trains with weak dependencies <doi:10.1007/s10827-016-0635-3>. Michael Messer, Stefan Albert and Gaby Schneider (2018), The multiple filter test for change point detection in time series <doi:10.1007/s00184-018-0672-1>. Michael Messer, Hendrik Backhaus, Albrecht Stroh and Gaby Schneider (2019+) Peak detection in time series.

Convenience functions for multivariate MCMC using univariate samplers including:
  slice sampler with stepout and shrinkage (Neal (2003) <DOI:10.1214/aos/1056562461>),
  adaptive rejection sampler (Gilks and Wild (1992) <DOI:10.2307/2347565>),
  adaptive rejection Metropolis (Gilks et al (1995) <DOI:10.2307/2986138>), and
  univariate Metropolis with Gaussian proposal.

Estimates probit, logit, Poisson, negative binomial, 
  and beta regression models, returning their marginal effects, odds ratios, 
  or incidence rate ratios as an output.
  Greene (2008, pp. 780-7) provides a textbook introduction to this topic.

The idea of a computational algorithm described in the article by Andronov M. et al. (2022) <https://link.springer.com/chapter/10.1007/978-3-030-92507-9_13>.
    The purpose of this package is to automate computations for a Markov-Modulated M/G/1 queuing system with alternating Poisson flow of arrivals. It offers a set of functions to calculate various mean indices of the system, including mean flow intensity, mean service busy and idle times, and the system's stationary probability.