seacode / gmacs

A generic size-structured stock assessment model
https://seacode.github.io/gmacs
18 stars 14 forks source link

Catchability priors code unfinished #133

Closed quantifish closed 9 years ago

quantifish commented 9 years ago
// ---Continue with catchability priors-----------------------
int iprior = ntheta + nGrwth + 1; 
for (int i=1; i<=nSurveys; i++)
{
    int itype = int(prior_qtype(i));
    switch(itype)
    {
        // Analytical soln, no prior (uniform, uniformative)
        case 0:
        break;
        // Prior on analytical soln, log-normal
        case 1:
            priorDensity(iprior) = dnorm(log(survey_q(i)),log(prior_qbar(i)),prior_qsd(i));
        break;
    }
    iprior++;
}
smartell commented 9 years ago

what is unfinished? Can add more options.

quantifish commented 9 years ago

I presume case 0 should be a uniform rather than just a constant prior? Or could add both as options...

quantifish commented 9 years ago

I've made this change from constant to uniform line 2128 of gmacs.tpl