R version 2.7.2 (2008-08-25) 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(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,1621,1539,1637,1523,1418,1819,1594,1359,1261,1722,1407,1380,1642,1681,1542,1704,1431) > par1 = '36' > #'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] 97 > (np <- floor(n / par1)) [1] 2 > 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] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 2 2 > arr [,1] [,2] [,3] [1,] 1671 1393 1765 [2,] 1385 1515 1607 [3,] 1632 1510 1577 [4,] 1313 1225 1493 [5,] 1300 1577 1615 [6,] 1431 1417 1700 [7,] 1398 1224 1335 [8,] 1198 1693 1523 [9,] 1292 1633 1621 [10,] 1434 1639 1539 [11,] 1660 1914 1637 [12,] 1837 1586 1523 [13,] 1455 1552 1418 [14,] 1315 2081 1819 [15,] 1642 1500 1594 [16,] 1069 1437 1359 [17,] 1209 1470 1261 [18,] 1586 1849 1722 [19,] 1122 1387 1407 [20,] 1063 1592 1380 [21,] 1125 1589 1642 [22,] 1414 1798 1681 [23,] 1347 1935 1542 [24,] 1403 1887 1704 [25,] 1299 2027 1431 [26,] 1547 2080 NA [27,] 1515 1556 NA [28,] 1247 1682 NA [29,] 1639 1785 NA [30,] 1296 1869 NA [31,] 1063 1781 NA [32,] 1282 2082 NA [33,] 1365 2570 NA [34,] 1268 1862 NA [35,] 1532 1936 NA [36,] 1455 1504 NA > darr [,1] [,2] [,3] [1,] -286 122 -158 [2,] 247 -5 -30 [3,] -319 -285 -84 [4,] -13 352 122 [5,] 131 -160 85 [6,] -33 -193 -365 [7,] -200 469 188 [8,] 94 -60 98 [9,] 142 6 -82 [10,] 226 275 98 [11,] 177 -328 -114 [12,] -382 -34 -105 [13,] -140 529 401 [14,] 327 -581 -225 [15,] -573 -63 -235 [16,] 140 33 -98 [17,] 377 379 461 [18,] -464 -462 -315 [19,] -59 205 -27 [20,] 62 -3 262 [21,] 289 209 39 [22,] -67 137 -139 [23,] 56 -48 162 [24,] -104 140 -273 [25,] 248 53 NA [26,] -32 -524 NA [27,] -268 126 NA [28,] 392 103 NA [29,] -343 84 NA [30,] -233 -88 NA [31,] 219 301 NA [32,] 83 488 NA [33,] -97 -708 NA [34,] 264 74 NA [35,] -77 -432 NA [36,] -62 261 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/1qeyc1225738506.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/2lpms1225738506.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/3dg971225738506.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/4cibi1225738506.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,] 1393 1385 1510.0 1225 1300.0 1417.0 1224.0 1198.0 1292.0 1434.0 1637.0 [2,] 1532 1450 1543.5 1269 1438.5 1424.0 1279.5 1360.5 1456.5 1486.5 1648.5 [3,] 1671 1515 1577.0 1313 1577.0 1431.0 1335.0 1523.0 1621.0 1539.0 1660.0 [4,] 1718 1561 1604.5 1403 1596.0 1565.5 1366.5 1608.0 1627.0 1589.0 1787.0 [5,] 1765 1607 1632.0 1493 1615.0 1700.0 1398.0 1693.0 1633.0 1639.0 1914.0 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [1,] 1523.0 1418.0 1315 1500 1069 1209.0 1586.0 1122.0 1063.0 1125.0 1414.0 [2,] 1554.5 1436.5 1567 1547 1214 1235.0 1654.0 1254.5 1221.5 1357.0 1547.5 [3,] 1586.0 1455.0 1819 1594 1359 1261.0 1722.0 1387.0 1380.0 1589.0 1681.0 [4,] 1711.5 1503.5 1950 1618 1398 1365.5 1785.5 1397.0 1486.0 1615.5 1739.5 [5,] 1837.0 1552.0 2081 1642 1437 1470.0 1849.0 1407.0 1592.0 1642.0 1798.0 [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [1,] 1347.0 1403.0 1299 1547.0 1515.0 1247.0 1639 1296.0 1063 1282 1365.0 [2,] 1444.5 1553.5 1365 1547.0 1515.0 1247.0 1639 1296.0 1063 1282 1365.0 [3,] 1542.0 1704.0 1431 1813.5 1535.5 1464.5 1712 1582.5 1422 1682 1967.5 [4,] 1738.5 1795.5 1729 2080.0 1556.0 1682.0 1785 1869.0 1781 2082 2570.0 [5,] 1935.0 1887.0 2027 2080.0 1556.0 1682.0 1785 1869.0 1781 2082 2570.0 [,34] [,35] [,36] [1,] 1268 1532 1455.0 [2,] 1268 1532 1455.0 [3,] 1565 1734 1479.5 [4,] 1862 1936 1504.0 [5,] 1862 1936 1504.0 $n [1] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 2 2 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1501.328 1413.744 1521.355 1190.763 1433.326 1301.922 1255.637 1297.227 [2,] 1840.672 1616.256 1632.645 1435.237 1720.674 1560.078 1414.363 1748.773 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 1465.468 1445.498 1533.658 1442.782 1393.882 1469.622 1529.233 1191.153 [2,] 1776.532 1632.502 1786.342 1729.218 1516.118 2168.378 1658.767 1526.847 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 1141.956 1602.044 1257.010 1138.720 1353.193 1505.855 1273.809 1483.244 [2,] 1380.044 1841.956 1516.990 1621.280 1824.807 1856.145 1810.191 1924.756 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 1098.954 1218.017 1489.694 978.5055 1548.885 942.328 619.8298 788.217 [2,] 1763.046 2408.983 1581.306 1950.4945 1875.115 2222.672 2224.1702 2575.783 [,33] [,34] [,35] [,36] [1,] 621.2394 901.3661 1282.640 1424.756 [2,] 3313.7606 2228.6339 2185.360 1534.244 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" Warning message: In bxp(list(stats = c(1393, 1532, 1671, 1718, 1765, 1385, 1450, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5hp1f1225738506.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,] -286 -30.0 -319.0 -13.0 -160.0 -365 -200.0 -60 -82 98.0 -328.0 -382.0 [2,] -222 -17.5 -302.0 54.5 -37.5 -279 -6.0 17 -38 162.0 -221.0 -243.5 [3,] -158 -5.0 -285.0 122.0 85.0 -193 188.0 94 6 226.0 -114.0 -105.0 [4,] -18 121.0 -184.5 237.0 108.0 -113 328.5 96 74 250.5 31.5 -69.5 [5,] 122 247.0 -84.0 352.0 131.0 -33 469.0 98 142 275.0 177.0 -34.0 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] -140.0 -581 -573 -98.0 377 -464.0 -59 -3.0 39 -139 -48 -273.0 [2,] 130.5 -403 -404 -32.5 378 -463.0 -43 29.5 124 -103 4 -188.5 [3,] 401.0 -225 -235 33.0 379 -462.0 -27 62.0 209 -67 56 -104.0 [4,] 465.0 51 -149 86.5 420 -388.5 89 162.0 249 35 109 18.0 [5,] 529.0 327 -63 140.0 461 -315.0 205 262.0 289 137 162 140.0 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [1,] 53.0 -524 -268 103.0 -343.0 -233.0 219 83.0 -708.0 74 -432.0 [2,] 53.0 -524 -268 103.0 -343.0 -233.0 219 83.0 -708.0 74 -432.0 [3,] 150.5 -278 -71 247.5 -129.5 -160.5 260 285.5 -402.5 169 -254.5 [4,] 248.0 -32 126 392.0 84.0 -88.0 301 488.0 -97.0 264 -77.0 [5,] 248.0 -32 126 392.0 84.0 -88.0 301 488.0 -97.0 264 -77.0 [,36] [1,] -62.0 [2,] -62.0 [3,] 99.5 [4,] 261.0 [5,] 261.0 $n [1] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 2 2 2 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -344.09154 -131.3416 -392.1851 -44.47895 -47.72705 -344.42743 -117.1354 [2,] 28.09154 121.3416 -177.8149 288.47895 217.72705 -41.57257 493.1354 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 21.93514 -96.1679 145.2691 -344.3339 -263.72514 95.8646 -639.1449 [2,] 166.06486 108.1679 306.7309 116.3339 53.72514 706.1354 189.1449 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] -467.614423 -75.5534 340.6870 -529.9599 -147.41217 -58.86828 94.97332 [2,] -2.385577 141.5534 417.3130 -394.0401 93.41217 182.86828 323.02668 [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] -192.88545 -39.78241 -292.37207 -67.3596 -827.6765 -511.1881 -75.3791 [2,] 58.88545 151.78241 84.37207 368.3596 271.6765 369.1881 570.3791 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [1,] -606.5567 -322.498164 168.3872 -166.9776 -1085.1267 -43.27346 -651.1162 [2,] 347.5567 1.498164 351.6128 737.9776 280.1267 381.27346 142.1162 [,36] [1,] -261.3649 [2,] 460.3649 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" Warning message: In bxp(list(stats = c(-286, -222, -158, -18, 122, -30, -17.5, -5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6kab01225738506.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] [1,] 1063.0 1224 1261 [2,] 1275.0 1507 1431 [3,] 1375.0 1636 1577 [4,] 1523.5 1878 1642 [5,] 1837.0 2082 1819 $n [1] 36 36 25 $conf [,1] [,2] [,3] [1,] 1309.562 1538.303 1510.324 [2,] 1440.438 1733.697 1643.676 $out [1] 2570 $group [1] 2 $names [1] "1" "2" NA Warning message: In bxp(list(stats = c(1063, 1275, 1375, 1523.5, 1837, 1224, 1507, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/77rsf1225738506.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,] 1288.333 1261.0 1300.250 [2,] 1467.917 1443.0 1474.750 [3,] 1551.167 1571.0 1544.375 [4,] 1656.667 1676.0 1659.000 [5,] 1813.500 1967.5 1813.500 $n [1] 36 36 36 $conf [,1] [,2] [,3] [1,] 1501.462 1509.643 1495.856 [2,] 1600.871 1632.357 1592.894 $out [1] 1967.5 1967.5 $group [1] 1 3 $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > system("convert tmp/1qeyc1225738506.ps tmp/1qeyc1225738506.png") > system("convert tmp/2lpms1225738506.ps tmp/2lpms1225738506.png") > system("convert tmp/3dg971225738506.ps tmp/3dg971225738506.png") > system("convert tmp/4cibi1225738506.ps tmp/4cibi1225738506.png") > system("convert tmp/5hp1f1225738506.ps tmp/5hp1f1225738506.png") > system("convert tmp/6kab01225738506.ps tmp/6kab01225738506.png") > system("convert tmp/77rsf1225738506.ps tmp/77rsf1225738506.png") > > > proc.time() user system elapsed 1.514 0.976 1.849