Cumulative Normal Percentage Probability (CPP) Plot
Displays a cumulative normal percentage probability (CPP) plot, equivalent to a Q-Q plot, as has been traditionally used by physical scientists and engineers.
cnpplt(xx, xlab = deparse(substitute(xx)), ylab = "% Cumulative Probability", log = FALSE, xlim = NULL, main = "", ifqs = FALSE, ifshape = FALSE, pch = 3, cex = 0.8, cexp = 1, cex.axis = 0.8, ...)
xx |
name of the variable to be plotted. |
xlab |
by default the character string for |
ylab |
a title for the y-axis, defaults to |
log |
to display the data with logarithmic (x-axis) scaling, set |
xlim |
when used in the |
main |
when used stand-alone a title may be added optionally above the plot by setting |
ifqs |
setting |
ifshape |
when used with function |
pch |
by default the plotting symbol is set to a plus, |
cex |
by default the size of the text for data set size, N, is set to 80%, i.e. |
cexp |
by default the size of the plotting symbol, |
cex.axis |
if overplotting occurs in the y-axis labelling the size of the y-axis labels may be reduced by setting |
... |
further arguments to be passed to methods. For example, the size of the axis titles by setting |
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 logarithmically scaled plot, i.e. log = TRUE
. If the defined limits lie within the observed data range a truncated plot will be displayed. If this occurs the number of data points omitted is displayed below the total number of observations.
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
.
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 NA
s in the data vector are removed prior to displaying the plot.
The available symbols are:
pch: 0 = square, 1 = circle, 2 = triangle, 3 = plus, 4 = X,
5 = diamond, 6 = upside-down triangle, 7 = square with X,
8 = asterisk, 9 = diamond with plus, 10 = circle with plus,
11 = double triangles, 12 = square with plus,
13 = circle with X, 14 = square with upside-down triangle.
Symbols 15 to 18 are solid in the colour specified:
15 = square, 16 = circle, 17 = triangle, 18 = diamond.
Robert G. Garrett
## Make test data available data(kola.o) attach(kola.o) ## A stand-alone cumulative normal percentage probability plot cnpplt(Cu) ## A more appropriately labelled and scaled cumulative normal percentage ## probability plot using a cross/x rather than a plus cnpplt(Cu, xlab = "Cu (mg/kg) in <2 mm O-horizon soil", log = TRUE, pch = 4) ## Detach test data detach(kola.o)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.