for语句，chisq分布的临界值

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

> z<-array(0,dim=c(20,10))
> for(i in 1:20)
+   {for(j in seq(0.95,0.05,-0.1))
+       {z[i,j]<-qchisq(1-j,i)}
+    }
> z
      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,]    0    0    0    0    0    0    0    0    0     0
[2,]    0    0    0    0    0    0    0    0    0     0
[3,]    0    0    0    0    0    0    0    0    0     0
[4,]    0    0    0    0    0    0    0    0    0     0
[5,]    0    0    0    0    0    0    0    0    0     0
[6,]    0    0    0    0    0    0    0    0    0     0
[7,]    0    0    0    0    0    0    0    0    0     0
[8,]    0    0    0    0    0    0    0    0    0     0
[9,]    0    0    0    0    0    0    0    0    0     0
[10,]    0    0    0    0    0    0    0    0    0     0
[11,]    0    0    0    0    0    0    0    0    0     0
[12,]    0    0    0    0    0    0    0    0    0     0
[13,]    0    0    0    0    0    0    0    0    0     0
[14,]    0    0    0    0    0    0    0    0    0     0
[15,]    0    0    0    0    0    0    0    0    0     0
[16,]    0    0    0    0    0    0    0    0    0     0
[17,]    0    0    0    0    0    0    0    0    0     0
[18,]    0    0    0    0    0    0    0    0    0     0
[19,]    0    0    0    0    0    0    0    0    0     0
[20,]    0    0    0    0    0    0    0    0    0     0
想用for循环计算卡方临界值得分布表，但是总是出来空矩阵？

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：chisq R语句 For 临界值 for循环

你的j是小数啊

for(j in seq(0.95,0.05,-0.1))
       {z[i,j]<-qchisq(1-j,i)}
j是小数，z[i,j]里j必须整数

万人往LVR 发表于 2015-11-9 17:41
for(j in seq(0.95,0.05,-0.1))
       {z

被自己蠢哭了，然后我把它改成整数之后，卡方分布临界值还是出不来，刚学R，请赐教

1、参考代码如下：
z<-array(0,dim=c(20,10))
for(i in 1:20){
    k=1;
    for(j in seq(0.95,0.05,-0.1)){
        z[i,k]<-qchisq(1-j,i);
        k=k+1
    }
}

2、结果
> z
             [,1]        [,2]       [,3]       [,4]       [,5]       [,6]       [,7]
[1,]  0.00393214  0.03576578  0.1015310  0.2059001  0.3573172  0.5706519  0.8734571
[2,]  0.10258659  0.32503786  0.5753641  0.8615658  1.1956740  1.5970154  2.0996442
[3,]  0.35184632  0.79777144  1.2125329  1.6415756  2.1094665  2.6430053  3.2831125
[4,]  0.71072302  1.36647723  1.9225575  2.4700874  3.0469464  3.6871338  4.4376892
[5,]  1.14547623  1.99381635  2.6746028  3.3251074  3.9959445  4.7277566  5.5730700
[6,]  1.63538289  2.66127318  3.4545988  4.1972695  4.9518766  5.7651993  6.6947609
[7,]  2.16734991  3.35828438  4.2548522  5.0816470  5.9125232  6.7999702  7.8061229
[8,]  2.73263679  4.07819910  5.0706404  5.9752891  6.8766311  7.8325090  8.9093589
[9,]  3.32511284  4.81652384  5.8988259  6.8762679  7.8434163  8.8631658 10.0059962
[10,]  3.94029914  5.57005944  6.7372008  7.7832430  8.8123518  9.8922157 11.0971423
[11,]  4.57481308  6.33643471  7.5841428  8.6952384  9.7830635 10.9198769 12.1836284
[12,]  5.22602949  7.11383528  8.4384188  9.6115174 10.7552746 11.9463244 13.2660971
[13,]  5.89186434  7.90083665  9.2990655 10.5315076 11.7287740 12.9717003 14.3450565
[14,]  6.57063138  8.69629635 10.1653138 11.4547541 12.7033958 13.9961222 15.4209165
[15,]  7.26094393  9.49928160 11.0365377 12.3808878 13.6790071 15.0196876 16.4940135
[16,]  7.96164557 10.30901918 11.9122197 13.3096047 14.6554995 16.0424793 17.5646273
[17,]  8.67176020 11.12485971 12.7919264 14.2406507 15.6327830 17.0645673 18.6329938
[18,]  9.39045508 11.94625154 13.6752904 15.1738110 16.6107822 18.0860121 19.6993136
[19,] 10.11701306 12.77272138 14.5619967 16.1089018 17.5894329 19.1068661 20.7637596
[20,] 10.85081139 13.60385954 15.4517735 17.0457646 18.5686799 20.1271750 21.8264814
           [,8]      [,9]     [,10]
[1,]  1.323304  2.072251  3.841459
[2,]  2.772589  3.794240  5.991465
[3,]  4.108345  5.317048  7.814728
[4,]  5.385269  6.744883  9.487729
[5,]  6.625680  8.115199 11.070498
[6,]  7.840804  9.446103 12.591587
[7,]  9.037148 10.747895 14.067140
[8,] 10.218855 12.027074 15.507313
[9,] 11.388751 13.288040 16.918978
[10,] 12.548861 14.533936 18.307038
[11,] 13.700693 15.767095 19.675138
[12,] 14.845404 16.989307 21.026070
[13,] 15.983906 18.201977 22.362032
[14,] 17.116934 19.406236 23.684791
[15,] 18.245086 20.603008 24.995790
[16,] 19.368860 21.793057 26.296228
[17,] 20.488676 22.977029 27.587112
[18,] 21.604890 24.155471 28.869299
[19,] 22.717807 25.328852 30.143527
[20,] 23.827692 26.497580 31.410433

出来的z是什么东西看不懂

hongweigg 发表于 2015-11-10 17:27
1、参考代码如下：
z

棒棒哒

万人往LVR 发表于 2015-11-10 17:41
出来的z是什么东西看不懂

卡方分布临界值表