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':
predict(object, newdata, ...)

Arguments

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.

Value

A vector of predictions.

Details

This is a wrapper around the predict method for object$model.

Examples

Run this code

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)

Run the code above in your browser using DataLab

Last chance! 50% off unlimited learning