Hi all, I have some some problem with regard to finding the integral of a function containing an indicator function. please see the code below:
func1 <- function(x, mu){ (mu^2)*dnorm(x, mean = mu, sd = 1)*dgamma(x, shape=2)} m1star <- function(x){ integrate(func1, lower = 0, upper = Inf,x)$val} T <- function(x){ 0.3*dnorm(x)/(0.3*dnorm(x)+0.7*m1star(x))} func2 <- function(x,c){(T(x) <=c)*0.3*dnorm(x)} func3 <- function(x,c){(T(x) <= c)*(0.3*dnorm(x)+0.7*m1star(x))} numer <- function(c){ integrate(func2, -Inf, Inf, c)$val} denom <- function(c){ integrate(func3, lower, Inf,c)$val} The error message is as below : > numer(0.5) Error in integrate(func1, lower = 0, upper = Inf, x) : the integral is probably divergent Thank you very much! Hannah [[alternative HTML version deleted]] ______________________________________________ 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.