# NOT RUN {
# Yaw rate, rolling, pitching and rudder angle of a ship
data(HAKUSAN)
y <- as.matrix(HAKUSAN[, c(1,2,4)]) # Yaw rate, Rolling, Rudder angle
z <- marlsq(y)
z
marspc(z$arcoef, v = z$v)
# }
Run the code above in your browser using DataLab