Yes, the BLC parameter is optional, and you opted to specify it by using the CON argument. But the CON argument is not the correct format. If you don't want a constrained optimization, then omit the ...
See more...
Yes, the BLC parameter is optional, and you opted to specify it by using the CON argument. But the CON argument is not the correct format. If you don't want a constrained optimization, then omit the CON argument. Otherwise, specify it correctly.
I don't know what distribution you are trying to fit, but I want to point out that your input data set only has three unique pieces of information:
When (x1,x2)=(0,0), the value P(x1,x2)=2.259E-6
When (x1,x2)=(1,1), the value P(x1,x2)=0.999
When (x1,x2)=(1,0) or (0,1), the value P(x1,x2)=0.00033
I also do not think you can perform an optimization in this function because it is unbounded. For example, look at the cross section for theta1=3 as you let theta2 vary:
/* visualize the LL when theta1=3 */
theta1 = 3;
theta2 = T( do(0,8,0.1) );
LL = j(nrow(theta2), 1);
do i = 1 to nrow(theta2);
LL[i] = LogLik( theta1 || theta2[i] );
end;
title "LL(3, theta2)";
call series(theta2, LL) grid={x y};
As you can see, the objective function does not have a maximum. I suggest you check the computations that you used to construct the LogLik function.