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(219.20,232.50,235.60,171.00,165.90,187.60,218.20,249.80,256.50,224.90,200.00,182.50,230.30,252.80,270.60,196.90,184.70,202.50,258.20,283.10,268.50,283.80,231.10,212.10,238.50,262.80,245.50,198.20,167.20,184.20,254.90,246.40,264.50,242.40,186.70,254.70,230.10,253.60,228.00,183.80,150.00,178.50,228.40,228.70,236.70,218.20,173.50,189.10,194.60,213.70,216.30,173.90,156.90,182.90,216.40,234.00,257.30,225.70,201.70,189.20) > 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,] 219.2 230.3 238.5 230.1 194.6 NA [2,] 232.5 252.8 262.8 253.6 213.7 NA [3,] 235.6 270.6 245.5 228.0 216.3 NA [4,] 171.0 196.9 198.2 183.8 173.9 NA [5,] 165.9 184.7 167.2 150.0 156.9 NA [6,] 187.6 202.5 184.2 178.5 182.9 NA [7,] 218.2 258.2 254.9 228.4 216.4 NA [8,] 249.8 283.1 246.4 228.7 234.0 NA [9,] 256.5 268.5 264.5 236.7 257.3 NA [10,] 224.9 283.8 242.4 218.2 225.7 NA [11,] 200.0 231.1 186.7 173.5 201.7 NA [12,] 182.5 212.1 254.7 189.1 189.2 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 13.3 22.5 24.3 23.5 19.1 NA [2,] 3.1 17.8 -17.3 -25.6 2.6 NA [3,] -64.6 -73.7 -47.3 -44.2 -42.4 NA [4,] -5.1 -12.2 -31.0 -33.8 -17.0 NA [5,] 21.7 17.8 17.0 28.5 26.0 NA [6,] 30.6 55.7 70.7 49.9 33.5 NA [7,] 31.6 24.9 -8.5 0.3 17.6 NA [8,] 6.7 -14.6 18.1 8.0 23.3 NA [9,] -31.6 15.3 -22.1 -18.5 -31.6 NA [10,] -24.9 -52.7 -55.7 -44.7 -24.0 NA [11,] -17.5 -19.0 68.0 15.6 -12.5 NA [12,] 47.8 26.4 -24.6 5.5 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/wessaorg/rcomp/tmp/1qjp71331398140.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/2ck891331398140.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/3qzo51331398140.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/42zry1331398140.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] [,12] [1,] 219.2 213.7 216.3 171.0 150.0 178.5 216.4 228.7 256.5 218.2 173.5 182.5 [2,] 219.2 232.5 228.0 173.9 156.9 182.9 218.2 234.0 256.5 224.9 186.7 189.1 [3,] 230.1 252.8 235.6 183.8 165.9 184.2 228.4 246.4 257.3 225.7 200.0 189.2 [4,] 230.3 253.6 245.5 196.9 167.2 187.6 254.9 249.8 264.5 242.4 201.7 212.1 [5,] 238.5 262.8 270.6 198.2 167.2 187.6 258.2 249.8 268.5 242.4 201.7 212.1 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 222.2568 237.8908 223.2345 167.5483 158.622 180.879 202.4679 235.2358 [2,] 237.9432 267.7092 247.9655 200.0517 173.178 187.521 254.3321 257.5642 [,9] [,10] [,11] [,12] [1,] 251.6472 213.3345 189.401 172.9483 [2,] 262.9528 238.0655 210.599 205.4517 $out [1] 194.6 184.7 202.5 283.1 236.7 283.8 231.1 254.7 $group [1] 1 5 6 8 9 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(219.2, 219.2, 230.1, 230.3, 238.5, 213.7, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5dog41331398140.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] [,12] [1,] 13.3 -25.6 -73.7 -33.8 17.0 30.6 -8.5 6.7 -31.6 -55.7 -19.0 -24.60 [2,] 19.1 -17.3 -64.6 -31.0 17.8 33.5 0.3 6.7 -31.6 -52.7 -17.5 -9.55 [3,] 22.5 2.6 -47.3 -17.0 21.7 49.9 17.6 8.0 -22.1 -44.7 -12.5 15.95 [4,] 23.5 3.1 -44.2 -12.2 26.0 55.7 24.9 18.1 -18.5 -24.9 15.6 37.10 [5,] 24.3 17.8 -42.4 -5.1 28.5 70.7 31.6 23.3 -18.5 -24.0 15.6 47.80 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 19.39097 -11.81459 -61.71459 -30.284033 15.9059 34.21354 0.217702 [2,] 25.60903 17.01459 -32.88541 -3.715967 27.4941 65.58646 34.982298 [,8] [,9] [,10] [,11] [,12] [1,] -0.05521128 -31.35643 -64.34341 -35.88838 -20.9035 [2,] 16.05521128 -12.84357 -25.05659 10.88838 52.8035 $out [1] -14.6 15.3 68.0 $group [1] 8 9 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(13.3, 19.1, 22.5, 23.5, 24.3, -25.6, -17.3, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6mtyo1331398140.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] [1,] 165.90 184.70 167.20 150.00 156.90 NA [2,] 185.05 207.30 192.45 181.15 186.05 NA [3,] 218.70 241.95 243.95 223.10 207.70 NA [4,] 234.05 269.55 254.80 229.40 221.05 NA [5,] 256.50 283.80 264.50 253.60 257.30 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 196.3508 213.5574 215.5117 201.0929 191.7363 NA [2,] 241.0492 270.3426 272.3883 245.1071 223.6637 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(165.9, 185.05, 218.7, 234.05, 256.5, 184.7, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7kqhe1331398140.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,] 164.94 165.90 162.050 [2,] 192.87 186.70 189.800 [3,] 228.88 227.05 229.200 [4,] 241.14 241.00 239.325 [5,] 256.70 257.30 260.500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 206.8637 202.2834 206.6113 [2,] 250.8963 251.8166 251.7887 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(164.94, 192.87, 228.88, 241.14, 256.7, 165.9, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1qjp71331398140.ps tmp/1qjp71331398140.png",intern=TRUE)) character(0) > try(system("convert tmp/2ck891331398140.ps tmp/2ck891331398140.png",intern=TRUE)) character(0) > try(system("convert tmp/3qzo51331398140.ps tmp/3qzo51331398140.png",intern=TRUE)) character(0) > try(system("convert tmp/42zry1331398140.ps tmp/42zry1331398140.png",intern=TRUE)) character(0) > try(system("convert tmp/5dog41331398140.ps tmp/5dog41331398140.png",intern=TRUE)) character(0) > try(system("convert tmp/6mtyo1331398140.ps tmp/6mtyo1331398140.png",intern=TRUE)) character(0) > try(system("convert tmp/7kqhe1331398140.ps tmp/7kqhe1331398140.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.589 0.393 1.991