R统计分布绘图

Posted 2021-04-18

我在统计学上的作业遇到了一些麻烦。

我需要以图形方式表示两个内联图中的密度和分布函数，用于我选择的一组参数（必须至少为4），用于Student，Fisher和ChiS重新分区。

我们只举例说明学生重新分配。

根据我在互联网上搜索的内容，我来了：

首先，我需要生成一些随机值。

x <- rnorm( 20, 0, 1 )

问题1：我需要生成4个这个？

然后我必须用以下方法绘制这些值：

plot(dt( x, df = 1))
plot(pt( x, df = 1))

但是，如何为四组参数执行此操作？它们应该在同一个图中表示。

这是我到目前为止的好方法吗？请告诉我，如果我错了。

答案

要绘制某个分布的几个密度，您必须首先得到一个支持向量，在本例中为x。 然后使用您选择的参数计算密度值。 然后绘制它们。

在接下来的代码中，我将绘制4个Sudent-t pdf，其自由度为1到4。

x <- seq(-5, 5, by = 0.01)    # The support vector
y <- sapply(1:4, function(d) dt(x, df = d))

# Open an empty plot first
plot(1, type = "n", xlim = c(-5, 5), ylim = c(0, 0.5))
for(i in 1:4){
    lines(x, y[, i], col = i)
}

然后，您可以通过添加主标题，更改轴标题等使图形更漂亮。

如果你想要其他发行版，比如F或卡方，你将使用x严格为正，例如x <- seq(0.0001, 10, by = 0.01)。