R version 2.12.2 (2011-02-25) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-redhat-linux-gnu (64-bit) 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(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/www/wessaorg/rcomp/tmp/17i6a1304967106.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/www/wessaorg/rcomp/tmp/28dof1304967106.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/www/wessaorg/rcomp/tmp/3vegt1304967106.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/www/wessaorg/rcomp/tmp/4h30q1304967106.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/www/wessaorg/rcomp/tmp/5s36u1304967106.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/www/wessaorg/rcomp/tmp/6khs31304967106.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/17i6a1304967106.ps tmp/17i6a1304967106.png",intern=TRUE)) character(0) > try(system("convert tmp/28dof1304967106.ps tmp/28dof1304967106.png",intern=TRUE)) character(0) > try(system("convert tmp/3vegt1304967106.ps tmp/3vegt1304967106.png",intern=TRUE)) character(0) > try(system("convert tmp/4h30q1304967106.ps tmp/4h30q1304967106.png",intern=TRUE)) character(0) > try(system("convert tmp/5s36u1304967106.ps tmp/5s36u1304967106.png",intern=TRUE)) character(0) > try(system("convert tmp/6khs31304967106.ps tmp/6khs31304967106.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 4.010 0.410 4.515