R version 3.3.2 (2016-10-31) -- "Sincere Pumpkin Patch"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- c(4998251
+ ,5003638
+ ,5009320
+ ,5016883
+ ,5022453
+ ,5029346
+ ,5035808
+ ,5044609
+ ,5053333
+ ,5061445
+ ,5068136
+ ,5077099
+ ,5083759
+ ,5089916
+ ,5094249
+ ,5098873
+ ,5103289
+ ,5108045
+ ,5113586
+ ,5120164
+ ,5124961
+ ,5131173
+ ,5135907
+ ,5140142
+ ,5145772
+ ,5149446
+ ,5153500
+ ,5157738
+ ,5161326
+ ,5165704
+ ,5169562
+ ,5174322
+ ,5179859
+ ,5184968
+ ,5189502
+ ,5193158
+ ,5198033
+ ,5202300
+ ,5205826
+ ,5209858
+ ,5214325
+ ,5219325
+ ,5224464
+ ,5229458
+ ,5235289
+ ,5241026
+ ,5245792
+ ,5250588
+ ,5256108
+ ,5260122
+ ,5264764
+ ,5269258
+ ,5274105
+ ,5279421
+ ,5285329
+ ,5291269
+ ,5297320
+ ,5303155
+ ,5308602
+ ,5314221
+ ,5320759
+ ,5325517
+ ,5330531
+ ,5335960
+ ,5341703
+ ,5347862
+ ,5354070
+ ,5360839
+ ,5368096
+ ,5373789
+ ,5378892
+ ,5384250
+ ,5390751
+ ,5395677
+ ,5400639
+ ,5406138
+ ,5411535
+ ,5418271
+ ,5424410
+ ,5430256
+ ,5436765
+ ,5442767
+ ,5448135
+ ,5453327
+ ,5459578
+ ,5463642
+ ,5467841
+ ,5473170
+ ,5478905
+ ,5485402
+ ,5491730
+ ,5497742
+ ,5504218
+ ,5509942
+ ,5515659
+ ,5521800
+ ,5529065
+ ,5536374
+ ,5544361
+ ,5552806
+ ,5559276
+ ,5566782
+ ,5573459
+ ,5579345
+ ,5585850
+ ,5592246
+ ,5597749
+ ,5608254
+ ,5618540
+ ,5623505
+ ,5628726
+ ,5634071
+ ,5639470
+ ,5645747
+ ,5651625
+ ,5658077
+ ,5664201
+ ,5670233
+ ,5675532
+ ,5680696
+ ,5686152
+ ,5690973
+ ,5695251
+ ,5700230
+ ,5705336
+ ,5711019
+ ,5717009
+ ,5722743
+ ,5729432
+ ,5736349
+ ,5742667
+ ,5748162
+ ,5754229
+ ,5759749
+ ,5765715
+ ,5773280
+ ,5780631
+ ,5788252
+ ,5795551
+ ,5802529
+ ,5809982
+ ,5817405
+ ,5824530
+ ,5830936
+ ,5838368
+ ,5844289
+ ,5850301
+ ,5856373
+ ,5862425
+ ,5870063
+ ,5876855
+ ,5883791
+ ,5891806
+ ,5900621
+ ,5909153
+ ,5916920
+ ,5923811
+ ,5929343
+ ,5934355
+ ,5941185
+ ,5948326
+ ,5955997
+ ,5963128
+ ,5965395
+ ,5972699
+ ,5979236
+ ,5988217
+ ,5994664
+ ,6001376
+ ,6007398
+ ,6014003
+ ,6020980
+ ,6028883
+ ,6036838
+ ,6044315
+ ,6051848
+ ,6059199
+ ,6067500
+ ,6075436
+ ,6082476
+ ,6090529
+ ,6096378
+ ,6103018
+ ,6110289
+ ,6117855
+ ,6125668
+ ,6132968
+ ,6140802
+ ,6148178
+ ,6156148
+ ,6163139
+ ,6170003
+ ,6178283
+ ,6183865
+ ,6189544
+ ,6196260
+ ,6203691
+ ,6211761
+ ,6219028
+ ,6226687
+ ,6234498
+ ,6242595
+ ,6249692
+ ,6255524
+ ,6262052
+ ,6265163
+ ,6269186
+ ,6273813
+ ,6278669
+ ,6284014
+ ,6288870
+ ,6293425
+ ,6299089
+ ,6303822
+ ,6308440
+ ,6312739
+ ,6317025
+ ,6320023
+ ,6323941
+ ,6329519
+ ,6335055
+ ,6341317
+ ,6347485
+ ,6353486
+ ,6359869
+ ,6366546
+ ,6372165
+ ,6377357
+ ,6382960
+ ,6387765
+ ,6393026
+ ,6399611
+ ,6406939
+ ,6414767
+ ,6421632
+ ,6428999
+ ,6436955
+ ,6444694
+ ,6451896
+ ,6458331
+ ,6465081
+ ,6469425
+ ,6475394
+ ,6482235
+ ,6489759
+ ,6497267
+ ,6504287
+ ,6511602
+ ,6519176
+ ,6527511
+ ,6534732
+ ,6541556
+ ,6549696
+ ,6554493
+ ,6560234
+ ,6566772
+ ,6573430
+ ,6580327
+ ,6586595
+ ,6593091
+ ,6599959
+ ,6606596
+ ,6612326
+ ,6617607
+ ,6624210
+ ,6626698
+ ,6630622
+ ,6635693
+ ,6640927
+ ,6646065
+ ,6650526
+ ,6654668
+ ,6660142
+ ,6664417
+ ,6668643
+ ,6671891
+ ,6675531
+ ,6677500
+ ,6680676
+ ,6685001
+ ,6689345
+ ,6694213
+ ,6698408
+ ,6702566
+ ,6707215
+ ,6711634
+ ,6715311
+ ,6718863
+ ,6722096
+ ,6724545
+ ,6727567
+ ,6731115
+ ,6735403
+ ,6740211
+ ,6744302
+ ,6748748
+ ,6754000
+ ,6759471
+ ,6764308
+ ,6767928
+ ,6771613
+ ,6774802
+ ,6777915
+ ,6781924
+ ,6786926
+ ,6792619
+ ,6797827
+ ,6804797
+ ,6811494
+ ,6819245
+ ,6826377
+ ,6832912
+ ,6836875
+ ,6839328
+ ,6841605
+ ,6843807
+ ,6846217
+ ,6849397
+ ,6852167
+ ,6855609
+ ,6859834
+ ,6863597
+ ,6866921
+ ,6869566
+ ,6871547
+ ,6873441
+ ,6875711
+ ,6878545
+ ,6881123
+ ,6884437
+ ,6887801
+ ,6891378
+ ,6895674
+ ,6899202
+ ,6902272
+ ,6904886
+ ,6907216
+ ,6908830
+ ,6910103
+ ,6913296
+ ,6916712
+ ,6920576
+ ,6923880
+ ,6927498
+ ,6931958
+ ,6935870
+ ,6939141
+ ,6942733
+ ,6945442
+ ,6948319
+ ,6951045
+ ,6954314
+ ,6958044
+ ,6962174
+ ,6965763
+ ,6969974
+ ,6975150
+ ,6980183
+ ,6985916
+ ,6990649
+ ,6994280
+ ,6998216
+ ,7001695
+ ,7005501
+ ,7009834
+ ,7014339
+ ,7018458
+ ,7023545
+ ,7029430
+ ,7035913
+ ,7040815
+ ,7044756
+ ,7048474
+ ,7050785
+ ,7053385
+ ,7056098
+ ,7059113
+ ,7061977
+ ,7064453
+ ,7067642
+ ,7071008
+ ,7074430
+ ,7077063
+ ,7079442
+ ,7081566
+ ,7082757
+ ,7083916
+ ,7085386
+ ,7087290
+ ,7089240
+ ,7090953
+ ,7092800
+ ,7095383
+ ,7098158
+ ,7099909
+ ,7101273
+ ,7102598
+ ,7103151
+ ,7104262
+ ,7105641
+ ,7107536
+ ,7109144
+ ,7111059
+ ,7113153
+ ,7116139
+ ,7118675
+ ,7120622
+ ,7122394
+ ,7124152
+ ,7125698
+ ,7127274
+ ,7128715
+ ,7130356
+ ,7132625
+ ,7134972
+ ,7137715
+ ,7141000
+ ,7143928
+ ,7146246
+ ,7148302
+ ,7149620
+ ,7151061
+ ,7153065
+ ,7155098
+ ,7157836
+ ,7160907
+ ,7163556
+ ,7167284
+ ,7171777
+ ,7176146
+ ,7179890
+ ,7182462
+ ,7184538
+ ,7186261
+ ,7187760
+ ,7189697
+ ,7193421
+ ,7197055
+ ,7200225
+ ,7204328
+ ,7208975
+ ,7213626
+ ,7218168
+ ,7221063
+ ,7224323
+ ,7227592
+ ,7230954
+ ,7234928
+ ,7238879
+ ,7243218
+ ,7247345
+ ,7252076
+ ,7257372
+ ,7263090
+ ,7267429
+ ,7270582
+ ,7273611
+ ,7276980
+ ,7280264
+ ,7283913
+ ,7287325
+ ,7290758
+ ,7294247
+ ,7297687
+ ,7302749
+ ,7307804
+ ,7310947
+ ,7314118
+ ,7316590
+ ,7318927
+ ,7321089
+ ,7324129
+ ,7327626
+ ,7331125
+ ,7334324
+ ,7338334
+ ,7343700
+ ,7348571
+ ,7352525
+ ,7356433
+ ,7358482
+ ,7361911
+ ,7365870
+ ,7370596
+ ,7375100
+ ,7380482
+ ,7385235
+ ,7390710
+ ,7397473
+ ,7404111
+ ,7410271
+ ,7415460
+ ,7419501
+ ,7424168
+ ,7428191
+ ,7432789
+ ,7437234
+ ,7442263
+ ,7447314
+ ,7452727
+ ,7459611
+ ,7466251
+ ,7472177
+ ,7476635
+ ,7480422
+ ,7485033
+ ,7488669
+ ,7492766
+ ,7496579
+ ,7500777
+ ,7504798
+ ,7509683
+ ,7515290
+ ,7521531
+ ,7526647
+ ,7530625
+ ,7535268
+ ,7538600
+ ,7542228
+ ,7545926
+ ,7550077
+ ,7554598
+ ,7558978
+ ,7564670
+ ,7571143
+ ,7577108
+ ,7580944
+ ,7583477
+ ,7585887
+ ,7588559
+ ,7591214
+ ,7595207
+ ,7599164
+ ,7602596
+ ,7605835
+ ,7608298
+ ,7612572
+ ,7616928
+ ,7620390
+ ,7624558
+ ,7627482
+ ,7630502
+ ,7633494
+ ,7636360
+ ,7638508
+ ,7641578
+ ,7644806
+ ,7647289
+ ,7651132
+ ,7654750
+ ,7658292
+ ,7660970
+ ,7662289
+ ,7663984
+ ,7665441
+ ,7667100
+ ,7668938
+ ,7671672
+ ,7674587
+ ,7678101
+ ,7682904
+ ,7686986
+ ,7691383
+ ,7694849
+ ,7696803
+ ,7699332
+ ,7701907
+ ,7704693
+ ,7708286
+ ,7711976
+ ,7715648
+ ,7720023
+ ,7724353
+ ,7729836
+ ,7734046
+ ,7737380
+ ,7740074
+ ,7744788
+ ,7748121
+ ,7751042
+ ,7753878
+ ,7757491
+ ,7761290
+ ,7765877
+ ,7772444
+ ,7779463
+ ,7785438
+ ,7790248
+ ,7793271
+ ,7796751
+ ,7800010
+ ,7803554
+ ,7807762
+ ,7811726
+ ,7816076
+ ,7821014
+ ,7827255
+ ,7833731
+ ,7838330
+ ,7842800
+ ,7846317
+ ,7849498
+ ,7854534
+ ,7859444
+ ,7863737
+ ,7869194
+ ,7874109
+ ,7879691
+ ,7887162
+ ,7894275
+ ,7900310
+ ,7906193
+ ,7909855
+ ,7914160
+ ,7918742
+ ,7922820
+ ,7927200
+ ,7932429
+ ,7937106
+ ,7942595
+ ,7949467
+ ,7956241
+ ,7962056
+ ,7966928
+ ,7971967
+ ,7975761
+ ,7978782
+ ,7981367
+ ,7984421
+ ,7987694
+ ,7990541
+ ,7994469
+ ,7999834
+ ,8005331
+ ,8009706
+ ,8012422
+ ,8015471
+ ,8017853
+ ,8020041
+ ,8021677
+ ,8023645
+ ,8026355
+ ,8028301
+ ,8030053
+ ,8033515
+ ,8037945
+ ,8040861
+ ,8042317
+ ,8045914
+ ,8046294
+ ,8046948
+ ,8047395
+ ,8047837
+ ,8049443
+ ,8051361
+ ,8052736
+ ,8055118
+ ,8058952
+ ,8060720
+ ,8062433
+ ,8065979
+ ,8066467
+ ,8066688
+ ,8066257
+ ,8066714
+ ,8067547
+ ,8068098
+ ,8069174
+ ,8071058
+ ,8074816
+ ,8076190
+ ,8076510
+ ,8077407
+ ,8077602
+ ,8077916
+ ,8077185
+ ,8076868
+ ,8077328
+ ,8078300
+ ,8078531
+ ,8080989
+ ,8084190
+ ,8085995
+ ,8086086
+ ,8088514
+ ,8089088
+ ,8090979
+ ,8092041
+ ,8093302
+ ,8095132
+ ,8096756
+ ,8098149
+ ,8101094
+ ,8106329
+ ,8108958
+ ,8109892
+ ,8112073
+ ,8114285
+ ,8117237
+ ,8119470
+ ,8122461
+ ,8125487
+ ,8128957
+ ,8132066
+ ,8136720
+ ,8144900
+ ,8149659
+ ,8152801
+ ,8156396
+ ,8158251
+ ,8161557
+ ,8164994
+ ,8167762
+ ,8171387
+ ,8175119
+ ,8179193
+ ,8184457
+ ,8192647
+ ,8196713
+ ,8200018
+ ,8203476
+ ,8205586
+ ,8209197
+ ,8210899
+ ,8213427
+ ,8216286
+ ,8219951
+ ,8222904
+ ,8228652
+ ,8235192
+ ,8239730
+ ,8242556
+ ,8243482
+ ,8246483
+ ,8250539
+ ,8253158
+ ,8254803
+ ,8257440
+ ,8260524
+ ,8262642
+ ,8269045
+ ,8276541
+ ,8279969
+ ,8282336
+ ,8282871
+ ,8283273
+ ,8285849
+ ,8286974
+ ,8288019
+ ,8289646
+ ,8290866
+ ,8291787
+ ,8297340
+ ,8304056
+ ,8306959
+ ,8308098
+ ,8307339
+ ,8308772
+ ,8310421
+ ,8311363
+ ,8311926
+ ,8313781
+ ,8315731
+ ,8317691
+ ,8323013
+ ,8329897
+ ,8332639
+ ,8334485
+ ,8334385
+ ,8337055
+ ,8340482
+ ,8343066
+ ,8345517
+ ,8347964
+ ,8350453
+ ,8351952
+ ,8357487
+ ,8365236
+ ,8369371
+ ,8372420
+ ,8372858
+ ,8374861
+ ,8377729
+ ,8379637
+ ,8381851
+ ,8384251
+ ,8386614
+ ,8388680
+ ,8395572
+ ,8404511
+ ,8409838
+ ,8414637
+ ,8417135
+ ,8420985
+ ,8426542
+ ,8431202
+ ,8435176
+ ,8440119
+ ,8445299)
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa P., (2012), Variability (v1.0.5) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_variability.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> num <- 50
> res <- array(NA,dim=c(num,3))
> q1 <- function(data,n,p,i,f) {
+ np <- n*p;
+ i <<- floor(np)
+ f <<- np - i
+ qvalue <- (1-f)*data[i] + f*data[i+1]
+ }
> q2 <- function(data,n,p,i,f) {
+ np <- (n+1)*p
+ i <<- floor(np)
+ f <<- np - i
+ qvalue <- (1-f)*data[i] + f*data[i+1]
+ }
> q3 <- function(data,n,p,i,f) {
+ np <- n*p
+ i <<- floor(np)
+ f <<- np - i
+ if (f==0) {
+ qvalue <- data[i]
+ } else {
+ qvalue <- data[i+1]
+ }
+ }
> q4 <- function(data,n,p,i,f) {
+ np <- n*p
+ i <<- floor(np)
+ f <<- np - i
+ if (f==0) {
+ qvalue <- (data[i]+data[i+1])/2
+ } else {
+ qvalue <- data[i+1]
+ }
+ }
> q5 <- function(data,n,p,i,f) {
+ np <- (n-1)*p
+ i <<- floor(np)
+ f <<- np - i
+ if (f==0) {
+ qvalue <- data[i+1]
+ } else {
+ qvalue <- data[i+1] + f*(data[i+2]-data[i+1])
+ }
+ }
> q6 <- function(data,n,p,i,f) {
+ np <- n*p+0.5
+ i <<- floor(np)
+ f <<- np - i
+ qvalue <- data[i]
+ }
> q7 <- function(data,n,p,i,f) {
+ np <- (n+1)*p
+ i <<- floor(np)
+ f <<- np - i
+ if (f==0) {
+ qvalue <- data[i]
+ } else {
+ qvalue <- f*data[i] + (1-f)*data[i+1]
+ }
+ }
> q8 <- function(data,n,p,i,f) {
+ np <- (n+1)*p
+ i <<- floor(np)
+ f <<- np - i
+ if (f==0) {
+ qvalue <- data[i]
+ } else {
+ if (f == 0.5) {
+ qvalue <- (data[i]+data[i+1])/2
+ } else {
+ if (f < 0.5) {
+ qvalue <- data[i]
+ } else {
+ qvalue <- data[i+1]
+ }
+ }
+ }
+ }
> iqd <- function(x,def) {
+ x <-sort(x[!is.na(x)])
+ n<-length(x)
+ if (def==1) {
+ qvalue1 <- q1(x,n,0.25,i,f)
+ qvalue3 <- q1(x,n,0.75,i,f)
+ }
+ if (def==2) {
+ qvalue1 <- q2(x,n,0.25,i,f)
+ qvalue3 <- q2(x,n,0.75,i,f)
+ }
+ if (def==3) {
+ qvalue1 <- q3(x,n,0.25,i,f)
+ qvalue3 <- q3(x,n,0.75,i,f)
+ }
+ if (def==4) {
+ qvalue1 <- q4(x,n,0.25,i,f)
+ qvalue3 <- q4(x,n,0.75,i,f)
+ }
+ if (def==5) {
+ qvalue1 <- q5(x,n,0.25,i,f)
+ qvalue3 <- q5(x,n,0.75,i,f)
+ }
+ if (def==6) {
+ qvalue1 <- q6(x,n,0.25,i,f)
+ qvalue3 <- q6(x,n,0.75,i,f)
+ }
+ if (def==7) {
+ qvalue1 <- q7(x,n,0.25,i,f)
+ qvalue3 <- q7(x,n,0.75,i,f)
+ }
+ if (def==8) {
+ qvalue1 <- q8(x,n,0.25,i,f)
+ qvalue3 <- q8(x,n,0.75,i,f)
+ }
+ iqdiff <- qvalue3 - qvalue1
+ return(c(iqdiff,iqdiff/2,iqdiff/(qvalue3 + qvalue1)))
+ }
> range <- max(x) - min(x)
> lx <- length(x)
> biasf <- (lx-1)/lx
> varx <- var(x)
> bvarx <- varx*biasf
> sdx <- sqrt(varx)
> mx <- mean(x)
> bsdx <- sqrt(bvarx)
> x2 <- x*x
> mse0 <- sum(x2)/lx
> xmm <- x-mx
> xmm2 <- xmm*xmm
> msem <- sum(xmm2)/lx
> axmm <- abs(x - mx)
> medx <- median(x)
> axmmed <- abs(x - medx)
> xmmed <- x - medx
> xmmed2 <- xmmed*xmmed
> msemed <- sum(xmmed2)/lx
> qarr <- array(NA,dim=c(8,3))
> for (j in 1:8) {
+ qarr[j,] <- iqd(x,j)
+ }
> sdpo <- 0
> adpo <- 0
> for (i in 1:(lx-1)) {
+ for (j in (i+1):lx) {
+ ldi <- x[i]-x[j]
+ aldi <- abs(ldi)
+ sdpo = sdpo + ldi * ldi
+ adpo = adpo + aldi
+ }
+ }
> denom <- (lx*(lx-1)/2)
> sdpo = sdpo / denom
> adpo = adpo / denom
> gmd <- 0
> for (i in 1:lx) {
+ for (j in 1:lx) {
+ ldi <- abs(x[i]-x[j])
+ gmd = gmd + ldi
+ }
+ }
> gmd <- gmd / (lx*(lx-1))
> sumx <- sum(x)
> pk <- x / sumx
> ck <- cumsum(pk)
> dk <- array(NA,dim=lx)
> for (i in 1:lx) {
+ if (ck[i] <= 0.5) dk[i] <- ck[i] else dk[i] <- 1 - ck[i]
+ }
> bigd <- sum(dk) * 2 / (lx-1)
> iod <- 1 - sum(pk*pk)
> res[1,] <- c('Absolute range','http://www.xycoon.com/absolute.htm', range)
> res[2,] <- c('Relative range (unbiased)','http://www.xycoon.com/relative.htm', range/sd(x))
> res[3,] <- c('Relative range (biased)','http://www.xycoon.com/relative.htm', range/sqrt(varx*biasf))
> res[4,] <- c('Variance (unbiased)','http://www.xycoon.com/unbiased.htm', varx)
> res[5,] <- c('Variance (biased)','http://www.xycoon.com/biased.htm', bvarx)
> res[6,] <- c('Standard Deviation (unbiased)','http://www.xycoon.com/unbiased1.htm', sdx)
> res[7,] <- c('Standard Deviation (biased)','http://www.xycoon.com/biased1.htm', bsdx)
> res[8,] <- c('Coefficient of Variation (unbiased)','http://www.xycoon.com/variation.htm', sdx/mx)
> res[9,] <- c('Coefficient of Variation (biased)','http://www.xycoon.com/variation.htm', bsdx/mx)
> res[10,] <- c('Mean Squared Error (MSE versus 0)','http://www.xycoon.com/mse.htm', mse0)
> res[11,] <- c('Mean Squared Error (MSE versus Mean)','http://www.xycoon.com/mse.htm', msem)
> res[12,] <- c('Mean Absolute Deviation from Mean (MAD Mean)', 'http://www.xycoon.com/mean2.htm', sum(axmm)/lx)
> res[13,] <- c('Mean Absolute Deviation from Median (MAD Median)', 'http://www.xycoon.com/median1.htm', sum(axmmed)/lx)
> res[14,] <- c('Median Absolute Deviation from Mean', 'http://www.xycoon.com/mean3.htm', median(axmm))
> res[15,] <- c('Median Absolute Deviation from Median', 'http://www.xycoon.com/median2.htm', median(axmmed))
> res[16,] <- c('Mean Squared Deviation from Mean', 'http://www.xycoon.com/mean1.htm', msem)
> res[17,] <- c('Mean Squared Deviation from Median', 'http://www.xycoon.com/median.htm', msemed)
>
> #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/rcomp/createtable")
>
> mylink1 <- hyperlink('http://www.xycoon.com/difference.htm','Interquartile Difference','')
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_1.htm','(Weighted Average at Xnp)',''),sep=' ')
> res[18,] <- c('', mylink2, qarr[1,1])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_2.htm','(Weighted Average at X(n+1)p)',''),sep=' ')
> res[19,] <- c('', mylink2, qarr[2,1])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_3.htm','(Empirical Distribution Function)',''),sep=' ')
> res[20,] <- c('', mylink2, qarr[3,1])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_4.htm','(Empirical Distribution Function - Averaging)',''),sep=' ')
> res[21,] <- c('', mylink2, qarr[4,1])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_5.htm','(Empirical Distribution Function - Interpolation)',''),sep=' ')
> res[22,] <- c('', mylink2, qarr[5,1])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_6.htm','(Closest Observation)',''),sep=' ')
> res[23,] <- c('', mylink2, qarr[6,1])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_7.htm','(True Basic - Statistics Graphics Toolkit)',''),sep=' ')
> res[24,] <- c('', mylink2, qarr[7,1])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_8.htm','(MS Excel (old versions))',''),sep=' ')
> res[25,] <- c('', mylink2, qarr[8,1])
> mylink1 <- hyperlink('http://www.xycoon.com/deviation.htm','Semi Interquartile Difference','')
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_1.htm','(Weighted Average at Xnp)',''),sep=' ')
> res[26,] <- c('', mylink2, qarr[1,2])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_2.htm','(Weighted Average at X(n+1)p)',''),sep=' ')
> res[27,] <- c('', mylink2, qarr[2,2])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_3.htm','(Empirical Distribution Function)',''),sep=' ')
> res[28,] <- c('', mylink2, qarr[3,2])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_4.htm','(Empirical Distribution Function - Averaging)',''),sep=' ')
> res[29,] <- c('', mylink2, qarr[4,2])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_5.htm','(Empirical Distribution Function - Interpolation)',''),sep=' ')
> res[30,] <- c('', mylink2, qarr[5,2])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_6.htm','(Closest Observation)',''),sep=' ')
> res[31,] <- c('', mylink2, qarr[6,2])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_7.htm','(True Basic - Statistics Graphics Toolkit)',''),sep=' ')
> res[32,] <- c('', mylink2, qarr[7,2])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_8.htm','(MS Excel (old versions))',''),sep=' ')
> res[33,] <- c('', mylink2, qarr[8,2])
> mylink1 <- hyperlink('http://www.xycoon.com/variation1.htm','Coefficient of Quartile Variation','')
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_1.htm','(Weighted Average at Xnp)',''),sep=' ')
> res[34,] <- c('', mylink2, qarr[1,3])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_2.htm','(Weighted Average at X(n+1)p)',''),sep=' ')
> res[35,] <- c('', mylink2, qarr[2,3])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_3.htm','(Empirical Distribution Function)',''),sep=' ')
> res[36,] <- c('', mylink2, qarr[3,3])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_4.htm','(Empirical Distribution Function - Averaging)',''),sep=' ')
> res[37,] <- c('', mylink2, qarr[4,3])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_5.htm','(Empirical Distribution Function - Interpolation)',''),sep=' ')
> res[38,] <- c('', mylink2, qarr[5,3])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_6.htm','(Closest Observation)',''),sep=' ')
> res[39,] <- c('', mylink2, qarr[6,3])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_7.htm','(True Basic - Statistics Graphics Toolkit)',''),sep=' ')
> res[40,] <- c('', mylink2, qarr[7,3])
> mylink2 <- paste(mylink1,hyperlink('http://www.xycoon.com/method_8.htm','(MS Excel (old versions))',''),sep=' ')
> res[41,] <- c('', mylink2, qarr[8,3])
> res[42,] <- c('Number of all Pairs of Observations', 'http://www.xycoon.com/pair_numbers.htm', lx*(lx-1)/2)
> res[43,] <- c('Squared Differences between all Pairs of Observations', 'http://www.xycoon.com/squared_differences.htm', sdpo)
> res[44,] <- c('Mean Absolute Differences between all Pairs of Observations', 'http://www.xycoon.com/mean_abs_differences.htm', adpo)
> res[45,] <- c('Gini Mean Difference', 'http://www.xycoon.com/gini_mean_difference.htm', gmd)
> res[46,] <- c('Leik Measure of Dispersion', 'http://www.xycoon.com/leiks_d.htm', bigd)
> res[47,] <- c('Index of Diversity', 'http://www.xycoon.com/diversity.htm', iod)
> res[48,] <- c('Index of Qualitative Variation', 'http://www.xycoon.com/qualitative_variation.htm', iod*lx/(lx-1))
> res[49,] <- c('Coefficient of Dispersion', 'http://www.xycoon.com/dispersion.htm', sum(axmm)/lx/medx)
> res[50,] <- c('Observations', '', lx)
> res
[,1]
[1,] "Absolute range"
[2,] "Relative range (unbiased)"
[3,] "Relative range (biased)"
[4,] "Variance (unbiased)"
[5,] "Variance (biased)"
[6,] "Standard Deviation (unbiased)"
[7,] "Standard Deviation (biased)"
[8,] "Coefficient of Variation (unbiased)"
[9,] "Coefficient of Variation (biased)"
[10,] "Mean Squared Error (MSE versus 0)"
[11,] "Mean Squared Error (MSE versus Mean)"
[12,] "Mean Absolute Deviation from Mean (MAD Mean)"
[13,] "Mean Absolute Deviation from Median (MAD Median)"
[14,] "Median Absolute Deviation from Mean"
[15,] "Median Absolute Deviation from Median"
[16,] "Mean Squared Deviation from Mean"
[17,] "Mean Squared Deviation from Median"
[18,] ""
[19,] ""
[20,] ""
[21,] ""
[22,] ""
[23,] ""
[24,] ""
[25,] ""
[26,] ""
[27,] ""
[28,] ""
[29,] ""
[30,] ""
[31,] ""
[32,] ""
[33,] ""
[34,] ""
[35,] ""
[36,] ""
[37,] ""
[38,] ""
[39,] ""
[40,] ""
[41,] ""
[42,] "Number of all Pairs of Observations"
[43,] "Squared Differences between all Pairs of Observations"
[44,] "Mean Absolute Differences between all Pairs of Observations"
[45,] "Gini Mean Difference"
[46,] "Leik Measure of Dispersion"
[47,] "Index of Diversity"
[48,] "Index of Qualitative Variation"
[49,] "Coefficient of Dispersion"
[50,] "Observations"
[,2]
[1,] "http://www.xycoon.com/absolute.htm"
[2,] "http://www.xycoon.com/relative.htm"
[3,] "http://www.xycoon.com/relative.htm"
[4,] "http://www.xycoon.com/unbiased.htm"
[5,] "http://www.xycoon.com/biased.htm"
[6,] "http://www.xycoon.com/unbiased1.htm"
[7,] "http://www.xycoon.com/biased1.htm"
[8,] "http://www.xycoon.com/variation.htm"
[9,] "http://www.xycoon.com/variation.htm"
[10,] "http://www.xycoon.com/mse.htm"
[11,] "http://www.xycoon.com/mse.htm"
[12,] "http://www.xycoon.com/mean2.htm"
[13,] "http://www.xycoon.com/median1.htm"
[14,] "http://www.xycoon.com/mean3.htm"
[15,] "http://www.xycoon.com/median2.htm"
[16,] "http://www.xycoon.com/mean1.htm"
[17,] "http://www.xycoon.com/median.htm"
[18,] "Interquartile Difference (Weighted Average at Xnp)"
[19,] "Interquartile Difference (Weighted Average at X(n+1)p)"
[20,] "Interquartile Difference (Empirical Distribution Function)"
[21,] "Interquartile Difference (Empirical Distribution Function - Averaging)"
[22,] "Interquartile Difference (Empirical Distribution Function - Interpolation)"
[23,] "Interquartile Difference (Closest Observation)"
[24,] "Interquartile Difference (True Basic - Statistics Graphics Toolkit)"
[25,] "Interquartile Difference (MS Excel (old versions))"
[26,] "Semi Interquartile Difference (Weighted Average at Xnp)"
[27,] "Semi Interquartile Difference (Weighted Average at X(n+1)p)"
[28,] "Semi Interquartile Difference (Empirical Distribution Function)"
[29,] "Semi Interquartile Difference (Empirical Distribution Function - Averaging)"
[30,] "Semi Interquartile Difference (Empirical Distribution Function - Interpolation)"
[31,] "Semi Interquartile Difference (Closest Observation)"
[32,] "Semi Interquartile Difference (True Basic - Statistics Graphics Toolkit)"
[33,] "Semi Interquartile Difference (MS Excel (old versions))"
[34,] "Coefficient of Quartile Variation (Weighted Average at Xnp)"
[35,] "Coefficient of Quartile Variation (Weighted Average at X(n+1)p)"
[36,] "Coefficient of Quartile Variation (Empirical Distribution Function)"
[37,] "Coefficient of Quartile Variation (Empirical Distribution Function - Averaging)"
[38,] "Coefficient of Quartile Variation (Empirical Distribution Function - Interpolation)"
[39,] "Coefficient of Quartile Variation (Closest Observation)"
[40,] "Coefficient of Quartile Variation (True Basic - Statistics Graphics Toolkit)"
[41,] "Coefficient of Quartile Variation (MS Excel (old versions))"
[42,] "http://www.xycoon.com/pair_numbers.htm"
[43,] "http://www.xycoon.com/squared_differences.htm"
[44,] "http://www.xycoon.com/mean_abs_differences.htm"
[45,] "http://www.xycoon.com/gini_mean_difference.htm"
[46,] "http://www.xycoon.com/leiks_d.htm"
[47,] "http://www.xycoon.com/diversity.htm"
[48,] "http://www.xycoon.com/qualitative_variation.htm"
[49,] "http://www.xycoon.com/dispersion.htm"
[50,] ""
[,3]
[1,] "3447048"
[2,] "3.49396245472426"
[3,] "3.49615096884723"
[4,] "973325707827.299"
[5,] "972107527967.691"
[6,] "986572.707826088"
[7,] "985955.134865523"
[8,] "0.141378467516415"
[9,] "0.141289967684571"
[10,] "49667924309234.5"
[11,] "972107527967.691"
[12,] "831128.21628099"
[13,] "822833.564455569"
[14,] "818512.197747184"
[15,] "804214"
[16,] "972107527967.691"
[17,] "988338847958.503"
[18,] "1614675.25"
[19,] "1616113"
[20,] "1616113"
[21,] "1616113"
[22,] "1609972.5"
[23,] "1611643"
[24,] "1616113"
[25,] "1616113"
[26,] "807337.625"
[27,] "808056.5"
[28,] "808056.5"
[29,] "808056.5"
[30,] "804986.25"
[31,] "805821.5"
[32,] "808056.5"
[33,] "808056.5"
[34,] "0.114807310942365"
[35,] "0.114866519298109"
[36,] "0.114866519298109"
[37,] "0.114866519298109"
[38,] "0.114416493454785"
[39,] "0.114585215218723"
[40,] "0.114866519298109"
[41,] "0.114866519298109"
[42,] "318801"
[43,] "1946651415656.52"
[44,] "1127545.42012102"
[45,] "1127545.42012102"
[46,] "0.488514747487185"
[47,] "0.99872345074472"
[48,] "0.999974983891017"
[49,] "0.116967380744537"
[50,] "799"
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Variability - Ungrouped Data',2,TRUE)
> a<-table.row.end(a)
> for (i in 1:num) {
+ a<-table.row.start(a)
+ if (res[i,1] != '') {
+ a<-table.element(a,hyperlink(res[i,2],res[i,1],''),header=TRUE)
+ } else {
+ a<-table.element(a,res[i,2],header=TRUE)
+ }
+ a<-table.element(a,res[i,3])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/1l9jl1479581549.tab")
>
>
>
> proc.time()
user system elapsed
2.063 0.029 2.091