R help es - Oct 2010 - what´s wrong with this code?

Hello,  I  want  to  maximize  a  likelihood  function expressed as an
integral  that  can not be symbolically evaluated. I expose my problem
in a reduced form.

g<- function(x){
integrand<-function(y) {exp(-x^2)*y}
g<-integrate(integrand,0,1)
}                         
h<-function(x) log((g(x)))

g  is an object of the class function, but g(2) is a integrate object,
I can print(g(2)) an get a result, but
if I try h(2) R says is a nonnumeric argument for a mathematical
function. My goal is to maximize h.
what´s wrong?

thanks

José,

escoge más bien la función

h<-function(x) log((g(x)$value))

¡En todo caso, el foro es en español!

Un saludo. Olivier
--  
____________________________________

Olivier G. Nuñez
Email: onunez en iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es

____________________________________




El 29/10/2010, a las 11:20, José Manuel Gavilán Ruiz escribió:
> g<- function(x){
> integrand<-function(y) {exp(-x^2)*y}
> g<-integrate(integrand,0,1)
> }
> h<-function(x) log((g(x)))

Por otro lado, el calculo de la integral de exp(-x^2)*y sobre y=0..1,  
parece más bien trivial!
¿no habrá un error en la formulación de tu función g?
Un saludo. Olivier
--  
____________________________________

Olivier G. Nuñez
Email: onunez en iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es

____________________________________




El 29/10/2010, a las 11:20, José Manuel Gavilán Ruiz escribió:
> Hello,  I  want  to  maximize  a  likelihood  function expressed as an
> integral  that  can not be symbolically evaluated. I expose my problem
> in a reduced form.
>
> g<- function(x){
> integrand<-function(y) {exp(-x^2)*y}
> g<-integrate(integrand,0,1)
> }
> h<-function(x) log((g(x)))
>
> g  is an object of the class function, but g(2) is a integrate object,
> I can print(g(2)) an get a result, but
> if I try h(2) R says is a nonnumeric argument for a mathematical
> function. My goal is to maximize h.
> what´s wrong?
>
> thanks
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es