Learn R Programming
oce (version 0.9-13)

predict.tidem: Predict a time series from a tidal model

Description

Predict a time series from a tidal model.

Usage

## S3 method for class 'tidem}(object, newdata, ...)':
predictundefined
object{a tidem object.}
  newdata{optional vector of POSIXt times at which to make the
    prediction.  If not present, predict.tidem uses the times
    that were provided in the original call to tidem.}
  ...{optional arguments passed on to children.}
This is a wrapper around the predict method for
  object$model.
A vector of predictions.
tidem fits a tidal model.
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(sealevelHalifax)
m <- tidem(sealevelHalifax)
## Check the fit over 2 days (interpolating to a finer timescale)
look <- 1:48
time <- sealevelHalifax[["time"]]
elevation <- sealevelHalifax[["elevation"]]
oce.plot.ts(time[look], elevation[look])
t <- seq(from=time[1], to=time[max(look)], by=360) # 360s = 10 minute timescale
lines(t, predict(m,newdata=t), col='red')
legend("topright", col=c("black","red"),legend=c("data","model"),lwd=1)
[object Object]
misc

Arguments