plot.inla.climate: Plot-function for the 'inla.climate' model

In eirikmn/INLA.climate.test: Full Bayesian analysis of radiative forcing climate model

Description

Plots results obtained from the INLA.climate fit and the TCR sampling procedure.

Usage

## S3 method for class 'inla.climate'
plot(x,  plot.random=TRUE,
                                plot.ar=TRUE,
                                plot.hyperpars=TRUE,
                                plot.tcr.samples=TRUE,
                                plot.mu = TRUE,
                                postscript = FALSE,
                                pdf = FALSE,
                                prefix = "inla.climate.plots/figure-",
                                ...)

Arguments

x	inla.climate object returned by inla.climate function. Also works if the TCR sampling procedure has been skipped.
plot.random	Boolean variable indicating if marginal posterior means and 95% credible intervals for the latent GMRF corresponding to the AR(1) mixture should be plotted.
plot.ar	Boolean variable indicating if marginal posterior means and 95% credible intervals for the underlying AR(1) components should be plotted.
plot.hyperpars	Boolean indicating if the marginal posterior distributions of the hyperparameters should be plotted.
plot.tcr.samples	Boolean variable indicating if a histogram for the TCR samples should be plotted (should there be any).
plot.mu	Boolean variable indicating if the forcing response that can be interpreted as the mean of the fGn should be plotted. The observations are also included in the plot.
postscript	Boolean variable indicating if postscript files should be produced instead.
pdf	Boolean variable indicating if pdf files should be produced instead.
prefix	The prefix for created files. Additional numbering is added.
...	Additional arguments to postscripts(), pdf() or dev.new().

Author(s)

Eirik Myrvoll-Nilsen [e-mail]

See Also

inla.climate

Examples

if(require("INLA",quietly=TRUE)){
  data(GISS_E2_R)
  result.climate <- inla.climate(data=GISS_E2_R$Temperature,forcing=GISS_E2_R$Forcing)
  plot(result.climate)
}

eirikmn/INLA.climate.test documentation built on Feb. 2, 2020, 9:17 a.m.