Estimate a stochastic model based on the MLE and the Hector implementation.

estimate_hector(
  x,
  n_seasonal = 1,
  model_string,
  likelihood_method = "AmmarGrag",
  cleanup = TRUE
)

Arguments

x

A gnssts object

n_seasonal

An integer specifying the number of seasonal component in the time series.

model_string

A string specifying the model to be estimated.

likelihood_method

A string taking either value "FullCov" or "AmmarGrag" that specify the method for the Likelihood computation.

cleanup

A boolean specifying if the files created by the estimation procedure should be cleaned.

Value

A gnsstsmodel object.

Examples

if (FALSE) {
cola = PBO_get_station(station_name = "COLA", column = "dE", time_range = c(51130, 52000))
fit_mle = estimate_hector(x = cola,
                          n_seasonal = 1, 
                          model_string = "wn+matern")

}