r - プロット密度関数

Question

確率変数 Y は、次の式で与えられる密度 f の分布を持ちます。

f(y)= y-1/6 if 1<=y<=4

  -(y-5/6) if 4<=y<=5

   0  else

R を使用して、 Y の密度関数のプロットを作成します。これは私の質問です。Rプログラムを使用してプロットするにはどうすればよいですか?

Answer 1

使用するifelse

  (y <- c(1:6)) 
[1] 1 2 3 4 5 6
(fy <-  ifelse(1<=y & y<4,  y-1/6, ifelse(4<=y & y<=5, -(y-5/6), 0)))
[1]  0.8333333  1.8333333  2.8333333 -3.1666667 -4.1666667  0.0000000
plot(fy~y, type='l') 

私は次のように仮定しました: if 1 <= y < 4 then y-1/6 4 <= y <=5 then -(y-5/6) y>5 then 0

これにより、1<=y<=4 と 4<=y<=5 の重複が回避されます