问题描述
我想绘制一张带有自由度$ df \ in {1,3,5,10,100} $的标准卡方分布的cdfs的图形。
我知道标准化的卡方分布由
给出
$$ T = \ frac {\ chi ^ {2} _ {df} -df} {\ sqrt {2df}} $$
我尝试使用lines()
在一个图上绘制卡方分布,
z_score <- seq(-3,3,by=0.1)
plot((pchisq(z_score,10)-10)/(sqrt(2*10)),type = "l")
lines((pchisq(z_score,5)-5)/(sqrt(2*5)))
lines((pchisq(z_score,100)-100)/(sqrt(2*100)))
lines((pchisq(z_score,3)-3)/(sqrt(2*3)))
lines((pchisq(z_score,1)-1)/(sqrt(2*1)))
但是这些线并没有加入。我认为不是在轴上有不同限制的情况。你能帮我吗?
解决方法
是的,是轴上的限制不同。查看绘图上的y轴。从-2.236到-2.232。现在,用第一个lines
语句查看您要绘制的数据。
summary((pchisq(z_score,5)-5)/(sqrt(2*5)))
Min. 1st Qu. Median Mean 3rd Qu. Max.
-1.581 -1.581 -1.581 -1.564 -1.554 -1.486
该线的y值都不会出现在您的第一幅图中所示的范围内。一切都会在屏幕外显示。
,这确实是一件大事:我随机尝试了此命令
plot((pchisq(z_score,10)-10)/(sqrt(2*10)),type = "l",ylim=c(-3,0))
并获得多行。然后,您将需要进行一些摆弄才能找到最佳选择。