R version 2.12.0 (2010-10-15) Copyright (C) 2010 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(9700,9081,9084,9743,8587,9731,9563,9998,9437,10038,9918,9252,9737,9035,9133,9487,8700,9627,8947,9283,8829,9947,9628,9318,9605,8640,9214,9567,8547,9185,9470,9123,9278,10170,9434,9655,9429,8739,9552,9687,9019,9672,9206,9069,9788,10312,10105,9863,9656,9295,9946,9701,9049,10190,9706,9765,9893,9994,10433,10073,10112,9266,9820,10097,9115,10411,9678,10408,10153,10368,10581,10597,10680,9738,9556) > 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] 75 > (np <- floor(n / par1)) [1] 6 > 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] 7 7 7 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 9700 9737 9605 9429 9656 10112 10680 [2,] 9081 9035 8640 8739 9295 9266 9738 [3,] 9084 9133 9214 9552 9946 9820 9556 [4,] 9743 9487 9567 9687 9701 10097 NA [5,] 8587 8700 8547 9019 9049 9115 NA [6,] 9731 9627 9185 9672 10190 10411 NA [7,] 9563 8947 9470 9206 9706 9678 NA [8,] 9998 9283 9123 9069 9765 10408 NA [9,] 9437 8829 9278 9788 9893 10153 NA [10,] 10038 9947 10170 10312 9994 10368 NA [11,] 9918 9628 9434 10105 10433 10581 NA [12,] 9252 9318 9655 9863 10073 10597 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -619 -702 -965 -690 -361 -846 -942 [2,] 3 98 574 813 651 554 -182 [3,] 659 354 353 135 -245 277 NA [4,] -1156 -787 -1020 -668 -652 -982 NA [5,] 1144 927 638 653 1141 1296 NA [6,] -168 -680 285 -466 -484 -733 NA [7,] 435 336 -347 -137 59 730 NA [8,] -561 -454 155 719 128 -255 NA [9,] 601 1118 892 524 101 215 NA [10,] -120 -319 -736 -207 439 213 NA [11,] -666 -310 221 -242 -360 16 NA [12,] 485 287 -226 -207 39 83 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/rcomp/tmp/1vyor1322658331.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/www/rcomp/tmp/2k9ex1322658331.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/www/rcomp/tmp/3uo0x1322658331.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/www/rcomp/tmp/4rs651322658331.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] [1,] 9429.0 8640.0 9084.0 9487 8547.0 9185.0 8947.0 9069 8829.0 9947 [2,] 9630.5 8887.0 9173.5 9567 8587.0 9627.0 9206.0 9123 9278.0 9994 [3,] 9700.0 9081.0 9552.0 9694 8859.5 9701.5 9516.5 9524 9612.5 10104 [4,] 9924.5 9280.5 9688.0 9743 9049.0 10190.0 9678.0 9998 9893.0 10312 [5,] 10112.0 9738.0 9946.0 9743 9115.0 10411.0 9706.0 10408 10153.0 10368 [,11] [,12] [1,] 9434.0 9252 [2,] 9628.0 9318 [3,] 10011.5 9759 [4,] 10433.0 10073 [5,] 10581.0 10597 $n [1] 7 7 7 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 9524.428 8846.008 9244.749 9580.474 8561.495 9338.347 9212.045 8959.597 [2,] 9875.572 9315.992 9859.251 9807.526 9157.505 10064.653 9820.955 10088.403 [,9] [,10] [,11] [,12] [1,] 9215.805 9898.88 9492.249 9272 [2,] 10009.195 10309.12 10530.751 10246 $out [1] 10680 10097 $group [1] 1 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(9429, 9630.5, 9700, 9924.5, 10112, 8640, 8887, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5p4ev1322658331.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,] -965.0 -182.0 135 -1156.0 638 -733 -347.0 -561.0 101.0 -736.0 -666 [2,] -894.0 50.5 135 -1020.0 653 -680 -137.0 -454.0 215.0 -319.0 -360 [3,] -702.0 554.0 315 -884.5 1034 -475 197.5 -63.5 562.5 -163.5 -276 [4,] -654.5 612.5 354 -668.0 1144 -168 435.0 155.0 892.0 213.0 16 [5,] -361.0 813.0 659 -652.0 1296 285 730.0 719.0 1118.0 439.0 221 [,12] [1,] -226 [2,] -207 [3,] 61 [4,] 287 [5,] 485 $n [1] 7 7 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -845.0255 218.3827 173.7379 -1111.5514 717.2891 -805.2565 -171.4585 [2,] -558.9745 889.6173 456.2621 -657.4486 1350.7109 -144.7435 566.4585 [,8] [,9] [,10] [,11] [,12] [1,] -456.3247 125.8131 -506.6572 -518.53214 -257.6460 [2,] 329.3247 999.1869 179.6572 -33.46786 379.6460 $out [1] -245 $group [1] 3 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-965, -894, -702, -654.5, -361, -182, 50.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6xo2w1322658331.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] [1,] 8587.0 8700.0 8547 8739.0 9295.0 9115.0 9556 [2,] 9168.0 8991.0 9154 9137.5 9678.5 9749.0 9647 [3,] 9631.5 9300.5 9356 9612.0 9829.0 10132.5 9738 [4,] 9830.5 9627.5 9586 9825.5 10033.5 10409.5 10209 [5,] 10038.0 9947.0 10170 10312.0 10433.0 10597.0 10680 $n [1] 12 12 12 12 12 12 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 9329.33 9010.188 9158.962 9298.199 9667.082 9831.242 9225.336 [2,] 9933.67 9590.812 9553.038 9925.801 9990.918 10433.758 10250.664 $out [1] 9049 $group [1] 5 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(8587, 9168, 9631.5, 9830.5, 10038, 8700, 8991, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/78vef1322658331.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,] 9113.429 9516.50 9083.750 [2,] 9450.238 9520.25 9445.688 [3,] 9660.667 9653.25 9628.500 [4,] 9824.119 9730.25 9817.125 [5,] 10138.167 10011.50 10140.750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 9490.137 9557.468 9459.085 [2,] 9831.196 9749.032 9797.915 $out [1] 8836.167 9081.000 8859.500 10104.000 8828.375 $group [1] 1 2 2 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(9113.42857142857, 9450.2380952381, 9660.66666666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1vyor1322658331.ps tmp/1vyor1322658331.png",intern=TRUE)) character(0) > try(system("convert tmp/2k9ex1322658331.ps tmp/2k9ex1322658331.png",intern=TRUE)) character(0) > try(system("convert tmp/3uo0x1322658331.ps tmp/3uo0x1322658331.png",intern=TRUE)) character(0) > try(system("convert tmp/4rs651322658331.ps tmp/4rs651322658331.png",intern=TRUE)) character(0) > try(system("convert tmp/5p4ev1322658331.ps tmp/5p4ev1322658331.png",intern=TRUE)) character(0) > try(system("convert tmp/6xo2w1322658331.ps tmp/6xo2w1322658331.png",intern=TRUE)) character(0) > try(system("convert tmp/78vef1322658331.ps tmp/78vef1322658331.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.680 0.260 2.946