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)+1000000000*sum(((t(u)%*%x)-c)^2) result } thankssss! -- View this message in context: http://r.789695.n4.nabble.com/BFGS-and-Neldear-Mead-tp3388017p3393167.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.