Important: Read this before posting to this forum

  1. This forum is for questions related to the use of Apollo. We will answer some general choice modelling questions too, where appropriate, and time permitting. We cannot answer questions about how to estimate choice models with other software packages.
  2. There is a very detailed manual for Apollo available at http://www.ApolloChoiceModelling.com/manual.html. This contains detailed descriptions of the various Apollo functions, and numerous examples are available at http://www.ApolloChoiceModelling.com/examples.html. In addition, help files are available for all functions, using e.g. ?apollo_mnl
  3. Before asking a question on the forum, users are kindly requested to follow these steps:
    1. Check that the same issue has not already been addressed in the forum - there is a search tool.
    2. Ensure that the correct syntax has been used. For any function, detailed instructions are available directly in Apollo, e.g. by using ?apollo_mnl for apollo_mnl
    3. Check the frequently asked questions section on the Apollo website, which discusses some common issues/failures. Please see http://www.apollochoicemodelling.com/faq.html
    4. Make sure that R is using the latest official release of Apollo.
  4. If the above steps do not resolve the issue, then users should follow these steps when posting a question:
    1. provide full details on the issue, including the entire code and output, including any error messages
    2. posts will not immediately appear on the forum, but will be checked by a moderator first. This may take a day or two at busy times. There is no need to submit the post multiple times.

Errors in Hybrid Choice Model with Bayesian Estimator

Ask questions about the results reported after estimation. If the output includes errors, please include your model code if possible.
Post Reply
janak12_jp
Posts: 6
Joined: 08 Sep 2021, 16:52

Errors in Hybrid Choice Model with Bayesian Estimator

Post by janak12_jp »

Hello,

I have tried several specifications of Hybrid Choice Model with Bayesian Estimator to solve my choice problem. My data contains several mode specific LVs in addition to other observable attributes like trip characteristics. I used initial 100k iterations as burn-in and next 100k for estimation. But after completing the estimation, I come across the same error each time:

Code: Select all

WARNING: RSGHB has censored the probabilities. Please note that in at least some iterations RSGHB has avoided numerical
issues by left censoring the probabilities. This has the side effect of zero or negative probabilities not leading to
failures!
Warning messages:
1: In log(test2_LL) : NaNs produced
2: In log(test1_LL) : NaNs produced
The problem as I seen in my results is most probably related to estimation of covariance matrix (though not 100% sure).
It shows the initial loglikelihood value to -Inf though it was somewhere around -56000 when I started estimation.

I produces correlation matrices but NaNs in covariance matrices, what is the reason for that? I read somewhere that correlation matrix is unstandardized version of covariance matrix. Can you please suggest on this also?

I am attaching the code as well as my results.
Attachments
ICLV.zip
Code and Results
(28.29 KiB) Downloaded 367 times
dpalma
Posts: 190
Joined: 24 Apr 2020, 17:54

Re: Errors in Hybrid Choice Model with Bayesian Estimator

Post by dpalma »

Hi,

It's difficult to diagnose the problem without looking at your data. Could you share your database? If you don;t want to share it in the forum, you can email it to D.Palma [at] leeds.ac.uk

Cheers
David
janak12_jp
Posts: 6
Joined: 08 Sep 2021, 16:52

Re: Errors in Hybrid Choice Model with Bayesian Estimator

Post by janak12_jp »

Dear Dr Palma,

I have mailed the data file to the stated email. Thank you.

Regards,
Janak
dpalma
Posts: 190
Joined: 24 Apr 2020, 17:54

Re: Errors in Hybrid Choice Model with Bayesian Estimator

Post by dpalma »

Hi Janak,

Sorry for the slow response.

First of all, you are not getting errors, but only warnings. The first warning (related to the censoring of probabilities) means that during the estimation process (i.e. for some values in your chain) the likelihood values you obtained were so small that they were indistinguishable from zero. To avoid numerical issues Apollo replaced those values by a small (but bigger than zero) value. If the chain looks like it converged correctly, then this was probably due to the chain going through these problematic values only for a short time, and then moved on to better values.

The second warning (NaN in log(test2_LL)) has to do with your starting values leading to very a small likelihood, i.e. too close to zero. Considering the previous warning, my guess is your starting values are not very good. I would recommend using your estimated values as starting values, and performing the estimation again. You will probably run into fewer warnings then.

Finally, concerning the NaNs in the "Covariances of random parameters", this was to be expected as you set the option gFULLCV= FALSE inside apollo_HB. This means you are forcing the chains of the random parameters to be uncorrelated, that is why the covariance matrix only has values for the diagonal (e.g. asc_bus_asc_bus, asc_air_asc_air), but not for off-diagonal elements (e.g. asc_air_asc_bus).

Finally, note that the number of draws to keep in the chain is defined in the setting "gNEREP", not "GNEREP". R is case sensitive, so using upper or lower case does matter.

Cheers
David
Post Reply