想得出卡方临界值表,程序如下:
df<-seq(1,50,1)
p<-seq(0.95,0.05,-0.1)
dd<-matrix(pchisq(p,df),nrow=50,ncol=10,byrow=TRUE,dimnames=list(df,p))
dd
得出的结果:
0.95 0.85 0.75 0.65 0.55 0.45 0.35 0.25 0.15
1 6.702807e-01 3.462302e-01 1.386149e-01 4.265126e-02 9.822070e-03 1.605190e-03 1.683014e-04 9.206414e-06 1.557030e-07
2 3.880745e-05 5.693874e-06 6.581850e-07 5.721815e-08 3.488211e-09 1.334106e-10 2.637966e-12 1.834834e-14 1.695448e-17
3 2.196457e-11 1.387277e-12 6.538425e-14 2.148289e-15 4.446509e-17 4.936338e-19 2.222357e-21 2.321012e-24 1.358584e-28
4 1.202148e-18 3.630845e-20 7.680696e-22 1.041618e-23 7.937817e-26 2.756758e-28 3.028690e-31 5.065600e-35 1.994429e-40
5 1.351986e-26 2.049807e-28 2.037126e-30 1.190415e-32 3.478515e-35 3.927140e-38 1.091859e-41 3.027083e-46 8.283304e-53
6 6.702807e-01 3.462302e-01 1.386149e-01 4.265126e-02 9.822070e-03 1.605190e-03 1.683014e-04 9.206414e-06 1.557030e-07
7 3.880745e-05 5.693874e-06 6.581850e-07 5.721815e-08 3.488211e-09 1.334106e-10 2.637966e-12 1.834834e-14 1.695448e-17
8 2.196457e-11 1.387277e-12 6.538425e-14 2.148289e-15 4.446509e-17 4.936338e-19 2.222357e-21 2.321012e-24 1.358584e-28
9 1.202148e-18 3.630845e-20 7.680696e-22 1.041618e-23 7.937817e-26 2.756758e-28 3.028690e-31 5.065600e-35 1.994429e-40
10 1.351986e-26 2.049807e-28 2.037126e-30 1.190415e-32 3.478515e-35 3.927140e-38 1.091859e-41 3.027083e-46 8.283304e-53
11 6.702807e-01 3.462302e-01 1.386149e-01 4.265126e-02 9.822070e-03 1.605190e-03 1.683014e-04 9.206414e-06 1.557030e-07
12 3.880745e-05 5.693874e-06 6.581850e-07 5.721815e-08 3.488211e-09 1.334106e-10 2.637966e-12 1.834834e-14 1.695448e-17
13 2.196457e-11 1.387277e-12 6.538425e-14 2.148289e-15 4.446509e-17 4.936338e-19 2.222357e-21 2.321012e-24 1.358584e-28
14 1.202148e-18 3.630845e-20 7.680696e-22 1.041618e-23 7.937817e-26 2.756758e-28 3.028690e-31 5.065600e-35 1.994429e-40
15 1.351986e-26 2.049807e-28 2.037126e-30 1.190415e-32 3.478515e-35 3.927140e-38 1.091859e-41 3.027083e-46 8.283304e-53
临界值随着自由度的增大反而变小,并且出现5个一循环???