如何在 r 中绘制 t 分布
为了绘制 R 中 at 分布的概率密度函数,我们可以使用以下函数:
- dt(x, df)创建概率密度函数
- curve(function, from = NULL, to = NULL)绘制概率密度函数
为了绘制概率密度函数,我们需要在dt()函数中指定df (自由度)以及curve()函数中的from和to值。
例如,以下代码说明了如何绘制自由度为 10 的分布的概率密度函数,其中绘图的 x 轴介于 -4 和 4 之间:
curve(dt(x, df=10), from=-4, to=4)
与正态分布类似,t 分布围绕均值 0 对称。
我们可以添加标题、更改 Y 轴标签、增加行宽,甚至更改行颜色,以使绘图更加美观:
curve(dt(x, df=10), from=-4, to=4, main = 't Distribution (df = 10)', #add title ylab = 'Density', #change y-axis label lwd = 2, #increase line width to 2 col = 'steelblue') #change line color to steelblue
我们还可以在图中添加多条曲线来比较不同自由度的 t 分布。例如,以下代码创建 df = 6、df = 10 和 df = 30 的 t 分布图:
curve(dt(x, df=6), from=-4, to=4, col='blue') curve(dt(x, df=10), from=-4, to=4, col='red', add=TRUE) curve(dt(x, df=30), from=-4, to=4, col='green', add=TRUE)
我们可以使用legend()函数向图中添加图例,该函数采用以下语法:
图例(x,y = NULL,图例,填充,col,bg,lty,cex)
金子:
- x, y:用于定位图例的 x 和 y 坐标
- 图例:要放入图例中的文本
- fill:图例内的填充颜色
- col:用于图例内线条的颜色列表
- bg:图例的背景颜色
- lty:线条样式
- cex:图例中文本的大小
在我们的示例中,我们将使用以下语法来创建图例:
#create density plots curve(dt(x, df=6), from=-4, to=4, col='blue') curve(dt(x, df=10), from=-4, to=4, col='red', add=TRUE) curve(dt(x, df=30), from=-4, to=4, col='green', add=TRUE) #add legend legend(-4, .3, legend=c("df=6", "df=10", "df=30"), col=c("blue", "red", "green"), lty=1, cex=1.2)