MCMC routine in MARK crashes

[Morten Frederiksen]

In specific circumstances, the MCMC routine in MARK (v 4.1) crashes with the following error message:

forrtl: severe(161): Program Exception - array bounds exceeded

This happens when I try run to use a design matrix with more than one column to model the mean of a hyperdistribution for a specific set of parameters, which includes a couple of unidentifiable parameters (cohorts with no recoveries in a recovery only model). There's no problem if I just want to model the mean of this hyperdistribution, or if I model a different set of parameters with more than one column in the hyperdistribution design matrix. It's not just an issue with the unidentifiable parameters, as excluding these from the hyperdistribution doesn't help.

Any suggestions for why this happens, and what I can do to avoid it?

Morten

[gwhite]

Morton:
Can you send me the input file that is causing this error, so that I can fix the code? The easiest way to create the file is to use the "Save Structure" button, and then open up the saved model with the full output option into notepad. Save this file, and then send it to me via email. Thanks.

Gary

[egc]

In specific circumstances, the MCMC routine in MARK (v 4.1) crashes with the following error message:

As per the MARK FAQ, step (1) should be to see if the problem *goes away* if you use the most recent build of MARK. Version 4.1 is pretty out of date (current version is 4.3, and is date stamped 3 March, 2007).