# pocket.plot: graphs the probability or standardized variance in the... In geospt: Geostatistical Analysis and Design of Optimal Spatial Sampling Networks

## Description

The pocket-plot (so named because of its use in detecting pockets of non-stationarity) is a technique necessary to identify a localized area that is atypical with respect to the stationarity model. It is built to exploit the spatial nature of the data through the coordinates of rows and columns (east "X" and north "Y", respectively).

## Usage

 `1` ```pocket.plot(data, graph, X, Y, Z, Iden, ...) ```

## Arguments

 `data` data frame should contain the dependent variable and coordinates X and Y, data must be gridded `graph` type of graph associated with the probability or standardized variance plot pocket in the directions north-south or east-west; Probabilities PocketPlot by rows, ie horizontal "south-north" (PPR), Probabilities PocketPlot by columns, ie vertical "east-west" (PPC), PocketPlot of variance by rows, ie horizontal "south-north" (PVR) and PocketPlot of variance by columns, ie vertical "east-west" (PVC) `X` defined by the spatial coordinates `Y` defined by the spatial coordinates `Z` regionalized variable with which you construct the statistics associated with the probability or standardized variance, these are plotted in the so-called pocket plot `Iden` logical. The users can identify the points by themselves, TRUE or FALSE `...` arguments to be passed to ...

## Details

For identifying outliers, this function uses a modification of the boxplot.with.outlier.label function, available at http://www.r-statistics.com/wp-content/uploads/2011/01/boxplot-with-outlier-label-r.txt

## Value

returns (or plots) the pocket plot

## References

Cressie, N.A.C. 1993. Statistics for Spatial Data. Wiley.

Gomez, M., Hazen, K. 1970. Evaluating sulfur and ash distribution in coal seems by statistical response surface regression analysis. U.S. Bureau of Mines Report RI 7377.

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36``` ```# Core measurements (in % coal ash) at reoriented locations. # Units on the vertical axis are % coal ash. # These data was found in mining samples originally reported by # Gomez and Hazen (1970), and later used by Cressie (1993). # These data are available in the sp and gstat packages library(gstat) data(coalash) plot(coalash[,1:2], type="n", xlab="x", ylab="y") text(coalash\$x,coalash\$y,coalash\$coalash,cex=0.6) # Pocket plot in the north-south direction. # Units on the vertical axis are root (% coal ash) # Plot generated with the function pocket.plot # Clearly rows 2, 6, and 8 are atypical # This serves as verification that these rows are potentially problematic # Analysis of local stationarity in probabilities of the coal in south-north direction pocket.plot(coalash, "PPR", coalash\$x, coalash\$y, coalash\$coalash, FALSE) # Analysis of local stationarity in variance of the coal in south-north direction pocket.plot(coalash, "PVR", coalash\$x, coalash\$y, coalash\$coalash, FALSE) # Analysis of local stationarity in probabilities of the coal in east-west direction pocket.plot(coalash, "PPC", coalash\$x, coalash\$y, coalash\$coalash, FALSE) # Analysis of local stationarity in variance of the coal in east-west direction pocket.plot(coalash, "PVC", coalash\$x, coalash\$y, coalash\$coalash, FALSE) ```

### Example output

```Loading required package: gstat