R version 2.10.1 (2009-12-14) 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(1686,1591,2304,1712,1471,1377,1966,2453,1984,2596,4087,5179,1530,1523,1633,1976,1170,1480,1781,2472,1981,2273,3857,4551,1510,1329,1518,1790,1537,1449,1954,1897,1706,2514,3593,4524,1609,1638,2030,1375,1320,1245,1600,2298,2191,2511,3440,4923,1609,1435,2061,1789,1567,1404,1597,3159,1759,2504,4273,5274,1771,1682,1846,1589,1896,1379,1645,2512,1771,3727,4388,5434,1606,1523,1577,1605,1765,1403,2584,3318,1562,2349,3987,5891,1389,1442,1548,1935,1518,1250,1847,1930,2638,3114,4405,7242,1853,1779,2108,2336,1728,1661,2230,1645,2421,3740,4988,6757,1757,1394,1982,1650,1654,1406,1971,1968,2608,3845,4514,6694,1720,1321,1859,1628,1615,1457,1899,1605,2424,3116,4286,6047,1902,2049,1874,1279,1432,1540,2214,1857,2408,3252,3627,6153,1577,1667,1993,1997,1783,1625,2076,1773,2377,3088,4096,6119,1494,1564,1898,2121,1831,1515,2048,2795,1749,3339,4227,6410,1197,1968,1720,1725,1674,1693,2031,1495,2968,3385,3729,5999,1070,1402,1897,1862,1670,1688,2031) > 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,] 1686 1530 1510 1609 1609 1771 1606 1389 1853 1757 1720 1902 1577 [2,] 1591 1523 1329 1638 1435 1682 1523 1442 1779 1394 1321 2049 1667 [3,] 2304 1633 1518 2030 2061 1846 1577 1548 2108 1982 1859 1874 1993 [4,] 1712 1976 1790 1375 1789 1589 1605 1935 2336 1650 1628 1279 1997 [5,] 1471 1170 1537 1320 1567 1896 1765 1518 1728 1654 1615 1432 1783 [6,] 1377 1480 1449 1245 1404 1379 1403 1250 1661 1406 1457 1540 1625 [7,] 1966 1781 1954 1600 1597 1645 2584 1847 2230 1971 1899 2214 2076 [8,] 2453 2472 1897 2298 3159 2512 3318 1930 1645 1968 1605 1857 1773 [9,] 1984 1981 1706 2191 1759 1771 1562 2638 2421 2608 2424 2408 2377 [10,] 2596 2273 2514 2511 2504 3727 2349 3114 3740 3845 3116 3252 3088 [11,] 4087 3857 3593 3440 4273 4388 3987 4405 4988 4514 4286 3627 4096 [12,] 5179 4551 4524 4923 5274 5434 5891 7242 6757 6694 6047 6153 6119 [,14] [,15] [,16] [1,] 1494 1197 1070 [2,] 1564 1968 1402 [3,] 1898 1720 1897 [4,] 2121 1725 1862 [5,] 1831 1674 1670 [6,] 1515 1693 1688 [7,] 2048 2031 2031 [8,] 2795 1495 NA [9,] 1749 2968 NA [10,] 3339 3385 NA [11,] 4227 3729 NA [12,] 6410 5999 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] -95 -7 -181 29 -174 -89 -83 53 -74 -363 -399 147 [2,] 713 110 189 392 626 164 54 106 329 588 538 -175 [3,] -592 343 272 -655 -272 -257 28 387 228 -332 -231 -595 [4,] -241 -806 -253 -55 -222 307 160 -417 -608 4 -13 153 [5,] -94 310 -88 -75 -163 -517 -362 -268 -67 -248 -158 108 [6,] 589 301 505 355 193 266 1181 597 569 565 442 674 [7,] 487 691 -57 698 1562 867 734 83 -585 -3 -294 -357 [8,] -469 -491 -191 -107 -1400 -741 -1756 708 776 640 819 551 [9,] 612 292 808 320 745 1956 787 476 1319 1237 692 844 [10,] 1491 1584 1079 929 1769 661 1638 1291 1248 669 1170 375 [11,] 1092 694 931 1483 1001 1046 1904 2837 1769 2180 1761 2526 [12,] -3649 -3041 -2915 -3314 -3503 -3828 -4502 -5389 -5000 -4974 -4145 -4576 [,13] [,14] [,15] [,16] [1,] 90 70 771 332 [2,] 326 334 -248 495 [3,] 4 223 5 -35 [4,] -214 -290 -51 -192 [5,] -158 -316 19 18 [6,] 451 533 338 343 [7,] -303 747 -536 NA [8,] 604 -1046 1473 NA [9,] 711 1590 417 NA [10,] 1008 888 344 NA [11,] 2023 2183 2270 NA [12,] -4625 -5213 -4929 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/yougetitorg/rcomp/tmp/1kkzj1301575837.ps",horizontal=F,onefile=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/yougetitorg/rcomp/tmp/2ef7i1301575837.ps",horizontal=F,onefile=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/yougetitorg/rcomp/tmp/3pyfp1301575837.ps",horizontal=F,onefile=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/yougetitorg/rcomp/tmp/4jhs61301575837.ps",horizontal=F,onefile=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,] 1197.0 1321.0 1518.0 1279.0 1170.0 1245.0 1597.0 1495 1562.0 2273.0 3440 [2,] 1502.0 1418.5 1676.5 1616.5 1494.5 1391.0 1814.0 1815 1765.0 2512.5 3793 [3,] 1607.5 1543.5 1885.5 1757.0 1634.5 1453.0 1968.5 1968 2191.0 3114.0 4096 [4,] 1738.5 1674.5 2011.5 1955.5 1746.5 1582.5 2062.0 2492 2422.5 3362.0 4337 [5,] 1902.0 2049.0 2304.0 2336.0 1896.0 1693.0 2230.0 3318 2968.0 3845.0 4988 [,12] [1,] 4524.0 [2,] 5226.5 [3,] 5999.0 [4,] 6281.5 [5,] 7242.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,] 1514.082 1442.38 1753.175 1623.095 1534.96 1377.358 1870.54 1691.815 [2,] 1700.918 1644.62 2017.825 1890.905 1734.04 1528.642 2066.46 2244.185 [,9] [,10] [,11] [,12] [1,] 1922.77 2767.443 3874.073 5568.608 [2,] 2459.23 3460.557 4317.927 6429.392 $out [1] 1070 2584 $group [1] 1 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(1197, 1502, 1607.5, 1738.5, 1902, 1321, 1418.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/5phhk1301575837.ps",horizontal=F,onefile=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,] -399.0 -248.0 -655.0 -608.0 -517.0 193.0 -585.0 -1756 292.0 344.0 694.0 [2,] -134.5 108.0 -302.0 -271.5 -258.0 340.5 -298.5 -616 544.0 778.5 1069.0 [3,] -40.5 327.5 -15.5 -203.0 -126.0 478.0 83.0 -107 745.0 1079.0 1769.0 [4,] 80.0 516.5 225.5 -4.5 -24.5 579.0 716.0 674 1040.5 1391.0 2181.5 [5,] 332.0 713.0 387.0 307.0 310.0 674.0 1562.0 1473 1590.0 1769.0 2837.0 [,12] [1,] -5389.0 [2,] -4951.5 [3,] -4502.0 [4,] -3576.0 [5,] -2915.0 $n [1] 16 16 16 16 16 16 15 15 15 15 15 15 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -125.2275 166.1425 -223.8625 -308.465 -218.2325 383.7925 -330.8696 [2,] 44.2275 488.8575 192.8625 -97.535 -33.7675 572.2075 496.8696 [,8] [,9] [,10] [,11] [,12] [1,] -633.261 542.4507 829.128 1315.151 -5063.141 [2,] 419.261 947.5493 1328.872 2222.849 -3940.859 $out [1] 771 -806 1181 1956 $group [1] 1 4 6 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-399, -134.5, -40.5, 80, 332, -248, 108, 327.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/6hpvz1301575837.ps",horizontal=F,onefile=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,] 1377.0 1170.0 1329 1245.0 1404.0 1379.0 1403.0 1250.0 1645.0 1394.0 1321.0 [2,] 1638.5 1526.5 1514 1487.5 1582.0 1663.5 1569.5 1480.0 1753.5 1652.0 1610.0 [3,] 1975.0 1878.5 1748 1834.0 1774.0 1808.5 1685.5 1888.5 2169.0 1969.5 1789.5 [4,] 2524.5 2372.5 2234 2404.5 2831.5 3119.5 2951.0 2876.0 3080.5 3226.5 2770.0 [5,] 2596.0 2472.0 2514 3440.0 4273.0 4388.0 3987.0 4405.0 4988.0 4514.0 4286.0 [,12] [,13] [,14] [,15] [,16] [1,] 1279.0 1577.0 1494.0 1197.0 1070.0 [2,] 1698.5 1720.0 1656.5 1683.5 1536.0 [3,] 1975.5 1995.0 1973.0 1846.5 1688.0 [4,] 2830.0 2732.5 3067.0 3176.5 1879.5 [5,] 3627.0 4096.0 4227.0 3729.0 2031.0 $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,] 1570.889 1492.634 1419.603 1415.750 1204.095 1144.409 1055.389 1251.775 [2,] 2379.111 2264.366 2076.397 2252.250 2343.905 2472.591 2315.611 2525.225 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 1563.746 1251.36 1260.416 1459.415 1533.192 1329.661 1165.533 1482.867 [2,] 2774.254 2687.64 2318.584 2491.585 2456.808 2616.339 2527.467 1893.133 $out [1] 4087 5179 3857 4551 3593 4524 4923 5274 5434 5891 7242 6757 6694 6047 6153 [16] 6119 6410 5999 $group [1] 1 1 2 2 3 3 4 5 6 7 8 9 10 11 12 13 14 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(1377, 1638.5, 1975, 2524.5, 2596, 1170, 1526.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/72j0l1301575837.ps",horizontal=F,onefile=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,] 1473.250 1453.00 1479.125 [2,] 1591.812 1621.00 1619.688 [3,] 1916.312 1926.75 1896.500 [4,] 2617.667 2652.50 2545.375 [5,] 4099.800 4096.00 2937.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1448.414 1456.276 1474.288 [2,] 2384.211 2397.224 2318.712 $out [1] 5813.133 5999.000 4065.000 5754.000 $group [1] 1 2 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1473.25, 1591.8125, 1916.3125, 2617.66666666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1kkzj1301575837.ps tmp/1kkzj1301575837.png",intern=TRUE)) character(0) > try(system("convert tmp/2ef7i1301575837.ps tmp/2ef7i1301575837.png",intern=TRUE)) character(0) > try(system("convert tmp/3pyfp1301575837.ps tmp/3pyfp1301575837.png",intern=TRUE)) character(0) > try(system("convert tmp/4jhs61301575837.ps tmp/4jhs61301575837.png",intern=TRUE)) character(0) > try(system("convert tmp/5phhk1301575837.ps tmp/5phhk1301575837.png",intern=TRUE)) character(0) > try(system("convert tmp/6hpvz1301575837.ps tmp/6hpvz1301575837.png",intern=TRUE)) character(0) > try(system("convert tmp/72j0l1301575837.ps tmp/72j0l1301575837.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.540 1.200 2.059