The initial state SDs are assumed to follow half-Cauchy priors, C+(0,A), where the SDs may be common or distinct among the states.
sampleEvol0(mu0, evolParams0, commonSD = FALSE, A = 1)
List of relevant components:
the p x 1
evolution error SD sigma_w0
and the p x 1
parameter-expanded RV's px_sigma_w0
p x 1
vector of initial values (undifferenced)
list of relevant components (see below)
logical; if TRUE, use common SDs (otherwise distict)
prior scale parameter from the half-Cauchy prior, C+(0,A)
This function samples the parameters for a PX-Gibbs sampler.