genasis (version 1.0)

genstatistic: Computes a set of statistics on input data

Description

Function genstatistic computes a set of statistic according to demands on evaluation of Genasis and GMP data.

Usage

genstatistic(x, y=NA, input="openair", pollutant=NA, by.years=FALSE)

Arguments

x

vector of concentration values or data frame of "genasis"/"openair" type. See 'Details' for more detailed description of both data types.

y

vector of measurement dates in the case of vector input only.

input

type of data frame in the case of data frame input. The allowed values are "openair" (default) and "genasis". In case of vector input, this argument is meaningless.

pollutant

name of the pollutant(s), which will be included into resulting data frame.

by.years

logical. Should the statistics be computed for each year separately?

Value

res

data frame containing 18 (19 in case of by.years=TRUE) as described in paragraph 'Details'.

Details

The genstatistic function provides a set of descriptive statistics, commonly used when evaluating Genasis datasets. The resulting table contains 18 columns (19 in the variant with by.years=TRUE): pollutant name of the compound as taken from x. year the year, for which the summary statistics are computed (only if by.years=TRUE). n number of valid results used for calculations of the summary statistics. mean arithmetic mean of concentration values. sd standard deviation of concentration values. geom. mean geometric mean of concentration values. geom. sd geometric standard deviation. min minimal value. median median of concentration values. max maximal value. Pearson Pearson correlation (product moment) to time, available if at least 3 concentration values are present. Pp p value of Pearson trend test, available if at least 3 concentration values are present. Daniels Daniels trend test, available if at least 3 concentration values are present. Dp p value of Daniels trend test, available if at least 3 concentration values are present. Mann-Kendall Mann-Kendall tau, see MannKendall, available if at least 4 concentration values are present. MKp p-value of Mann-Kendall trend test, for metodic, see MannKendall, available if at least 3 concentration values are present. LS slope linear regression line (least squares fitting) slope. TS slope Theil-Sen linear regression line slope. delta difference between the final and the initial concentration value in the period.

See Also

genloq, genoutlier, genhistogram, genpastoact, genanaggr, genplot, gentransform, genwhisker

Examples

## Vector input.
genstatistic(c(0.123,0.158,0.087,0.252,0.211,0.154),
             c("2012-01-10","2012-08-17","2012-12-12",
               "2013-04-09","2013-08-08","2013-12-10"),
             pollutant=c("pp-DDE"),by.years=TRUE)

## Use of example data from the package:
data(kosetice.pas.genasis)
genstatistic(kosetice.pas.genasis,input="genasis")
data(kosetice.pas.openair)
genstatistic(genpastoact(genoutlier(kosetice.pas.openair[,c(1:4,20:26)],
                         plot=FALSE)$res,method="auto"))