R version 2.12.2 (2011-02-25) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-redhat-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(814,1150,1225,1691,1759,1754,2100,2062,2012,1897,1964,2186,966,1549,1538,1612,2078,2137,2907,2249,1883,1739,1828,1868,1138,1430,1809,1763,2200,2067,2503,2141,2103,1972,2181,2344,970,1199,1718,1683,2025,2051,2439,2353,2230,1852,2147,2286,1007,1665,1642,1518,1831,2207,2822,2393,2306,1785,2047,2171,1212,1335,2011,1860,1954,2152,2835,2224,2182,1992,2389,2724,891,1247,2017,2257,2255,2255,3057,3330,1896,2096,2374,2535,1041,1728,2201,2455,2204,2660,3670,2665,2639,2226,2586,2684,1185,1749,2459,2618,2585,3310,3923) > 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] 103 > (np <- floor(n / par1)) [1] 8 > 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] 9 9 9 9 9 9 9 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 814 966 1138 970 1007 1212 891 1041 1185 [2,] 1150 1549 1430 1199 1665 1335 1247 1728 1749 [3,] 1225 1538 1809 1718 1642 2011 2017 2201 2459 [4,] 1691 1612 1763 1683 1518 1860 2257 2455 2618 [5,] 1759 2078 2200 2025 1831 1954 2255 2204 2585 [6,] 1754 2137 2067 2051 2207 2152 2255 2660 3310 [7,] 2100 2907 2503 2439 2822 2835 3057 3670 3923 [8,] 2062 2249 2141 2353 2393 2224 3330 2665 NA [9,] 2012 1883 2103 2230 2306 2182 1896 2639 NA [10,] 1897 1739 1972 1852 1785 1992 2096 2226 NA [11,] 1964 1828 2181 2147 2047 2389 2374 2586 NA [12,] 2186 1868 2344 2286 2171 2724 2535 2684 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 336 583 292 229 658 123 356 687 564 [2,] 75 -11 379 519 -23 676 770 473 710 [3,] 466 74 -46 -35 -124 -151 240 254 159 [4,] 68 466 437 342 313 94 -2 -251 -33 [5,] -5 59 -133 26 376 198 0 456 725 [6,] 346 770 436 388 615 683 802 1010 613 [7,] -38 -658 -362 -86 -429 -611 273 -1005 NA [8,] -50 -366 -38 -123 -87 -42 -1434 -26 NA [9,] -115 -144 -131 -378 -521 -190 200 -413 NA [10,] 67 89 209 295 262 397 278 360 NA [11,] 222 40 163 139 124 335 161 98 NA [12,] -1220 -730 -1374 -1279 -959 -1833 -1494 -1499 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/wessaorg/rcomp/tmp/1dj6i1301754953.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/wessaorg/rcomp/tmp/2cwly1301754953.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/wessaorg/rcomp/tmp/3r0hy1301754953.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/wessaorg/rcomp/tmp/4oyv11301754953.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] [,10] [,11] [,12] [1,] 814 1150 1225 1518 1759 2051 2100 2062.0 1883.0 1739.0 1828.0 1868.0 [2,] 966 1247 1642 1683 1954 2067 2503 2182.5 1954.0 1818.5 2005.5 2178.5 [3,] 1007 1430 1809 1763 2078 2152 2835 2301.0 2142.5 1934.5 2164.0 2315.0 [4,] 1138 1665 2017 2257 2204 2255 3057 2529.0 2268.0 2044.0 2381.5 2609.5 [5,] 1212 1749 2459 2618 2255 2255 3670 2665.0 2639.0 2226.0 2586.0 2724.0 $n [1] 9 9 9 9 9 9 9 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 916.4133 1209.853 1611.5 1460.693 1946.333 2052.987 2543.227 2107.44 [2,] 1097.5867 1650.147 2006.5 2065.307 2209.667 2251.013 3126.773 2494.56 [,9] [,10] [,11] [,12] [1,] 1967.095 1808.532 1953.961 2074.237 [2,] 2317.905 2060.468 2374.039 2555.763 $out [1] 2585 1754 2660 3310 3923 3330 $group [1] 5 6 6 6 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(814, 966, 1007, 1138, 1212, 1150, 1247, 1430, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/5vtcd1301754953.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] [,11] [,12] [1,] 123 -23 -151 -251 -133 346 -1005.0 -366.0 -521.0 67.0 40.0 -1833.0 [2,] 292 75 -46 -2 0 436 -634.5 -244.5 -395.5 149.0 111.0 -1496.5 [3,] 356 473 74 94 59 615 -395.5 -68.5 -167.0 270.0 150.0 -1326.5 [4,] 583 676 240 342 376 770 -62.0 -40.0 -123.0 327.5 192.5 -1089.5 [5,] 687 770 466 466 725 1010 273.0 -26.0 200.0 397.0 222.0 -730.0 $n [1] 9 9 9 9 9 9 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 202.74 156.4733 -76.62667 -87.17333 -139.0267 439.0933 -715.30672 [2,] 509.26 789.5267 224.62667 275.17333 257.0267 790.9067 -75.69328 [,8] [,9] [,10] [,11] [,12] [1,] -182.73664 -319.22241 170.2873 104.4729 -1553.856 [2,] 45.73664 -14.77759 369.7127 195.5271 -1099.144 $out [1] -1434 335 $group [1] 8 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(123, 292, 356, 583, 687, -23, 75, 473, 676, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/6leqj1301754953.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,] 814 966.0 1430.0 970.0 1007.0 1335.0 1247.0 1728.0 1185 [2,] 1458 1580.5 1786.0 1700.5 1653.5 1907.0 1956.5 2202.5 2104 [3,] 1828 1848.0 2085.0 2038.0 1939.0 2081.5 2255.0 2520.5 2585 [4,] 2037 2107.5 2190.5 2258.0 2256.5 2306.5 2454.5 2662.5 2964 [5,] 2186 2249.0 2503.0 2439.0 2822.0 2835.0 3057.0 2684.0 3923 $n [1] 12 12 12 12 12 12 12 12 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1563.914 1607.632 1900.505 1783.721 1663.968 1899.285 2027.859 2310.691 [2,] 2092.086 2088.368 2269.495 2292.279 2214.032 2263.715 2482.141 2730.309 [,9] [1,] 2071.422 [2,] 3098.578 $out [1] 2907 1138 1212 891 3330 1041 3670 $group [1] 2 3 6 7 7 8 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(814, 1458, 1828, 2037, 2186, 966, 1580.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/7u7ms1301754953.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,] 1450.222 1430.00 1456.000 [2,] 1893.167 1786.00 1878.062 [3,] 2127.688 2110.25 2097.500 [4,] 2318.931 2232.50 2267.062 [5,] 2917.333 2835.00 2780.000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1933.494 1906.598 1920.074 [2,] 2321.881 2313.902 2274.926 $out [1] 1024.889 1007.000 1052.000 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1450.22222222222, 1893.16666666667, 2127.6875, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1dj6i1301754953.ps tmp/1dj6i1301754953.png",intern=TRUE)) character(0) > try(system("convert tmp/2cwly1301754953.ps tmp/2cwly1301754953.png",intern=TRUE)) character(0) > try(system("convert tmp/3r0hy1301754953.ps tmp/3r0hy1301754953.png",intern=TRUE)) character(0) > try(system("convert tmp/4oyv11301754953.ps tmp/4oyv11301754953.png",intern=TRUE)) character(0) > try(system("convert tmp/5vtcd1301754953.ps tmp/5vtcd1301754953.png",intern=TRUE)) character(0) > try(system("convert tmp/6leqj1301754953.ps tmp/6leqj1301754953.png",intern=TRUE)) character(0) > try(system("convert tmp/7u7ms1301754953.ps tmp/7u7ms1301754953.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.490 0.280 1.884