Re: [R] BFGS and Neldear-Mead

2011-03-21 Thread LC-Bea
Ravi please look my last mail. I tried again, and in my last mail i described the results, i reach the convergence but the values are still very diferent. thank you so much! Lucía -- View this message in context: http://r.789695.n4.nabble.com/BFGS-and-Neldear-Mead-tp3388017p3394334.html Sent f

Re: [R] BFGS and Neldear-Mead

2011-03-21 Thread LC-Bea
I´m triying to give you more information -other function (the problem its the same with both) pen.wls<-function(x){ ax<-x[32:40] indice<-x[1:31] bx<-x[41:49] awlsbind1<-matrix(rep(ax,31),nrow=9,ncol=31) indice<-as.vector(indice-mean(indice)) bx<-as.vector(bx/sum(bx)) me<-(awlsbind1+bx%*%t(indice)

Re: [R] BFGS and Neldear-Mead

2011-03-21 Thread LC-Bea
Sorry for posting too much! i found that in Nelder-Mead the result is "1" i run this algorith too many times with several data maybe the problem is that i have to give a bigger number of iterations Lucía -- View this message in context: http://r.789695.n4.nabble.com/BFGS-and-Neldear-Mead-tp3

Re: [R] BFGS and Neldear-Mead

2011-03-21 Thread LC-Bea
i´ve tried with optimx these are the results pleae help me with this... i´m triying to understand what they mean... > lcmleQN

Re: [R] BFGS and Neldear-Mead

2011-03-21 Thread LC-Bea
Thank you for your answer! the message that gives as result is 0, successful completion so that point wouldnt be the problem sorry if my english its not good (i speak spanish!) thanks! -- View this message in context: http://r.789695.n4.nabble.com/BFGS-and-Neldear-Mead-tp3388017p3393173.

Re: [R] BFGS and Neldear-Mead

2011-03-21 Thread LC-Bea
hello! this is one of the two functions i use pen.mle<-function(x){ ax<-x[32:40] indice<-x[1:31] bx<-x[41:49] amlebind1<-matrix(rep(ax,31),nrow=9,ncol=31) mv<-(-(defunciones*log((exp(amlebind1+bx%*%t(indice)))*poblacion)-((exp(amlebind1+bx%*%t(indice)))*poblacion))) result<-sum(mv)+10*su

[R] BFGS and Neldear-Mead

2011-03-18 Thread LC-Bea
Hello! I´m doing this: A<-optim(vv,pen.wlsv,method="BFGS") B<-optim(vv,pen.wlsv,method="Nelder-Mead") the function is the same and the initial values too the estimate values are different, very different, at the first decimal is there an error? or maybe its because of my data, or funtion to mi