sdpdth (version 0.2)

msdpd: M-estimator for spatial dynamic panel data model

Description

Estimating the spatial dynamic panel data model with M-estimator

Usage

msdpd(
  y,
  x,
  w1,
  correction = TRUE,
  hessian_er = FALSE,
  true_range = FALSE,
  max_try = 5,
  w2 = w1,
  w3 = w1,
  no_tf = FALSE,
  model = "full",
  rcpp = TRUE,
  cma_pop_multi = 1
)

Arguments

y

matrix, containing regional index (first column), time index (second column, numeric) and dependent variable (third column, numeric).

x

matrix, containing regional index (first column), time index (second column, numeric) and regressors (numeric).

w1

matrix, the spatial weight matrix. If w2 and w3 are supplied, the spatial weight matrix for spatial lag.

correction

logical, whether to use adjusted score function. Default value is TRUE.

hessian_er

logical, whether to output hessian based se. Ignored if correction is set to False. Default value is FALSE.

true_range

logical, whether to used the accurate stationary check. Default value is FALSE due to performance reasons.

max_try

integer, maximum attempt for the solver. Default value is 5.

w2

matrix, the spatial weight matrix for spatio-temporal lag. Default value is the same as w1.

w3

matrix, the spatial weight matrix for spatial error. Default value is the same as w1.

no_tf

logical, whether to account for time effect. Default value is TRUE.

model

character, indicates the model used for estimation, can be "full", "slm", "sem", "sltl". See Details.

rcpp

logical, whether to use the rcpp implementation to calculate the score function. Default value is TRUE.

cma_pop_multi

integer, multiplier for the population size used in CMA-ES. Default value is 1.

Value

A list of estimation results of S3 class "msdpd"

• "coefficient" list, coefficients and standard errors

• "model" character, model used for estimation

• "vc_mat" matrix, variance-covariance matrix

• "hessian" matrix, optional, hessian matrix

Details

Estimating the spatial dynamic panel data model with Yang(2018)'s M-estimator $$ y_{ti} = \mu_{i}+\alpha_t + x_{ti}\beta + \rho y_{t-1,i} + \lambda_1 \sum_{j =1}^{n}w_{1,ij}y_{tj} + \lambda_2 \sum_{j =1}^{n}w_{2,ij}y_{t-1,j} + u_{ti},\\ u_{ti} = \lambda_3\sum_{j =1}^{n}w_{3,ij}u_{tj} + v_{ti}, i=1,\ldots,n,t=1,\ldots,T$$ The minimum number of time-periods is 4. Make sure the rows and columns of w1, w2, and w3 are lined up with the regional index. Sub-models can be specified by argument "model"

• "full" Full model

• "slm" \(\lambda_2 = \lambda_3 = 0\)

• "sem" \(\lambda_1 = \lambda_2 = 0\)

• "sltl" \(\lambda_3 = 0\)

Some suggestions when the optimizer fails:

• Increase max_try

• Increase cma_pop_multi

• try a different submodel

References

Yang, Z. (2018). Unified M-estimation of fixed-effects spatial dynamic models with short panels. Journal of Econometrics, 205(2), 423-447.

Examples

# NOT RUN {
data(data_n, data_nw)
result <- msdpd(y = data_n$y, x = data_n$x, w1 = data_nw)
# }
# NOT RUN {
# }