GenerateMDCEVData: GenerateMDCEVData

In rmdcev: Kuhn-Tucker and Multiple Discrete-Continuous Extreme Value Models

GenerateMDCEVData

Description

Simulate data for KT models

Usage

GenerateMDCEVData(
  model,
  nobs = 1000,
  nalts = 10,
  income = stats::runif(nobs, 1e+05, 150000),
  price = matrix(stats::runif(nobs * nalts, 100, 500), nobs, nalts),
  alpha_parms = 0.5,
  scale_parms = 1,
  gamma_parms = stats::runif(nalts, 1, 10),
  psi_i_parms = c(-1.5, 2, -1),
  psi_j_parms = c(-5, 0.5, 2),
  phi_parms = c(-5, 0.5, 2),
  dat_psi_i = matrix(2 * stats::runif(nobs * length(psi_i_parms)), nobs,
    length(psi_i_parms)),
  dat_psi_j = cbind(matrix(stats::runif(nalts * (length(psi_j_parms)), 0, 1), nrow =
    nalts)),
  dat_phi = cbind(matrix(stats::runif(nalts * (length(phi_parms)), 0, 1), nrow =
    nalts)),
  nerrs = 1,
  tol = 1e-20,
  max_loop = 999
)

Arguments

model	A string indicating which model specification is estimated. The options are "alpha", "gamma", "hybrid" and "hybrid0" for the MDCEV model and "kt_ee" for the environmental economics Kuhn-Tucker specification.
nobs	Number of individuals
nalts	Number of non-numeraire alts
income	Vector of individual income
price	Matrix of prices for non-numeraire alternatives.
alpha_parms	Parameter value for alpha term
scale_parms	Parameter value for scale term
gamma_parms	Parameter value for gamma terms
psi_i_parms	Parameter value for psi terms that vary by individual
psi_j_parms	Parameter value for psi terms that vary by alt (all models except kt_ee)
phi_parms	Parameter value for phi terms that vary by alt (kt_ee model only)
dat_psi_i	(nobs X # psi_i_parms) matrix with individual-specific characteristics
dat_psi_j	(nalts X # psi_j_parms) matrix with alternative-specific variables (all models except kt_ee)
dat_phi	(nalts X # phi_parms) matrix with alternative-specific variables (kt_ee model only)
nerrs	Number of error draws for demand simulation
tol	Tolerance level for simulations if using general approach
max_loop	maximum number of loops for simulations if using general approach

Value

A 'mdcev.data' object, which is a 'data.frame' in long format. Also includes parms_true with parameter values

Examples

data <- GenerateMDCEVData(model = "gamma")

rmdcev documentation built on March 31, 2023, 6:49 p.m.