我有以下任务:从t分布和受某些分位数限制的阴影区域绘制密度曲线。我想在ggplot2中做到这一点。
使用以下代码绘制曲线:
ggplot(data = data.frame(x = c(-5, 5)), aes(x)) +
stat_function(fun = dt, args = list(df = 70))+
geom_segment(aes(x=qt(.95,70),xend=qt(.95,70),y=0,yend=dt(qt(.95,70),70)))
ylab("f(t)")+xlab("t")
然后我要对大于qt(.95,70)
但小于密度的数字段的区域进行阴影处理
我希望有一种解决此问题的方法,而不必模拟密度数据来绘制并遮盖所需的区域,如此处Shaded area under density curve in GGPLOT2所示,因为密度不是由任何数据产生的,而是由方程。