plotDurationUncertainty: plotDurationUncertainty

Description Usage Arguments Details Value Author(s) Examples

View source: R/duration.R

Description

plot the increase of uncertainty with decreaseing measurement duration

Usage

1
2
3
4
5
plotDurationUncertainty(ds, colTime = "TIMESTAMP", 
    fRegress = c(exp = regressFluxExp, lin = regressFluxLinear, 
        tanh = regressFluxTanh), ..., durations = seq(max(65, 
        resFit0$tLag), max(times0), length.out = nDur + 
        1), nDur = 20, maxSdFlux = 1)

Arguments

ds
colTime

column name of time [s]

fRegress

list of functions to yield << a single flux estimate, see details of calcClosedChamberFlux

...

further arguments to calcClosedChamberFlux

durations

durations to check. Default is equally spaced between tLag and maximum duration

nDur

number of durations to check

maxSdFlux

maxium allowed standard deviation of flux in [mumol / s]

Details

Produces a plot with standard deviation of the flux estimate versus the duration of the measurment. The lines correspond to the given maxium acceptable standard deviation and the duration that matches this criterion.

Value

tibble result of calcClosedChamberFlux for the minimum duration, with additional component

statAll

tibble: each row a fit for a given duration

Author(s)

Thomas Wutzler, Oscar Perez Priego

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
#data(chamberLoggerEx2)
ds <- subset(chamberLoggerEx2, iChunk==99)	# very strong (and therefore precise) uptake
#plot( CO2_dry ~ TIMESTAMP, ds )
resDur <- plotDurationUncertainty( ds, colTemp="AirTemp", volume = 0.6*0.6*0.6
	, maxSdFlux = 0.8
	, nDur=10
	, durations = c(100,120,150)
)
resDur$duration
#plot( flux ~ duration, resDur$statAll[[1]] )
#plot( sdFlux ~ duration, resDur$statAll[[1]] )

RespChamberProc documentation built on May 2, 2019, 5:53 p.m.