2017-12-04 4 views
1

逆変換法を使用してパレート分布をシミュレートします。私は比較のためにパレート(2,2)の密度を に重ねてサンプルの濃度ヒストグラムを作成しました。グラフは私が望んでいたようには見えません。R - 軸の値と位取りを変更する

私のy軸は規模と私のx軸(0120)から行くの私が好きな0.02

の規模で(0.00,0.08)から行くの私がしたいです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") 

enter image description here

答えて

1

ソリューションへの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") 

enter image description here

+0

XLIMであり、一般的ではなくxticksまたはyticksを使用しようとするよりも、Rに移動するための方法をYLimプロパティ?私はオプション2が良いbtwが好きです。それはハードコードの代わりに私がやっていることに変わるからです。ありがとうございました。 – DudeWah

関連する問題