SMC Procedure
STATE Statement
STATE state-variable ~ distribution </ options> ;
The STATE statement specifies the transition distributions of the state variable given the parameters, the lagged state variables, the lagged and current independent variables; that is, ,
. For more information about the transition equation in the state space model, see the section Sequential Monte Carlo Methods and State Space Models. The STATE statement is ignored if you specify the SUBMITMODEL statement. The STATE statement has the following requirements if you omit the SUBMITMODEL statement:
The STATE statement is required. Each state variable must have a corresponding STATE statement.
Multiple STATE statements are permitted. When you specify multiple STATE statements, be careful about the order if there is any dependency between state variables.
The STATE statement must have the following components in the order shown:
State variable: A single state variable that is specified in the STATEVAR statement.
Tilde: A tilde (~) symbol between the state variable and the distribution.
-
Distribution: A distribution that is specified in the form distribution-name(distribution-arguments). The available distribution-name and distribution-arguments in the STATE statements are displayed in Table 6. For the description of each distribution, see the section Standard Distributions.
Table 6: Distribution Names and Arguments
The distribution-arguments are specified by expressions that are enclosed in parentheses and separated by commas. Such expressions must follow any SAS programming statements. To avoid the possible computation errors, you should make some adjustments to the expressions. For more information about the adjustments, see the section Distribution Argument Restriction. The number of expressions is equal to the number of arguments of the given distribution. For the STATE statement, the variables that can appear in the expressions must be selected from the following list:
the model parameters
the lagged state variables
the lagged dependent variables
the lagged and current independent variables
the state variables that have been specified in the preceding STATE statements
the actual time index,
_time
For more information about the actual time index, see the section Time Index. For more information about the parameters, the state variables, and the dependent/independent variables, see the sections PARAMETERS Statement, STATEVAR Statement, and VAR Statement, respectively.
You can specify the following option after a forward slash (/):
Moreover, the STATE statements must satisfy the following conditions:
The maximum lag order of all the state variables in the distribution-arguments of all STATE statements must be less than or equal to r, where r is defined in the section Method 1: Maximum Lag Order of State Variables, r.
The maximum lag order of all the dependent variables in the distribution-arguments of all the STATE statements, including the PROPOSAL= option, must be less than or equal to p, where p is defined in the section Method 1: Maximum Lag Order of Dependent Variables, p.
The maximum lag order of all the independent variables in the distribution-arguments of all the STATE statements, including the PROPOSAL= option, must be less than or equal to s, where s is defined in the section Method 1: Maximum Lag Order of Independent Variables, s.
The state variables appear as the variable-name in all STATE statements must not have duplicates; that is, each state variable should appear only once as the variable-name in all STATE statements.
For example, suppose the transition equations are
Then the following statements define these initialization equations:
state x1 ~ normal(x1.l1, gamma);
state x2 ~ normal(phi*x2.l1, beta);
where phi, gamma, and beta are parameters and x1 and x2 are current state variables. The first-order lagged state variables x1.l1 and x2.l1 are constructed by appending the suffix ".l1" to the current state variables x1 and x2 (see the section Method 1: Lagged Variables).
If has a proposal transition equation
, the transition equations can be defined by
state x1 ~ normal(x1.l, gamma) / proposal = normal(y, gamma);
state x2 ~ normal(phi*x2.l, beta);
where the state variable x1 has a proposal transition distribution that incorporates the current dependent variable y and the state variable x2 uses the original transition as its proposal transition distribution.