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.
______________________________________________
[email protected] 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.
