R version 2.8.0 (2008-10-20) Copyright (C) 2008 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 = '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] 187 > (np <- floor(n / par1)) [1] 15 > 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] 16 16 16 16 16 16 16 15 15 15 15 15 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 464 530 544 615 699 809 779 814 966 1138 970 1007 1212 [2,] 675 883 635 722 830 997 1005 1150 1549 1430 1199 1665 1335 [3,] 703 894 804 832 996 1164 1193 1225 1538 1809 1718 1642 2011 [4,] 887 1045 980 977 1124 1205 1522 1691 1612 1763 1683 1518 1860 [5,] 1139 1199 1018 1270 1458 1538 1539 1759 2078 2200 2025 1831 1954 [6,] 1077 1287 1064 1437 1270 1513 1546 1754 2137 2067 2051 2207 2152 [7,] 1318 1565 1404 1520 1753 1378 2116 2100 2907 2503 2439 2822 2835 [8,] 1260 1577 1286 1708 2258 2083 2326 2062 2249 2141 2353 2393 2224 [9,] 1120 1076 1104 1151 1208 1357 1596 2012 1883 2103 2230 2306 2182 [10,] 963 918 999 934 1241 1536 1356 1897 1739 1972 1852 1785 1992 [11,] 996 1008 996 1159 1265 1526 1553 1964 1828 2181 2147 2047 2389 [12,] 960 1063 1015 1209 1828 1376 1613 2186 1868 2344 2286 2171 2724 [,14] [,15] [,16] [1,] 891 1041 1185 [2,] 1247 1728 1749 [3,] 2017 2201 2459 [4,] 2257 2455 2618 [5,] 2255 2204 2585 [6,] 2255 2660 3310 [7,] 3057 3670 3923 [8,] 3330 2665 NA [9,] 1896 2639 NA [10,] 2096 2226 NA [11,] 2374 2586 NA [12,] 2535 2684 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 211 353 91 107 131 188 226 336 583 292 229 658 123 [2,] 28 11 169 110 166 167 188 75 -11 379 519 -23 676 [3,] 184 151 176 145 128 41 329 466 74 -46 -35 -124 -151 [4,] 252 154 38 293 334 333 17 68 466 437 342 313 94 [5,] -62 88 46 167 -188 -25 7 -5 59 -133 26 376 198 [6,] 241 278 340 83 483 -135 570 346 770 436 388 615 683 [7,] -58 12 -118 188 505 705 210 -38 -658 -362 -86 -429 -611 [8,] -140 -501 -182 -557 -1050 -726 -730 -50 -366 -38 -123 -87 -42 [9,] -157 -158 -105 -217 33 179 -240 -115 -144 -131 -378 -521 -190 [10,] 33 90 -3 225 24 -10 197 67 89 209 295 262 397 [11,] -36 55 19 50 563 -150 60 222 40 163 139 124 335 [12,] -430 -519 -400 -510 -1019 -597 -799 -1220 -730 -1374 -1279 -959 -1833 [,14] [,15] [,16] [1,] 356 687 564 [2,] 770 473 710 [3,] 240 254 159 [4,] -2 -251 -33 [5,] 0 456 725 [6,] 802 1010 613 [7,] 273 -1005 NA [8,] -1434 -26 NA [9,] 200 -413 NA [10,] 278 360 NA [11,] 161 98 NA [12,] -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/html/freestat/rcomp/tmp/12t6l1272310816.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/html/freestat/rcomp/tmp/22t6l1272310816.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/html/freestat/rcomp/tmp/32t6l1272310816.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/html/freestat/rcomp/tmp/4vkno1272310816.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] [1,] 464.0 635.0 703.0 887.0 1018 1064.0 1318.0 1260.0 1076.0 918.0 996 [2,] 657.0 856.5 945.0 1084.5 1364 1362.0 1542.5 1885.0 1179.5 1120.0 1212 [3,] 852.5 1174.5 1381.5 1567.0 1795 1902.5 2277.5 2224.0 1883.0 1739.0 1828 [4,] 1024.0 1489.5 1910.0 1811.5 2139 2179.5 2871.0 2339.5 2142.5 1934.5 2164 [5,] 1212.0 1749.0 2459.0 2618.0 2585 3310.0 3923.0 2665.0 2639.0 2226.0 2586 [,12] [1,] 960.0 [2,] 1292.5 [3,] 1868.0 [4,] 2315.0 [5,] 2724.0 $n [1] 16 16 16 16 16 16 16 15 15 15 15 15 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 707.535 924.465 1000.325 1279.835 1488.875 1579.588 1752.742 2038.585 [2,] 997.465 1424.535 1762.675 1854.165 2101.125 2225.412 2802.258 2409.415 [,9] [,10] [,11] [,12] [1,] 1490.14 1406.721 1439.628 1450.867 [2,] 2275.86 2071.279 2216.372 2285.133 $out [1] 3330 $group [1] 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(464, 657, 852.5, 1024, 1212, 635, 856.5, 1174.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5vkno1272310816.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] [1,] 91.0 -23.0 -151 -251.0 -188.0 -135.0 -1005.0 -1434.0 -378.0 -10 -36 [2,] 159.5 51.5 3 27.5 -15.0 309.0 -395.5 -641.5 -228.5 50 45 [3,] 260.5 168.0 148 203.0 36.0 459.5 -58.0 -182.0 -157.0 197 98 [4,] 460.0 496.0 212 333.5 182.5 649.0 199.0 -68.5 -110.0 270 162 [5,] 687.0 770.0 466 466.0 456.0 1010.0 705.0 -26.0 33.0 397 335 [,12] [1,] -1833.0 [2,] -1326.5 [3,] -959.0 [4,] -558.0 [5,] -400.0 $n [1] 16 16 16 16 16 16 15 15 15 15 15 15 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 141.8025 -7.5775 65.445 82.13 -42.0125 325.2 -300.5288 -415.75778 [2,] 379.1975 343.5775 230.555 323.87 114.0125 593.8 184.5288 51.75778 [,9] [,10] [,11] [,12] [1,] -205.3426 107.2501 50.26935 -1272.5128 [2,] -108.6574 286.7499 145.73065 -645.4872 $out [1] 725 179 -521 200 -413 563 -150 $group [1] 5 9 9 9 9 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(91, 159.5, 260.5, 460, 687, -23, 51.5, 168, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/66u5r1272310816.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 530 635 615.0 699.0 809.0 1005.0 814 966.0 1430.0 970.0 [2,] 795.0 906 892 883.0 1060.0 1184.5 1274.5 1458 1580.5 1786.0 1700.5 [3,] 979.5 1054 1007 1155.0 1253.0 1377.0 1542.5 1828 1848.0 2085.0 2038.0 [4,] 1129.5 1243 1084 1353.5 1605.5 1531.0 1604.5 2037 2107.5 2190.5 2258.0 [5,] 1318.0 1577 1286 1708.0 2258.0 1538.0 1613.0 2186 2249.0 2503.0 2439.0 [,12] [,13] [,14] [,15] [,16] [1,] 1007.0 1335.0 1247.0 1728.0 1185 [2,] 1653.5 1907.0 1956.5 2202.5 2104 [3,] 1939.0 2081.5 2255.0 2520.5 2585 [4,] 2256.5 2306.5 2454.5 2662.5 2964 [5,] 2822.0 2835.0 3057.0 2684.0 3923 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 826.9323 900.292 919.4275 940.4018 1004.194 1218.959 1391.985 1563.914 [2,] 1132.0677 1207.708 1094.5725 1369.5982 1501.806 1535.041 1693.015 2092.086 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 1607.632 1900.505 1783.721 1663.968 1899.285 2027.859 2310.691 2071.422 [2,] 2088.368 2269.495 2292.279 2214.032 2263.715 2482.141 2730.309 3098.578 $out [1] 544 1404 2083 779 2116 2326 2907 1138 1212 891 3330 1041 3670 $group [1] 3 3 6 7 7 7 9 10 13 14 14 15 15 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] NA Warning message: In bxp(list(stats = c(464, 795, 979.5, 1129.5, 1318, 530, 906, 1054, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/76u5r1272310816.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,] 1174.938 852.50 1164.750 [2,] 1508.721 1474.25 1430.375 [3,] 1729.400 1811.50 1674.500 [4,] 1859.577 1892.75 1793.188 [5,] 2331.875 2277.50 2203.375 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1569.372 1620.619 1509.019 [2,] 1889.428 2002.381 1839.981 $out [1] 854.00 846.75 $group [1] 1 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1174.9375, 1508.72083333333, 1729.4, 1859.57708333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/12t6l1272310816.ps tmp/12t6l1272310816.png",intern=TRUE)) character(0) > try(system("convert tmp/22t6l1272310816.ps tmp/22t6l1272310816.png",intern=TRUE)) character(0) > try(system("convert tmp/32t6l1272310816.ps tmp/32t6l1272310816.png",intern=TRUE)) character(0) > try(system("convert tmp/4vkno1272310816.ps tmp/4vkno1272310816.png",intern=TRUE)) character(0) > try(system("convert tmp/5vkno1272310816.ps tmp/5vkno1272310816.png",intern=TRUE)) character(0) > try(system("convert tmp/66u5r1272310816.ps tmp/66u5r1272310816.png",intern=TRUE)) character(0) > try(system("convert tmp/76u5r1272310816.ps tmp/76u5r1272310816.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.897 1.501 2.210