|
DI plot.probCreating a Probability-of-Exceedence plotWe begin with help on the command line side and then present a demonstration of how to plot probabilities of exceeding a given value using the GUI display dialog.Plot Probabilities for off-network values. DESCRIPTION: Creates image plot for spatial probabilities. USAGE: plot.prob(nw.obj=new.network, row.name=NULL, mean.obj=0, sd.obj=1, add.map=T, add.county=T, ...) REQUIRED ARGUMENTS: nw.obj A network object usually created by make.network.obj or edit.network OPTIONAL ARGUMENTS: K Threshold for exceedence. row.name Row name (usually a date) of the measureements matrix (part of nw.obj) to be used for kriging and determining probabilities for a given day. mean.obj Mean. sd.obj Standard Deviation. add.map Overlays US map if True. add.county Overlays US state county boundaries if true. Also maps state boundaries regardless of the value of add.map ... Other parameters for plotting or for prob.cov VALUE: Invisible. SIDE EFFECTS: Creates a plot. SEE ALSO: network.plotter, make.network.obj, edit.network EXAMPLES: net <- make.network.obj( locs = example.locs, cov.obj = example.cov.obj ) plot.prob( net ) |
Plotting Probabilities of Exceedence
Select the fourth DI toolbar button or Plot Probabilities in the DI Menu to begin.
Network Object
Here we will select our NC.net network object created in our
example) from the Network Object.
Date
(or Row Name of Measurements Matrix)
For the date, we chose August 5, 1998 because it had some values in exceedence
of 85 ppb for ozone content.
Exceedence
Choose the value, K, that you wish to plot probabilities of--P(Z > K).
Plots
We can also select which plots we want to see; here select "all" to
see all of them. The diagnostic plot shows the results of the kriging. We can
also select whether we want to have a US map with or without county boundaries.
Other Optional Arguments
There are also some choices for the Fields function,
Krig: whether or not to use GCV criterion or to use a predetermined
number of degrees of freedom. Click on the Advanced tab for more Krig
options.
Because I selected "all" for type of plot, I can click on the different pages (you may have to re-size your plot as I have done in order to click on the different pages) to see different plots. For example, page 2 shows a larger plot of the means found from kriging (as shown below).
This is software for statistical research and should not be used for commercial purposes. The authors do not guarantee the correctness of any function or program in this package.