On 13. Juni 2012 13:15, Michel Juillard wrote: What about lb = [lb; min(max(estim_params_.corrx(:,4),-1),1)]; Yes you are right, this needs to be changed for both ncx and ncn. -------- Johannes Pfeifer Haußerstr. 29 72076 Tübingen Tel.: +49-(0)7071-6396184 Mobil: +49-(0)170-6936820 jpfeifer@gmx.de Von: dev-bounces@dynare.org [mailto:dev-bounces@dynare.org] Im Auftrag von Michel Juillard Gesendet: Mittwoch, 13. Juni 2012 13:15 An: List for Dynare developers Betreff: Re: [DynareDev] set_prior.m On 06/11/2012 11:15 PM, Johannes Pfeifer wrote: Could someone have a look at the lines for setting the priors for correlation in set_prior.m (occurring in at least two places): ub = [ub; max(min(estim_params_.corrx(:,5),1),-1)]; lb = [lb; max(min(estim_params_.corrx(:,4),1),-1)]; In case estim_params_.corrx(:,4) is NaN, the lower bound defaults to +1 instead of -1. Good catch! My conjecture is it should be lb = [lb; max(min(estim_params_.corrx(:,4),-1),-1)]; No, this would always be -1? What about lb = [lb; min(max(estim_params_.corrx(:,4),-1),1)]; Best Michel -- Johannes Pfeifer Haußerstr. 29 72076 Tübingen Tel.: +49-(0)7071-6396184 Mobil.: +49-(0)170-6936820 Germany _______________________________________________ Dev mailing list Dev@dynare.org https://www.dynare.org/cgi-bin/mailman/listinfo/dev