Prepare a Concentration-Number (C-N) Plot
Display a log-log Concentration-Number (C-N) Plot to inspect a data set for multifractality.
gx.mf(xx, xlab = deparse(substitute(xx)),
ylab = "Cumulative Percentage of Data", main = "Multifractality Plot",
ifrev = FALSE, xlim = range(xx, na.rm = TRUE), ...)xx |
name of the variable to be plotted. |
xlab |
by default the character string for |
ylab |
by default the y-axia is labelled ‘Cumulative Percentage of Data’. This replaced with no title, |
main |
by default a plot title of ‘Multifractality Plot’ is provided. This replaced with no title, |
ifrev |
by default the empirical C-N function is plotted from highest value to lowest, |
xlim |
by default this is the range of the data with any NAs removed, |
... |
further arguments to be passed to methods. For example, the size of the axis scale annotation can be changed by setting |
Any multi-fractality in the data, will be revealed as 'lines' of data points with different slopes.
Any less than detection limit values represented by negative values, or zeros or other numeric codes representing blanks in the data, must be removed prior to executing this function, see ltdl.fix.df.
Any NAs in the data vector are removed prior to displaying the plot.
If the default selection for xlim is inappropriate it can be set, e.g., xlim = c(0, 200) or c(2, 200), the latter being appropriate for a logarithmcally scaled plot, i.e. log = TRUE. If the defined limits lie within the observed data range truncated plots will be displayed. If this occurs the number of data points omitted is displayed below the total number of observations in the various panels.
If it is desired to prepare a display of data falling within a defined part of the actual data range, then either a data subset can be prepared externally using the appropriate R syntax, or xx may be defined in the function call as, for example, Cu[Cu < some.value] which would remove the influence of one or more outliers having values greater than some.value. In this case the number of data values displayed will be the number that are <some.value.
In some R installations the generation of multi-panel displays and the use of function eqscplot from package MASS causes warning messages related to graphics parameters to be displayed on the current device. These may be suppressed by entering options(warn = -1) on the R command line, or that line may be included in a ‘first’ function prepared by the user that loads the ‘rgr’ package, etc.
For summary statistics displays to complement the graphics see, gx.summary1, gx.summary2 and inset.
Robert G. Garrett
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.
## 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)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.