PlotEval | R Documentation |
This function creates an evaluation plot for the Monte Carlo simulation result.
PlotEval(eval, ts, gof1, namePlot, pos1, pos2, pos3)
eval |
A |
ts |
An |
gof1 |
A matrix with the output of |
namePlot |
A character string defining the name of the plot to be created. |
pos1 |
Location to place the legend on the inside of the first sub-plot frame. Can be one of "bottomright", "bottom", "bottomleft", "left", "topleft", "top", "topright", "right" and "center". |
pos2 |
Location to place the legend on the inside of the second sub-plot frame. Can be one of "bottomright", "bottom", "bottomleft", "left", "topleft", "top", "topright", "right" and "center". |
pos3 |
Location to place the legend on the inside of the third sub-plot frame. Can be one of "bottomright", "bottom", "bottomleft", "left", "topleft", "top", "topright", "right" and "center". |
The function creates a plot in the current working directory with the goodness-of-fit between simulations and observations. The plot is provided in pdf format.
J.A. Torres-Matallana
time <- seq(from = as.POSIXct("2017-11-09"), by = 60*60*24, length.out = 230) # the time vector
data <- cbind.data.frame(time, NA) # a NA vector
data[,3] <- rnorm(230, .25, .1) # random normal distributed data, obs
data[,4] <- data[,3]*1.2 # positive correlated data, sim
data[,5] <- data[,4] - data[,3] # difference sim and obs
data[,6] <- 0 # driving force
data[,7] <- NA # a NA vector
colnames(data) <- c("time", "var1", "obs", "sim", "difference", "Rainfall", "Rainfall2")
head(data)
ts <- IsReg.ts(data, "%Y-%m-%d", "ECT")
ts <- ts[[2]]
gof.new <- GoF(data, 4, 3, "")
gof.new
## not run
## creating the plot (uncomment to run)
#PlotEval(data, ts, gof.new, "ExamplePlot", "topright", "topright", "topright")
#
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.