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