我有以下pdf,
f(x)=2(1-x)^2; 2=>x>1
现在我必须在cdf
中绘制随机变量的r
。
#cumulative distribution function
Fx <- function(x){
(2/3) * (x^3-3*x^2+3*x-1)
}
R <- runif(100,1,2)
cum_prob <- Fx(R) # determining cumulative probability
但我不知道应该使用哪个命令来绘制cdf
。
hist()
&amp; barplot()
不适合cdf。因为据我所知绘制一个cdf,它需要X轴上的随机变量值和Y轴上的累积概率。
我的变量也是连续的。
答案 0 :(得分:2)
回应评论中的代码:
fx <- function(x)2*(x-1)^2
integrate(fx,1,2)
#------------
0.6666667 with absolute error < 7.4e-15
plot(seq(1,2,by=0.01), cumsum(fx(seq(1,2,by=0.01))*0.01 ), type="l" )
答案 1 :(得分:2)
您可以使用:
定义密度:
f <- function(x) (2*(1-x)^2)*(1<x & x<=2)
绘图密度:
plot(f,0,3)
绘制CDF:
plot(Vectorize(function(X)integrate(f,0,X)$value),0,3,add=TRUE)
注意:0和3是图中x轴的限制;你可以改变它们。另请注意,您的密度未标准化。