R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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(56190,54300,51362,49802,48088,46696,56586,64148,56449,52538,49359,49583,51050,49610,48321,47692,46243,46248,56381,62329,60673,58393,55742,57135,57961,56571,55615,53494,52623,52820,66825,70695,65660,63238,61741,63642,65521,64006,62728,62438,61109,63422,78094,82030,75892,72431,69194,71171,72545,71503,69624,67407,66103,67466,81088,86781,79964,80407,76589,78083,78000,76431,75461,73739,71988,72929,85785,89261,84012,80924,76588,77546,73054,73430,71093,72202,70872,70452,80506,80400,77613,69056,65321,64018,64767,61099,58329,56396,54656,55259,66912,66631,59907,56274,54045,55792,55499,53216,52259,51257,48150,51125,61046,61022,56742,54485,53862,58228,61951,62874,64013,62937,61897,65267,75228,76161,71480,69070,68293,74685,72664,71965,69238,67738,65187,66170,77309,77134,70957,67749,65081) > par1 = '12' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > par1 <- as.numeric(par1) > (n <- length(x)) [1] 131 > (np <- floor(n / par1)) [1] 10 > arr <- array(NA,dim=c(par1,np+1)) > darr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > dx <- diff(x) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + darr[j,ari[j]] <- dx[i] + if (j == par1) j = 0 + } > ari [1] 11 11 11 11 11 11 11 11 11 11 11 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 56190 51050 57961 65521 72545 78000 73054 64767 55499 61951 72664 [2,] 54300 49610 56571 64006 71503 76431 73430 61099 53216 62874 71965 [3,] 51362 48321 55615 62728 69624 75461 71093 58329 52259 64013 69238 [4,] 49802 47692 53494 62438 67407 73739 72202 56396 51257 62937 67738 [5,] 48088 46243 52623 61109 66103 71988 70872 54656 48150 61897 65187 [6,] 46696 46248 52820 63422 67466 72929 70452 55259 51125 65267 66170 [7,] 56586 56381 66825 78094 81088 85785 80506 66912 61046 75228 77309 [8,] 64148 62329 70695 82030 86781 89261 80400 66631 61022 76161 77134 [9,] 56449 60673 65660 75892 79964 84012 77613 59907 56742 71480 70957 [10,] 52538 58393 63238 72431 80407 80924 69056 56274 54485 69070 67749 [11,] 49359 55742 61741 69194 76589 76588 65321 54045 53862 68293 65081 [12,] 49583 57135 63642 71171 78083 77546 64018 55792 58228 74685 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -1890 -1440 -1390 -1515 -1042 -1569 376 -3668 -2283 923 -699 [2,] -2938 -1289 -956 -1278 -1879 -970 -2337 -2770 -957 1139 -2727 [3,] -1560 -629 -2121 -290 -2217 -1722 1109 -1933 -1002 -1076 -1500 [4,] -1714 -1449 -871 -1329 -1304 -1751 -1330 -1740 -3107 -1040 -2551 [5,] -1392 5 197 2313 1363 941 -420 603 2975 3370 983 [6,] 9890 10133 14005 14672 13622 12856 10054 11653 9921 9961 11139 [7,] 7562 5948 3870 3936 5693 3476 -106 -281 -24 933 -175 [8,] -7699 -1656 -5035 -6138 -6817 -5249 -2787 -6724 -4280 -4681 -6177 [9,] -3911 -2280 -2422 -3461 443 -3088 -8557 -3633 -2257 -2410 -3208 [10,] -3179 -2651 -1497 -3237 -3818 -4336 -3735 -2229 -623 -777 -2668 [11,] 224 1393 1901 1977 1494 958 -1303 1747 4366 6392 NA [12,] 1467 826 1879 1374 -83 -4492 749 -293 3723 -2021 NA > arr.mean <- array(NA,dim=par1) > arr.median <- array(NA,dim=par1) > arr.midrange <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.mean[j] <- mean(arr[j,],na.rm=TRUE) + arr.median[j] <- median(arr[j,],na.rm=TRUE) + arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2 + } > overall.mean <- mean(x) > overall.median <- median(x) > overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2 > postscript(file="/var/www/html/rcomp/tmp/1afex1293185635.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.mean,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/237vi1293185635.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.median,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/337vi1293185635.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.midrange,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4vgul1293185635.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(arr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 51050.0 49610.0 48321 47692.0 46243.0 46248.0 56381.0 61022.0 56449.0 [2,] 57075.5 55435.5 53937 52375.5 50386.5 51972.5 63935.5 65389.5 60290.0 [3,] 64767.0 62874.0 62728 62438.0 61109.0 63422.0 75228.0 76161.0 70957.0 [4,] 72604.5 71734.0 69431 67572.5 65645.0 66818.0 79300.0 81215.0 76752.5 [5,] 78000.0 76431.0 75461 73739.0 71988.0 72929.0 85785.0 89261.0 84012.0 [,10] [,11] [,12] [1,] 52538.0 49359.0 49583 [2,] 57333.5 54893.5 57135 [3,] 67749.0 65081.0 63830 [4,] 70750.5 68743.5 74685 [5,] 80924.0 76589.0 78083 $n [1] 11 11 11 11 11 11 11 11 11 11 11 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 57369.17 55109.59 55346.85 55198.33 53840.03 56349.78 67908.54 68621.92 [2,] 72164.83 70638.41 70109.15 69677.67 68377.97 70494.22 82547.46 83700.08 [,9] [,10] [,11] [,12] [1,] 63114.46 61357.3 58483.03 55061.32 [2,] 78799.54 74140.7 71678.97 72598.68 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(51050, 57075.5, 64767, 72604.5, 78000, 49610, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5vgul1293185635.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(darr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] -2283.0 -2938.0 -2217.0 -1751.0 -1392 9890.0 -281.0 -7699.0 -3911 -4336 [2,] -1729.5 -2532.0 -1827.5 -1745.5 101 10007.5 -65.0 -6450.5 -3547 -3486 [3,] -1440.0 -1289.0 -1500.0 -1449.0 941 11139.0 3476.0 -5249.0 -3088 -2668 [4,] -870.5 -963.5 -815.5 -1316.5 1838 13239.0 4814.5 -4480.5 -2345 -1863 [5,] 376.0 1139.0 -290.0 -871.0 3370 14672.0 7562.0 -1656.0 -2257 -623 [,11] [,12] [1,] 224.0 -2021.0 [2,] 958.0 -293.0 [3,] 1620.5 787.5 [4,] 1977.0 1467.0 [5,] 1977.0 3723.0 $n [1] 11 11 11 11 11 11 11 11 11 11 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1849.217 -2036.2145 -1982.105 -1653.370 113.5142 9599.552 1151.465 [2,] -1030.783 -541.7855 -1017.895 -1244.630 1768.4858 12678.448 5800.535 [,8] [,9] [,10] [,11] [,12] [1,] -6187.484 -3660.618 -3441.178 1111.367 -91.86617 [2,] -4310.516 -2515.382 -1894.822 2129.633 1666.86617 $out [1] -3668 923 1109 -3107 -2551 443 -8557 -1303 4366 6392 -4492 $group [1] 1 1 3 4 4 9 9 11 11 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-2283, -1729.5, -1440, -870.5, 376, -2938, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6vgul1293185635.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(arr) > names(z) <- c(1:np) > (boxplot(z,notch=TRUE,col='grey',xlab='Block Index',ylab='Value',main='Notched Box Plots - Sequential Blocks')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 46696.0 46243.0 52623.0 61109.0 66103.0 71988 64018.0 54045.0 48150.0 [2,] 49471.0 48006.5 54554.5 63075.0 68545.0 74600 69754.0 55525.5 51758.0 [3,] 51950.0 53396.0 59851.0 67357.5 74567.0 77067 71647.5 57362.5 54173.5 [4,] 56319.5 57764.0 64651.0 74161.5 80185.5 82468 75521.5 62933.0 57485.0 [5,] 64148.0 62329.0 70695.0 82030.0 86781.0 89261 80506.0 66912.0 61046.0 [,10] [,11] [1,] 61897.0 65081.0 [2,] 62905.5 66954.0 [3,] 66780.0 69238.0 [4,] 73082.5 72314.5 [5,] 76161.0 77309.0 $n [1] 12 12 12 12 12 12 12 12 12 12 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 48826.35 48945.54 55245.92 62300.87 69257.69 73478.35 69016.9 53983.89 [2,] 55073.65 57846.46 64456.08 72414.13 79876.31 80655.65 74278.1 60741.11 [,9] [,10] [,11] [1,] 51561.38 62138.2 66684.32 [2,] 56785.62 71421.8 71791.68 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(46696, 49471, 51950, 56319.5, 64148, 46243, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/76pco1293185635.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.mean,arr.median,arr.midrange)) > names(z) <- list('mean','median','midrange') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] 58810.55 61109.0 58015.75 [2,] 61052.05 62801.0 60829.00 [3,] 63864.41 64298.5 63813.38 [4,] 67450.64 69353.0 67064.31 [5,] 74235.64 76161.0 73302.25 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 60945.97 61310.09 60969.41 [2,] 66782.85 67286.91 66657.34 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(58810.5454545455, 61052.0454545455, 63864.4090909091, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1afex1293185635.ps tmp/1afex1293185635.png",intern=TRUE)) character(0) > try(system("convert tmp/237vi1293185635.ps tmp/237vi1293185635.png",intern=TRUE)) character(0) > try(system("convert tmp/337vi1293185635.ps tmp/337vi1293185635.png",intern=TRUE)) character(0) > try(system("convert tmp/4vgul1293185635.ps tmp/4vgul1293185635.png",intern=TRUE)) character(0) > try(system("convert tmp/5vgul1293185635.ps tmp/5vgul1293185635.png",intern=TRUE)) character(0) > try(system("convert tmp/6vgul1293185635.ps tmp/6vgul1293185635.png",intern=TRUE)) character(0) > try(system("convert tmp/76pco1293185635.ps tmp/76pco1293185635.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.205 0.950 2.907