R version 2.9.0 (2009-04-17) 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(1687,1508,1507,1385,1632,1511,1559,1630,1579,1653,2152,2148,1752,1765,1717,1558,1575,1520,1805,1800,1719,2008,2242,2478,2030,1655,1693,1623,1805,1746,1795,1926,1619,1992,2233,2192,2080,1768,1835,1569,1976,1853,1965,1689,1778,1976,2397,2654,2097,1963,1677,1941,2003,1813,2012,1912,2084,2080,2118,2150,1608,1503,1548,1382,1731,1798,1779,1887,2004,2077,2092,2051,1577,1356,1652,1382,1519,1421,1442,1543,1656,1561,1905,2199,1473,1655,1407,1395,1530,1309,1526,1327,1627,1748,1958,2274,1648,1401,1411,1403,1394,1520,1528,1643,1515,1685,2000,2215,1956,1462,1563,1459,1446,1622,1657,1638,1643,1683,2050,2262,1813,1445,1762,1461,1556,1431,1427,1554,1645,1653,2016,2207,1665,1361,1506,1360,1453,1522,1460,1552,1548,1827,1737,1941,1474,1458,1542,1404,1522,1385,1641,1510,1681,1938,1868,1726,1456,1445,1456,1365,1487,1558,1488,1684,1594,1850,1998,2079,1494,1057,1218,1168,1236,1076,1174,1139,1427,1487,1483,1513,1357,1165,1282,1110,1297,1185,1222,1284,1444,1575,1737,1763) > par20 = '' > par19 = '' > par18 = '' > par17 = '' > par16 = '' > par15 = '' > par14 = '' > par13 = '' > par12 = '' > par11 = '' > par10 = '' > par9 = '' > par8 = '' > par7 = '' > par6 = '' > par5 = '' > par4 = '' > par3 = '' > par2 = '' > par1 = '12' > ylab = '' > xlab = '' > main = '' > #'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] 192 > (np <- floor(n / par1)) [1] 16 > 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 16 16 16 16 16 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 1687 1752 2030 2080 2097 1608 1577 1473 1648 1956 1813 1665 1474 [2,] 1508 1765 1655 1768 1963 1503 1356 1655 1401 1462 1445 1361 1458 [3,] 1507 1717 1693 1835 1677 1548 1652 1407 1411 1563 1762 1506 1542 [4,] 1385 1558 1623 1569 1941 1382 1382 1395 1403 1459 1461 1360 1404 [5,] 1632 1575 1805 1976 2003 1731 1519 1530 1394 1446 1556 1453 1522 [6,] 1511 1520 1746 1853 1813 1798 1421 1309 1520 1622 1431 1522 1385 [7,] 1559 1805 1795 1965 2012 1779 1442 1526 1528 1657 1427 1460 1641 [8,] 1630 1800 1926 1689 1912 1887 1543 1327 1643 1638 1554 1552 1510 [9,] 1579 1719 1619 1778 2084 2004 1656 1627 1515 1643 1645 1548 1681 [10,] 1653 2008 1992 1976 2080 2077 1561 1748 1685 1683 1653 1827 1938 [11,] 2152 2242 2233 2397 2118 2092 1905 1958 2000 2050 2016 1737 1868 [12,] 2148 2478 2192 2654 2150 2051 2199 2274 2215 2262 2207 1941 1726 [,14] [,15] [,16] [,17] [1,] 1456 1494 1357 NA [2,] 1445 1057 1165 NA [3,] 1456 1218 1282 NA [4,] 1365 1168 1110 NA [5,] 1487 1236 1297 NA [6,] 1558 1076 1185 NA [7,] 1488 1174 1222 NA [8,] 1684 1139 1284 NA [9,] 1594 1427 1444 NA [10,] 1850 1487 1575 NA [11,] 1998 1483 1737 NA [12,] 2079 1513 1763 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] -179 13 -375 -312 -134 -105 -221 182 -247 -494 -368 -304 -16 [2,] -1 -48 38 67 -286 45 296 -248 10 101 317 145 84 [3,] -122 -159 -70 -266 264 -166 -270 -12 -8 -104 -301 -146 -138 [4,] 247 17 182 407 62 349 137 135 -9 -13 95 93 118 [5,] -121 -55 -59 -123 -190 67 -98 -221 126 176 -125 69 -137 [6,] 48 285 49 112 199 -19 21 217 8 35 -4 -62 256 [7,] 71 -5 131 -276 -100 108 101 -199 115 -19 127 92 -131 [8,] -51 -81 -307 89 172 117 113 300 -128 5 91 -4 171 [9,] 74 289 373 198 -4 73 -95 121 170 40 8 279 257 [10,] 499 234 241 421 38 15 344 210 315 367 363 -90 -70 [11,] -4 236 -41 257 32 -41 294 316 215 212 191 204 -142 [12,] -396 -448 -112 -557 -542 -474 -726 -626 -259 -449 -542 -467 -270 [,14] [,15] [,16] [,17] [1,] -11 -437 -192 NA [2,] 11 161 117 NA [3,] -91 -50 -172 NA [4,] 122 68 187 NA [5,] 71 -160 -112 NA [6,] -70 98 37 NA [7,] 196 -35 62 NA [8,] -90 288 160 NA [9,] 256 60 131 NA [10,] 148 -4 162 NA [11,] 81 30 26 NA [12,] -585 -156 NA 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/1tzj41290722814.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/2tzj41290722814.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/38tpw1290722814.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/48tpw1290722814.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,] 1357.0 1057 1218.0 1360.0 1236.0 1076 1174.0 1284.0 1427.0 1487.0 1737.0 [2,] 1484.0 1381 1433.5 1373.5 1449.5 1403 1451.0 1526.5 1563.5 1653.0 1886.5 [3,] 1656.5 1460 1545.0 1399.0 1526.0 1520 1543.5 1634.0 1635.0 1787.5 2008.0 [4,] 1884.5 1655 1685.0 1509.5 1681.5 1684 1787.0 1744.5 1700.0 1984.0 2135.0 [5,] 2097.0 1963 1835.0 1623.0 2003.0 1853 2012.0 1926.0 1778.0 2080.0 2397.0 [,12] [1,] 1726.0 [2,] 1996.0 [3,] 2171.0 [4,] 2238.5 [5,] 2478.0 $n [1] 16 16 16 16 16 16 16 16 16 16 16 16 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1498.303 1351.77 1445.658 1345.28 1434.36 1409.005 1410.78 1547.89 [2,] 1814.697 1568.23 1644.342 1452.72 1617.64 1630.995 1676.22 1720.11 [,9] [,10] [,11] [,12] [1,] 1581.082 1656.755 1909.842 2075.213 [2,] 1688.918 1918.245 2106.157 2266.787 $out [1] 1941 1168 1110 1139 2084 2004 1483 2654 1513 $group [1] 4 4 4 8 9 9 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(1357, 1484, 1656.5, 1884.5, 2097, 1057, 1381, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/58tpw1290722814.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,] -494.0 -48.0 -301 -13.0 -221 -70.0 -276.0 -307.0 -95.0 -90.0 -142.0 -726.0 [2,] -340.0 4.5 -169 65.0 -131 2.0 -67.5 -66.0 50.0 26.5 11.0 -549.5 [3,] -206.5 56.0 -130 120.0 -105 42.5 66.5 90.0 126.0 222.0 136.0 -467.0 [4,] -60.5 131.0 -60 184.5 68 155.5 111.5 165.5 256.5 353.5 225.5 -333.0 [5,] 182.0 317.0 -8 349.0 176 285.0 196.0 300.0 373.0 499.0 316.0 -112.0 $n [1] 16 16 16 16 16 16 16 16 16 16 16 15 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -316.9025 6.0325 -173.055 72.7975 -183.605 -18.1325 -4.205 -1.4425 [2,] -96.0975 105.9675 -86.945 167.2025 -26.395 103.1325 137.205 181.4425 [,9] [,10] [,11] [,12] [1,] 44.4325 92.835 51.2725 -555.3221 [2,] 207.5675 351.165 220.7275 -378.6779 $out [1] -286 -248 264 407 $group [1] 2 2 3 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-494, -340, -206.5, -60.5, 182, -48, 4.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/61kph1290722814.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,] 1385.0 1520.0 1619 1569 1813.0 1382.0 1356.0 1309.0 1394.0 1446.0 1427.0 [2,] 1509.5 1646.0 1674 1773 1926.5 1578.0 1431.5 1401.0 1407.0 1512.5 1453.0 [3,] 1604.5 1758.5 1800 1909 2007.5 1788.5 1552.0 1528.0 1524.0 1640.5 1600.5 [4,] 1670.0 1906.5 2011 2028 2090.5 2027.5 1654.0 1701.5 1666.5 1819.5 1787.5 [5,] 1687.0 2242.0 2233 2397 2150.0 2092.0 1905.0 1958.0 2000.0 2262.0 2207.0 [,12] [,13] [,14] [,15] [,16] [,17] [1,] 1360.0 1385.0 1365 1057.0 1110.0 NA [2,] 1456.5 1466.0 1456 1153.5 1203.5 NA [3,] 1535.0 1532.0 1523 1227.0 1290.5 NA [4,] 1701.0 1703.5 1767 1485.0 1509.5 NA [5,] 1941.0 1938.0 2079 1513.0 1763.0 NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1531.295 1639.684 1646.292 1792.693 1932.698 1583.48 1450.516 1390.94 [2,] 1677.705 1877.316 1953.708 2025.307 2082.302 1993.52 1653.484 1665.06 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 1405.640 1500.475 1447.932 1423.482 1423.675 1381.151 1075.801 1150.931 [2,] 1642.360 1780.525 1753.068 1646.518 1640.325 1664.849 1378.199 1430.069 [,17] [1,] NA [2,] NA $out [1] 2152 2148 2478 2654 1677 2199 2274 2215 $group [1] 1 1 2 4 5 7 8 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" NA Warning message: In bxp(list(stats = c(1385, 1509.5, 1604.5, 1670, 1687, 1520, 1646, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/71kph1290722814.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,] 1435.312 1399.0 1431.500 [2,] 1532.688 1523.0 1543.250 [3,] 1599.938 1589.5 1622.500 [4,] 1748.750 1722.0 1742.688 [5,] 1999.125 2008.0 2011.125 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1501.390 1498.735 1531.535 [2,] 1698.485 1680.265 1713.465 $out [1] 2115.750 2171.000 2125.125 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1435.3125, 1532.6875, 1599.9375, 1748.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1tzj41290722814.ps tmp/1tzj41290722814.png",intern=TRUE)) character(0) > try(system("convert tmp/2tzj41290722814.ps tmp/2tzj41290722814.png",intern=TRUE)) character(0) > try(system("convert tmp/38tpw1290722814.ps tmp/38tpw1290722814.png",intern=TRUE)) character(0) > try(system("convert tmp/48tpw1290722814.ps tmp/48tpw1290722814.png",intern=TRUE)) character(0) > try(system("convert tmp/58tpw1290722814.ps tmp/58tpw1290722814.png",intern=TRUE)) character(0) > try(system("convert tmp/61kph1290722814.ps tmp/61kph1290722814.png",intern=TRUE)) character(0) > try(system("convert tmp/71kph1290722814.ps tmp/71kph1290722814.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.234 0.981 3.388