Dear list,
I am using Stata 11
Stata/IC 11.0 for Windows (32-bit)
Born 21 Oct 2009
After my last update I have problems with -optimize-. The following example
does not work anymore:
#delimit;
mata ;
x=(0.35 \ 0.29 \ 0.3\ 0.3\ 0.65\ 0.56\ 0.37\ 0.16\ 0.26\ 0.19)
void lnbetaden04(todo, p, x, lnf, S, H)
{
a=p[1]
b=p[2]
lnf= lngamma(a+b) :- lngamma(a) :- lngamma(b) :+ (a-1)*log(x) :+
(b-1)*log(1:-x)
}
mata
S=optimize_init()
optimize_init_which(S, "max")
optimize_init_evaluator(S, &lnbetaden04())
optimize_init_evaluatortype(S, "v0")
optimize_init_params(S, (1,1))
optimize_init_argument(S, 1, x)
optimize_init_conv_ptol(S, 1e-8)
optimize_init_tracelevel(S, "params")
optimize_init_tracelevel(S, "step")
optimize_init_tracelevel(S, "gradient")
optimize_init_nmsimplexdeltas(S, (1e-3, 1e-3))
optimize_init_technique(S,"nr")
p=optimize(S)
end;
After "iteration 0" I get the following message:
Iteration 1:
Parameter vector:
c1 c2
r1 4.084429 7.868267
Changed elements of parameter vector:
select(): 3200 conformability error
opt__looputil_step(): - function returned error
opt__loop_nr(): - function returned error
optimize(): - function returned error
<istmt>: - function returned error
r(3200);
I found that the technique "nm" does work, but the techniques "nr", "bhhh"
and "dfp" are not working anymore. They do work in the old Stata 10 version
on my old computer (Stata/IC 10.1 for Windows
Born 02 Feb 2009).
What could be the problem?
Best regards,
Marike
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
