R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(278.1,283.3,287,287.9,279.7,273.2,272.6,274.8,277.2,277.5,278.7,279,280.1,280.8,280.1,270.5,260,258.2,260.1,262.2,263.4,264.7,267.4,266.8,265.2,262.5,255,252.8,249.6,254.2,259.4,260.8,261.6,257.2,249.8,235.2,227.4,223.1,221.6,218.5,220.1,222.1,222.9,222.8,222.1,209,201.7,204.5,208.7,213,218.6,221.9,224.6,226.8,229.6,231.4,229.1,232.7,236.6,242.2,251.5,261.5,266.4,280.9,294.4,308.8,328.8,348.2,365.9,387.5,395.9,395.8) > 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] 72 > (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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 278.1 280.1 265.2 227.4 208.7 251.5 NA [2,] 283.3 280.8 262.5 223.1 213.0 261.5 NA [3,] 287.0 280.1 255.0 221.6 218.6 266.4 NA [4,] 287.9 270.5 252.8 218.5 221.9 280.9 NA [5,] 279.7 260.0 249.6 220.1 224.6 294.4 NA [6,] 273.2 258.2 254.2 222.1 226.8 308.8 NA [7,] 272.6 260.1 259.4 222.9 229.6 328.8 NA [8,] 274.8 262.2 260.8 222.8 231.4 348.2 NA [9,] 277.2 263.4 261.6 222.1 229.1 365.9 NA [10,] 277.5 264.7 257.2 209.0 232.7 387.5 NA [11,] 278.7 267.4 249.8 201.7 236.6 395.9 NA [12,] 279.0 266.8 235.2 204.5 242.2 395.8 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 5.2 0.7 -2.7 -4.3 4.3 10.0 NA [2,] 3.7 -0.7 -7.5 -1.5 5.6 4.9 NA [3,] 0.9 -9.6 -2.2 -3.1 3.3 14.5 NA [4,] -8.2 -10.5 -3.2 1.6 2.7 13.5 NA [5,] -6.5 -1.8 4.6 2.0 2.2 14.4 NA [6,] -0.6 1.9 5.2 0.8 2.8 20.0 NA [7,] 2.2 2.1 1.4 -0.1 1.8 19.4 NA [8,] 2.4 1.2 0.8 -0.7 -2.3 17.7 NA [9,] 0.3 1.3 -4.4 -13.1 3.6 21.6 NA [10,] 1.2 2.7 -7.4 -7.3 3.9 8.4 NA [11,] 0.3 -0.6 -14.6 2.8 5.6 -0.1 NA [12,] 1.1 -1.6 -7.8 4.2 9.3 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/rcomp/tmp/1jmbm1275814897.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/rcomp/tmp/2uvap1275814897.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/rcomp/tmp/3uvap1275814897.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/rcomp/tmp/4nn9s1275814897.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,] 208.70 213.0 218.6 218.50 220.1 222.1 222.90 222.8 222.1 209.00 201.7 [2,] 227.40 223.1 221.6 221.90 224.6 226.8 229.60 231.4 229.1 232.70 236.6 [3,] 258.35 262.0 260.7 261.65 254.8 256.2 259.75 261.5 262.5 260.95 258.6 [4,] 278.10 280.8 280.1 280.90 279.7 273.2 272.60 274.8 277.2 277.50 278.7 [5,] 280.10 283.3 287.0 287.90 294.4 308.8 328.80 274.8 277.2 277.50 278.7 [,12] [1,] 204.5 [2,] 235.2 [3,] 254.5 [4,] 279.0 [5,] 279.0 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 225.6469 224.7816 222.9656 223.5931 219.2587 226.2705 232.0136 233.5056 [2,] 291.0531 299.2184 298.4344 299.7069 290.3413 286.1295 287.4864 289.4944 [,9] [,10] [,11] [,12] [1,] 231.4739 232.0526 231.4441 226.2476 [2,] 293.5261 289.8474 285.7559 282.7524 $out [1] 348.2 365.9 387.5 395.9 395.8 $group [1] 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(208.7, 227.4, 258.35, 278.1, 280.1, 213, 223.1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5nn9s1275814897.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,] -4.3 -7.5 -9.60 -10.5 -6.5 -0.60 1.40 -2.3 -13.1 -7.40 -0.6 -7.8 [2,] -2.7 -1.5 -3.10 -8.2 -1.8 0.80 1.40 -0.7 -4.4 -7.30 -0.6 -1.6 [3,] 2.5 1.5 -0.65 -0.8 2.1 2.35 1.95 1.0 0.8 1.95 0.1 1.1 [4,] 5.2 4.9 3.30 2.7 4.6 5.20 2.20 2.4 3.6 3.90 2.8 4.2 [5,] 10.0 5.6 3.30 13.5 4.6 5.20 2.20 2.4 3.6 8.40 5.6 9.3 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2.595755 -2.628207 -4.778207 -7.830852 -2.028207 -0.4881421 1.433974 [2,] 7.595755 5.628207 3.478207 6.230852 6.228207 5.1881421 2.466026 [,8] [,9] [,10] [,11] [,12] [1,] -0.9996001 -4.360258 -5.274362 -2.09311 -2.998265 [2,] 2.9996001 5.960258 9.174362 2.29311 5.198265 $out [1] 14.5 14.4 20.0 -0.1 19.4 17.7 21.6 -14.6 $group [1] 3 5 6 7 7 8 9 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-4.30000000000001, -2.69999999999999, 2.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6nn9s1275814897.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] [1,] 272.6 258.20 249.6 201.70 208.70 251.50 NA [2,] 276.0 261.15 251.3 213.75 220.25 273.65 NA [3,] 278.4 265.75 256.1 221.85 227.95 318.80 NA [4,] 281.5 275.30 261.2 222.85 232.05 376.70 NA [5,] 287.9 280.80 265.2 227.40 242.20 395.90 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 275.8914 259.2961 251.5845 217.6994 222.5679 271.7982 NA [2,] 280.9086 272.2039 260.6155 226.0006 233.3321 365.8018 NA $out [1] 235.2 $group [1] 3 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(272.6, 276, 278.4, 281.5, 287.9, 258.2, 261.15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7ywqd1275814897.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,] 251.8333 254.500 251.5500 [2,] 254.7583 257.275 253.2875 [3,] 259.7250 260.225 254.3063 [4,] 270.2333 261.575 255.9688 [5,] 271.6833 262.500 257.8875 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 252.6667 258.2637 253.0833 [2,] 266.7833 262.1863 255.5292 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(251.833333333333, 254.758333333333, 259.725, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1jmbm1275814897.ps tmp/1jmbm1275814897.png",intern=TRUE)) character(0) > try(system("convert tmp/2uvap1275814897.ps tmp/2uvap1275814897.png",intern=TRUE)) character(0) > try(system("convert tmp/3uvap1275814897.ps tmp/3uvap1275814897.png",intern=TRUE)) character(0) > try(system("convert tmp/4nn9s1275814897.ps tmp/4nn9s1275814897.png",intern=TRUE)) character(0) > try(system("convert tmp/5nn9s1275814897.ps tmp/5nn9s1275814897.png",intern=TRUE)) character(0) > try(system("convert tmp/6nn9s1275814897.ps tmp/6nn9s1275814897.png",intern=TRUE)) character(0) > try(system("convert tmp/7ywqd1275814897.ps tmp/7ywqd1275814897.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.179 0.943 1.536