11.07.2015 Views

2DkcTXceO

2DkcTXceO

2DkcTXceO

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

N. Cressie 443is the set of all pixels in the Olentangy River watershed. Then squared-errorloss implies that E[1{Y (s C ) >k}|Z] =Pr{Y (s C ) >k|Z} given by (38.21) isoptimal for estimating 1{Y (s C ) >k}. Critically, other loss functions wouldyield different optimal estimates, since (38.7) depends on the loss functionL. A policy-maker’s question translated into a loss function yields a tailoredanswer to that question. Quite naturally in statistical decision theory, differentquestions are treated differently and result in different answers.More critically, the average (climate change over 30 years) Y can be replacedwith an extreme quantile, say the .95 quantile, which I denote here asg (.95) (Y ); this hidden process corresponds to temperature change that couldcause extreme stress to agricultural production in, for example, the HunterValley, NSW, Australia. Such projections for farmers in the “stressed” regionswould be invaluable for planning crop varieties that are more conducive tohigher temperature/lower rainfall conditions. That is, I propose making inferencedirectly on extremal processes, and decisions should be made withloss functions that are tailor-made to the typical “what if” queries made bydecision-makers.Furthermore, output could have non-Gaussian distributions; for example,quantiles of temperature or rainfall would be skewed, for which spatial generalisedlinear models (Diggle et al., 1998; Sengupta and Cressie, 2013) wouldbe well suited: In this framework, (38.20) is replaced with the data model,[Z(s) =z|Y,θ]=EF[z;E{Z(s)|Y (s),θ}], (38.22)which are conditionally independent for pixels s ∈ North America. In (38.22),EF denotes the one-parameter exponential family of probability distributions.Now consider a link function l that satisfies, l[E{Z(s)|Y (s),θ}] =Y (s); on thistransformed scale, climate change Y is modelled as the spatial process givenby (38.19). In Sengupta and Cressie (2013) and Sengupta et al. (2012), wehave developed spatial-statistical methodology for very large remote sensingdatasets based on (38.22), that could be adapted to RCM projections. Thatmethodology gives the predictive distribution, (38.10), which is summarised bymapping the predictive means, the predictive standard deviations (a measureof uncertainty), and the predictive extreme quantiles. Other data models couldalso be used in place of (38.20), such as the extreme-value distributions.Increases in temperature generally lead to decreases in water availability,due to an increase in evaporation. By developing conditional-probability distributionsof [Temperature] and [Rainfall | Temperature], we can infer thejoint behaviour of [Rainfall, Temperature]. This is in contrast to the bivariateanalysis in Sain et al. (2011), and it is a further example of the utility of aconditional-probability modelling approach, here embedded in a multivariatehierarchical statistical model.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!