Setter or modifier methods for objects of class
1 2 3 4 5 6 7 8 9 10
## S4 method for signature 'stsm' set.cpar(x, value, check = TRUE, inplace = FALSE) ## S4 method for signature 'stsm' set.nopars(x, v, check = TRUE, inplace = FALSE) ## S4 method for signature 'stsm' set.pars(x, v, check = TRUE, inplace = FALSE) ## S4 method for signature 'stsm' set.sgfc(x, inplace = FALSE) ## S4 method for signature 'stsm' set.xreg(x, xreg, coefs = NULL)
an object of class
a numeric value.
a numeric vector.
a matrix or numeric vector of external regressors. The number of
rows or length of the vector must be equal to the length of
an optional vector containing the value of the coefficients
related to the regressors
Models parameterized with non-null
If the model is parameterized according to a non-null value
of the slot
transPars, the argument
must contain the values of the auxiliary set of parameters θ
rather than the actual parameters (variances and autoregressive coefficients).
For example, with
x@transPars = "square" the variances are θ^2.
Although this design may seem to disagree with the
the relevant input for the setter methods is actually the auxiliary values θ.
Be aware that if
transPars is not null the parameters are
get.pars according to the selected parameterization.
v must be referred to the non-transformed parameters.
The previous comment does not apply to the argument
is not affected by
Setter methods are safer.
For those users that are not familiar with the design and internal structure of the
stsm, it is safer to use setter methods
rather than modifying the contents of the slots through the
See the examples below.
Modifying the input object in-place.
Instead of returning the whole object and create a new one or overwrite
the original, it is possible to modify just the desired slot in the original
object that is passed as input. In the former case the
stsm object returned by
the method must be assigned to another object using the usual operator
In the latter approach, the
stsm object that is passed as argument is modified
in-place. See the example below. The solution to modify an object in-place is taken from
This option is not a customary solution in R, however, it seems suitable in this context.
The real benefit of this approach would depend on how R deals with objects that
are returned from functions. If assigning the output to a new object involves making
copies of all the slots, then modifying the object in-place would most likely be more efficient
since the desired slot is directly modified avoiding copying the whole object.
After R version 3.1 this issue may become less critical. One of the new features reported in the release of R 3.1 states: Avoid duplicating the right hand side values in complex assignments when possible. This reduces copying of replacement values in expressions such as Z$a <- a0. A related discussion for S4 classes can be found in this post.
Constant terms in the spectral generating function.
In pure variance models, some elements of the spectral generating function (s.g.f.)
do not depend on the parameters and can be stored as constants.
set.sgfc computes and stores those elements as a matrix
in the slot
sgfc. This is useful for example when working with maximum likelihood methods
in the frequency domain. In that context, the spectral generating function
has to be updated several times for different parameter values. Having the information about
the constant terms in the slot
sgfc saves several computations whenever
the s.g.f. is requested. For details about the s.g.f see
Further setter methods.
Future versions may include additional setter methods, for example to
change the slot
model or to modify the time series
The latter would also require updating the slots
if requested. Additional methods are not available in the current version because
defining a new object by means of
often be better than modifying one of those slots that do not
have a setter method.
If the slot is modified in place,
nothing is returned, the corresponding slot of the object
passed as argument is modified in place.
inplace=FALSE, a new
stsm object is returned.
It contains the same information as the input object
except for the slot that has been modified.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86
# sample models with arbitrary parameter values m <- stsm.model(model = "llm+seas", y = JohnsonJohnson, pars = c("var1" = 2, "var2" = 15, "var3" = 30)) get.pars(m) # correct modification m1 <- set.pars(m, c(1, 2, 3)) get.pars(m1) m1 <- set.pars(m, c(var1 = 11)) get.pars(m1) # correct but error prone m1@pars <- c(4, 22, 33) get.pars(m1) m1@pars <- c(var1 = 1, var2 = 2, var3 = 3) get.pars(m1) # inconsistent assignment (error returned) # 'var4' is not a parameter of model 'llm+seas' try(m1 <- set.pars(m, c(var4 = 4))) # inconsistent assignment (no error returned) # the error is not noticed at this point # unless 'validObject' is called m1 <- m m1@pars["var4"] <- 4 get.pars(m1) try(validObject(m1)) # modify only one element m1 <- set.pars(m, v=c(var1=22)) get.pars(m1) # wrong assignment, the whole vector in the slot is overwritten # no error returned at the time of doing the assignment m1@pars <- c(var1 = 1) get.pars(m1) try(validObject(m1)) # consistent assignment but maybe not really intended # all the elements are set equal to 12 m1 <- m m1@pars <- 12 get.pars(m1) # warning returned by 'set.pars' m2 <- set.pars(m, 12) get.pars(m2) # wrong value unnoticed (negative variance) m1 <- m m1@pars <- c(-11, 22, 33) get.pars(m1) # negative sign detected by 'set.pars' try(m1 <- set.pars(m, c(-11, 22, 33))) # inplace = FALSE # the whole object 'm' is assigned to a new object, # which will probably involve making a copy of all the slots m <- set.pars(m, c(1,2,3), inplace = FALSE) get.pars(m) # inplace = TRUE # the output is not assigned to a new object # the only operation is the modification of the slot 'pars' # no apparent additional internal operations such as copying unmodified slots get.pars(m) set.pars(m, c(11,22,33), inplace = TRUE) get.pars(m) # set a matrix of regressors xreg <- cbind(xreg1 = seq_len(84), xreg2 = c(rep(0, 40), rep(1, 44))) m <- stsm.model(model = "llm+seas", y = JohnsonJohnson, xreg = xreg) m # set a new matrix of regressors to an existing xreg3 <- seq(length(m@y)) m2 <- set.xreg(m, xreg3) m2 # remove the external regressors m3 <- set.xreg(m, NULL) m3 m3@xreg # initialize the coefficients to some values m <- stsm.model(model = "llm+seas", y = JohnsonJohnson, pars = c("xreg1" = 10), xreg = xreg) m m <- stsm.model(model = "llm+seas", y = JohnsonJohnson, pars = c("xreg2" = 20, "xreg1" = 10), xreg = xreg) m