uGMAR (version 3.2.6)

parameterChecks: Check the parameter vector is specified correctly

Description

parameterChecks checks dimension, restrictions, and stationarity of the given parameter of a GMAR, StMAR, or G-StMAR model. Throws an error if any check fails. Does NOT consider the identifiability condition!

Usage

parameterChecks(
  p,
  M,
  params,
  model = c("GMAR", "StMAR", "G-StMAR"),
  restricted = FALSE,
  constraints = NULL
)

Arguments

p

a positive integer specifying the autoregressive order of the model.

M

For GMAR and StMAR models:

a positive integer specifying the number of mixture components.

For G-StMAR models:

a size (2x1) integer vector specifying the number of GMAR type components M1 in the first element and StMAR type components M2 in the second element. The total number of mixture components is M=M1+M2.

params

a real valued parameter vector specifying the model.

For GMAR model:

Size $(M(p+3)-1x1)$ vector $\theta$$=$(${\upsilon }_1$,...,${\upsilon }_M$, ${\alpha }_1,...,{\alpha }_{M-1}$), where ${\upsilon }_m$$=({\varphi }_{m,0},$${\varphi }_m$$,{\sigma }_m^2)$ and ${\varphi }_m$=$({\varphi }_{m,1},...,{\varphi }_{m,p}),m=1,...,M$.

For StMAR model:

Size $(M(p+4)-1x1)$ vector ($\theta ,\nu$)$=$(${\upsilon }_1$,...,${\upsilon }_M$, ${\alpha }_1,...,{\alpha }_{M-1},{\nu }_1,...,{\nu }_M$).

For G-StMAR model:

Size $(M(p+3)+M2-1x1)$ vector ($\theta ,\nu$)$=$(${\upsilon }_1$,...,${\upsilon }_M$, ${\alpha }_1,...,{\alpha }_{M-1},{\nu }_{M1+1},...,{\nu }_M$).

model

is "GMAR", "StMAR", or "G-StMAR" model considered? In the G-StMAR model, the first M1 components are GMAR type and the rest M2 components are StMAR type.

restricted

a logical argument stating whether the AR coefficients ${\varphi }_{m,1},...,{\varphi }_{m,p}$ are restricted to be the same for all regimes.

constraints

specifies linear constraints applied to the autoregressive parameters.

For non-restricted models:

a list of size $(px{q}_m)$ constraint matrices $C_m$ of full column rank satisfying ${\varphi }_m$$=$$C_m{\psi }_m$ for all $m=1,...,M$, where ${\varphi }_m$$=({\varphi }_{m,1},...,{\varphi }_{m,p})$ and ${\psi }_m$$=({\psi }_{m,1},...,{\psi }_{m,q_m})$.

For restricted models:

a size $(pxq)$ constraint matrix $C$ of full column rank satisfying $\varphi$$=$$C\psi$, where $\varphi$$=({\varphi }_1,...,{\varphi }_p)$ and $\psi$$={\psi }_1,...,{\psi }_q$.

Symbol $\varphi$ denotes an AR coefficient. Note that regardless of any constraints, the nominal autoregressive order is always p for all regimes. Ignore or set to NULL if applying linear constraints is not desired.

Value

Throws an informative error if any check fails. Does not return anything.