genpastoact
provides two different tools for recalculating POPs concentration results from passive air samplers (in mass per filter units) to concentrations (in mass per air volume) comparable with results from active air samples.
genpastoact(x, date_start=NA, date_end=NA, input="openair", output=NA, method="auto", days=28, pollutant=NA, temp=NA, wind=NA)
"character"
, "Date"
, "POSIXct"
or "POSIXlt"
, as described in gendate. Only meaningful for vector input."character"
, "Date"
, "POSIXct"
or "POSIXlt"
, as described in gendate. Only meaningful for vector input.input
argument, both data frames "openair" and "genasis" are available, with the default value equal to input
."const"
, "temp"
, "tempwind"
, "auto"
, "gaps4"
and "gaps7"
. For detailed description see 'Details'.date_start
and date_end
aguments. If date_start
is specified and date_end
is not, the date of end of the measurement is computed and add.pollutant
argument in a case of numeric input.output
argument, containing recalculated values of POPs air concentrations (in mass unit per cubic meter). genpastoact
function enables to use two different models for recalculating values of POPs concentration (in mass per filter units) from passive air samplers (of Genasis design) to real air concentrations (in mass per air volume units) as obtained from active air samplers. The function recognises three different input formats: Option input="openair"
uses "openair" format of data frame with first column of name 'date' and type "Date"
, optional columns of names "date_end"
, "temp"
, "wind"
and "note"
and other columns of type "numeric"
containing concentration values and named by names of the compounds. input="genasis"
is used for the data frame with six columns "valu"
, "comp"
, "date_start"
, "date_end"
, "temp"
and "wind"
where the first, fifth and sixth are of type "numeric"
, second of type "character"
and third and fourth columns could be both "character"
or "Date"
type. The names of columns in input="genasis"
are not rigid, only their order is assumed. There is also a possibility to specify x
and y
as two vectors of equal lenght, first of type "numeric"
containing concentration values, second of type "character"
or "Date"
containing measurement dates.
There are six possibilities of model choice. Three similar variants method="const"
, method="temp"
and method="tempwind"
leads to the Genasis model taking into account gradually only the concentration, concentration and temperature and finally concentration, temperature and wind speed. If temperature and/or wind are necessary for the selected model, they need to be specified in appropriate column of data frame or as a temp
and/or wind
arguments in case of the vector input. The option method="auto"
combines above mentioned three approaches, using as many variables, as possible (in such a cases NA represents not-known temperature/wind.)
Remaining two options method="gaps4"
and method="gaps7"
lead to use of Global Air Passive Sampling (GAPS) model with two different effective gas-phase sampling rates (4 for usual and 7 for significantly windy sites). In all cases this model uses temperature (either as temp
argument in the vector input mode, or as column temp
in the data frame input mode) and does not use wind speed for the recalculation.
The argument pollutant
could take any value from the following set: "naphthalene", "acenaphtylene", "acenaphtene", "fluorene", "phenantrene", "anthracene", "fluoranthene", "pyrene", "benzo(a)anthracene", "chrysene", "benzo(b)fluoranthene", "benzo(k)fluoranthene", "benzo(a)pyrene", "indeno(123cd)pyrene", "benzo(ghi)perylene", "PCB 28", "PCB 52", "PCB 101", "PCB 118", "PCB 138", "PCB 153", "PCB 180", "alpha-HCH", "beta-HCH", "gamma-HCH", "pp-DDE", "pp-DDD", "pp-DDT", "PeCB", "HCB"
. For other pollutants no model constant are known and thus the concentration can not be racalculated.
genloq, genoutlier, genhistogram, genanaggr,
genplot, genstatistic, gentransform, genwhisker
## Vector input
genpastoact(c(1,2,1.2,1.1,0.1),days=28,pollutant=c("pp-DDE"))
## Use of example data from the package:
data(kosetice.pas.genasis)
genpastoact(kosetice.pas.genasis,input="genasis",
method="gaps7",pollutant=c("PCB 28"))
data(kosetice.pas.openair)
genpastoact(genoutlier(kosetice.pas.openair[,c(1:4,20)],plot=FALSE)$res,
method="auto")
Run the code above in your browser using DataLab