R version 2.10.1 (2009-12-14) 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(464,675,703,887,1139,1077,1318,1260,1120,963,996,960,530,883,894,1045,1199,1287,1565,1577,1076,918,1008,1063,544,635,804,980,1018,1064,1404,1286,1104,999,996,1015,615,722,832,977,1270,1437,1520,1708,1151,934,1159,1209,699,830,996,1124,1458,1270,1753,2258,1208,1241,1265,1828,809,997,1164,1205,1538,1513,1378,2083,1357,1536,1526,1376,779,1005,1193,1522,1539,1546,2116,2326,1596,1356,1553,1613,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 = '15' > #'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] 187 > (np <- floor(n / par1)) [1] 12 > 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] 13 13 13 13 13 13 13 12 12 12 12 12 12 12 12 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 464 1045 1404 934 809 1522 2100 1739 970 1518 2835 2096 1185 [2,] 675 1199 1286 1159 997 1539 2062 1828 1199 1831 2224 2374 1749 [3,] 703 1287 1104 1209 1164 1546 2012 1868 1718 2207 2182 2535 2459 [4,] 887 1565 999 699 1205 2116 1897 1138 1683 2822 1992 1041 2618 [5,] 1139 1577 996 830 1538 2326 1964 1430 2025 2393 2389 1728 2585 [6,] 1077 1076 1015 996 1513 1596 2186 1809 2051 2306 2724 2201 3310 [7,] 1318 918 615 1124 1378 1356 966 1763 2439 1785 891 2455 3923 [8,] 1260 1008 722 1458 2083 1553 1549 2200 2353 2047 1247 2204 NA [9,] 1120 1063 832 1270 1357 1613 1538 2067 2230 2171 2017 2660 NA [10,] 963 544 977 1753 1536 814 1612 2503 1852 1212 2257 3670 NA [11,] 996 635 1270 2258 1526 1150 2078 2141 2147 1335 2255 2665 NA [12,] 960 804 1437 1208 1376 1225 2137 2103 2286 2011 2255 2639 NA [13,] 530 980 1520 1241 779 1691 2907 1972 1007 1860 3057 2226 NA [14,] 883 1018 1708 1265 1005 1759 2249 2181 1665 1954 3330 2586 NA [15,] 894 1064 1151 1828 1193 1754 1883 2344 1642 2152 1896 2684 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 211 154 -118 225 188 17 -38 89 229 313 -611 278 564 [2,] 28 88 -182 50 167 7 -50 40 519 376 -42 161 710 [3,] 184 278 -105 -510 41 570 -115 -730 -35 615 -190 -1494 159 [4,] 252 12 -3 131 333 210 67 292 342 -429 397 687 -33 [5,] -62 -501 19 166 -25 -730 222 379 26 -87 335 473 725 [6,] 241 -158 -400 128 -135 -240 -1220 -46 388 -521 -1833 254 613 [7,] -58 90 107 334 705 197 583 437 -86 262 356 -251 NA [8,] -140 55 110 -188 -726 60 -11 -133 -123 124 770 456 NA [9,] -157 -519 145 483 179 -799 74 436 -378 -959 240 1010 NA [10,] 33 91 293 505 -10 336 466 -362 295 123 -2 -1005 NA [11,] -36 169 167 -1050 -150 75 59 -38 139 676 0 -26 NA [12,] -430 176 83 33 -597 466 770 -131 -1279 -151 802 -413 NA [13,] 353 38 188 24 226 68 -658 209 658 94 273 360 NA [14,] 11 46 -557 563 188 -5 -366 163 -23 198 -1434 98 NA [15,] 151 340 -217 -1019 329 346 -144 -1374 -124 683 200 -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/rcomp/tmp/1qoo61272314828.ps",horizontal=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/rcomp/tmp/2qoo61272314828.ps",horizontal=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/rcomp/tmp/3qoo61272314828.ps",horizontal=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/rcomp/tmp/4b7841272314829.ps",horizontal=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,] 464 675 703 699 830 996 615 722.0 832.0 544.0 635 804.0 [2,] 970 1199 1209 1041 1430 1077 966 1253.5 1195.0 970.0 1210 1216.5 [3,] 1404 1539 1718 1565 1728 1809 1356 1551.0 1575.5 1574.0 1802 1724.0 [4,] 1739 1831 2182 1992 2326 2201 1785 2141.5 2119.0 2054.5 2201 2196.0 [5,] 2835 2374 2535 2822 2585 3310 2455 2353.0 2660.0 3670.0 2665 2639.0 [,13] [,14] [,15] [1,] 530.0 883.0 894 [2,] 993.5 1141.5 1172 [3,] 1605.5 1733.5 1791 [4,] 2099.0 2215.0 2024 [5,] 3057.0 3330.0 2684 $n [1] 13 13 13 13 13 13 13 12 12 12 12 12 12 12 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1067.014 1262.049 1291.619 1148.259 1335.361 1316.448 997.1034 1145.977 [2,] 1740.986 1815.951 2144.381 1981.741 2120.639 2301.552 1714.8966 1956.023 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [1,] 1154.057 1079.352 1349.998 1277.243 1101.274 1243.869 1402.397 [2,] 1996.943 2068.648 2254.002 2170.757 2109.726 2223.131 2179.603 $out [1] 3923 $group [1] 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" Warning message: In bxp(list(stats = c(464, 970, 1404, 1739, 2835, 675, 1199, 1539, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5b7841272314829.ps",horizontal=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,] -118 -182 -730 -429 -501 -1220 -251.0 -188.0 -959.0 -362.0 -150.0 -1279.0 [2,] 17 7 -190 12 -62 -400 16.0 -136.5 -448.5 -6.0 -37.0 -421.5 [3,] 188 50 -35 210 26 -135 229.5 22.0 109.5 107.0 29.5 -49.0 [4,] 229 167 184 333 335 241 396.5 117.0 338.0 315.5 153.0 321.0 [5,] 313 376 615 687 725 613 705.0 456.0 1010.0 505.0 169.0 802.0 [,13] [,14] [,15] [1,] 24.0 -557.0 -1499.0 [2,] 53.0 -194.5 -618.0 [3,] 198.5 28.5 13.5 [4,] 313.0 175.5 334.5 [5,] 658.0 563.0 683.0 $n [1] 13 13 13 13 13 13 12 12 12 12 12 12 12 12 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 95.09881 -20.11410 -198.8917 69.33358 -147.9706 -415.8946 55.9514 [2,] 280.90119 120.11410 128.8917 350.66642 199.9706 145.8946 403.0486 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] -93.62305 -249.2279 -39.63831 -57.16028 -387.6592 79.91225 -140.2595 [2,] 137.62305 468.2279 253.63831 116.16028 289.6592 317.08775 197.2595 [,15] [1,] -420.9416 [2,] 447.9416 $out [1] -611 564 519 710 -1494 -730 -1833 -726 770 -1005 -1050 676 [13] -658 -1434 $group [1] 1 1 2 2 3 5 6 8 8 10 11 11 13 14 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" Warning message: In bxp(list(stats = c(-118, 17, 188, 229, 313, -182, 7, 50, 167, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6mgpp1272314829.ps",horizontal=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] [,10] [,11] [1,] 464.0 804.0 615.0 699 779.0 1150.0 1538.0 1430 970.0 1335 1247.0 [2,] 793.0 949.0 986.5 1060 1084.5 1439.0 1747.5 1786 1653.5 1808 2004.5 [3,] 960.0 1045.0 1104.0 1209 1357.0 1553.0 2012.0 1972 1852.0 2011 2255.0 [4,] 1098.5 1137.5 1345.0 1364 1519.5 1722.5 2118.5 2161 2188.5 2189 2556.5 [5,] 1318.0 1287.0 1708.0 1753 2083.0 2116.0 2249.0 2503 2439.0 2393 3330.0 [,12] [,13] [1,] 1728.0 1185 [2,] 2202.5 2104 [3,] 2455.0 2585 [4,] 2649.5 2964 [5,] 2684.0 3923 $n [1] 15 15 15 15 15 15 15 15 15 15 15 15 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 835.37 968.1006 957.7484 1084.982 1179.54 1437.345 1860.649 1819.017 [2,] 1084.63 1121.8994 1250.2516 1333.018 1534.46 1668.655 2163.351 2124.983 [,9] [,10] [,11] [,12] [,13] [1,] 1633.744 1855.569 2029.809 2272.644 2071.422 [2,] 2070.256 2166.431 2480.191 2637.356 3098.578 $out [1] 1565 1577 544 635 2258 1828 2326 814 966 2907 1138 2822 1212 891 1041 [16] 3670 $group [1] 2 2 2 2 4 4 6 6 7 7 8 10 10 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" NA Warning message: In bxp(list(stats = c(464, 793, 960, 1098.5, 1318, 804, 949, 1045, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7mgpp1272314829.ps",horizontal=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,] 1547.846 1356.00 1354.500 [2,] 1625.205 1558.00 1515.750 [3,] 1661.500 1605.50 1639.000 [4,] 1705.875 1730.75 1694.562 [5,] 1800.250 1809.00 1878.000 $n [1] 15 15 15 $conf [,1] [,2] [,3] [1,] 1628.590 1535.026 1566.053 [2,] 1694.410 1675.974 1711.947 $out [1] 1432.385 1835.385 $group [1] 1 1 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1547.84615384615, 1625.20512820513, 1661.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1qoo61272314828.ps tmp/1qoo61272314828.png",intern=TRUE)) character(0) > try(system("convert tmp/2qoo61272314828.ps tmp/2qoo61272314828.png",intern=TRUE)) character(0) > try(system("convert tmp/3qoo61272314828.ps tmp/3qoo61272314828.png",intern=TRUE)) character(0) > try(system("convert tmp/4b7841272314829.ps tmp/4b7841272314829.png",intern=TRUE)) character(0) > try(system("convert tmp/5b7841272314829.ps tmp/5b7841272314829.png",intern=TRUE)) character(0) > try(system("convert tmp/6mgpp1272314829.ps tmp/6mgpp1272314829.png",intern=TRUE)) character(0) > try(system("convert tmp/7mgpp1272314829.ps tmp/7mgpp1272314829.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.580 1.470 2.799