Displaying 1 result from an estimated 1 matches for "skellamreg_ll".
2009 Jul 02
2
constrained optimisation in R.
i want to estimate parameters with maximum likelihood method with contraints (contant numbers).
for example
sum(Ai)=0 and sum(Bi)=0
i have done it without the constraints but i realised that i have to use the contraints.
Without constraints(just a part-not complete):
skellamreg_LL=function(parameters,z,design)
{
n=length(z);
mu=parameters[1];
H=parameters[2];
Apar=parameters[3:10];
Dpar=parameters[11:18];
res=optim(par,skellamreg_LL,method="Nelder-Mead", hessian = FALSE, control=list(trace=1000,maxit=100000),z=z,design=x,)
_________________________...