R version 2.7.2 (2008-08-25) 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. 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(1671,1385,1632,1313,1300,1431,1398,1198,1292,1434,1660,1837,1455,1315,1642,1069,1209,1586,1122,1063,1125,1414,1347,1403,1299,1547,1515,1247,1639,1296,1063,1282,1365,1268,1532,1455,1393,1515,1510,1225,1577,1417,1224,1693,1633,1639,1914,1586,1552,2081,1500,1437,1470,1849,1387,1592,1589,1798,1935,1887,2027,2080,1556,1682,1785,1869,1781,2082,2570,1862,1936,1504,1765,1607,1577,1493,1615,1700,1335,1523,1621,1539,1637,1523,1418,1819,1594,1359,1261,1722,1407,1380,1642,1681,1542,1704,1431) > 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] 97 > (np <- floor(n / par1)) [1] 8 > 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] 9 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 1671 1455 1299 1393 1552 2027 1765 1418 1431 [2,] 1385 1315 1547 1515 2081 2080 1607 1819 NA [3,] 1632 1642 1515 1510 1500 1556 1577 1594 NA [4,] 1313 1069 1247 1225 1437 1682 1493 1359 NA [5,] 1300 1209 1639 1577 1470 1785 1615 1261 NA [6,] 1431 1586 1296 1417 1849 1869 1700 1722 NA [7,] 1398 1122 1063 1224 1387 1781 1335 1407 NA [8,] 1198 1063 1282 1693 1592 2082 1523 1380 NA [9,] 1292 1125 1365 1633 1589 2570 1621 1642 NA [10,] 1434 1414 1268 1639 1798 1862 1539 1681 NA [11,] 1660 1347 1532 1914 1935 1936 1637 1542 NA [12,] 1837 1403 1455 1586 1887 1504 1523 1704 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -286 -140 248 122 529 53 -158 401 NA [2,] 247 327 -32 -5 -581 -524 -30 -225 NA [3,] -319 -573 -268 -285 -63 126 -84 -235 NA [4,] -13 140 392 352 33 103 122 -98 NA [5,] 131 377 -343 -160 379 84 85 461 NA [6,] -33 -464 -233 -193 -462 -88 -365 -315 NA [7,] -200 -59 219 469 205 301 188 -27 NA [8,] 94 62 83 -60 -3 488 98 262 NA [9,] 142 289 -97 6 209 -708 -82 39 NA [10,] 226 -67 264 275 137 74 98 -139 NA [11,] 177 56 -77 -328 -48 -432 -114 162 NA [12,] -382 -104 -62 -34 140 261 -105 -273 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/1ls4u1226313061.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/2eylj1226313061.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/3wr231226313061.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/4u2ru1226313061.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,] 1299 1315.0 1500.0 1069 1209.0 1296.0 1063.0 1063.0 1125.0 1268.0 1347.0 [2,] 1418 1450.0 1512.5 1236 1280.5 1424.0 1173.0 1240.0 1328.5 1424.0 1537.0 [3,] 1455 1577.0 1566.5 1336 1523.5 1643.0 1361.0 1451.5 1605.0 1589.0 1648.5 [4,] 1671 1949.5 1613.0 1465 1627.0 1785.5 1402.5 1642.5 1637.5 1739.5 1924.5 [5,] 2027 2081.0 1642.0 1682 1785.0 1869.0 1407.0 2082.0 1642.0 1862.0 1936.0 [,12] [1,] 1403.0 [2,] 1479.5 [3,] 1554.5 [4,] 1770.5 [5,] 1887.0 $n [1] 9 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1321.753 1297.972 1510.359 1208.077 1329.94 1441.061 1232.798 1226.658 [2,] 1588.247 1856.028 1622.641 1463.923 1717.06 1844.939 1489.202 1676.342 [,9] [,10] [,11] [,12] [1,] 1432.388 1412.757 1432.037 1391.943 [2,] 1777.612 1765.243 1864.963 1717.057 $out [1] 1781 2570 $group [1] 7 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(1299, 1418, 1455, 1671, 2027, 1315, 1450, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/50i231226313061.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] [1,] -286.0 -581.0 -573.0 -98.0 -343 -464.0 -200.0 -60.0 -97.0 -139.0 -432.0 [2,] -149.0 -374.5 -302.0 10.0 -38 -413.5 -43.0 29.5 -89.5 3.5 -221.0 [3,] 87.5 -31.0 -251.5 112.5 108 -274.0 196.5 88.5 22.5 117.5 -62.5 [4,] 324.5 121.0 -73.5 246.0 378 -140.5 260.0 180.0 175.5 245.0 109.0 [5,] 529.0 327.0 126.0 392.0 461 -33.0 469.0 262.0 289.0 275.0 177.0 [,12] [1,] -382 [2,] -189 [3,] -83 [4,] 53 [5,] 261 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -177.0039 -307.7934 -379.1434 -19.33299 -124.3836 -426.5017 27.23985 [2,] 352.0039 245.7934 -123.8566 244.33299 340.3836 -121.4983 365.76015 [,8] [,9] [,10] [,11] [,12] [1,] 4.428539 -125.5328 -17.40537 -246.8427 -218.18467 [2,] 172.571461 170.5328 252.40537 121.8427 52.18467 $out [1] 488 -708 $group [1] 8 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-286, -149, 87.5, 324.5, 529, -581, -374.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6hfee1226313061.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] [1,] 1198.0 1063.0 1063.0 1224 1387.0 1504.0 1493 1261.0 1431 [2,] 1306.5 1123.5 1275.0 1405 1485.0 1731.5 1523 1393.5 1431 [3,] 1414.5 1331.0 1332.0 1546 1590.5 1865.5 1592 1568.0 1431 [4,] 1646.0 1434.5 1523.5 1636 1868.0 2053.5 1629 1692.5 1431 [5,] 1837.0 1642.0 1639.0 1914 2081.0 2082.0 1765 1819.0 1431 $n [1] 12 12 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1259.652 1189.151 1218.657 1440.639 1415.811 1718.634 1543.653 1431.624 [2,] 1569.348 1472.849 1445.343 1651.361 1765.189 2012.366 1640.347 1704.376 [,9] [1,] 1431 [2,] 1431 $out [1] 2570 1335 $group [1] 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(1198, 1306.5, 1414.5, 1646, 1837, 1063, 1123.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/74kj21226313061.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,] 1339.625 1336.00 1299.375 [2,] 1479.312 1453.25 1447.375 [3,] 1572.562 1560.50 1551.562 [4,] 1610.562 1597.00 1602.562 [5,] 1687.875 1648.50 1729.375 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1512.698 1494.935 1480.780 [2,] 1632.427 1626.065 1622.345 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1339.625, 1479.3125, 1572.5625, 1610.5625, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1ls4u1226313061.ps tmp/1ls4u1226313061.png") > system("convert tmp/2eylj1226313061.ps tmp/2eylj1226313061.png") > system("convert tmp/3wr231226313061.ps tmp/3wr231226313061.png") > system("convert tmp/4u2ru1226313061.ps tmp/4u2ru1226313061.png") > system("convert tmp/50i231226313061.ps tmp/50i231226313061.png") > system("convert tmp/6hfee1226313061.ps tmp/6hfee1226313061.png") > system("convert tmp/74kj21226313061.ps tmp/74kj21226313061.png") > > > proc.time() user system elapsed 2.182 1.489 3.335