Re: [R] Optim function & Hessian matrix

2013-11-19 Thread David Winsemius
This is the code that was attached. I arrived in the copy sent to my email address: #-- L1=function(X){ B1=X[1] C1=X[2] B2=X[3] C2=X[4] mu=X[5] S=-(B1*(C1^x)*((C1^v)-1)/log(C1)) +(d2*log(B1*(C1^(x+v -(B2*(C2^y)*((C2^v)-1)/log(C2)) +(d1*log(B2*(C2^(y+v -(v*mu) +

Re: [R] Optim function & Hessian matrix

2013-11-19 Thread David Winsemius
On Nov 19, 2013, at 12:09 PM, sofeira taajobian wrote: > Dear R Users > Hi, > > > I have very emergency problems in my programming about finding MLE > with optim command. I reproduced it with real data. I guess that my > function object in optim is very sensitive because it has power > func