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(310.95,312.97,315.07,315.43,315.73,315.77,315.77,315.77,313.99,314.57,314.63,314.65,314.65,314.93,315.27,316.26,316.98,317.01,317.07,317.07,317,317.08,317.04,317,317.05,321.59,325.59,326.23,326.28,326.35,326.35,326.35,326.39,326.74,326.9,326.9,326.91,336.93,348.5,349.43,349.26,349.26,349.28,349.61,349.66,349.68,349.91,349.91,350.89,355.52,356.36,357.04,360.28,360.63,360.79,360.97,361,361.01,361,361,361.58,363.19,363.61,364.14,365.51,365.51,365.5,365.5,364.59,364.63,364.54,363.67,365.22,369.05,370.45,370.46,370.46,370.58,370.58,370.22,370.21,370.29,370.29,370.2,370.2,372.55,374.51,375.58,375.75,375.75,375.75,375.69,375.76,377.5,377.51,377.74) > 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] 96 > (np <- floor(n / par1)) [1] 8 > 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] 8 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 310.95 314.65 317.05 326.91 350.89 361.58 365.22 370.20 NA [2,] 312.97 314.93 321.59 336.93 355.52 363.19 369.05 372.55 NA [3,] 315.07 315.27 325.59 348.50 356.36 363.61 370.45 374.51 NA [4,] 315.43 316.26 326.23 349.43 357.04 364.14 370.46 375.58 NA [5,] 315.73 316.98 326.28 349.26 360.28 365.51 370.46 375.75 NA [6,] 315.77 317.01 326.35 349.26 360.63 365.51 370.58 375.75 NA [7,] 315.77 317.07 326.35 349.28 360.79 365.50 370.58 375.75 NA [8,] 315.77 317.07 326.35 349.61 360.97 365.50 370.22 375.69 NA [9,] 313.99 317.00 326.39 349.66 361.00 364.59 370.21 375.76 NA [10,] 314.57 317.08 326.74 349.68 361.01 364.63 370.29 377.50 NA [11,] 314.63 317.04 326.90 349.91 361.00 364.54 370.29 377.51 NA [12,] 314.65 317.00 326.90 349.91 361.00 363.67 370.20 377.74 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 2.02 0.28 4.54 10.02 4.63 1.61 3.83 2.35 NA [2,] 2.10 0.34 4.00 11.57 0.84 0.42 1.40 1.96 NA [3,] 0.36 0.99 0.64 0.93 0.68 0.53 0.01 1.07 NA [4,] 0.30 0.72 0.05 -0.17 3.24 1.37 0.00 0.17 NA [5,] 0.04 0.03 0.07 0.00 0.35 0.00 0.12 0.00 NA [6,] 0.00 0.06 0.00 0.02 0.16 -0.01 0.00 0.00 NA [7,] 0.00 0.00 0.00 0.33 0.18 0.00 -0.36 -0.06 NA [8,] -1.78 -0.07 0.04 0.05 0.03 -0.91 -0.01 0.07 NA [9,] 0.58 0.08 0.35 0.02 0.01 0.04 0.08 1.74 NA [10,] 0.06 -0.04 0.16 0.23 -0.01 -0.09 0.00 0.01 NA [11,] 0.02 -0.04 0.00 0.00 0.00 -0.87 -0.09 0.23 NA [12,] 0.00 0.05 0.01 0.98 0.58 1.55 0.00 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/11met1243248063.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/2884h1243248063.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/3mrs91243248063.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/4atjx1243248063.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] [1,] 310.95 312.970 315.07 315.430 315.730 315.770 315.770 315.77 313.990 [2,] 315.85 318.260 320.43 321.245 321.630 321.680 321.710 321.71 321.695 [3,] 338.90 346.225 352.43 353.235 354.770 354.945 355.035 355.29 355.330 [4,] 363.40 366.120 367.03 367.300 367.985 368.045 368.040 367.86 367.400 [5,] 370.20 372.550 374.51 375.580 375.750 375.750 375.750 375.69 375.760 [,10] [,11] [,12] [1,] 314.570 314.630 314.650 [2,] 321.910 321.970 321.950 [3,] 355.345 355.455 355.455 [4,] 367.460 367.415 366.935 [5,] 377.500 377.510 377.740 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 312.3379 319.4897 326.3986 327.508 328.8754 329.0448 329.1544 329.5099 [2,] 365.4621 372.9603 378.4614 378.962 380.6646 380.8452 380.9156 381.0701 [,9] [,10] [,11] [,12] [1,] 329.7985 329.9001 330.0688 330.3257 [2,] 380.8615 380.7899 380.8412 380.5843 $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(310.95, 315.85, 338.9, 363.4, 370.2, 312.97, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5wbo21243248063.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,] 0.280 0.34 0.010 -0.170 0.000 -0.01 -0.06 -0.910 0.010 -0.090 -0.090 0.000 [2,] 1.815 0.63 0.445 0.025 0.000 0.00 -0.03 -0.490 0.030 -0.025 -0.065 0.005 [3,] 3.090 1.68 0.660 0.235 0.035 0.00 0.00 0.010 0.080 0.005 0.000 0.050 [4,] 4.585 3.05 0.960 1.045 0.095 0.04 0.09 0.045 0.465 0.110 0.010 0.780 [5,] 4.630 4.00 1.070 1.370 0.120 0.06 0.18 0.070 0.580 0.230 0.020 1.550 $n [1] 8 8 8 8 8 8 8 8 8 8 8 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1.542638 0.3281533 0.3723136 -0.3347866 -0.01806836 -0.02234457 [2,] 4.637362 3.0318467 0.9476864 0.8047866 0.08806836 0.02234457 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.06703372 -0.2888587 -0.1629972 -0.07041294 -0.04189608 -0.4128175 [2,] 0.06703372 0.3088587 0.3229972 0.08041294 0.04189608 0.5128175 $out [1] 10.02 11.57 3.24 0.35 0.16 0.33 -0.36 -1.78 1.74 -0.87 0.23 $group [1] 1 2 4 5 6 7 7 8 9 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(0.280000000000030, 1.81500000000003, 3.09000000000000, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6esn91243248063.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] [1,] 312.97 314.650 325.590 348.500 350.89 361.580 370.200 374.510 NA [2,] 314.28 315.765 325.910 348.880 356.70 363.640 370.205 375.045 NA [3,] 314.86 317.000 326.350 349.355 360.71 364.565 370.290 375.750 NA [4,] 315.75 317.055 326.565 349.670 361.00 365.500 370.460 376.630 NA [5,] 315.77 317.080 326.900 349.910 361.01 365.510 370.580 377.740 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 314.1895 316.4116 326.0513 348.9947 358.7487 363.7166 370.1737 375.0271 [2,] 315.5305 317.5884 326.6487 349.7153 362.6713 365.4134 370.4063 376.4729 [,9] [1,] NA [2,] NA $out [1] 310.95 317.05 321.59 326.91 336.93 365.22 369.05 370.20 372.55 $group [1] 1 3 3 4 4 7 7 8 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(312.97, 314.28, 314.86, 315.75, 315.77, 314.65, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7mnqn1243248063.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,] 346.1700 352.4300 344.1650 [2,] 346.4956 352.8325 344.4469 [3,] 347.5694 354.9900 345.1019 [4,] 347.6419 355.3375 345.3531 [5,] 347.7275 355.4550 345.4000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 347.0466 353.8475 344.6885 [2,] 348.0922 356.1325 345.5152 $out [1] 339.6812 343.3413 338.9000 346.2250 339.4700 342.2900 $group [1] 1 1 2 2 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(346.17, 346.495625, 347.569375, 347.641875, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/11met1243248063.ps tmp/11met1243248063.png") > system("convert tmp/2884h1243248063.ps tmp/2884h1243248063.png") > system("convert tmp/3mrs91243248063.ps tmp/3mrs91243248063.png") > system("convert tmp/4atjx1243248063.ps tmp/4atjx1243248063.png") > system("convert tmp/5wbo21243248063.ps tmp/5wbo21243248063.png") > system("convert tmp/6esn91243248063.ps tmp/6esn91243248063.png") > system("convert tmp/7mnqn1243248063.ps tmp/7mnqn1243248063.png") > > > proc.time() user system elapsed 1.342 1.062 4.042