R中的卡方密度图

Question

如何在R？中绘制卡方密度图？

我收到了以下代码，但我不确定如何操纵它们：

curve( dchisq(x, df=28), col='red', main = "Chi-Square Density Graph",
          from=0,to=60)
xvec <- seq(7.5,60,length=101)
pvec <- dchisq(xvec,df=28)
polygon(c(xvec,rev(xvec)),c(pvec,rep(0,length(pvec))),
        col=adjustcolor("black",alpha=0.3))

有人可以解释代码的含义吗？

Answer 1

包ggplot2提供了一种简单的方法来绘制卡方分布。您只需指定一个stat_function dchisq作为您的函数，然后指定一个列表args来表示自由度。

例如，下面是4自由度的卡方分布的示例代码：

library(ggplot2)

ggplot(data.frame(x = c(0, 20)), aes(x = x)) +
     stat_function(fun = dchisq, args = list(df = 4))