plot.tuts_poisBFS: Plots and visual diagnostics of tuts_BFS objects

In tuts: Time Uncertain Time Series Analysis

Description

plot.tuts_poisBFS generates plots and visual diagnostics of tuts_BFS objects.

Usage

## S3 method for class 'tuts_poisBFS'
plot(x, type, ...)

Arguments

x	A tuts_BFS objects.
type	plot type with the following options: - 'periodogram' plots estimates of power spectrum. - 'predTUTS' plots one step predictions of the model. - 'GR' plots Gelman-Rubin diagnostics. - 'cv' plots 5-fold cross validation. - 'mcmc' plots diagnostics of MCMC/JAGS objects. - 'lambda' plots lambda realizations.
...	list of optional parameters: - burn: burn-in parameter ranging from 0 to 0.7 with default value set to 0. - CI: credible interval ranging from 0.3 to 1 with default value set to 0.95.

Examples

# Note: Most of models included in tuts package are computationally intensive. In the example
# below I set parameters to meet CRAN<e2><80><99>s testing requirement of maximum 5 sec per example.
# A more practical example would contain N=50 in the first line of the code and n.sim=10000.

#1. Import or simulate the data (simulation is chosen for illustrative purposes):
DATA=simtuts(N=7,Harmonics=c(2,0,0), sin.ampl=c(10,0, 0), cos.ampl=c(0,0,0),
            trend=0,y.sd=2, ti.sd=0.2)
y=DATA$observed$y.obs
y=round(y-min(y))
ti.mu=DATA$observed$ti.obs.tnorm
ti.sd= rep(0.2, length(ti.mu))

#2. Fit the model:
n.sim=10
TUPOIS=tupoisbsf(y=y,ti.mu=ti.mu,ti.sd=ti.sd,freqs='internal',n.sim=n.sim,n.chains=2,
                 CV=TRUE,n.cores=2)

#3. Generate plots and diagnostics of the model (optional parameters are listed in brackets):
plot(TUPOIS,type='periodogram')               # spectral analysis (CI, burn).
plot(TUPOIS,type='predTUTS', CI=0.99)         # One step predictions (CI, burn).
plot(TUPOIS,type='cv')                        # 5 fold cross validation (CI, burn).
plot(TUPOIS,type='GR')                        # Gelman-Rubin diagnostics (CI, burn).
plot(TUPOIS,type='mcmc')                      # MCMC diagnostics.
plot(TUPOIS,type='lambda')                    # Realizaitons of lambda.

tuts documentation built on May 1, 2019, 7:56 p.m.