我在统计数据上的作业遇到了一些麻烦。
我需要以图形方式表示两个内联图中的密度和分布函数,以获得我选择的一组参数(必须至少为4),用于Student,Fisher和ChiS重新分区。
我们只举例说明学生重新分配。
根据我在互联网上搜索的内容,我来了:
首先,我需要生成一些随机值。
x <- rnorm( 20, 0, 1 )
问题1:我需要生成4个吗?
然后我必须用以下方式绘制这些值:
plot(dt( x, df = 1))
plot(pt( x, df = 1))
但是,如何为四组参数执行此操作?它们应该在同一个图中表示。
这对我到目前为止来说是一个好方法吗? 请告诉我,如果我错了。
答案 0 :(得分:2)
要绘制某个分布的多个密度,您必须先得到一个支持向量,在这种情况下为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)
。