R: Create a gnssts object

create.gnssts

R Documentation

Create a gnssts object

Description

Create a gnssts object

Usage

create.gnssts(t, y, jumps = NULL, sampling_period = 1)

Arguments

`t`	A `vector` specifying the time of each observation of the time series.
`y`	A `vector` specifying the values of each observation of the time series.
`jumps`	A `vector` specifying the time values for which there is a jump.
`sampling_period`	An `integer` specifying the sampling period.

Value

A gnssts object.

Examples

phase <- 0.45
amplitude <- 2.5
sigma2_wn <- 15
bias <- 0
trend <- 5 / 365.25
cosU <- amplitude * cos(phase)
sinU <- amplitude * sin(phase)
year <- 5
n <- year * 365
jump_vec <- c(200, 300, 500)
jump_height <- c(10, 15, 20)
nbr_sin <- 1 
A <- create_A_matrix(1:n, jump_vec, n_seasonal = nbr_sin)
x_0 <- c(bias, trend, cosU, sinU, jump_height)
eps <- rnorm(n = n, sd = sqrt(sigma2_wn))
yy <- A %*% x_0 + eps
gnssts_obj <- create.gnssts(t = 1:length(yy), y = yy, jumps = jump_vec)
str(gnssts_obj)