concarea: Plot Concentration Area

Description Usage Arguments Details Value Author(s) References See Also Examples

Description

Displays a concentration-area plot (see also concareaExampleKola). This function is preferable since it can be applied to non-Kola data!

Usage

1
2
3
4
5
concarea(x, y, z, zname = deparse(substitute(z)),
caname = deparse(substitute(z)), borders=NULL, logx = FALSE, ifjit = FALSE,
ifrev = FALSE, ngrid = 100, ncp = 0, xlim = NULL, xcoord = "Easting",
ycoord = "Northing", ifbw = FALSE, x.logfinetick = c(2, 5, 10),
y.logfinetick = c(2, 5, 10))

Arguments

x

name of the x-axis spatial coordinate, the eastings

y

name of the y-axis spatial coordinate, the northings

z

name of the variable to be processed and plotted

zname

a title for the x-axes of the qp-plot and concentration area plot.

caname

a title for the image of interpolated data.

borders

either NULL or character string with the name of the list with list elements x and y for x- and y-coordinates of map borders

logx

if it is required to make a logarithmis data transformation for the interpolation

ifrev

if FALSE the empirical function ist plotted from highest value to lowest

ngrid

default value is 100

xlim

the range for the x-axis

xcoord

a title for the x-axis, defaults to "Easting"

ycoord

a title for the y-axis, defaults to "Northing"

ifbw

if the plot is drawn in black and white

x.logfinetick

how fine are the tick marks on log-scale on x-axis

y.logfinetick

how fine are the tick marks on log-scale on y-axis

ifjit

default value is FALSE

ncp

default value is 0

Details

The function assumes that the area is proportional to the count of grid points. To be a reasonable model the data points should be 'evenly' spread over the plane. The interpolated grid size ist computed as (max(x) - min(x))/ngrid, with a default value of 100 for ngrid. Akima's interpolation function is used to obtain a linear interpolation between the spatial data values.

Value

The concentration area plot, in both directions, is created.

Author(s)

Peter Filzmoser <P.Filzmoser@tuwien.ac.at> http://cstat.tuwien.ac.at/filz/

References

C. Reimann, P. Filzmoser, R.G. Garrett, and R. Dutter: Statistical Data Analysis Explained. Applied Environmental Statistics with R. John Wiley and Sons, Chichester, 2008.

See Also

caplot, concareaExampleKola

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
data(ohorizon)
data(kola.background)
data(bordersKola)

Cu=ohorizon[,"Cu"]
X=ohorizon[,"XCOO"]
Y=ohorizon[,"YCOO"]

par(mfrow=c(1,2),mar=c(4,4,2,2))
concarea(X,Y,Cu,log=TRUE,zname="Cu in O-horizon [mg/kg]",borders="bordersKola", ifrev=FALSE,
         x.logfinetick=c(2,5,10),y.logfinetick=c(10))
	

Example output

Loading required package: geoR
--------------------------------------------------------------
 Analysis of Geostatistical Data
 For an Introduction to geoR go to http://www.leg.ufpr.br/geoR
 geoR version 1.7-5.2.1 (built on 2016-05-02) is now loaded
--------------------------------------------------------------

Loading required package: sgeostat
Warning messages:
1: no DISPLAY variable so Tk is not available 
2: In rgl.init(initValue, onlyNULL) : RGL: unable to open X11 display
3: 'rgl_init' failed, running with rgl.useNULL = TRUE 
4: .onUnload failed in unloadNamespace() for 'rgl', details:
  call: fun(...)
  error: object 'rgl_quit' not found 

StatDA documentation built on March 13, 2020, 2:42 a.m.

Related to concarea in StatDA...