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(2560 + ,2491 + ,2380 + ,2291 + ,2079 + ,1929 + ,1851 + ,1607 + ,1661 + ,2259 + ,1668 + ,2011 + ,1944 + ,1958 + ,1844 + ,1868 + ,1701 + ,2338 + ,2018 + ,1302 + ,2168 + ,2139 + ,1560 + ,2093 + ,1973 + ,2090 + ,2811 + ,1984 + ,1849 + ,2433 + ,2071 + ,1855 + ,1756 + ,1898 + ,1770 + ,1969 + ,1769 + ,2139 + ,3013 + ,2061 + ,2132 + ,2973 + ,2081 + ,2257 + ,2075 + ,2084 + ,1747 + ,2092 + ,1919 + ,2551 + ,2643 + ,2153 + ,2496 + ,2645 + ,2035 + ,2294 + ,2205 + ,2044 + ,1762 + ,1897 + ,1821 + ,1905 + ,2111 + ,1643 + ,1956 + ,1977 + ,1685 + ,1393 + ,1574 + ,1793 + ,1562 + ,1510 + ,1675 + ,1965 + ,2173 + ,2395 + ,2197 + ,2257 + ,2885 + ,1594 + ,1950 + ,1772 + ,1280 + ,1724 + ,1473 + ,1461 + ,1576 + ,1900 + ,1618 + ,2303 + ,1994 + ,1575 + ,1893 + ,1788 + ,1817 + ,3233 + ,727 + ,1121 + ,1665 + ,1401 + ,1415 + ,2058 + ,1544 + ,1379 + ,1402 + ,1313 + ,1296 + ,1398 + ,1288 + ,1563 + ,1972 + ,1496 + ,1481 + ,1819 + ,1479 + ,1635 + ,1511 + ,1547 + ,1388 + ,1958 + ,1390 + ,1597 + ,1842 + ,1396 + ,1671 + ,1385 + ,1632 + ,1313 + ,1300 + ,1431 + ,1398 + ,1198 + ,1292 + ,1434 + ,1660 + ,1837 + ,1455 + ,1315 + ,1642 + ,1069 + ,1209 + ,1586 + ,1122 + ,1063 + ,1125 + ,1414 + ,1347 + ,1403 + ,1299 + ,1547 + ,1515 + ,1247 + ,1639 + ,1296 + ,1063 + ,1282 + ,1365 + ,1268 + ,1532 + ,1455 + ,1393 + ,1515 + ,1510 + ,1225 + ,1577 + ,1417 + ,1224 + ,1693 + ,1633 + ,1639 + ,1914 + ,1586 + ,1552 + ,2081 + ,1500 + ,1437 + ,1470 + ,1849 + ,1387 + ,1592 + ,1589 + ,1798 + ,1935 + ,1887 + ,2027 + ,2080 + ,1556 + ,1682 + ,1785 + ,1869 + ,1781 + ,2082 + ,2570 + ,1862 + ,1936 + ,1504 + ,1765 + ,1607 + ,1577 + ,1493 + ,1615 + ,1700 + ,1335 + ,1523 + ,1623 + ,1540 + ,1637 + ,1524 + ,1419 + ,1821 + ,1593 + ,1357 + ,1263 + ,1750 + ,1405 + ,1393 + ,1639 + ,1679 + ,1551 + ,1744 + ,1429 + ,1784) > 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] 222 > (np <- floor(n / par1)) [1] 18 > 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] 19 19 19 19 19 19 18 18 18 18 18 18 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 2560 1944 1973 1769 1919 1821 1675 1473 727 1288 1390 1292 1125 [2,] 2491 1958 2090 2139 2551 1905 1965 1461 1121 1563 1597 1434 1414 [3,] 2380 1844 2811 3013 2643 2111 2173 1576 1665 1972 1842 1660 1347 [4,] 2291 1868 1984 2061 2153 1643 2395 1900 1401 1496 1396 1837 1403 [5,] 2079 1701 1849 2132 2496 1956 2197 1618 1415 1481 1671 1455 1299 [6,] 1929 2338 2433 2973 2645 1977 2257 2303 2058 1819 1385 1315 1547 [7,] 1851 2018 2071 2081 2035 1685 2885 1994 1544 1479 1632 1642 1515 [8,] 1607 1302 1855 2257 2294 1393 1594 1575 1379 1635 1313 1069 1247 [9,] 1661 2168 1756 2075 2205 1574 1950 1893 1402 1511 1300 1209 1639 [10,] 2259 2139 1898 2084 2044 1793 1772 1788 1313 1547 1431 1586 1296 [11,] 1668 1560 1770 1747 1762 1562 1280 1817 1296 1388 1398 1122 1063 [12,] 2011 2093 1969 2092 1897 1510 1724 3233 1398 1958 1198 1063 1282 [,14] [,15] [,16] [,17] [,18] [,19] [1,] 1365 1633 1589 2570 1623 1639 [2,] 1268 1639 1798 1862 1540 1679 [3,] 1532 1914 1935 1936 1637 1551 [4,] 1455 1586 1887 1504 1524 1744 [5,] 1393 1552 2027 1765 1419 1429 [6,] 1515 2081 2080 1607 1821 1784 [7,] 1510 1500 1556 1577 1593 NA [8,] 1225 1437 1682 1493 1357 NA [9,] 1577 1470 1785 1615 1263 NA [10,] 1417 1849 1869 1700 1750 NA [11,] 1224 1387 1781 1335 1405 NA [12,] 1693 1592 2082 1523 1393 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] -69 14 117 370 632 84 290 -12 394 275 207 142 289 [2,] -111 -114 721 874 92 206 208 115 544 409 245 226 -67 [3,] -89 24 -827 -952 -490 -468 222 324 -264 -476 -446 177 56 [4,] -212 -167 -135 71 343 313 -198 -282 14 -15 275 -382 -104 [5,] -150 637 584 841 149 21 60 685 643 338 -286 -140 248 [6,] -78 -320 -362 -892 -610 -292 628 -309 -514 -340 247 327 -32 [7,] -244 -716 -216 176 259 -292 -1291 -419 -165 156 -319 -573 -268 [8,] 54 866 -99 -182 -89 181 356 318 23 -124 -13 140 392 [9,] 598 -29 142 9 -161 219 -178 -105 -89 36 131 377 -343 [10,] -591 -579 -128 -337 -282 -231 -492 29 -17 -159 -33 -464 -233 [11,] 343 533 199 345 135 -52 444 1416 102 570 -200 -59 219 [12,] -67 -120 -200 -173 -76 165 -251 -2506 -110 -568 94 62 83 [,14] [,15] [,16] [,17] [,18] [,19] [1,] -97 6 209 -708 -83 40 [2,] 264 275 137 74 97 -128 [3,] -77 -328 -48 -432 -113 193 [4,] -62 -34 140 261 -105 -315 [5,] 122 529 53 -158 402 355 [6,] -5 -581 -524 -30 -228 NA [7,] -285 -63 126 -84 -236 NA [8,] 352 33 103 122 -94 NA [9,] -160 379 84 85 487 NA [10,] -193 -462 -88 -365 -345 NA [11,] 469 205 301 188 -12 NA [12,] -60 -3 488 100 246 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/1ga4d1229336903.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/rcomp/tmp/2sc8i1229336903.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/rcomp/tmp/3i4yj1229336903.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/rcomp/tmp/4emaw1229336903.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,] 727.0 1121.0 1347.0 1396 1299.0 1315.0 1479 1069 1209 1296 1063.0 1063.0 [2,] 1377.5 1500.5 1648.5 1500 1442.0 1695.5 1544 1313 1470 1547 1296.0 1398.0 [3,] 1633.0 1679.0 1914.0 1744 1671.0 1977.0 1637 1465 1627 1780 1401.5 1708.5 [4,] 1870.0 1961.5 2142.0 1942 1991.5 2280.0 2018 1635 1893 1898 1747.0 2011.0 [5,] 2570.0 2551.0 2811.0 2395 2496.0 2973.0 2081 1855 2205 2259 1817.0 2093.0 $n [1] 19 19 19 19 19 19 18 18 18 18 18 18 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1454.480 1511.898 1735.118 1583.785 1471.819 1765.132 1460.478 1345.084 [2,] 1811.520 1846.102 2092.882 1904.215 1870.181 2188.868 1813.522 1584.916 [,9] [,10] [,11] [,12] [1,] 1469.471 1649.284 1233.543 1480.213 [2,] 1784.529 1910.716 1569.457 1936.787 $out [1] 3013 2885 2257 2294 3233 $group [1] 3 7 8 8 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(727, 1377.5, 1633, 1870, 2570, 1121, 1500.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5tntb1229336903.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,] -97 -128.0 -952 -382.0 -286.0 -892.0 -573 -182.0 -343 -591.0 -200 -568.0 [2,] -3 83.0 -457 -182.5 37.0 -514.0 -319 -89.0 -105 -462.0 102 -173.0 [3,] 117 206.0 -113 -62.0 248.0 -300.5 -240 78.5 60 -257.5 212 -63.5 [4,] 282 269.5 40 105.5 556.5 -30.0 -63 318.0 219 -128.0 444 94.0 [5,] 632 544.0 324 343.0 841.0 628.0 259 866.0 598 29.0 570 488.0 $n [1] 19 19 19 19 19 18 18 18 18 18 18 18 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 13.69410 138.3981 -293.151 -166.39334 59.69327 -480.7462 -335.3369 [2,] 220.30590 273.6019 67.151 42.39334 436.30673 -120.2538 -144.6631 [,8] [,9] [,10] [,11] [,12] [1,] -73.0707 -60.6607 -381.8848 84.63593 -162.93336 [2,] 230.0707 180.6607 -133.1152 339.36407 35.93336 $out [1] -708 721 874 -716 -1291 1416 -2506 $group [1] 1 2 2 7 7 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-97, -3, 117, 282, 632, -128, 83, 206, 269.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6zkxu1229336903.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,] 1607.0 1302.0 1756.0 2061 1762.0 1393.0 1280.0 1461.0 1121.0 1288 1198.0 [2,] 1759.5 1772.5 1852.0 2068 1977.0 1568.0 1699.5 1575.5 1304.5 1480 1349.0 [3,] 2045.0 1951.0 1971.0 2088 2179.0 1739.0 1957.5 1802.5 1399.5 1529 1397.0 [4,] 2335.5 2116.0 2080.5 2198 2523.5 1930.5 2227.0 1947.0 1479.5 1727 1614.5 [5,] 2560.0 2338.0 2090.0 2257 2645.0 2111.0 2885.0 2303.0 1665.0 1972 1842.0 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [1,] 1063.0 1063.0 1224.0 1387 1556.0 1335.0 1263 1429 [2,] 1165.5 1264.5 1316.5 1485 1731.5 1513.5 1399 1551 [3,] 1374.5 1323.0 1436.0 1589 1833.5 1611.0 1532 1659 [4,] 1614.0 1464.5 1523.5 1744 1981.0 1813.5 1630 1744 [5,] 1837.0 1639.0 1693.0 2081 2082.0 1936.0 1821 1784 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1782.283 1794.327 1866.780 2028.706 1929.738 1573.661 1716.904 1633.056 [2,] 2307.717 2107.673 2075.220 2147.294 2428.262 1904.339 2198.096 1971.944 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 1319.681 1416.342 1275.904 1169.936 1231.779 1341.586 1470.868 1719.701 [2,] 1479.319 1641.658 1518.096 1579.064 1414.221 1530.414 1707.132 1947.299 [,17] [,18] [,19] [1,] 1474.168 1426.639 1534.509 [2,] 1747.832 1637.361 1783.491 $out [1] 2811 2433 1769 3013 2973 1747 3233 727 2058 2570 $group [1] 3 3 4 4 4 4 8 9 9 17 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" NA Warning message: In bxp(list(stats = c(1607, 1759.5, 2045, 2335.5, 2560, 1302, 1772.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7iodg1229336903.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,] 1539.667 1465 1516.500 [2,] 1660.463 1630 1648.438 [3,] 1756.833 1675 1718.875 [4,] 1775.871 1762 1755.250 [5,] 1787.111 1914 1895.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1704.195 1614.794 1670.157 [2,] 1809.472 1735.206 1767.593 $out [1] 1975.895 1993.000 1475.833 1977.000 1401.500 1987.750 1476.000 $group [1] 1 1 1 2 2 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1539.66666666667, 1660.46345029240, 1756.83333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1ga4d1229336903.ps tmp/1ga4d1229336903.png") > system("convert tmp/2sc8i1229336903.ps tmp/2sc8i1229336903.png") > system("convert tmp/3i4yj1229336903.ps tmp/3i4yj1229336903.png") > system("convert tmp/4emaw1229336903.ps tmp/4emaw1229336903.png") > system("convert tmp/5tntb1229336903.ps tmp/5tntb1229336903.png") > system("convert tmp/6zkxu1229336903.ps tmp/6zkxu1229336903.png") > system("convert tmp/7iodg1229336903.ps tmp/7iodg1229336903.png") > > > proc.time() user system elapsed 1.243 0.950 1.614