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 = '6' > #'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] 60 > 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] 60 60 60 60 60 60 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 255.0 389.2 346.5 308.9 274.5 284.7 321.4 321.3 274.4 246.1 248.0 256.6 [2,] 280.2 381.7 342.9 294.0 270.4 311.3 327.0 318.2 258.8 244.5 259.4 257.9 [3,] 299.9 375.2 340.3 285.6 263.4 322.1 320.0 307.2 253.3 243.6 260.5 256.5 [4,] 339.2 369.0 328.3 281.2 259.9 327.0 314.7 301.3 251.0 244.0 260.8 254.2 [5,] 374.2 357.4 322.9 280.3 258.0 331.3 316.7 287.5 248.4 240.8 261.3 253.3 [6,] 393.5 352.1 314.3 278.8 262.7 333.3 314.4 277.7 249.5 249.8 259.5 253.8 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 255.5 250.7 255.6 267.1 268.5 270.1 276.2 359.3 324.4 287.4 283.0 287.4 [2,] 257.1 252.2 261.1 266.4 270.5 264.9 360.1 345.8 316.5 277.7 283.3 283.1 [3,] 257.3 250.0 258.9 267.7 270.9 263.7 380.5 326.2 310.9 281.6 276.8 284.0 [4,] 253.2 251.0 259.9 268.6 270.1 264.8 373.7 324.5 301.5 288.0 284.5 285.5 [5,] 252.8 253.4 261.2 267.5 269.3 263.7 369.8 328.1 291.7 276.0 282.7 289.2 [6,] 252.0 251.2 264.7 268.5 269.8 255.9 366.6 327.5 290.4 272.9 281.2 292.5 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 296.4 322.5 291.7 305.4 294.5 280.5 266.8 265.2 235.2 254.8 253.0 334.1 [2,] 305.2 317.1 296.5 302.7 297.1 272.4 268.8 261.3 245.2 250.0 252.9 448.0 [3,] 303.9 309.8 289.1 302.5 289.4 269.2 263.1 253.7 246.2 249.3 251.5 445.8 [4,] 311.5 303.8 288.5 303.0 292.4 270.6 261.2 249.2 247.7 241.5 251.6 445.0 [5,] 316.3 290.3 293.8 294.5 287.9 267.3 266.0 239.1 251.4 243.3 253.5 448.2 [6,] 316.7 293.7 297.7 294.1 286.6 262.5 262.5 236.4 253.3 248.0 259.8 438.2 [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 439.8 402.7 357.7 338.0 330.0 462.8 402.5 377.1 343.5 342.0 354.0 344.6 [2,] 423.4 405.1 359.0 339.0 331.6 466.9 397.3 370.4 342.8 342.8 367.7 344.4 [3,] 410.8 399.6 355.0 333.3 351.2 461.7 403.3 362.0 347.6 339.3 363.3 333.9 [4,] 408.4 386.5 352.7 334.4 389.4 439.2 395.9 350.3 346.6 348.2 358.4 331.7 [5,] 406.7 381.4 344.4 328.3 410.9 430.3 387.8 348.2 349.5 333.7 353.1 324.3 [6,] 405.9 375.2 343.8 330.7 442.8 416.1 378.6 344.6 342.1 334.7 343.1 321.2 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [1,] 322.4 309.7 293.6 297.7 299.9 294.4 289.2 287.8 304.9 333.4 323.2 315.4 [2,] 321.7 304.6 294.6 292.9 292.4 292.1 285.6 287.8 294.0 348.2 317.4 319.6 [3,] 320.5 302.5 285.9 292.1 293.3 291.9 293.2 287.4 300.9 344.7 330.1 317.3 [4,] 312.8 301.5 297.6 287.2 300.8 282.5 290.8 284.0 324.0 344.7 329.2 313.8 [5,] 309.7 298.8 301.1 288.2 293.7 277.9 283.1 277.8 332.9 329.3 334.9 315.8 [6,] 315.6 291.3 293.8 283.8 293.1 287.5 275.0 277.6 341.6 323.5 315.8 311.3 [,61] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 25.2 -7.5 -3.6 -14.9 -4.1 26.6 5.6 -3.1 -15.6 -1.6 11.4 1.3 1.6 [2,] 19.7 -6.5 -2.6 -8.4 -7.0 10.8 -7.0 -11.0 -5.5 -0.9 1.1 -1.4 0.2 [3,] 39.3 -6.2 -12.0 -4.4 -3.5 4.9 -5.3 -5.9 -2.3 0.4 0.3 -2.3 -4.1 [4,] 35.0 -11.6 -5.4 -0.9 -1.9 4.3 2.0 -13.8 -2.6 -3.2 0.5 -0.9 -0.4 [5,] 19.3 -5.3 -8.6 -1.5 4.7 2.0 -2.3 -9.8 1.1 9.0 -1.8 0.5 -0.8 [6,] -4.3 -5.6 -5.4 -4.3 22.0 -11.9 6.9 -3.3 -3.4 -1.8 -2.9 1.7 -1.3 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 1.5 5.5 -0.7 2.0 -5.2 83.9 -13.5 -7.9 -9.7 0.3 -4.3 8.8 [2,] -2.2 -2.2 1.3 0.4 -1.2 20.4 -19.6 -5.6 3.9 -6.5 0.9 -1.3 [3,] 1.0 1.0 0.9 -0.8 1.1 -6.8 -1.7 -9.4 6.4 7.7 1.5 7.6 [4,] 2.4 1.3 -1.1 -0.8 -1.1 -3.9 3.6 -9.8 -12.0 -1.8 3.7 4.8 [5,] -2.2 3.5 1.0 0.5 -7.8 -3.2 -0.6 -1.3 -3.1 -1.5 3.3 0.4 [6,] 4.4 2.4 0.0 0.3 20.3 -7.3 -3.1 -3.0 10.1 6.2 3.9 5.8 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] -5.4 4.8 -2.7 2.6 -8.1 2.0 -3.9 10.0 -4.8 -0.1 113.9 -16.4 [2,] -7.3 -7.4 -0.2 -7.7 -3.2 -5.7 -7.6 1.0 -0.7 -1.4 -2.2 -12.6 [3,] -6.0 -0.6 0.5 3.0 1.4 -1.9 -4.5 1.5 -7.8 0.1 -0.8 -2.4 [4,] -13.5 5.3 -8.5 -4.5 -3.3 4.8 -10.1 3.7 1.8 1.9 3.2 -1.7 [5,] 3.4 3.9 -0.4 -1.3 -4.8 -3.5 -2.7 1.9 4.7 6.3 -10.0 -0.8 [6,] -2.0 7.7 0.4 -6.1 4.3 2.7 -1.2 1.5 5.0 74.3 1.6 -3.2 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 2.4 1.3 1.0 1.6 4.1 -5.2 -6.7 -0.7 0.8 13.7 -0.2 -0.7 [2,] -5.5 -4.0 -5.7 19.6 -5.2 6.0 -8.4 4.8 -3.5 -4.4 -10.5 -1.2 [3,] -13.1 -2.3 1.1 38.2 -22.5 -7.4 -11.7 -1.0 8.9 -4.9 -2.2 -7.7 [4,] -5.1 -8.3 -6.1 21.5 -8.9 -8.1 -2.1 2.9 -14.5 -5.3 -7.4 -3.1 [5,] -6.2 -0.6 2.4 31.9 -14.2 -9.2 -3.6 -7.4 1.0 -10.0 -3.1 5.9 [6,] -17.5 -5.8 -0.7 20.0 -13.6 -1.5 -1.1 -0.1 19.3 1.5 1.2 -5.9 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] -5.1 1.0 -4.8 -7.5 -2.3 -3.6 0.0 -10.9 14.8 -5.8 4.2 NA [2,] -2.1 -8.7 -0.8 0.9 -0.2 7.6 -0.4 6.9 -3.5 12.7 -2.3 NA [3,] -1.0 11.7 -4.9 7.5 -9.4 -2.4 -3.4 23.1 0.0 -0.9 -3.5 NA [4,] -2.7 3.5 1.0 -7.1 -4.6 -7.7 -6.2 8.9 -15.4 5.7 2.0 NA [5,] -7.5 -7.3 -4.4 -0.6 9.6 -8.1 -0.2 8.7 -5.8 -19.1 -4.5 NA [6,] 2.3 3.9 16.1 1.3 1.7 12.8 27.3 -8.2 -0.3 -0.4 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/1c2w91289653705.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/2c2w91289653705.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/3nbwt1289653705.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/4nbwt1289653705.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] [1,] 235.20 244.50 243.60 241.50 239.10 236.4 [2,] 272.25 270.45 268.45 269.35 267.40 266.6 [3,] 297.05 296.80 300.40 301.40 293.75 293.4 [4,] 333.75 342.80 336.60 336.80 333.30 334.0 [5,] 402.70 448.00 410.80 408.40 430.30 416.1 $n [1] 60 60 60 60 60 60 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 284.5054 282.0423 286.4990 287.6417 280.3079 279.6519 [2,] 309.5946 311.5577 314.3010 315.1583 307.1921 307.1481 $out [1] 439.8 462.8 466.9 445.8 461.7 445.0 439.2 448.2 438.2 442.8 $group [1] 1 1 2 3 3 4 4 5 6 6 $names [1] "1" "2" "3" "4" "5" "6" > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5nbwt1289653705.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] [1,] -16.40 -12.60 -13.10 -15.40 -14.20 -13.60 [2,] -5.15 -6.10 -4.90 -6.65 -5.05 -3.15 [3,] -0.45 -2.20 -1.35 -1.85 -1.30 0.30 [4,] 2.50 0.65 1.10 2.20 1.95 4.35 [5,] 13.70 7.60 8.90 8.90 9.60 12.80 $n [1] 60 60 60 60 60 59 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -2.010425 -3.5768456 -2.5738627 -3.65519754 -2.7278399 -1.242739 [2,] 1.110425 -0.8231544 -0.1261373 -0.04480246 0.1278399 1.842739 $out [1] 25.2 26.6 83.9 113.9 14.8 19.7 10.8 20.4 -19.6 19.6 12.7 39.3 [13] 38.2 -22.5 11.7 23.1 35.0 21.5 19.3 31.9 -19.1 22.0 20.3 74.3 [25] -17.5 20.0 19.3 16.1 27.3 $group [1] 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 4 4 5 5 5 6 6 6 6 6 6 6 6 $names [1] "1" "2" "3" "4" "5" "6" > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6xkve1289653705.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.00 352.1 314.3 278.8 258.00 284.70 314.40 277.70 248.40 240.80 259.4 [2,] 280.20 357.4 322.9 280.3 259.90 311.30 314.70 287.50 249.50 243.60 259.4 [3,] 319.55 372.1 334.3 283.4 263.05 324.55 318.35 304.25 252.15 244.25 260.0 [4,] 374.20 381.7 342.9 294.0 270.40 331.30 321.40 318.20 258.80 246.10 260.8 [5,] 393.50 389.2 346.5 308.9 274.50 333.30 327.00 321.30 258.80 249.80 261.3 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [1,] 253.30 252.00 250.0 255.6 266.4 268.50 263.70 360.1 324.5 290.4 272.90 [2,] 253.80 252.80 250.7 258.9 267.1 269.30 263.70 360.1 326.2 291.7 276.00 [3,] 255.35 254.35 251.1 260.5 267.6 269.95 264.25 368.2 327.8 306.2 279.65 [4,] 256.60 257.10 252.2 261.2 268.5 270.50 264.90 373.7 345.8 316.5 287.40 [5,] 257.90 257.30 253.4 261.2 268.6 270.90 264.90 380.5 359.3 324.4 288.00 [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [1,] 281.20 283.10 296.40 290.3 288.50 294.1 286.6 262.5 261.20 236.40 245.20 [2,] 281.20 284.00 303.90 293.7 289.10 294.5 287.9 267.3 262.50 239.10 245.20 [3,] 282.85 286.45 308.35 306.8 292.75 302.6 290.9 269.9 264.55 251.45 246.95 [4,] 283.30 289.20 316.30 317.1 296.50 303.0 294.5 272.4 266.80 261.30 251.40 [5,] 284.50 292.50 316.70 322.5 297.70 305.4 297.1 272.4 268.80 265.20 253.30 [,34] [,35] [,36] [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [1,] 241.50 251.50 438.2 405.9 375.20 343.80 328.30 330.0 416.10 378.6 344.60 [2,] 243.30 251.60 438.2 406.7 381.40 344.40 330.70 331.6 430.30 387.8 348.20 [3,] 248.65 252.95 445.4 409.6 393.05 353.85 333.85 370.3 450.45 396.6 356.15 [4,] 250.00 253.50 448.0 423.4 402.70 357.70 338.00 410.9 462.80 402.5 370.40 [5,] 254.80 253.50 448.2 439.8 405.10 359.00 339.00 442.8 466.90 403.3 377.10 [,45] [,46] [,47] [,48] [,49] [,50] [,51] [,52] [,53] [,54] [,55] [1,] 342.10 333.70 343.1 321.2 309.70 291.3 293.6 283.80 292.4 277.9 275.0 [2,] 342.80 334.70 353.1 324.3 312.80 298.8 293.6 287.20 293.1 282.5 283.1 [3,] 345.05 340.65 356.2 332.8 318.05 302.0 294.2 290.15 293.5 289.7 287.4 [4,] 347.60 342.80 363.3 344.4 321.70 304.6 297.6 292.90 299.9 292.1 290.8 [5,] 349.50 348.20 367.7 344.6 322.40 309.7 301.1 297.70 300.8 294.4 293.2 [,56] [,57] [,58] [,59] [,60] [,61] [1,] 277.6 294.00 323.50 315.8 311.3 NA [2,] 277.8 300.90 329.30 317.4 313.8 NA [3,] 285.7 314.45 339.05 326.2 315.6 NA [4,] 287.8 332.90 344.70 330.1 317.3 NA [5,] 287.8 341.60 348.20 334.9 319.6 NA $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 [39] 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 258.9170 356.4257 321.3994 274.5631 256.2772 311.6494 314.0283 284.4475 [2,] 380.1830 387.7743 347.2006 292.2369 269.8228 337.4506 322.6717 324.0525 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 246.1512 242.6374 259.0970 253.5439 251.5764 250.1325 259.0164 266.6970 [2,] 258.1488 245.8626 260.9030 257.1561 257.1236 252.0675 261.9836 268.5030 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 269.1760 263.4760 359.4276 315.1574 290.2032 272.2966 281.4954 283.0958 [2,] 270.7240 265.0240 376.9724 340.4426 322.1968 287.0034 284.2046 289.8042 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 300.3516 291.7062 287.9768 297.1172 286.6428 266.6103 261.7764 237.1303 [2,] 316.3484 321.8938 297.5232 308.0828 295.1572 273.1897 267.3236 265.7697 [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 242.9508 244.3283 251.7244 439.0787 398.828 379.3108 345.2711 329.1413 [2,] 250.9492 252.9717 254.1756 451.7213 420.372 406.7892 362.4289 338.5587 [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 319.1489 429.4865 387.118 341.8303 341.9538 335.4252 349.6207 319.8349 [2,] 421.4511 471.4135 406.082 370.4697 348.1462 345.8748 362.7793 345.7651 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [1,] 312.3092 298.2588 291.6199 286.4733 289.1138 283.5077 282.4333 279.2497 [2,] 323.7908 305.7412 296.7801 293.8267 297.8862 295.8923 292.3667 292.1503 [,57] [,58] [,59] [,60] [,61] [1,] 293.8090 329.1165 318.0081 313.3424 NA [2,] 335.0910 348.9835 334.3919 317.8576 NA $out [1] 274.4 248.0 264.7 270.1 255.9 276.2 276.8 280.5 235.2 259.8 334.1 285.9 $group [1] 9 11 15 18 18 19 23 30 33 35 36 51 $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" "37" "38" "39" "40" "41" "42" "43" "44" "45" [46] "46" "47" "48" "49" "50" "51" "52" "53" "54" "55" "56" "57" "58" "59" "60" [61] NA Warning message: In bxp(list(stats = c(255, 280.2, 319.55, 374.2, 393.5, 352.1, 357.4, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7xkve1289653705.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,] 306.1800 293.400 300.2750 [2,] 307.1817 293.750 300.6000 [3,] 308.4225 296.925 302.3500 [4,] 309.1733 300.400 303.4500 [5,] 310.7250 301.400 306.6375 $n [1] 6 6 6 $conf [,1] [,2] [,3] [1,] 307.1378 292.6355 300.5117 [2,] 309.7072 301.2145 304.1883 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(306.18, 307.181666666667, 308.4225, 309.173333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1c2w91289653705.ps tmp/1c2w91289653705.png",intern=TRUE)) character(0) > try(system("convert tmp/2c2w91289653705.ps tmp/2c2w91289653705.png",intern=TRUE)) character(0) > try(system("convert tmp/3nbwt1289653705.ps tmp/3nbwt1289653705.png",intern=TRUE)) character(0) > try(system("convert tmp/4nbwt1289653705.ps tmp/4nbwt1289653705.png",intern=TRUE)) character(0) > try(system("convert tmp/5nbwt1289653705.ps tmp/5nbwt1289653705.png",intern=TRUE)) character(0) > try(system("convert tmp/6xkve1289653705.ps tmp/6xkve1289653705.png",intern=TRUE)) character(0) > try(system("convert tmp/7xkve1289653705.ps tmp/7xkve1289653705.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.036 1.482 2.269