load_french_creek_std_mle: Generate outputs using Bob's code for comparison

View source: R/load_french_creek_std.R

load_french_creek_std_mleR Documentation

Generate outputs using Bob's code for comparison

Description

Bob's code includes MLE and nighttime regression models. This function generates the output from those models, keeping the code as much intact as possible. The exception is that we're using solar.time rather than local.time, for consistency with streamMetabolizer's recommendations

Usage

load_french_creek_std_mle(
  french,
  K = 35,
  estimate = c("PRK", "K", "PR"),
  start = c(dates = "08/23/12", times = "22:00:00"),
  end = c(dates = "08/25/12", times = "06:00:00"),
  plot = FALSE
)

Arguments

french

the French Creek dataset

K

optional. If specified, a number for the K600 to assume (units of 1/d)

estimate

character indicating the type of model to fit

start

a character vector specifying the time at which the 'day' (the time period to use in producing an estimate for a single date) begins. The vector should have 2 elements, dates and times, to pass to chron()

end

a character vector specifying the time at which the 'day' ends. The vector should have 2 elements, dates and times, to pass to chron()

plot

logical - should plots be produced?

Details

This function requires the chron package, which is only suggested rather than required for the streamMetabolizer package. If you wish to run this function, ensure that chron is installed or install it with install.packages('chron').


USGS-R/streamMetabolizer documentation built on Aug. 15, 2023, 7:50 a.m.