inla.climate.mu: Temperature response estimation
In eirikmn/INLA.climate: Full Bayesian analysis of radiative forcing climate model

inla.climate.mu

R Documentation

Temperature response estimation

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 = NULL, 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 eirik.myrvoll-nilsen@uit.no

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