【发布时间】：2013-11-28 03:22:24
【问题描述】：
我想计算以下密度的中位数。但它不起作用。
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
【问题讨论】：
-
你需要解释你期望它做什么。
-
@joran 我已经编辑了我的问题。请检查一下。
-
所以对于这个密度，你可以直接写下答案。为什么要进行计算？这只是一个玩具例子吗？还是作业问题？
-
@Aaron 是否禁止在
stack overflow中提出与家庭作业有关的问题，如果一个人表现出他/她的努力？ -
不，不禁止。但背景很重要。在这种情况下，您不应该以尝试的方式得到答案，而应该简单地整合。除非有一些你没有告诉我们的背景。