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. Natural language support but running in an English locale 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) > 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] 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/freestat/rcomp/tmp/12f6g1288702005.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/2uon11288702005.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/3uon11288702005.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/4uon11288702005.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/freestat/rcomp/tmp/55fm41288702005.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/freestat/rcomp/tmp/65fm41288702005.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/freestat/rcomp/tmp/7g6mp1288702005.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/12f6g1288702005.ps tmp/12f6g1288702005.png",intern=TRUE)) character(0) > try(system("convert tmp/2uon11288702005.ps tmp/2uon11288702005.png",intern=TRUE)) character(0) > try(system("convert tmp/3uon11288702005.ps tmp/3uon11288702005.png",intern=TRUE)) character(0) > try(system("convert tmp/4uon11288702005.ps tmp/4uon11288702005.png",intern=TRUE)) character(0) > try(system("convert tmp/55fm41288702005.ps tmp/55fm41288702005.png",intern=TRUE)) character(0) > try(system("convert tmp/65fm41288702005.ps tmp/65fm41288702005.png",intern=TRUE)) character(0) > try(system("convert tmp/7g6mp1288702005.ps tmp/7g6mp1288702005.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.955 1.474 2.190