View source: R/convertToEGRET.R
convertToEGRET | R Documentation |
Convert a loadflex object into an EGRET object for plotting.
convertToEGRET(load.model = NULL, newdata = NULL, data = NULL,
meta = NULL)
load.model |
a load model (loadReg2, loadComp, loadInterp, loadLm, etc.) whose data and predictions are to be converted to EGRET format |
newdata |
data.frame of data used to generate predictions from an already-fitted model |
data |
data.frame of data used to fit a model. only required if load.model is omitted |
meta |
loadflex metadata object; it must include constituent, flow, dates, conc.units, site.id, and consti.name. only required if load.model is omitted |
data(lamprey_nitrate)
fitdat <- lamprey_nitrate
data(lamprey_discharge)
estdat <- subset(lamprey_discharge, DATE < as.POSIXct("2012-10-01 00:00:00", tz="EST5EDT"))
estdat <- estdat[seq(1, nrow(estdat), by=96/4),] # only keep 4 observations per day
meta <- metadata(constituent="NO3", flow="DISCHARGE", dates="DATE",
conc.units="mg L^-1", flow.units="cfs", load.units="kg", load.rate.units="kg d^-1",
site.name="Lamprey River, NH", site.id='NWIS 01073500', consti.name = "nitrate")
no3_lm <- loadLm(formula=log(NO3) ~ log(DISCHARGE), pred.format="conc",
data=fitdat, metadata=meta, retrans=exp)
minimalEGRET <- loadflex:::convertToEGRET(data=fitdat, meta=meta)
maximalEGRET <- loadflex:::convertToEGRET(load.model=no3_lm, newdata=estdat)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.