R软件做出得卡方临界值表总是不对，求高手看看

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

想得出卡方临界值表，程序如下：
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个一循环？？？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：临界值表 求高手 临界值 r软件 R软件 卡方分布 统计

1. df <- seq(1, 50, 1)
   
2. p <- seq(0.95, 0.05, -0.1)
   
3. df2 <- rep(df, length(p))
   
4. p2 <- rep(p, each = length(df))
   
5. dd <- matrix(qchisq(p2, df2), nrow=50, ncol=10, byrow=FALSE, dimnames=list(df, p))

复制代码

应该是qchisq，而不是pchisq

qoiqpwqr 发表于 2015-10-25 19:41
应该是qchisq，而不是pchisq

真棒