inst/doc/estimate_geodetic_time_series_from_ngl.R

In gmwmx2: Estimate Functional and Stochastic Parameters of Linear Models with Correlated Residuals and Missing Data

## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(
  fig.width = 10, # Set default plot width (adjust as needed)
  fig.height = 8, # Set default plot height (adjust as needed)
  fig.align = "center" # Center align all plots
)

knitr::opts_chunk$set(eval = FALSE)

## -----------------------------------------------------------------------------
# library(gmwmx2)

## -----------------------------------------------------------------------------
# station_data <- download_station_ngl("1LSU")

## -----------------------------------------------------------------------------
# plot(station_data)

## -----------------------------------------------------------------------------
# plot(station_data, component = "N")

## -----------------------------------------------------------------------------
# fit1 <- gmwmx2(station_data, n_seasonal = 2, component = "N", model = wn())
# fit1

## -----------------------------------------------------------------------------
# plot(fit1)

## -----------------------------------------------------------------------------
# fit2 <- gmwmx2(station_data, n_seasonal = 2, component = "N", model = wn()+ pl())
# fit2

## -----------------------------------------------------------------------------
# plot(fit2)

## -----------------------------------------------------------------------------
# fit3 = gmwmx2(station_data, n_seasonal = 2, component = "N", model = wn() + flicker())
# fit3

## -----------------------------------------------------------------------------
# plot(fit3)

## -----------------------------------------------------------------------------
# fit4 = gmwmx2(station_data, n_seasonal = 2, component = "N", model = wn() + ar1())
# fit4

## -----------------------------------------------------------------------------
# plot(fit4)