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(698,697,696,694,714,713,698,688,689,689,690,692,688,679,677,673,694,690,673,659,657,654,644,643,638,626,621,615,640,633,620,610,601,595,585,584,580,574,560,550,580,569,551,536,535,526,517,512,510,501,496,491,524,514,495,479,479,467,451,459,461,460,452,449,483,470,442,419,419,406,393,396,390,389,373,371,407,391,357,327,321,317,300,304,296,296,283,279,319,295,255,227,228,233,210,219,212,209,201,198,245,216,173,144,143,152,127,141,129,127,113,117,174,143,103,81,92,104,81,89) > 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] 120 > (np <- floor(n / par1)) [1] 10 > 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] 10 10 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 698 688 638 580 510 461 390 296 212 129 NA [2,] 697 679 626 574 501 460 389 296 209 127 NA [3,] 696 677 621 560 496 452 373 283 201 113 NA [4,] 694 673 615 550 491 449 371 279 198 117 NA [5,] 714 694 640 580 524 483 407 319 245 174 NA [6,] 713 690 633 569 514 470 391 295 216 143 NA [7,] 698 673 620 551 495 442 357 255 173 103 NA [8,] 688 659 610 536 479 419 327 227 144 81 NA [9,] 689 657 601 535 479 419 321 228 143 92 NA [10,] 689 654 595 526 467 406 317 233 152 104 NA [11,] 690 644 585 517 451 393 300 210 127 81 NA [12,] 692 643 584 512 459 396 304 219 141 89 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -1 -9 -12 -6 -9 -1 -1 0 -3 -2 NA [2,] -1 -2 -5 -14 -5 -8 -16 -13 -8 -14 NA [3,] -2 -4 -6 -10 -5 -3 -2 -4 -3 4 NA [4,] 20 21 25 30 33 34 36 40 47 57 NA [5,] -1 -4 -7 -11 -10 -13 -16 -24 -29 -31 NA [6,] -15 -17 -13 -18 -19 -28 -34 -40 -43 -40 NA [7,] -10 -14 -10 -15 -16 -23 -30 -28 -29 -22 NA [8,] 1 -2 -9 -1 0 0 -6 1 -1 11 NA [9,] 0 -3 -6 -9 -12 -13 -4 5 9 12 NA [10,] 1 -10 -10 -9 -16 -13 -17 -23 -25 -23 NA [11,] 2 -1 -1 -5 8 3 4 9 14 8 NA [12,] -4 -5 -4 -2 2 -6 -8 -7 -12 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/1w4lg1280157064.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/2w4lg1280157064.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/3w4lg1280157064.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/4pe311280157064.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,] 129.0 127.0 113 117 174.0 143 103.0 81 92 104.0 81 89.0 [2,] 296.0 296.0 283 279 319.0 295 255.0 227 228 233.0 210 219.0 [3,] 485.5 480.5 474 470 503.5 492 468.5 449 449 436.5 422 427.5 [4,] 638.0 626.0 621 615 640.0 633 620.0 610 601 595.0 585 584.0 [5,] 698.0 697.0 696 694 714.0 713 698.0 688 689 689.0 690 692.0 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 314.6232 315.6188 305.1217 302.121 343.1156 323.1217 286.1314 257.6379 [2,] 656.3768 645.3812 642.8783 637.879 663.8844 660.8783 650.8686 640.3621 [,9] [,10] [,11] [,12] [1,] 262.6343 255.6304 234.6350 245.1314 [2,] 635.3657 617.3696 609.3650 609.8686 $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(129, 296, 485.5, 638, 698, 127, 296, 480.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5pe311280157064.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,] -12.0 -16 -6.0 20.0 -31 -43.0 -30 -6.0 -13.0 -25.0 -5.0 -8 [2,] -9.0 -14 -5.0 25.0 -24 -40.0 -28 -2.0 -9.0 -23.0 -1.0 -7 [3,] -2.5 -8 -3.5 33.5 -12 -23.5 -19 -0.5 -3.5 -14.5 3.5 -5 [4,] -1.0 -5 -2.0 40.0 -7 -17.0 -14 1.0 5.0 -10.0 8.0 -4 [5,] 0.0 -1 -2.0 57.0 -1 -13.0 -10 1.0 12.0 1.0 14.0 -2 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -6.497119 -12.496759 -4.998920 26.0054 -20.493878 -34.99172 -25.99496 [2,] 1.497119 -3.503241 -2.001080 40.9946 -3.506122 -12.00828 -12.00504 [,8] [,9] [,10] [,11] [,12] [1,] -1.9989196 -10.494958 -20.995318 -0.9967588 -6.58 [2,] 0.9989196 3.494958 -8.004682 7.9967588 -3.42 $out [1] -10 4 -9 11 2 -12 $group [1] 3 3 8 8 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-12, -9, -2.5, -1, 0, -16, -14, -8, -5, -1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6hnkm1280157064.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,] 688.0 643.0 584.0 512.0 451.0 393.0 300.0 210.0 127.0 81.0 NA [2,] 689.5 655.5 598.0 530.5 473.0 412.5 319.0 227.5 143.5 90.5 NA [3,] 695.0 673.0 617.5 550.5 493.0 445.5 364.0 267.0 185.5 108.5 NA [4,] 698.0 683.5 629.5 571.5 505.5 460.5 389.5 295.5 210.5 128.0 NA [5,] 698.0 694.0 640.0 580.0 524.0 483.0 407.0 319.0 245.0 174.0 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 691.1231 660.229 603.1326 531.7996 478.1765 423.6069 331.8445 235.9847 [2,] 698.8769 685.771 631.8674 569.2004 507.8235 467.3931 396.1555 298.0153 [,9] [,10] [,11] [1,] 154.9409 91.396 NA [2,] 216.0591 125.604 NA $out [1] 714 713 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(688, 689.5, 695, 698, 698, 643, 655.5, 673, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7hnkm1280157064.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,] 399.80 422.00 400.2500 [2,] 415.35 442.75 416.8750 [3,] 440.20 469.25 446.0000 [4,] 458.00 483.00 467.0625 [5,] 478.00 503.50 483.0000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 420.7470 450.8917 423.1091 [2,] 459.6530 487.6083 468.8909 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(399.8, 415.35, 440.2, 458, 478, 422, 442.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1w4lg1280157064.ps tmp/1w4lg1280157064.png",intern=TRUE)) character(0) > try(system("convert tmp/2w4lg1280157064.ps tmp/2w4lg1280157064.png",intern=TRUE)) character(0) > try(system("convert tmp/3w4lg1280157064.ps tmp/3w4lg1280157064.png",intern=TRUE)) character(0) > try(system("convert tmp/4pe311280157064.ps tmp/4pe311280157064.png",intern=TRUE)) character(0) > try(system("convert tmp/5pe311280157064.ps tmp/5pe311280157064.png",intern=TRUE)) character(0) > try(system("convert tmp/6hnkm1280157064.ps tmp/6hnkm1280157064.png",intern=TRUE)) character(0) > try(system("convert tmp/7hnkm1280157064.ps tmp/7hnkm1280157064.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.972 1.497 2.264