trendLevel(mydata, pollutant = "nox", x = "month",
y = "hour", type = "year", rotate.axis = c(90, 0),
n.levels = c(10, 10, 4), limits = c(0, 100),
cols = "default", auto.text = TRUE,
key.header = "use.stat.name", key.footer = pollutant,
key.position = "right", key = TRUE,
statistic = c("mean", "max", "frequency"),
stat.args = NULL, stat.safe.mode = TRUE,
drop.unused.types = TRUE, ...)trendLevel plot.mydata to sample to produce the trendLevel
plot.trendLevel x-axis. This is used with the y
and type options to bin the data before applying
statistic (see below). Allowed options currently
include trendLevel y-axis and for additional conditioning,
respectively. As x above.trendLevel x and y axes. The
default, c(90, 0), rotates the x axis by 90
degrees but does not rotate the y axis. (Note: If only
one value is supplied, this is apx,
y and type data into if numeric. The
default, c(10, 10, 4), cuts numeric x and
y data into ten levels and numeric type
trendLevel plot.trendLevel surface. cols is passed to
openColours for evaluation. See
?openColours for more details.auto.text = TRUE passes common lattice
labelling terms (e.g. xlab for the x-axis,
ylab for the y-axis and main for the title)
to the plot via <key.header = "",
key.footer = c("mean","nox") adds the addition text as a
scale footer. If enabled (auto.text = TRU"top", "right", "bottom" and
"left".drawOpenKey. See ?drawOpenKey for further
details.pollutant measurements
with. Three options are currently encoded: "mean"
(default), "max" and "frequency". (Note:
Functions can alsostatistic if this is a function. The extra options
should be supplied as a list of named parameters. (see
Details below.)statistic that most
users can ignore. This option returns NA instead
of running statistic on binned subsamples that are
empty. Many commontype
conditioning cases. Some conditioning options may
generate empty cases for some data sets, e.g. a hour of
the day when no measurements were taken. Empty x
and y cases generate 'holcutData for type handling and
levelplot in lattice for finer control of
the plot itself.trendLevel
also returns an object of class ``openair''. The object
includes three main components: call, the command
used to generate the plot; data, the data frame of
summarised information used to make the plot; and
plot, the plot itself. If retained, e.g. using
output <- trendLevel(mydata), this output can be
used to recover the data, reproduce or rework the
original plot or undertake further analysis.
An openair output can be manipulated using a number of
generic operations, including print, plot
and summary. See openair.generics
for further details.
Summary statistics can also be extracted directly using
results, e.g. results(object) for
output <- trendLevel(mydata).trendLevel allows the use of third party
summarising functions via the statistic option.
Any additional function arguments not included within a
function called using statistic should be supplied
as a list of named parameters and sent using
stat.args. For example, the encoded option
statistic = "mean" is equivalent to
statistic = mean, stat.args = list(na.rm = TRUE)
or the R command mean(x, na.rm= TRUE). Many R
functions and user`s own code could be applied in a
similar fashion, subject to the following restrictions:
the first argument sent to the function must be the data
series to be analysed; the name `x' cannot be used for
any of the extra options supplied in stat.args;
and the function should return the required answer as a
numeric or NA. Note: If the supplied function returns
more than one answer, currently only the first of these
is retained and used by trendLevel. All other
returned information will be ignored without warning. If
the function terminates with an error when it is sent an
empty data series, the option stat.safe.mode
should not be set to FALSE or trendLevel
may fail. Note: The stat.safe.mode = TRUE option
returns an NA without warning for empty data series.openColours and drawOpenKey
for more detailed plot control and
openair.generics for output handling.#basic use
#default statistic = "mean"
trendLevel(mydata, pollutant = "nox")
#applying same as 'own' statistic
my.mean <- function(x) mean(x, na.rm = TRUE)
trendLevel(mydata, pollutant = "nox", statistic = my.mean)
#alternative for 'third party' statistic
#trendLevel(mydata, pollutant = "nox", statistic = mean,
# stat.args = list(na.rm = TRUE))Run the code above in your browser using DataLab