Вы можете добавить область под кривой плотности с добавлением полигонов. Границы зависят от x
в ваших seq
(желаемый интервал в х), а также от m1$z2
значений, когда x
находится между выбранным интервалом x
.
polygon(c(.3, x[x >= .3 & x <= .4], .4),
c(0, m1$z2[x >= .3 & x <= .4], 0), col="red")
polygon(c(.2, x[x >= .2 & x <= .3], .3),
c(0, m1$z2[x >= .2 & x <= .3], 0), col="light blue")