Compare graphically two gnsstsmodel objects.

Compare graphically two gnsstsmodel objects.

compare_fits(fit_1, fit_2, main = NULL, y_unit = "mm", x_unit = "days")

Arguments

fit_1

A gnsstsmodel object.

fit_2

A gnsstsmodel object.

main

A string specifying the plot title.

y_unit

A string specifying the y axis label.

x_unit

A string specifying the x axis label.

Value

No return value. Produce a plot comparing two estimated models.

Examples

if (FALSE) {
data(cola)
fit_gmwmx_1 = estimate_gmwmx(x = cola,
                             theta_0 = c(0.1,0.1,0.1,0.1), 
                             n_seasonal = 1, 
                             model_string = "wn+matern")
fit_gmwmx_2 = estimate_gmwmx(x = cola,
                             theta_0 = c(0.1,0.1,0.1), 
                             n_seasonal = 1, 
                             model_string = "wn+powerlaw")
compare_fits(fit_gmwmx_1, fit_gmwmx_2)
}