gibbSampler: Gibbs sampler for mixture model.

Description Usage Arguments Examples

View source: R/gibbSampler.R

Description

Gibbs sampler for mixture model.

Usage

1
2
3
4
gibbSampler(priors = list(theta = list(m0 = 0, s0 = 1), sigma2 = list(a0 = 2,
  b0 = 1)), y, p.sites, season.sites, basis, ID, N.run = 10000,
  mixture = TRUE, meanFunc = NULL, debug = FALSE, print.res = FALSE,
  ...)

Arguments

y

data

p.sites

the composition of each site.

season.sites

season data for each site.

ID

ID vector of the sites.

N.run

Number of run of the sampler.

theta

parameters to be estimated

Examples

1
2
data("wq_analysis_week2")
SPTMData(wq.raw.obs, frequency = "quarter")

ick003/SpTMixture documentation built on May 18, 2019, 2:32 a.m.