rgr (version 1.1.15)

gx.mf: Prepare a Concentration-Number (C-N) Plot

Description

Display a log-log Concentration-Number (C-N) Plot to inspect a data set for multifractality.

Usage

gx.mf(xx, xlab = deparse(substitute(xx)), 
        ylab = "Cumulative Percentage of Data", main = "Multifractality Plot",
        ifrev = FALSE, xlim = range(xx, na.rm = TRUE), ...)

Arguments

xx

name of the variable to be plotted.

xlab

by default the character string for xx is used for the x-axis plot titles. An alternate title can be displayed with xlab = "text string", see Examples.

ylab

by default the y-axia is labelled ‘Cumulative Percentage of Data’. This replaced with no title, main = " ", or text of the user's choice, main = "text string".

main

by default a plot title of ‘Multifractality Plot’ is provided. This replaced with no title, main = " ", or text of the user's choice, main = "text string".

ifrev

by default the empirical C-N function is plotted from highest value to lowest, ifrev = FALSE. As the C-N plot is a log-log display this provides greater detail for the highest values. The direction of accumulation can be key in detecting multifractal patterns, it is usually informative to also prepare a plot with ifrev = TRUE, i.e. accumulation from lowest to highest values. To see a dramatic example of this, run the Examples below.

xlim

by default this is the range of the data with any NAs removed, xlim = range (xx, ra.rm = TRUE). When called from caplot or shape.alt an appropriate value for xlim is determined by the calling functions to ensure all four panels have the same x-axis scaling. xlim may be defined, see Note below.

…

further arguments to be passed to methods. For example, the size of the axis scale annotation can be changed by setting cex.axis, the size of the axis titles by setting cex.lab, and the size of the plot title by setting cex.main. For example, if it is required to make the plot title smaller, add cex.main = 0.9 to reduce the font size by 10%. By default individual data points are marked by a plus sign, pch = 3, if a cross or open circle is desired, then set pch = 4 or pch = 1, respectively. See display.marks for all available symbols.

Details

Any multi-fractality in the data, will be revealed as 'lines' of data points with different slopes.

References

Agterberg, F.P., 1995. Multifractal Modeling of the Sizes and Grades of Giant and Supergiant Deposits. International Geology Review, Vol. 37, No. 1, pp. 1-8.

See Also

gx.hist, cnpplt, gx.ecdf, gx.mf, remove.na, display.lty, display.marks, ltdl.fix.df, inset

Examples

# NOT RUN {
## Make test data available
data(kola.o)
attach(kola.o)

## Generates an initial display to have a first look at the data and 
## decide how best to proceed
gx.mf(Cu)

## Provides a more appropriate initial display and indicates the 
## quartiles
gx.mf(Cu, xlab = "Cu (mg/kg) in <2 mm O-horizon soil", ifqs = TRUE)

## Causes the C-N plot to be cumulated in reverse order.  This will reveal
## any multi-fractal properties of the data at lower concentrations
gx.mf(Cu, xlab = "Cu (mg/kg) in <2 mm O-horizon soil", ifrev = TRUE)

## Detach test data
detach(kola.o)
# }