inla.climate.mu: Temperature response estimation

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

Description

Computes Bayesian inference about temperature response of forcing given parameter information contained in result obtained from Monte Carlo simulations generated by sampling from the posterior marginal distributions obtained by INLA.

Usage

inla.climate.mu(result, forcing, quick=FALSE, T0.corr = 0, nsamples = 100000, seed = 1234, print.progress = FALSE, model="fgn")

Arguments

result	An inla.climate or an inla object which can be accessed using results$inla.result, where results is the returned object from the inla.climate function.
forcing	The forcing for which the temperature response will be obtained.
quick	Boolean variable. If FALSE, all Monte Carlo simulations will be stored. If TRUE, the quantiles will not be computed.
T0.corr	The temperature dataset analysed with INLA is shifted by T0 = mean(data[1:20]) and the results generated by inla.climate.mu be shifted accordingly. If this variable is set equal to T0, the results will be aligned with the original input.
nsamples	The number of Monte Carlo simulations used in obtaining transient climate response inference.
seed	Seed used for random number generator.
print.progress	Prints progression if TRUE.
model	String which determines which model to be used when computing the temperature response. "fgn" and "arfima" produces a long-range dependent response, whilst "ar1" computes it as a sum of exponentials.

Value

Returns a list of class inla.climate.mu which contains the following objects:

mu.mean	The mean value of the temperature response due to forcing.
mu.sd	The standard deviation of the temperature response due to forcing.
mu.quant0.025	The 2.5% quantile value of the temperature response due to forcing.
mu.quant0.5	The 50% quantile of the temperature response due to forcing.
mu.quant0.975	The 97.5% quantile of the temperature response due to forcing.
samples	List containing the Monte Carlo simulations of the TCR and all hyperparameters.
time	The time spent doing the Monte Carlo simulations.

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)
  result.mu <- inla.climate.mu(result.climate,forcing = GISS_E2_R$Forcing,quick = FALSE)
}

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