predict.tidem {oce} | R Documentation |
Predict a time series from a tidal model.
## S3 method for class 'tidem' predict(object, newdata, ...)
object |
a |
newdata |
optional vector of POSIXt times at which to make the
prediction. If not present, |
... |
optional arguments passed on to children. |
This is a wrapper around the predict method for
object$model
.
A vector of predictions.
Dan Kelley
tidem
fits a tidal model.
## Not run: library(oce) # 1. tidal anomaly data(sealevelTuktoyaktuk) time <- sealevelTuktoyaktuk[["time"]] elevation <- sealevelTuktoyaktuk[["elevation"]] oce.plot.ts(time, elevation, type='l', ylab="Height [m]", ylim=c(-2,6)) tide <- tidem(sealevelTuktoyaktuk) lines(time, elevation - predict(tide), col="red") abline(h=0, col="red") # 2. prediction at specified times data(sealevel) m <- tidem(sealevel) ## Check fit over 2 days (interpolating to finer timescale) look <- 1:48 time <- sealevel[["time"]] elevation <- sealevel[["elevation"]] oce.plot.ts(time[look], elevation[look]) # 360s = 10 minute timescale t <- seq(from=time[1], to=time[max(look)], by=360) lines(t, predict(m,newdata=t), col='red') legend("topright", col=c("black","red"), legend=c("data","model"),lwd=1) ## End(Not run)