R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-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. 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(426.113 + ,383.703 + ,232.444 + ,70.939 + ,226.731 + ,947.293 + ,611.281 + ,158.047 + ,33.999 + ,37.028 + ,388.3 + ,506.652 + ,392.25 + ,180.818 + ,198.296 + ,217.465 + ,275.562 + ,1030.944 + ,57.47 + ,136.452 + ,556.277 + ,213.361 + ,274.482 + ,220.553 + ,236.71 + ,260.642 + ,2763.544 + ,213.923 + ,169.861 + ,403.064 + ,449.594 + ,406.167 + ,206.893 + ,156.187 + ,257.102 + ,62.156 + ,662.883 + ,251.422 + ,171.328 + ,350.089 + ,221.588 + ,4.813 + ,183.186 + ,190.379 + ,223.166 + ,232.669 + ,356.725 + ,109.215 + ,475.834 + ,315.955 + ,694.87 + ,8.95 + ,278.741 + ,308.16 + ,207.533 + ,192.797 + ,601.162 + ,289.714 + ,293.671 + ,386.688 + ,699.645 + ,85.094 + ,131.812 + ,645.285 + ,197.549 + ,308.174 + ,86.58 + ,242.205 + ,238.502 + ,187.881 + ,140.321 + ,440.31 + ,421.403 + ,218.761 + ,1305.923 + ,137.55 + ,262.517 + ,348.821 + ,150.034 + ,64.016 + ,261.596 + ,259.7 + ,171.26 + ,203.077 + ,249.148 + ,211.655 + ,252.64 + ,438.555 + ,239.89 + ,401.915 + ,216.886 + ,184.641 + ,380.155 + ,653.641 + ,313.906 + ,366.936 + ,236.302 + ,229.641 + ,235.577 + ,103.898 + ,263.906 + ,241.171 + ,216.548 + ,295.281 + ,193.299 + ,204.386 + ,257.567 + ,136.813 + ,240.755 + ,59.609 + ,213.511 + ,380.531 + ,242.344 + ,250.407 + ,183.613 + ,191.835 + ,266.793 + ,246.542 + ,330.563 + ,403.556 + ,208.108 + ,324.04 + ,308.532 + ,199.297 + ,200.156 + ,262.875 + ,287.069 + ,190.157 + ,199.746 + ,265.777 + ,435.956 + ,72.844 + ,756.46 + ,206.771 + ,4202.361 + ,401.422 + ,216.046 + ,39.047 + ,441.437) > 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] 139 > (np <- floor(n / par1)) [1] 11 > 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] 12 12 12 12 12 12 12 11 11 11 11 11 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 426.113 392.250 236.710 662.883 475.834 699.645 421.403 249.148 [2,] 383.703 180.818 260.642 251.422 315.955 85.094 218.761 211.655 [3,] 232.444 198.296 2763.544 171.328 694.870 131.812 1305.923 252.640 [4,] 70.939 217.465 213.923 350.089 8.950 645.285 137.550 438.555 [5,] 226.731 275.562 169.861 221.588 278.741 197.549 262.517 239.890 [6,] 947.293 1030.944 403.064 4.813 308.160 308.174 348.821 401.915 [7,] 611.281 57.470 449.594 183.186 207.533 86.580 150.034 216.886 [8,] 158.047 136.452 406.167 190.379 192.797 242.205 64.016 184.641 [9,] 33.999 556.277 206.893 223.166 601.162 238.502 261.596 380.155 [10,] 37.028 213.361 156.187 232.669 289.714 187.881 259.700 653.641 [11,] 388.300 274.482 257.102 356.725 293.671 140.321 171.260 313.906 [12,] 506.652 220.553 62.156 109.215 386.688 440.310 203.077 366.936 [,9] [,10] [,11] [,12] [1,] 236.302 240.755 208.108 756.460 [2,] 229.641 59.609 324.040 206.771 [3,] 235.577 213.511 308.532 4202.361 [4,] 103.898 380.531 199.297 401.422 [5,] 263.906 242.344 200.156 216.046 [6,] 241.171 250.407 262.875 39.047 [7,] 216.548 183.613 287.069 441.437 [8,] 295.281 191.835 190.157 NA [9,] 193.299 266.793 199.746 NA [10,] 204.386 246.542 265.777 NA [11,] 257.567 330.563 435.956 NA [12,] 136.813 403.556 72.844 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -42.410 -211.432 23.932 -411.461 -159.879 -614.551 -202.642 -37.493 [2,] -151.259 17.478 2502.902 -80.094 378.915 46.718 1087.162 40.985 [3,] -161.505 19.169 -2549.621 178.761 -685.920 513.473 -1168.373 185.915 [4,] 155.792 58.097 -44.062 -128.501 269.791 -447.736 124.967 -198.665 [5,] 720.562 755.382 233.203 -216.775 29.419 110.625 86.304 162.025 [6,] -336.012 -973.474 46.530 178.373 -100.627 -221.594 -198.787 -185.029 [7,] -453.234 78.982 -43.427 7.193 -14.736 155.625 -86.018 -32.245 [8,] -124.048 419.825 -199.274 32.787 408.365 -3.703 197.580 195.514 [9,] 3.029 -342.916 -50.706 9.503 -311.448 -50.621 -1.896 273.486 [10,] 351.272 61.121 100.915 124.056 3.957 -47.560 -88.440 -339.735 [11,] 118.352 -53.929 -194.946 -247.510 93.017 299.989 31.817 53.030 [12,] -114.402 16.157 600.727 366.619 312.957 -18.907 46.071 -130.634 [,9] [,10] [,11] [,12] [1,] -6.661 -181.146 115.932 -549.689 [2,] 5.936 153.902 -15.508 3995.590 [3,] -131.679 167.020 -109.235 -3800.939 [4,] 160.008 -138.187 0.859 -185.376 [5,] -22.735 8.063 62.719 -176.999 [6,] -24.623 -66.794 24.194 402.390 [7,] 78.733 8.222 -96.912 NA [8,] -101.982 74.958 9.589 NA [9,] 11.087 -20.251 66.031 NA [10,] 53.181 84.021 170.179 NA [11,] -120.754 72.993 -363.112 NA [12,] 103.942 -195.448 683.616 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/wessaorg/rcomp/tmp/11ymu1322575900.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/wessaorg/rcomp/tmp/2d5dy1322575900.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/wessaorg/rcomp/tmp/321gg1322575900.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/wessaorg/rcomp/tmp/4f34h1322575900.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] [1,] 208.1080 59.6090 131.8120 8.9500 169.8610 39.0470 57.4700 136.452 [2,] 238.7325 193.7945 205.9035 120.7240 208.1010 245.7890 166.6100 171.344 [3,] 406.8265 224.2010 244.1085 215.6940 233.3105 308.1670 212.0405 190.379 [4,] 569.3585 288.2985 1000.3965 390.9765 263.2115 402.4895 364.2530 217.501 [5,] 756.4600 383.7030 1305.9230 645.2850 278.7410 403.0640 611.2810 242.205 [,9] [,10] [,11] [,12] [1,] 33.9990 156.1870 140.3210 62.156 [2,] 203.3195 196.1335 257.3345 123.014 [3,] 238.5020 232.6690 293.6710 220.553 [4,] 323.4740 262.7385 343.6440 395.122 [5,] 380.1550 289.7140 435.9560 506.652 $n [1] 12 12 12 12 12 12 12 11 11 11 11 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 256.0258 181.0971 -118.2651 92.43002 208.1742 236.6949 121.8942 168.3904 [2,] 557.6272 267.3049 606.4821 338.95798 258.4468 379.6391 302.1868 212.3676 [,9] [,10] [,11] [,12] [1,] 181.2618 200.9392 252.5542 90.92403 [2,] 295.7422 264.3988 334.7878 350.18197 $out [1] 2763.544 4202.361 947.293 1030.944 4.813 406.167 64.016 295.281 [9] 556.277 601.162 37.028 653.641 $group [1] 3 3 6 6 6 8 8 8 9 9 10 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(208.108, 238.7325, 406.8265, 569.3585, 756.46, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/54ksg1322575900.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] [1,] -614.5510 -151.2590 -2549.6210 -447.7360 -216.7750 -336.0120 -96.9120 [2,] -311.4465 -4.7860 -927.1465 -161.7815 -7.3360 -210.1905 -64.7225 [3,] -170.5125 43.8515 -120.4570 -21.6015 74.5115 -83.7105 -14.7360 [4,] -22.0770 733.0385 172.8905 140.3795 197.6140 35.3620 43.4775 [5,] 115.9320 1087.1620 513.4730 269.7910 233.2030 402.3900 155.6250 [,8] [,9] [,10] [,11] [,12] [1,] -199.2740 -50.7060 -88.4400 -363.112 -195.4480 [2,] -52.8425 -50.6635 -21.8015 -157.850 -66.6545 [3,] 32.7870 -1.8960 61.1210 31.817 46.0710 [4,] 196.5470 10.2950 112.4855 83.005 339.7880 [5,] 419.8250 66.0310 170.1790 299.989 683.6160 $n [1] 12 12 12 12 12 12 11 11 11 11 11 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -302.49587 -292.6752 -622.1913 -159.4192 -18.96757 -195.70864 -66.28117 [2,] -38.52913 380.3782 381.2773 116.2162 167.99057 28.28764 36.80917 [,8] [,9] [,10] [,11] [,12] [1,] -86.01915 -30.93589 -2.851705 -82.92341 -147.5533 [2,] 151.59315 27.14389 125.093705 146.55741 239.6953 $out [1] 2502.902 3995.590 -3800.939 720.562 755.382 -973.474 -453.234 [8] -342.916 -311.448 273.486 351.272 -339.735 $group [1] 2 2 3 5 5 6 7 9 9 9 10 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-614.551, -311.4465, -170.5125, -22.077, 115.932, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6ym8k1322575900.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] [1,] 33.9990 57.470 62.1560 4.8130 8.9500 85.0940 64.0160 184.641 [2,] 114.4930 189.557 188.3770 177.2570 243.1370 136.0665 160.6470 228.388 [3,] 308.0735 219.009 246.9060 222.3770 300.9155 218.0255 239.2305 283.273 [4,] 466.3825 333.906 404.6155 300.7555 431.2610 374.2420 305.6690 391.035 [5,] 947.2930 392.250 449.5940 356.7250 694.8700 699.6450 421.4030 438.555 [,9] [,10] [,11] [,12] [1,] 136.8130 59.609 72.8440 39.0470 [2,] 198.8425 202.673 199.5215 211.4085 [3,] 232.6090 244.443 235.4915 401.4220 [4,] 249.3690 298.678 297.8005 598.9485 [5,] 295.2810 403.556 435.9560 756.4600 $n [1] 12 12 12 12 12 12 12 12 12 12 12 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 147.5743 153.1705 148.2782 166.0485 215.1109 109.3921 173.085 209.0886 [2,] 468.5727 284.8475 345.5338 278.7055 386.7201 326.6589 305.376 357.4574 [,9] [,10] [,11] [,12] [1,] 209.5635 200.6545 190.6658 169.9894 [2,] 255.6545 288.2315 280.3172 632.8546 $out [1] 1030.944 556.277 2763.544 662.883 1305.923 653.641 103.898 4202.361 $group [1] 2 2 3 4 7 8 9 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" NA Warning message: In bxp(list(stats = c(33.999, 114.493, 308.0735, 466.3825, 947.293, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7c2y31322575900.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,] 204.7252 190.3790 194.4225 [2,] 241.3122 218.1235 237.4227 [3,] 264.2142 232.9898 258.2567 [4,] 335.8021 268.8897 312.8197 [5,] 417.1343 308.1670 381.1700 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 221.1167 209.8349 223.8676 [2,] 307.3116 256.1446 292.6458 $out [1] 892.5698 406.8265 528.6703 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(204.725181818182, 241.312246212121, 264.214181818182, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/11ymu1322575900.ps tmp/11ymu1322575900.png",intern=TRUE)) character(0) > try(system("convert tmp/2d5dy1322575900.ps tmp/2d5dy1322575900.png",intern=TRUE)) character(0) > try(system("convert tmp/321gg1322575900.ps tmp/321gg1322575900.png",intern=TRUE)) character(0) > try(system("convert tmp/4f34h1322575900.ps tmp/4f34h1322575900.png",intern=TRUE)) character(0) > try(system("convert tmp/54ksg1322575900.ps tmp/54ksg1322575900.png",intern=TRUE)) character(0) > try(system("convert tmp/6ym8k1322575900.ps tmp/6ym8k1322575900.png",intern=TRUE)) character(0) > try(system("convert tmp/7c2y31322575900.ps tmp/7c2y31322575900.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.650 0.338 2.068