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 中的重叠