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

View source: R/plot.inla.climate.R

plot.inla.climateR Documentation

Plot-function for the inla.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 eirik.myrvoll-nilsen@uit.no

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 documentation built on Feb. 6, 2023, 11:41 a.m.