Hola, qué función tiene R para optimizar funciones bidimensionales. Gracias.
Hola:
En R tienes muchas opciones para resolver problemas de optimización, te recomiendo que mires el "Task View" de optimización en
http://cran.es.r-project.org/web/views/Optimization.html
Saludos.
En R tienes muchas opciones para resolver problemas de optimización, te recomiendo que mires el "Task View" de optimización en
http://cran.es.r-project.org/web/views/Optimization.html
Saludos.
Hola, he encontrado la función optim pero no se muy bien como implementarla. La función que necesito maximizar es la siguiente:
L<-function(pe,rho){
sum(A*log(pe)+(B-A)*log(1-pe)+sum(log(1+(rho/(2*pe^2*(1-pe)^2))*((A-B*pe)^2)+A*(2*pe-1)-B*pe^2)))}
Donde pe y rho toman valores en (0,1) y A y B son vectores introducidos previemente.
Si hago:
optim(c(0.01,0.01),L)
me dice que rho está perdido. No entiendo cuál sería la forma correcta de hacerlo.
Gracias.
Si miras la ayuda de optim (con ?optim) verás que optim espera una función cuyo argumento sea un vector y no una lista de argumentos.
Si defines una función que dependa de x y sustituyes pe por x[1] y rho por x[2]. Por hacerlo rápido, te debería funcionar:
M<- function(x) L(x[1], x[2])
optim(c(.01, .01), M)
El intervalo lo tienes que especificar usando los argumentos lower y upper de optim.
Saludos.
Si defines una función que dependa de x y sustituyes pe por x[1] y rho por x[2]. Por hacerlo rápido, te debería funcionar:
M<- function(x) L(x[1], x[2])
optim(c(.01, .01), M)
El intervalo lo tienes que especificar usando los argumentos lower y upper de optim.
Saludos.
Hola! muchas gracias, esto es lo que necesitaba aunque de todos modos no me sirve porque se producen NaN's en mi función y me da error.
Hola, he estado utilizando el comando "optim" para maximizar una función de verosimilitud y me funciona bien pero tarda demasiado tiempo. Qué otras funciones hay en R que hagan lo mismo que "optim"? Gracias.