使用逆变换方法来模拟帕累托分布。我用帕累托(2,2)密度制作了样本的密度直方图 叠加用于比较。尽管图表看起来并不像我想要的那样。
我希望我的y轴从(0.00,0.08)开始,比例为.02
我希望我的x轴从(0,120)开始,比例为20。
我尝试使用一些命令,如xticks,yticks,但无法获得我想要的结果。有什么建议吗?
x <- runif(100, min = 0, max = 1)
a <- 2
b <- 2
gen_pareto <- b/(1-x)^(1/a)
plot(x, gen_pareto)
hist(gen_pareto, prob = TRUE)
lines(density(gen_pareto), # density plot
lwd = 3, # thickness of line
col = "blue")
答案 0 :(得分:1)
以下是解决方案的一种方法:
par(mfrow=c(1,2))
hist(gen_pareto, prob=T,xlim = c(0,120), ylim = c(0,.8))
lines(density(gen_pareto),
lwd = 3,
col = "red")
但是,我觉得下面的解决方案要好得多,但它并没有限定你所说的边界条件。
hist(gen_pareto, prob=T,xlim = c(0, max(density(gen_pareto)$x)), ylim = c(0, max(density(gen_pareto)$y)))
lines(density(gen_pareto),
lwd = 3,
col = "red")