biogas (version 1.0.1)

interp: Interpolate Biogas Composition or Cumulative Production

Description

interp interpolates (or extrapolates) biogas composition (methane concentration) or cumulative production data to a wanted time using one of several possible methods.

Usage

interp(times, y, time.out, method = "linear", extrap = FALSE)

Arguments

times

measurement times. Numeric vector or POSIX.

y

response variable at times, e.g., mole fraction of methane in biogas, or biogas cumulative production. Numeric vector.

time.out

time or times at which interpolated values are needed. Numeric vector or POSIX.

method

method used for interpolation. Default is "linear", which uses the approx function. For spline interpolation, use "ffm" or the other options available (see

extrap

should y be exptrapolated? Logical. Default is FALSE. See 'Details'.

Value

• Interpolated estimates of y at given times.

concept

biogas

Details

interp is really a wrapper for the interpolation functions approx and spline. For cumulative production, which (usually) must monotonically increase, method="hyman" is the best choice. Extrapolation behavior depends on method. For method="linear", extrap=TRUE simply returns the value of the closest y. See approx and spline for more information.

See Also

cumBg

Examples

# Fake composition data 
dat <- data.frame(time=c(1,7,14,28),xCH4=c(0.3,0.5,0.61,0.65))
interp(dat$time,dat$xCH4,time.out=10)
interp(dat$time,dat$xCH4,time.out=10,method="natural")
interp(dat$time,dat$xCH4,time.out=c(10,30))
interp(dat$time,dat$xCH4,time.out=c(10,30),method="natural")
interp(dat$time,dat$xCH4,time.out=c(10,30),extrap=TRUE)

# Actual data
data(comp)
# Work with one reactor
bgc <- subset(comp,id=="2_1")
# With numeric time, interpolate to 1, 7, and 30 days
interp(bgc$days,bgc$xCH4,time.out=c(1,7,30))
# Or POSIXct
interp(bgc$date.time,bgc$xCH4,time.out=as.POSIXct("2014-07-12 13:00:00"))

# For cumulative gas production
data(vol)
# Work with one reactor
bgv <- subset(vol,id=="2_1")
# Calculate cumulative volume
bgv <- cumBg(bgv,time.name="days")
# The interpolate of cumulative production to 1, 7, and 30 days
interp(bgv$days,bgv$cvBg,time.out=c(1,7,30),method="hyman")