Hi,
I used the commands below to make Hartley's table,
but some values are NA.
require(SuppDists)
trat = seq(2, 15, 1)
gl = seq(2, 40, 1)
har = matrix(0, nr=length(gl), nc=length(trat))
for(i in 1:length(gl))
for(j in 1:length(trat))
har[i,j] <- qmaxFratio(.95, df=gl[i], k=trat[j])
rownames(har) <- gl
colnames(har) <- trat
head(har)
The output (head):
2 3 4 5 6
7
2 39.000000 87.488567 142.502114 202.378010 266.16460
333.18651
3 15.439182 27.758294 39.503063 50.885084 46.75297
72.83358
4 9.604530 15.457596 20.559223 25.211423 29.54387
33.62982
5 7.146382 10.751785 13.723953 NA 10.14968
20.87902
6 5.819757 8.362843 10.380280 12.108103 13.64260
15.03555
7 4.994909 6.939901 8.439993 9.697305 10.80480
11.79569
8 9 10 11 12
13
2 403.07945 475.372389 549.84302 626.22781 704.41272
784.22483
3 83.47794 93.943236 104.24551 114.40008 124.41965
80.84930
4 37.51656 41.237201 44.81356 48.26765 51.61274
54.86066
5 NA 2.173131 26.64526 NA NA
31.62688
6 16.31888 17.514001 18.63614 19.69659 20.70403
21.66528
7 12.69999 13.535307 14.31403 15.04529 15.73600
16.39159
14 15
2 865.53893 948.24795
3 82.18823 83.40749
4 58.02112 61.10224
5 19.09581 34.64536
6 22.58582 23.46946
7 17.01639 17.61392
What's wrong?
Thanks.
--------------------------------------
Silvano Cesar da Costa
Departamento de EstatÃstica
Universidade Estadual de Londrina
Fone: 3371-4346
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
