计算中位数。错误在哪里？

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]