Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!
Get Started for Free
Documentation
R
apollo

apollo_estimateHB

Estimates model

Description

Estimates a model using the likelihood function defined by apollo_probabilities.

Usage

apollo_estimateHB(
  apollo_beta,
  apollo_fixed,
  apollo_probabilities,
  apollo_inputs,
  estimate_settings = NA
)

Arguments

apollo_beta

Named numeric vector. Names and values for parameters.
apollo_fixed

Character vector. Names (as defined in apollo_beta) of parameters whose value should not change during estimation.
apollo_probabilities

Function. Returns probabilities of the model to be estimated. Must receive three arguments:

  • apollo_beta: Named numeric vector. Names and values of model parameters.

  • apollo_inputs: List containing options of the model. See apollo_validateInputs.

  • functionality: Character. Can be either "estimate" (default), "prediction", "validate", "conditionals", "zero_LL", or "raw".

apollo_inputs

List grouping most common inputs. Created by function apollo_validateInputs.
estimate_settings

List. Options controlling the estimation process.

  • estimationRoutine: Character. Estimation method. Can take values "bfgs", "bhhh", or "nr". Used only if apollo_control$HB is FALSE. Default is "bfgs".

  • maxIterations: Numeric. Maximum number of iterations of the estimation routine before stopping. Used only if apollo_control$HB is FALSE. Default is 200.

  • writeIter: Boolean. Writes value of the parameters in each iteration to a csv file. Works only if estimation_routine="bfgs". Default is TRUE.

  • hessianRoutine: Character. Name of routine used to calculate the Hessian of the loglikelihood function after estimation. Valid values are "numDeriv" (default) and "maxLik" to use the routines in those packages, and "none" to avoid estimating the Hessian (and the covariance matrix). Only used if apollo_control$HB=FALSE.

  • printLevel: Higher values render more verbous outputs. Can take values 0, 1, 2 or 3. Ignored if apollo_control$HB is TRUE. Default is 3.

  • constraints: Constraints on parameters to estimate. Should ignore fixed parameters. See argument constraints in maxBFGS for more details.

  • scaling: Named vector. Names of elements should match those in apollo_beta. Optional scaling for parameters. If provided, for each parameter i, (apollo_beta[i]/scaling[i]) is optimised, but scaling[i]*(apollo_beta[i]/scaling[i]) is used during estimation. For example, if parameter b3=10, while b1 and b2 are close to 1, then setting scaling = c(b3=10) can help estimation, specially the calculation of the Hessian. Reports will still be based on the non-scaled parameters.

  • numDeriv_settings: List. Additional arguments to the Richardson method used by numDeriv to calculate the Hessian. See argument method.args in grad for more details.

  • bootstrapSE: Numeric. Number of bootstrap samples to calculate standard errors. Default is 0, meaning no bootstrap s.e. will be calculated. Number must zero or a positive integer. Only used if apollo_control$HB is FALSE.

  • bootstrapSeed: Numeric scalar (integer). Random number generator seed to generate the bootstrap samples. Only used if bootstrapSE>0. Default is 24.

  • silent: Boolean. If TRUE, no information is printed to the console during estimation. Default is FALSE.

Details

This is the main function of the Apollo package. The estimation process begins by checking the definition of apollo_probabilities by estimating it at the starting values. Then it runs the function with argument functionality="validate". If the user requested more than one core for estimation (i.e. apollo_control$nCores>1), and no bayesian estimation is used (i.e. apollo_control$HB=FALSE), then a cluster is created. Using a cluster at least doubles the requires RAM, as the database must be copied into the cluster. If all checks are passed, estimation begins. There is no limit to estimation time other than reaching the maximum number of iterations. If bayesian estimation is used, estimation will finish once the predefined number of iterations are completed. This functions does not save results into a file nor prints them into the console, so if users want to see and store estimation the results, they must make sure to call function apollo_modelOutput and/or apollo_saveOutput afterwards.

Value

model object

apollo

Tools for Choice Model Estimation and Application

v0.2.4
GPL-2
Authors
Stephane Hess [aut], David Palma [aut, cre]
Initial release

We don't support your browser anymore

Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.