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(255 + ,280.2 + ,299.9 + ,339.2 + ,374.2 + ,393.5 + ,389.2 + ,381.7 + ,375.2 + ,369 + ,357.4 + ,352.1 + ,346.5 + ,342.9 + ,340.3 + ,328.3 + ,322.9 + ,314.3 + ,308.9 + ,294 + ,285.6 + ,281.2 + ,280.3 + ,278.8 + ,274.5 + ,270.4 + ,263.4 + ,259.9 + ,258 + ,262.7 + ,284.7 + ,311.3 + ,322.1 + ,327 + ,331.3 + ,333.3 + ,321.4 + ,327 + ,320 + ,314.7 + ,316.7 + ,314.4 + ,321.3 + ,318.2 + ,307.2 + ,301.3 + ,287.5 + ,277.7 + ,274.4 + ,258.8 + ,253.3 + ,251 + ,248.4 + ,249.5 + ,246.1 + ,244.5 + ,243.6 + ,244 + ,240.8 + ,249.8 + ,248 + ,259.4 + ,260.5 + ,260.8 + ,261.3 + ,259.5 + ,256.6 + ,257.9 + ,256.5 + ,254.2 + ,253.3 + ,253.8 + ,255.5 + ,257.1 + ,257.3 + ,253.2 + ,252.8 + ,252 + ,250.7 + ,252.2 + ,250 + ,251 + ,253.4 + ,251.2 + ,255.6 + ,261.1 + ,258.9 + ,259.9 + ,261.2 + ,264.7 + ,267.1 + ,266.4 + ,267.7 + ,268.6 + ,267.5 + ,268.5 + ,268.5 + ,270.5 + ,270.9 + ,270.1 + ,269.3 + ,269.8 + ,270.1 + ,264.9 + ,263.7 + ,264.8 + ,263.7 + ,255.9 + ,276.2 + ,360.1 + ,380.5 + ,373.7 + ,369.8 + ,366.6 + ,359.3 + ,345.8 + ,326.2 + ,324.5 + ,328.1 + ,327.5 + ,324.4 + ,316.5 + ,310.9 + ,301.5 + ,291.7 + ,290.4 + ,287.4 + ,277.7 + ,281.6 + ,288 + ,276 + ,272.9 + ,283 + ,283.3 + ,276.8 + ,284.5 + ,282.7 + ,281.2 + ,287.4 + ,283.1 + ,284 + ,285.5 + ,289.2 + ,292.5 + ,296.4 + ,305.2 + ,303.9 + ,311.5 + ,316.3 + ,316.7 + ,322.5 + ,317.1 + ,309.8 + ,303.8 + ,290.3 + ,293.7 + ,291.7 + ,296.5 + ,289.1 + ,288.5 + ,293.8 + ,297.7 + ,305.4 + ,302.7 + ,302.5 + ,303 + ,294.5 + ,294.1 + ,294.5 + ,297.1 + ,289.4 + ,292.4 + ,287.9 + ,286.6 + ,280.5 + ,272.4 + ,269.2 + ,270.6 + ,267.3 + ,262.5 + ,266.8 + ,268.8 + ,263.1 + ,261.2 + ,266 + ,262.5 + ,265.2 + ,261.3 + ,253.7 + ,249.2 + ,239.1 + ,236.4 + ,235.2 + ,245.2 + ,246.2 + ,247.7 + ,251.4 + ,253.3 + ,254.8 + ,250 + ,249.3 + ,241.5 + ,243.3 + ,248 + ,253 + ,252.9 + ,251.5 + ,251.6 + ,253.5 + ,259.8 + ,334.1 + ,448 + ,445.8 + ,445 + ,448.2 + ,438.2 + ,439.8 + ,423.4 + ,410.8 + ,408.4 + ,406.7 + ,405.9 + ,402.7 + ,405.1 + ,399.6 + ,386.5 + ,381.4 + ,375.2 + ,357.7 + ,359 + ,355 + ,352.7 + ,344.4 + ,343.8 + ,338 + ,339 + ,333.3 + ,334.4 + ,328.3 + ,330.7 + ,330 + ,331.6 + ,351.2 + ,389.4 + ,410.9 + ,442.8 + ,462.8 + ,466.9 + ,461.7 + ,439.2 + ,430.3 + ,416.1 + ,402.5 + ,397.3 + ,403.3 + ,395.9 + ,387.8 + ,378.6 + ,377.1 + ,370.4 + ,362 + ,350.3 + ,348.2 + ,344.6 + ,343.5 + ,342.8 + ,347.6 + ,346.6 + ,349.5 + ,342.1 + ,342 + ,342.8 + ,339.3 + ,348.2 + ,333.7 + ,334.7 + ,354 + ,367.7 + ,363.3 + ,358.4 + ,353.1 + ,343.1 + ,344.6 + ,344.4 + ,333.9 + ,331.7 + ,324.3 + ,321.2 + ,322.4 + ,321.7 + ,320.5 + ,312.8 + ,309.7 + ,315.6 + ,309.7 + ,304.6 + ,302.5 + ,301.5 + ,298.8 + ,291.3 + ,293.6 + ,294.6 + ,285.9 + ,297.6 + ,301.1 + ,293.8 + ,297.7 + ,292.9 + ,292.1 + ,287.2 + ,288.2 + ,283.8 + ,299.9 + ,292.4 + ,293.3 + ,300.8 + ,293.7 + ,293.1 + ,294.4 + ,292.1 + ,291.9 + ,282.5 + ,277.9 + ,287.5 + ,289.2 + ,285.6 + ,293.2 + ,290.8 + ,283.1 + ,275 + ,287.8 + ,287.8 + ,287.4 + ,284 + ,277.8 + ,277.6 + ,304.9 + ,294 + ,300.9 + ,324 + ,332.9 + ,341.6 + ,333.4 + ,348.2 + ,344.7 + ,344.7 + ,329.3 + ,323.5 + ,323.2 + ,317.4 + ,330.1 + ,329.2 + ,334.9 + ,315.8 + ,315.4 + ,319.6 + ,317.3 + ,313.8 + ,315.8 + ,311.3) > par1 = '10' > #'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] 360 > (np <- floor(n / par1)) [1] 36 > 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] 36 36 36 36 36 36 36 36 36 36 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 255.0 357.4 285.6 284.7 316.7 253.3 248.0 253.3 250.0 267.1 269.3 380.5 [2,] 280.2 352.1 281.2 311.3 314.4 251.0 259.4 253.8 251.0 266.4 269.8 373.7 [3,] 299.9 346.5 280.3 322.1 321.3 248.4 260.5 255.5 253.4 267.7 270.1 369.8 [4,] 339.2 342.9 278.8 327.0 318.2 249.5 260.8 257.1 251.2 268.6 264.9 366.6 [5,] 374.2 340.3 274.5 331.3 307.2 246.1 261.3 257.3 255.6 267.5 263.7 359.3 [6,] 393.5 328.3 270.4 333.3 301.3 244.5 259.5 253.2 261.1 268.5 264.8 345.8 [7,] 389.2 322.9 263.4 321.4 287.5 243.6 256.6 252.8 258.9 268.5 263.7 326.2 [8,] 381.7 314.3 259.9 327.0 277.7 244.0 257.9 252.0 259.9 270.5 255.9 324.5 [9,] 375.2 308.9 258.0 320.0 274.4 240.8 256.5 250.7 261.2 270.9 276.2 328.1 [10,] 369.0 294.0 262.7 314.7 258.8 249.8 254.2 252.2 264.7 270.1 360.1 327.5 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 324.4 276.0 284.0 322.5 293.8 289.4 266.8 239.1 249.3 334.1 406.7 355.0 [2,] 316.5 272.9 285.5 317.1 297.7 292.4 268.8 236.4 241.5 448.0 405.9 352.7 [3,] 310.9 283.0 289.2 309.8 305.4 287.9 263.1 235.2 243.3 445.8 402.7 344.4 [4,] 301.5 283.3 292.5 303.8 302.7 286.6 261.2 245.2 248.0 445.0 405.1 343.8 [5,] 291.7 276.8 296.4 290.3 302.5 280.5 266.0 246.2 253.0 448.2 399.6 338.0 [6,] 290.4 284.5 305.2 293.7 303.0 272.4 262.5 247.7 252.9 438.2 386.5 339.0 [7,] 287.4 282.7 303.9 291.7 294.5 269.2 265.2 251.4 251.5 439.8 381.4 333.3 [8,] 277.7 281.2 311.5 296.5 294.1 270.6 261.3 253.3 251.6 423.4 375.2 334.4 [9,] 281.6 287.4 316.3 289.1 294.5 267.3 253.7 254.8 253.5 410.8 357.7 328.3 [10,] 288.0 283.1 316.7 288.5 297.1 262.5 249.2 250.0 259.8 408.4 359.0 330.7 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 330.0 430.3 362.0 342.0 353.1 320.5 293.6 288.2 291.9 287.8 332.9 330.1 [2,] 331.6 416.1 350.3 342.8 343.1 312.8 294.6 283.8 282.5 287.8 341.6 329.2 [3,] 351.2 402.5 348.2 339.3 344.6 309.7 285.9 299.9 277.9 287.4 333.4 334.9 [4,] 389.4 397.3 344.6 348.2 344.4 315.6 297.6 292.4 287.5 284.0 348.2 315.8 [5,] 410.9 403.3 343.5 333.7 333.9 309.7 301.1 293.3 289.2 277.8 344.7 315.4 [6,] 442.8 395.9 342.8 334.7 331.7 304.6 293.8 300.8 285.6 277.6 344.7 319.6 [7,] 462.8 387.8 347.6 354.0 324.3 302.5 297.7 293.7 293.2 304.9 329.3 317.3 [8,] 466.9 378.6 346.6 367.7 321.2 301.5 292.9 293.1 290.8 294.0 323.5 313.8 [9,] 461.7 377.1 349.5 363.3 322.4 298.8 292.1 294.4 283.1 300.9 323.2 315.8 [10,] 439.2 370.4 342.1 358.4 321.7 291.3 287.2 292.1 275.0 324.0 317.4 311.3 [,37] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 25.2 -5.3 -4.4 26.6 -2.3 -2.3 11.4 0.5 1.0 -0.7 0.5 -6.8 -7.9 [2,] 19.7 -5.6 -0.9 10.8 6.9 -2.6 1.1 1.7 2.4 1.3 0.3 -3.9 -5.6 [3,] 39.3 -3.6 -1.5 4.9 -3.1 1.1 0.3 1.6 -2.2 0.9 -5.2 -3.2 -9.4 [4,] 35.0 -2.6 -4.3 4.3 -11.0 -3.4 0.5 0.2 4.4 -1.1 -1.2 -7.3 -9.8 [5,] 19.3 -12.0 -4.1 2.0 -5.9 -1.6 -1.8 -4.1 5.5 1.0 1.1 -13.5 -1.3 [6,] -4.3 -5.4 -7.0 -11.9 -13.8 -0.9 -2.9 -0.4 -2.2 0.0 -1.1 -19.6 -3.0 [7,] -7.5 -8.6 -3.5 5.6 -9.8 0.4 1.3 -0.8 1.0 2.0 -7.8 -1.7 -9.7 [8,] -6.5 -5.4 -1.9 -7.0 -3.3 -3.2 -1.4 -1.3 1.3 0.4 20.3 3.6 3.9 [9,] -6.2 -14.9 4.7 -5.3 -15.6 9.0 -2.3 1.5 3.5 -0.8 83.9 -0.6 6.4 [10,] -11.6 -8.4 22.0 2.0 -5.5 -1.8 -0.9 -2.2 2.4 -0.8 20.4 -3.1 -12.0 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] -3.1 1.5 -5.4 3.9 3.0 2.0 -2.7 -7.8 113.9 -0.8 -2.3 1.6 [2,] 10.1 3.7 -7.3 7.7 -4.5 -5.7 -1.2 1.8 -2.2 -3.2 -8.3 19.6 [3,] 0.3 3.3 -6.0 -2.7 -1.3 -1.9 10.0 4.7 -0.8 2.4 -0.6 38.2 [4,] -6.5 3.9 -13.5 -0.2 -6.1 4.8 1.0 5.0 3.2 -5.5 -5.8 21.5 [5,] 7.7 8.8 3.4 0.5 -8.1 -3.5 1.5 -0.1 -10.0 -13.1 1.0 31.9 [6,] -1.8 -1.3 -2.0 -8.5 -3.2 2.7 3.7 -1.4 1.6 -5.1 -5.7 20.0 [7,] -1.5 7.6 4.8 -0.4 1.4 -3.9 1.9 0.1 -16.4 -6.2 1.1 4.1 [8,] 6.2 4.8 -7.4 0.4 -3.3 -7.6 1.5 1.9 -12.6 -17.5 -6.1 -5.2 [9,] -4.3 0.4 -0.6 2.6 -4.8 -4.5 -4.8 6.3 -2.4 1.3 2.4 -22.5 [10,] 0.9 5.8 5.3 -7.7 4.3 -10.1 -0.7 74.3 -1.7 -4.0 -0.7 -8.9 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] -14.2 -11.7 0.8 -10.0 -7.7 1.0 -4.4 -9.4 0.0 8.7 -0.9 NA [2,] -13.6 -2.1 -3.5 1.5 -3.1 -8.7 16.1 -4.6 -0.4 -8.2 5.7 NA [3,] -5.2 -3.6 8.9 -0.2 5.9 11.7 -7.5 9.6 -3.4 14.8 -19.1 NA [4,] 6.0 -1.1 -14.5 -10.5 -5.9 3.5 0.9 1.7 -6.2 -3.5 -0.4 NA [5,] -7.4 -0.7 1.0 -2.2 -5.1 -7.3 7.5 -3.6 -0.2 0.0 4.2 NA [6,] -8.1 4.8 19.3 -7.4 -2.1 3.9 -7.1 7.6 27.3 -15.4 -2.3 NA [7,] -9.2 -1.0 13.7 -3.1 -1.0 -4.8 -0.6 -2.4 -10.9 -5.8 -3.5 NA [8,] -1.5 2.9 -4.4 1.2 -2.7 -0.8 1.3 -7.7 6.9 -0.3 2.0 NA [9,] -6.7 -7.4 -4.9 -0.7 -7.5 -4.9 -2.3 -8.1 23.1 -5.8 -4.5 NA [10,] -8.4 -0.1 -5.3 -1.2 2.3 1.0 -0.2 12.8 8.9 12.7 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/1esei1289922714.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/2esei1289922714.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/3p1d31289922714.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/4p1d31289922714.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] [1,] 239.10 236.40 235.20 245.2 246.10 244.50 243.60 244.00 240.80 249.20 [2,] 268.20 271.35 274.00 273.7 271.00 269.45 266.85 265.90 269.10 263.70 [3,] 292.75 296.15 302.65 302.1 298.75 301.05 296.10 294.05 294.45 293.05 [4,] 333.50 342.20 341.85 344.1 339.15 336.85 327.75 325.75 325.65 329.10 [5,] 430.30 448.00 402.70 445.0 410.90 395.90 389.20 381.70 377.10 408.40 $n [1] 36 36 36 36 36 36 36 36 36 36 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 275.5543 277.4928 284.7828 283.5613 280.8038 283.3013 280.063 278.2895 [2,] 309.9457 314.8072 320.5172 320.6387 316.6962 318.7987 312.137 309.8105 [,9] [,10] [1,] 279.5585 275.828 [2,] 309.3415 310.272 $out [1] 445.8 448.2 438.2 442.8 439.8 462.8 423.4 466.9 410.8 461.7 439.2 $group [1] 3 5 6 6 7 7 8 8 9 9 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5p1d31289922714.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] [1,] -14.20 -13.60 -9.4 -14.50 -13.50 -15.40 -16.40 -12.60 -15.60 -12.00 [2,] -5.35 -4.55 -3.3 -6.00 -4.60 -6.35 -6.00 -5.30 -5.55 -4.65 [3,] -0.85 -1.05 -0.4 -1.10 -0.45 -2.15 -1.25 -1.35 -2.35 -0.70 [4,] 1.55 3.05 4.8 3.35 1.75 0.80 1.20 1.70 1.95 3.35 [5,] 11.40 10.80 14.8 6.00 8.80 7.60 7.60 6.90 9.00 12.80 $n [1] 36 36 36 36 36 36 36 36 36 35 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -2.667 -3.0513333 -2.533 -3.562167 -2.122167 -4.0328333 -3.146 -3.1933333 [2,] 0.967 0.9513333 1.733 1.362167 1.222167 -0.2671667 0.646 0.4933333 [,9] [,10] [1,] -4.325 -2.83655 [2,] -0.375 1.43655 $out [1] 25.2 26.6 113.9 19.7 19.6 16.1 39.3 38.2 -19.1 35.0 21.5 19.3 [13] 31.9 -19.6 20.0 19.3 27.3 13.7 20.3 -17.5 83.9 -22.5 23.1 22.0 [25] 20.4 74.3 $group [1] 1 1 1 2 2 2 3 3 3 4 4 5 5 6 6 6 6 7 8 8 9 9 9 10 10 [26] 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/60bdo1289922714.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] [,9] [,10] [,11] [1,] 255.0 294.0 258.00 311.30 258.80 240.80 254.20 250.70 250.00 266.4 255.9 [2,] 299.9 314.3 262.70 314.70 277.70 244.00 256.50 252.20 251.20 267.5 263.7 [3,] 371.6 334.3 272.45 321.75 304.25 247.25 258.65 253.25 257.25 268.5 267.1 [4,] 381.7 346.5 280.30 327.00 316.70 249.80 260.50 255.50 261.10 270.1 270.1 [5,] 393.5 357.4 285.60 333.30 321.30 253.30 261.30 257.30 264.70 270.9 276.2 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [1,] 324.50 277.70 272.90 284.00 288.5 293.8 262.50 261.2 235.20 241.50 408.4 [2,] 327.50 287.40 276.80 289.20 290.3 294.5 269.20 261.2 239.10 248.00 410.8 [3,] 352.55 291.05 282.85 300.15 295.1 297.4 276.45 262.8 246.95 251.55 439.0 [4,] 369.80 310.90 283.30 311.50 309.8 302.7 287.90 266.0 251.40 253.00 445.8 [5,] 380.50 324.40 287.40 316.70 322.5 305.4 292.40 268.8 254.80 259.80 448.2 [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [1,] 357.70 328.3 330.00 370.4 342.1 333.7 321.2 291.30 285.9 292.1 275.00 [2,] 375.20 333.3 351.20 378.6 343.5 339.3 322.4 301.50 292.1 292.1 282.50 [3,] 393.05 338.5 425.05 396.6 347.1 345.5 332.8 307.15 293.7 293.2 286.55 [4,] 405.10 344.4 461.70 403.3 349.5 358.4 344.4 312.80 297.6 294.4 290.80 [5,] 406.70 355.0 466.90 430.3 350.3 367.7 353.1 320.50 301.1 294.4 293.20 [,34] [,35] [,36] [,37] [1,] 277.6 317.40 311.30 NA [2,] 284.0 323.50 315.40 NA [3,] 287.8 333.15 316.55 NA [4,] 300.9 344.70 329.20 NA [5,] 324.0 348.20 334.90 NA $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 [26] 10 10 10 10 10 10 10 10 10 10 10 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 330.7295 318.2116 263.6563 315.6044 284.7640 244.3521 256.6514 251.6012 [2,] 412.4705 350.3884 281.2437 327.8956 323.7360 250.1479 260.6486 254.8988 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 252.3036 267.2009 263.9023 331.4152 279.3085 279.6023 289.0080 285.357 [2,] 262.1964 269.7991 270.2977 373.6848 302.7915 286.0977 311.2920 304.843 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 293.3030 267.1067 260.4017 240.8044 249.0518 421.5126 378.1108 332.954 [2,] 301.4970 285.7933 265.1983 253.0956 254.0482 456.4874 407.9892 344.046 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 369.8398 384.2589 344.1022 335.9569 321.8079 301.5041 290.952 292.0508 [2,] 480.2602 408.9411 350.0978 355.0431 343.7921 312.7959 296.448 294.3492 [,33] [,34] [,35] [,36] [,37] [1,] 282.403 279.3561 322.5576 309.655 NA [2,] 290.697 296.2439 343.7424 323.445 NA $out [1] 284.7 248.0 360.1 253.7 249.2 334.1 362.0 288.2 283.8 299.9 300.8 $group [1] 4 7 11 19 19 22 27 32 32 32 32 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" NA Warning message: In bxp(list(stats = c(255, 299.9, 371.6, 381.7, 393.5, 294, 314.3, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/70bdo1289922714.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,] 305.5056 292.750 296.2500 [2,] 306.2333 294.050 297.2125 [3,] 308.8569 296.125 301.9125 [4,] 310.4111 301.050 307.0125 [5,] 311.3472 302.650 310.1000 $n [1] 10 10 10 $conf [,1] [,2] [,3] [1,] 306.7696 292.6275 297.016 [2,] 310.9443 299.6225 306.809 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(305.505555555556, 306.233333333333, 308.856944444444, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1esei1289922714.ps tmp/1esei1289922714.png",intern=TRUE)) character(0) > try(system("convert tmp/2esei1289922714.ps tmp/2esei1289922714.png",intern=TRUE)) character(0) > try(system("convert tmp/3p1d31289922714.ps tmp/3p1d31289922714.png",intern=TRUE)) character(0) > try(system("convert tmp/4p1d31289922714.ps tmp/4p1d31289922714.png",intern=TRUE)) character(0) > try(system("convert tmp/5p1d31289922714.ps tmp/5p1d31289922714.png",intern=TRUE)) character(0) > try(system("convert tmp/60bdo1289922714.ps tmp/60bdo1289922714.png",intern=TRUE)) character(0) > try(system("convert tmp/70bdo1289922714.ps tmp/70bdo1289922714.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.043 1.493 2.583