rasterConcentrationResponse: Raster Concentration Response

In christiaanpauw/AQoffset: Makes calculation needed for air quality offsetting

Description

Excecute a concentration response function on rasterised input data

Usage

rasterConcentrationResponse(conc, popr, base.conc, fun.form = c("linear",
  "log-linear")[1], beta = NULL, RR = NULL, delta = 10, cases = NULL,
  incidence.rate = NULL, est = c("low estimate", "cental estimate",
  "high estimate"), outcomename = NULL, pollutantname = NULL,
  unit.cost = NULL, costing.var = NULL, risk.only = FALSE,
  out = c("all", "popr", "cases", "conc", "RR", "AF", "AM")[1],
  verbose = FALSE, debug = FALSE, ...)

Arguments

conc	A raster of concentration values
popr	A raster of population numbers rasterConcentrationResponse
base.conc	A scalar giving the base concentration of the pollutant
fun.form	The form of the concentration response function. Currently supports "linear" (default) and "log-linear"
beta	The Beta values from the original regression model from which the concentration response function (CRF) is derived. Defaults to NULL. If NULL, it is calculated from RR and delta (beta = log(RR)/delta)
RR	A single numeric value. The relative risk associated with a concentration increase of delta.
delta	A single numeric value. The concentration difference for which RR is defined. Defaults to 10.
cases	A raster containing the number of cases of the outcome of concern in the population.
incidence.rate	The incidence rate of the outcome of concern.
est	The type of estimate that the values of RR represent. e.g. c("low estimate","cental estimate","high estimate"). Must be the same length as RR
outcomename	User specified name of the outcome. Defaults to NULL.
pollutantname	User specified name of the pollutant. Defaults to NULL.
unit.cost	The cost of the outcome per case
costing.var	The name of the varible containing the unit cost
risk.only	Logical. Returns the RR raster if TRUE
out	Character vector. The desired output. Defaults to "all"
verbose	Do you want to be annoyed with a lot of messages, or not?
debug	Set to TRUE if the function misbehaves so you can see what's going on
conc	A raster object of concentration values for a particular pollutant.
popr	A raster object of population numbers.

Details

Excecute a concentration response function on rasterised input data

pollutant This is used to match the name of the pollutant in the concentration raster with that in the concentration response information (typically given as a sicklist)

Value

IF risk.only, the function returns a raster stack with the same number of layers as conc; if not...? # CAN ONLY WORK WITH *ONE* OUTCOME ASSOCIATED WITH *ONE* POLLUTANT AT A TIME

Examples

cases = people$total * rnorm(n = ncell(people), mean = 0.01, sd = 0.003)

christiaanpauw/AQoffset documentation built on May 13, 2019, 6:57 p.m.