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)) > ari <- array(0,dim=par1) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[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 > arr.sd <- array(NA,dim=par1) > arr.range <- array(NA,dim=par1) > arr.iqr <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE)) + arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE) + arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE) + } > overall.sd <- sqrt(var(x)) > overall.range <- max(x) - min(x) > overall.iqr <- quantile(x,0.75) - quantile(x,0.25) > postscript(file="/var/yougetitorg/rcomp/tmp/1auru1304522778.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.sd,0) > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/2lcw91304522778.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.range,0) > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/3qt9g1304522778.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.iqr,0) > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/44y0u1304522778.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/586s61304522778.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/6c6p81304522778.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.sd,arr.range,arr.iqr)) > names(z) <- list('S.D.','Range','IQR') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] 140.0288 448 164.25 [2,] 216.2518 757 227.75 [3,] 260.5194 1022 313.50 [4,] 475.1343 1560 667.25 [5,] 812.9664 2718 1055.00 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 142.4414 655.7463 113.0411 [2,] 378.5975 1388.2537 513.9589 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(140.028806560174, 216.251758258677, 260.519437613619, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1auru1304522778.ps tmp/1auru1304522778.png",intern=TRUE)) character(0) > try(system("convert tmp/2lcw91304522778.ps tmp/2lcw91304522778.png",intern=TRUE)) character(0) > try(system("convert tmp/3qt9g1304522778.ps tmp/3qt9g1304522778.png",intern=TRUE)) character(0) > try(system("convert tmp/44y0u1304522778.ps tmp/44y0u1304522778.png",intern=TRUE)) character(0) > try(system("convert tmp/586s61304522778.ps tmp/586s61304522778.png",intern=TRUE)) character(0) > try(system("convert tmp/6c6p81304522778.ps tmp/6c6p81304522778.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.330 1.020 1.891