Last chance! 50% off unlimited learning
Sale ends in
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, ...)
xlab ="Cu (mg/kg) in <2 mm="" o-horizon="" soil"<="" code="">.
log = TRUE
.main
, e.g., main = "Kola Ecogeochemistry Project, 1995"
.ifqs = TRUE
results in horizontal and vertical dotted lines being plotted at the three central quartiles and their values, respectively.pch = 3
, alternate plotting symbols may be chosen from those displayed by display.marks
.cex = 0.8
, and may be changed if required.pch
, is set to 100%, and may be changed if required.cex.axis
to a number smaller than the default of cex.axis = 0.8
.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 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 .
display.marks
, ltdl.fix.df
, remove.na
## 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)
Run the code above in your browser using DataLab