Вычислительная медиана. Где ошибка?

Question

Я хочу вычислить медиану следующей плотности. Но это не работает. не

f <- function(x)(3/7)*x^2 

    r <- rep(0,5000) 
    x=seq(1,2,length=5000) 

for(i in 1:5000){ 
    r[i]=integrate(f,lower=1,upper=x[i])$value 
    } 

d <- data.frame(x,r) 
med <- d$x[d$r==0.5] 
med 

Answer 1

Нет элементов d$r равной точно 0,5. Я подозреваю, что вам нужно что-то вроде:

med <- d$x[d$r>=0.5 & d$r <=0.501]