R version 2.8.0 (2008-10-20) 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. 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(26862,28261,31761,32624,32950,34711,34124,31860,32404,33027,28944,29702,29144,29715,35710,37324,37687,39071,38986,37702,36444,36423,35266,35461,33073,35721,44670,43461,42998,45039,41555,43850,42253,42036,39552,38572,36862,37883,47264,43698,48335,50368,45633,48708,44588,44021,41601,39761,40848,45488,51253,49495,53348,50286,50191,50596,46689,49249,43597,41862) > 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] 60 > (np <- floor(n / par1)) [1] 5 > 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] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 26862 29144 33073 36862 40848 NA [2,] 28261 29715 35721 37883 45488 NA [3,] 31761 35710 44670 47264 51253 NA [4,] 32624 37324 43461 43698 49495 NA [5,] 32950 37687 42998 48335 53348 NA [6,] 34711 39071 45039 50368 50286 NA [7,] 34124 38986 41555 45633 50191 NA [8,] 31860 37702 43850 48708 50596 NA [9,] 32404 36444 42253 44588 46689 NA [10,] 33027 36423 42036 44021 49249 NA [11,] 28944 35266 39552 41601 43597 NA [12,] 29702 35461 38572 39761 41862 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1399 571 2648 1021 4640 NA [2,] 3500 5995 8949 9381 5765 NA [3,] 863 1614 -1209 -3566 -1758 NA [4,] 326 363 -463 4637 3853 NA [5,] 1761 1384 2041 2033 -3062 NA [6,] -587 -85 -3484 -4735 -95 NA [7,] -2264 -1284 2295 3075 405 NA [8,] 544 -1258 -1597 -4120 -3907 NA [9,] 623 -21 -217 -567 2560 NA [10,] -4083 -1157 -2484 -2420 -5652 NA [11,] 758 195 -980 -1840 -1735 NA [12,] -558 -2388 -1710 1087 NA 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/11yey1289060543.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/21yey1289060543.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/31yey1289060543.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/4bpdj1289060543.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] [,12] [1,] 26862 28261 31761 32624 32950 34711 34124 31860 32404 33027 28944 29702 [2,] 29144 29715 35710 37324 37687 39071 38986 37702 36444 36423 35266 35461 [3,] 33073 35721 44670 43461 42998 45039 41555 43850 42253 42036 39552 38572 [4,] 36862 37883 47264 43698 48335 50286 45633 48708 44588 44021 41601 39761 [5,] 40848 45488 51253 49495 53348 50368 50191 50596 46689 49249 43597 41862 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 27619.48 29949.51 36505.97 38957.15 35474.15 37114.51 36858.25 36073.19 [2,] 38526.52 41492.49 52834.03 47964.85 50521.85 52963.49 46251.75 51626.81 [,9] [,10] [,11] [,12] [1,] 36498.47 36667.27 35075.70 35533.63 [2,] 48007.53 47404.73 44028.30 41610.37 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(26862, 29144, 33073, 36862, 40848, 28261, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5bpdj1289060543.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] [,12] [1,] 571 3500 -3566 -463 1384 -4735 -2264 -4120 -567 -5652 -1840 -2388.0 [2,] 1021 5765 -1758 326 1384 -3484 -1284 -3907 -217 -4083 -1735 -2049.0 [3,] 1399 5995 -1209 363 1761 -587 405 -1597 -21 -2484 -980 -1134.0 [4,] 2648 8949 863 3853 2033 -95 2295 -1258 623 -2420 195 264.5 [5,] 4640 9381 1614 4637 2041 -85 3075 544 623 -1157 758 1087.0 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 249.3659 3745.194 -3060.992 -2129.169 1302.418 -2981.659 -2123.912 [2,] 2548.6341 8244.806 642.992 2855.169 2219.582 1807.659 2933.912 [,8] [,9] [,10] [,11] [,12] [1,] -3468.7767 -614.5419 -3659.072 -2343.7331 -2961.665 [2,] 274.7767 572.5419 -1308.928 383.7331 693.665 $out [1] -3062 2560 $group [1] 5 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(571, 1021, 1399, 2648, 4640, 3500, 5765, 5995, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/64zd41289060543.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] [1,] 26862.0 35266.0 33073.0 36862.0 40848.0 NA [2,] 29323.0 35363.5 39062.0 40681.0 44542.5 NA [3,] 32132.0 36433.5 42144.5 44304.5 49372.0 NA [4,] 32988.5 37694.5 43655.5 47799.5 50441.0 NA [5,] 34711.0 39071.0 45039.0 50368.0 53348.0 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 30460.14 35370.32 40049.37 41057.70 46681.65 NA [2,] 33803.86 37496.68 44239.63 47551.30 52062.35 NA $out [1] 29144 29715 $group [1] 2 2 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(26862, 29323, 32132, 32988.5, 34711, 35266, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/74zd41289060543.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,] 33357.8 33073.0 33003.00 [2,] 37431.8 39062.0 38022.25 [3,] 41135.8 42144.5 40513.50 [4,] 42337.4 43655.5 42660.25 [5,] 43895.0 45039.0 44678.50 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 38898.32 40049.37 38398.08 [2,] 43373.28 44239.63 42628.92 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(33357.8, 37431.8, 41135.8, 42337.4, 43895, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > try(system("convert tmp/11yey1289060543.ps tmp/11yey1289060543.png",intern=TRUE)) character(0) > try(system("convert tmp/21yey1289060543.ps tmp/21yey1289060543.png",intern=TRUE)) character(0) > try(system("convert tmp/31yey1289060543.ps tmp/31yey1289060543.png",intern=TRUE)) character(0) > try(system("convert tmp/4bpdj1289060543.ps tmp/4bpdj1289060543.png",intern=TRUE)) character(0) > try(system("convert tmp/5bpdj1289060543.ps tmp/5bpdj1289060543.png",intern=TRUE)) character(0) > try(system("convert tmp/64zd41289060543.ps tmp/64zd41289060543.png",intern=TRUE)) character(0) > try(system("convert tmp/74zd41289060543.ps tmp/74zd41289060543.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.932 1.468 2.166