QUALYPSOSS | R Documentation |
QUALYPSOSS
QUALYPSOSS( ClimateProjections, scenAvail, vecYears = NULL, predCont = NULL, predContUnique = NULL, iCpredCont = NULL, iCpredContUnique = NULL, listOption = NULL, RK = NULL )
ClimateProjections |
matrix |
scenAvail |
matrix of scenario characteristics |
vecYears |
(optional) vector of years of length |
predCont |
(optional) matrix |
predContUnique |
(optional) vector of length |
iCpredCont |
(optional) index in |
iCpredContUnique |
(optional) index in |
listOption |
(optional) list of options
|
RK |
Reproducing kernels: list |
list with the following fields:
POINT: list containing the mean estimate of different quantities: RESIDUALVAR
(residual variability), INTERNALVAR
(internal variability), GRANDMEAN
(grand mean for all time
steps), MAINEFFECT
(list with one item per discrete predictor i
, containing matrices nT
x
nEffi
, where nEffi
is the number of possible values for the discrete predictor i
).
EFFECTVAR
, uncertainty related to the different main effect, TOTVAR
Total variability,
DECOMPVAR
, decomposition of the total variability (percentages) for the different components,
CONTRIB_EACH_EFFECT
, contribution of each individual effects (percentages) to the corr. effect uncertainty.
BAYES: list containing quantiles of different estimated quantities, listed in POINT.
MCMC: list containing the MCMC chains (not returned by default).
climateResponse: list containing different objects related to the extraction of the climate response.
phiStar (φ^*) is an array nQ
x nS
x nP
containing climate change responses, where nQ
is the
number of returned quantiles, nS
is the number of scenarios and nP
is the length of predContUnique
(e.g. number
of future years).
Similarly, etaStar (η^*) contains the deviation from the climate change response.
phi (φ) contains the climate responses and eta (η) contains the deviations from the climate responses.
listCR: list containing objects created during the extraction of the climate responses
ClimateProjections: argument of the call to the function, for records.
predCont: (optional) argument of the call to the function, for records.
predContUnique: (optional) argument of the call to the function, for records.
predDiscreteUnique: list of possible values taken by the discrete predictors given in scenAvail
.
listOption: list of options
listScenario: list of scenario characteristics (obtained from QUALYPSOSS.process.scenario
)
RK: list containing the reproducing kernels
Guillaume Evin
########################################################################## # SYNTHETIC SCENARIOS ########################################################################## # create nS=3 fictive climate scenarios with 2 GCMs and 2 RCMs, for a period of nY=20 years n=20 t=1:n/n # GCM effects (sums to 0 for each t) effGCM1 = t*2 effGCM2 = t*-2 # RCM effects (sums to 0 for each t) effRCM1 = t*1 effRCM2 = t*-1 # These climate scenarios are a sum of effects and a random gaussian noise scenGCM1RCM1 = effGCM1 + effRCM1 + rnorm(n=n,sd=0.5) scenGCM1RCM2 = effGCM1 + effRCM2 + rnorm(n=n,sd=0.5) scenGCM2RCM1 = effGCM2 + effRCM1 + rnorm(n=n,sd=0.5) ClimateProjections = cbind(scenGCM1RCM1,scenGCM1RCM2,scenGCM2RCM1) # Here, scenAvail indicates that the first scenario is obtained with the combination of the # GCM "GCM1" and RCM "RCM1", the second scenario is obtained with the combination of # the GCM "GCM1" and RCM "RCM2" and the third scenario is obtained with the combination # of the GCM "GCM2" and RCM "RCM1". scenAvail = data.frame(GCM=c('GCM1','GCM1','GCM2'),RCM=c('RCM1','RCM2','RCM1')) listOption = list(nBurn=20,nKeep=30,type.temporal.dep="iid",type.hetero="constant") QUALYPSOSSOUT = QUALYPSOSS(ClimateProjections=ClimateProjections,scenAvail=scenAvail, listOption=listOption) # QUALYPSOSSOUT output contains many different information about climate projections uncertainties, # which can be plotted using the following functions. # plotQUALYPSOSSClimateResponse draws the climate responses, for all simulation chains, # in comparison to the raw climate responses. plotQUALYPSOSSClimateResponse(QUALYPSOSSOUT) # plotQUALYPSOSSClimateChangeResponse draws the climate change responses, for all simulation chains. plotQUALYPSOSSClimateChangeResponse(QUALYPSOSSOUT) # plotQUALYPSOSSeffect draws the estimated effects, for a discrete predictor specified by iEff, # as a function of the continuous predictor. plotQUALYPSOSSeffect(QUALYPSOSSOUT, iEff = 1) plotQUALYPSOSSeffect(QUALYPSOSSOUT, iEff = 2) # plotQUALYPSOSSgrandmean draws the estimated grand mean, as a function of the continuous predictor. plotQUALYPSOSSgrandmean(QUALYPSOSSOUT) # plotQUALYPSOSSTotalVarianceDecomposition draws the decomposition of the total variance responses, # as a function of the continuous predictor. plotQUALYPSOSSTotalVarianceDecomposition(QUALYPSOSSOUT)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.