ivregimes {hspm}R Documentation

Estimation of regime models with endogenous variables

Description

The function ivregimes deals with the estimation of regime models. Most of the times the variable identifying the regimes reveals some spatial aspects of the data (e.g., administrative boundaries). The model includes exogenous as well as endogenous variables among the regressors.

Usage

ivregimes(formula, data, rgv = NULL, vc = c("homoskedastic", "robust", "OGMM"))

Arguments

formula

a symbolic description of the model of the form y ~ x_f | x_v | h_f | h_v where y is the dependent variable, x_f are the regressors that do not vary by regimes, x_v are the regressors that vary by regimes, h_f are the fixed instruments and h_v are the instruments that vary by regimes.

data

the data of class data.frame.

rgv

an object of class formula to identify the regime variables

vc

one of c("homoskedastic", "robust", "OGMM"). If "OGMM" an optimal weighted GMM is used to estimate the VC matrix.

Details

The basic (non spatial) model with endogenous variables can be written in a general way as:

y=[X100X2][β1β2]+Xβ+[Y100Y2][π1π2]+Yπ+ε y = \begin{bmatrix} X_1& 0 \\ 0 & X_2 \\ \end{bmatrix} \begin{bmatrix} \beta_1 \\ \beta_2 \\ \end{bmatrix} + X\beta + \begin{bmatrix} Y_1& 0 \\ 0 & Y_2 \\ \end{bmatrix} \begin{bmatrix} \pi_1 \\ \pi_2 \\ \end{bmatrix} + Y\pi + \varepsilon

where y=[y1,y2]y = [y_1^\prime,y_2^\prime]^\prime, and the n1×1n_1 \times 1 vector y1y_1 contains the observations on the dependent variable for the first regime, and the n2×1n_2 \times 1 vector y2y_2 (with n1+n2=nn_1 + n_2 = n) contains the observations on the dependent variable for the second regime. The n1×kn_1 \times k matrix X1X_1 and the n2×kn_2 \times k matrix X2X_2 are blocks of a block diagonal matrix, the vectors of parameters β1\beta_1 and β2\beta_2 have dimension k1×1k_1 \times 1 and k2×1k_2 \times 1, respectively, XX is the n×pn \times p matrix of regressors that do not vary by regime, β\beta is a p×1p\times 1 vector of parameters. The three matrices Y1Y_1 (n1×qn_1 \times q), Y2Y_2 (n2×qn_2 \times q) and YY (n×rn \times r) with corresponding vectors of parameters π1\pi_1, π2\pi_2 and π\pi, contain the endogenous variables. Finally, ε=[ε1,ε2]\varepsilon = [\varepsilon_1^\prime,\varepsilon_2^\prime]^\prime is the n×1n\times 1 vector of innovations. The model is estimated by two stage least square. In particular:

Value

An object of class ivregimes. A list of five elements. The first element of the list contains the estimation results. The other elements are needed for printing the results.

Author(s)

Gianfranco Piras and Mauricio Sarrias

Examples

data("natreg")
form   <- HR90  ~ 0 | MA90 + PS90 + RD90 + UE90 | 0 | MA90 + PS90 + RD90 + FH90 + FP89 + GI89
split  <- ~ REGIONS
mod <- ivregimes(formula = form, data = natreg, rgv = split, vc = "robust")
summary(mod)
mod1 <- ivregimes(formula = form, data = natreg, rgv = split, vc = "OGMM")
summary(mod1)
form1   <- HR90  ~ MA90 + PS90 |  RD90 + UE90 -1 | MA90 + PS90 | RD90 + FH90 + FP89 + GI89 -1
mod2 <- ivregimes(formula = form1, data = natreg, rgv = split, vc = "homoskedastic")
summary(mod2)


[Package hspm version 1.1 Index]