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(430,429,428,426,424,423,424,426,427,427,428,430,432,435,426,411,405,403,402,399,392,387,380,379,386,385,365,356,338,338,343,338,320,316,317,315,317,321,303,303,290,285,300,291,278,273,277,269,275,278,255,254,245,240,261,247,229,213,218,206,217,219,196,193,188,171,190,180,149,135,151,134,145,151,137,124,125,109,131,133,103,85,104,82) > 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] 84 > (np <- floor(n / par1)) [1] 7 > 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 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 430 432 386 317 275 217 145 NA [2,] 429 435 385 321 278 219 151 NA [3,] 428 426 365 303 255 196 137 NA [4,] 426 411 356 303 254 193 124 NA [5,] 424 405 338 290 245 188 125 NA [6,] 423 403 338 285 240 171 109 NA [7,] 424 402 343 300 261 190 131 NA [8,] 426 399 338 291 247 180 133 NA [9,] 427 392 320 278 229 149 103 NA [10,] 427 387 316 273 213 135 85 NA [11,] 428 380 317 277 218 151 104 NA [12,] 430 379 315 269 206 134 82 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 3 -1 4 3 2 6 NA [2,] -1 -9 -20 -18 -23 -23 -14 NA [3,] -2 -15 -9 0 -1 -3 -13 NA [4,] -2 -6 -18 -13 -9 -5 1 NA [5,] -1 -2 0 -5 -5 -17 -16 NA [6,] 1 -1 5 15 21 19 22 NA [7,] 2 -3 -5 -9 -14 -10 2 NA [8,] 1 -7 -18 -13 -18 -31 -30 NA [9,] 0 -5 -4 -5 -16 -14 -18 NA [10,] 1 -7 1 4 5 16 19 NA [11,] 2 -1 -2 -8 -12 -17 -22 NA [12,] 2 7 2 6 11 11 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/1xwva1281023871.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/2xwva1281023871.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/3q5cd1281023871.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/4q5cd1281023871.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,] 145 151.0 137.0 124.0 125.0 109.0 131.0 133.0 103 85.0 104.0 82 [2,] 246 248.5 225.5 223.5 216.5 205.5 225.5 213.5 189 174.0 184.5 170 [3,] 317 321.0 303.0 303.0 290.0 285.0 300.0 291.0 278 273.0 277.0 269 [4,] 408 407.0 395.5 383.5 371.5 370.5 372.5 368.5 356 351.5 348.5 347 [5,] 432 435.0 428.0 426.0 424.0 423.0 424.0 426.0 427 427.0 428.0 430 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 220.2562 226.3464 201.4787 207.4506 197.4365 186.4647 212.2140 198.4365 [2,] 413.7438 415.6536 404.5213 398.5494 382.5635 383.5353 387.7860 383.5635 [,9] [,10] [,11] [,12] [1,] 178.2703 166.9999 179.0618 163.2985 [2,] 377.7297 379.0001 374.9382 374.7015 $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(145, 246, 317, 408, 432, 151, 248.5, 321, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5q5cd1281023871.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,] -1.0 -23.0 -15.0 -18.0 -17.0 -1 -14.0 -31 -18.0 -7.0 -22.0 2.0 [2,] 0.5 -21.5 -11.0 -11.0 -10.5 3 -9.5 -24 -15.0 1.0 -14.5 2.0 [3,] 3.0 -18.0 -3.0 -6.0 -5.0 15 -5.0 -18 -5.0 4.0 -8.0 6.5 [4,] 3.5 -11.5 -1.5 -3.5 -1.5 20 -0.5 -10 -4.5 10.5 -1.5 11.0 [5,] 6.0 -1.0 0.0 1.0 0.0 22 2.0 1 0.0 19.0 2.0 11.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.208448 -23.97184 -8.673247 -10.478879 -10.3746548 4.847874 -10.3746548 [2,] 4.791552 -12.02816 2.673247 -1.521121 0.3746548 25.152126 0.3746548 [,8] [,9] [,10] [,11] [,12] [1,] -26.360574 -11.270431 -1.673247 -15.7633903 0.6947093 [2,] -9.639426 1.270431 9.673247 -0.2366097 12.3052907 $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(-1, 0.5, 3, 3.5, 6, -23, -21.5, -18, -11.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6ifug1281023871.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] [,8] [1,] 423.0 379.0 315.0 269.0 206.0 134.0 82.0 NA [2,] 425.0 389.5 318.5 277.5 223.5 150.0 103.5 NA [3,] 427.0 402.5 338.0 290.5 246.0 184.0 124.5 NA [4,] 428.5 418.5 360.5 303.0 258.0 194.5 135.0 NA [5,] 430.0 435.0 386.0 321.0 278.0 219.0 151.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 425.4036 389.2729 318.8435 278.8693 230.2643 163.7033 110.1326 NA [2,] 428.5964 415.7271 357.1565 302.1307 261.7357 204.2967 138.8674 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(423, 425, 427, 428.5, 430, 379, 389.5, 402.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7ifug1281023871.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,] 259.2857 269.0 258.50 [2,] 269.5000 277.5 269.50 [3,] 287.7857 290.5 292.50 [4,] 298.3571 303.0 307.00 [5,] 316.8571 321.0 327.75 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 274.6238 278.8693 275.396 [2,] 300.9477 302.1307 309.604 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(259.285714285714, 269.5, 287.785714285714, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1xwva1281023871.ps tmp/1xwva1281023871.png",intern=TRUE)) character(0) > try(system("convert tmp/2xwva1281023871.ps tmp/2xwva1281023871.png",intern=TRUE)) character(0) > try(system("convert tmp/3q5cd1281023871.ps tmp/3q5cd1281023871.png",intern=TRUE)) character(0) > try(system("convert tmp/4q5cd1281023871.ps tmp/4q5cd1281023871.png",intern=TRUE)) character(0) > try(system("convert tmp/5q5cd1281023871.ps tmp/5q5cd1281023871.png",intern=TRUE)) character(0) > try(system("convert tmp/6ifug1281023871.ps tmp/6ifug1281023871.png",intern=TRUE)) character(0) > try(system("convert tmp/7ifug1281023871.ps tmp/7ifug1281023871.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.989 1.503 2.223