R中卡方的密度图

时间:2014-05-26 00:57:40

标签: r graph plot chi-squared probability-density

我试图在R中绘制Chi-Squared分布的密度函数图,其中28 df是x高于7.5。

直到现在,我从我能够聚集的地方得到了这个:

x <- pchisq(7.5, 28, lower.tail=FALSE)
hist(x, prob=TRUE)
curve( dchisq(x, df=28), col='red', main = "Chi-Square Density Graph")

但是绘图似乎不起作用._。

1 个答案:

答案 0 :(得分:3)

我想你想要:

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))

在这种情况下,图表有点愚蠢,因为曲线下的几乎所有区域都有阴影 - x>7.5的概率

pchisq(7.5,df=28,lower.tail=FALSE)
## [1] 0.9999611