
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,)

_________________________________________________________________


	[[alternative HTML version deleted]]


