R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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(1069108,1059362,1049495,1029082,1231089,1220388,1069108,968521,978233,978233,989056,1008514,1069108,1049495,1079775,1129547,1412683,1412683,1352244,1291650,1341422,1401983,1412683,1442964,1533833,1473244,1473244,1564114,1816014,1836427,1785734,1664578,1755420,1755420,1765165,1816014,1856040,1876453,1876453,1937014,2169452,2229890,2239603,2088322,2169452,2139171,2078582,2209478,2239603,2188909,2199610,2269916,2532640,2663352,2663352,2602913,2693660,2602913,2552092,2744509,2774634,2703378,2884967,2956228,3168103,3308711,3289258,3278430,3359559,3349692,3228692,3410253,3470847,3410253,3662154,3783309,4065362,4176650,4146492,4085897,4136624,4197185,3995056,4156204,4257779,4216798,4479366,4570080,4953837,5024138,4933418,4984117,5014398,5044678,4852261,5033856,5134437,5033856,5326854,5417607,5811037,5871631,5891089,5992631,5992631,6032657,5851063,5941938,6002377,5891089,6214246,6274812,6678021,6749283,6849864,6940739,6950451,6961152,6779563,6961152) > 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] [1,] 1069108 1069108 1533833 1856040 2239603 2774634 3470847 4257779 5134437 [2,] 1059362 1049495 1473244 1876453 2188909 2703378 3410253 4216798 5033856 [3,] 1049495 1079775 1473244 1876453 2199610 2884967 3662154 4479366 5326854 [4,] 1029082 1129547 1564114 1937014 2269916 2956228 3783309 4570080 5417607 [5,] 1231089 1412683 1816014 2169452 2532640 3168103 4065362 4953837 5811037 [6,] 1220388 1412683 1836427 2229890 2663352 3308711 4176650 5024138 5871631 [7,] 1069108 1352244 1785734 2239603 2663352 3289258 4146492 4933418 5891089 [8,] 968521 1291650 1664578 2088322 2602913 3278430 4085897 4984117 5992631 [9,] 978233 1341422 1755420 2169452 2693660 3359559 4136624 5014398 5992631 [10,] 978233 1401983 1755420 2139171 2602913 3349692 4197185 5044678 6032657 [11,] 989056 1412683 1765165 2078582 2552092 3228692 3995056 4852261 5851063 [12,] 1008514 1442964 1816014 2209478 2744509 3410253 4156204 5033856 5941938 [,10] [,11] [1,] 6002377 NA [2,] 5891089 NA [3,] 6214246 NA [4,] 6274812 NA [5,] 6678021 NA [6,] 6749283 NA [7,] 6849864 NA [8,] 6940739 NA [9,] 6950451 NA [10,] 6961152 NA [11,] 6779563 NA [12,] 6961152 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -9746 -19613 -60589 20413 -50694 -71256 -60594 -40981 -100581 [2,] -9867 30280 0 0 10701 181589 251901 262568 292998 [3,] -20413 49772 90870 60561 70306 71261 121155 90714 90753 [4,] 202007 283136 251900 232438 262724 211875 282053 383757 393430 [5,] -10701 0 20413 60438 130712 140608 111288 70301 60594 [6,] -151280 -60439 -50693 9713 0 -19453 -30158 -90720 19458 [7,] -100587 -60594 -121156 -151281 -60439 -10828 -60595 50699 101542 [8,] 9712 49772 90842 81130 90747 81129 50727 30281 0 [9,] 0 60561 0 -30281 -90747 -9867 60561 30280 40026 [10,] 10823 10700 9745 -60589 -50821 -121000 -202129 -192417 -181594 [11,] 19458 30281 50849 130896 192417 181561 161148 181595 90875 [12,] 60594 90869 40026 30125 30125 60594 101575 100581 60439 [,10] [,11] [1,] -111288 NA [2,] 323157 NA [3,] 60566 NA [4,] 403209 NA [5,] 71262 NA [6,] 100581 NA [7,] 90875 NA [8,] 9712 NA [9,] 10701 NA [10,] -181589 NA [11,] 181589 NA [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/wessaorg/rcomp/tmp/1svut1310302060.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2q9qo1310302060.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3agmh1310302060.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4cqmg1310302060.ps",horizontal=F,onefile=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,] 1069108 1049495 1049495 1029082 1231089 1220388 1069108 968521 978233 [2,] 1533833 1473244 1473244 1564114 1816014 1836427 1785734 1664578 1755420 [3,] 2507118 2446144 2542288 2613072 2850372 2986032 2976305 2940672 3026610 [4,] 4257779 4216798 4479366 4570080 4953837 5024138 4933418 4984117 5014398 [5,] 6002377 5891089 6214246 6274812 6678021 6749283 6849864 6940739 6950451 [,10] [,11] [,12] [1,] 978233 989056 1008514 [2,] 1755420 1765165 1816014 [3,] 2976302 2890392 3077381 [4,] 5044678 4852261 5033856 [5,] 6961152 6779563 6961152 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 1146126 1075355 1040310 1111172 1282590 1393324 1403597 1282097 1398294 [2,] 3868111 3816932 4044267 4114972 4418153 4578739 4549013 4599246 4654925 [,10] [,11] [,12] [1,] 1332858 1347956 1469619 [2,] 4619747 4432828 4685143 $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(1069108, 1533833, 2507118.5, 4257779, 6002377, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/592c01310302060.ps",horizontal=F,onefile=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] [1,] -111288.0 -9867.0 49772.0 202007.0 -10701.0 -151280.0 -151281.0 0.0 [2,] -71256.0 0.0 60561.0 232438.0 20413.0 -60439.0 -100587.0 9712.0 [3,] -55641.5 105934.5 70783.5 272388.5 65447.5 -24805.5 -60516.5 50249.5 [4,] -19613.0 262568.0 90753.0 383757.0 111288.0 9713.0 50699.0 81130.0 [5,] 20413.0 323157.0 121155.0 403209.0 140608.0 100581.0 101542.0 90842.0 [,9] [,10] [,11] [,12] [1,] -30281.0 -202129.0 19458 30125 [2,] -9867.0 -181594.0 50849 40026 [3,] 5350.5 -90794.5 146022 60594 [4,] 40026.0 9745.0 181589 90869 [5,] 60561.0 10823.0 192417 101575 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -81444.4 -25254.94 55698.37 196783.5 20042.73 -59856.24 -136105.02 [2,] -29838.6 237123.94 85868.63 347993.5 110852.27 10245.24 15072.02 [,8] [,9] [,10] [,11] [,12] [1,] 14566.22 -19578.03 -186395.093 80699.08 33816.69 [2,] 85932.78 30279.03 4806.093 211344.92 87371.31 $out [1] -20413 -90747 $group [1] 3 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-111288, -71256, -55641.5, -19613, 20413, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6f68w1310302060.ps",horizontal=F,onefile=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,] 968521.0 1049495 1473244 1856040 2188909 2703378 3410253 4216798 5033856 [2,] 983644.5 1104661 1548974 1906734 2254760 2920598 3722732 4524723 5372230 [3,] 1039288.5 1346833 1755420 2113746 2577502 3253561 4075630 4943628 5861347 [4,] 1069108.0 1412683 1800874 2189465 2663352 3329202 4151348 5019268 5967284 [5,] 1069108.0 1442964 1836427 2239603 2744509 3410253 4197185 5044678 6032657 [,10] [,11] [1,] 5891089 NA [2,] 6244529 NA [3,] 6764423 NA [4,] 6945595 NA [5,] 6961152 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 1000308 1206342 1640526 1984791 2391141 3067194 3880135 4718062 5589939 [2,] 1078269 1487324 1870314 2242702 2763864 3439928 4271124 5169193 6132755 [,10] [,11] [1,] 6444662 NA [2,] 7084184 NA $out [1] 1231089 1220388 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(968521, 983644.5, 1039288.5, 1069108, 1069108, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7ergi1310302060.ps",horizontal=F,onefile=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,] 2890284 2446144 2794604 [2,] 3058894 2577680 2969959 [3,] 3386802 2915532 3291491 [4,] 3442747 2981168 3334511 [5,] 3472488 3077381 3373529 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 3211724 2731498 3125216 [2,] 3561880 3099565 3457766 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(2890283.7, 3058893.65, 3386801.8, 3442746.7, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1svut1310302060.ps tmp/1svut1310302060.png",intern=TRUE)) character(0) > try(system("convert tmp/2q9qo1310302060.ps tmp/2q9qo1310302060.png",intern=TRUE)) character(0) > try(system("convert tmp/3agmh1310302060.ps tmp/3agmh1310302060.png",intern=TRUE)) character(0) > try(system("convert tmp/4cqmg1310302060.ps tmp/4cqmg1310302060.png",intern=TRUE)) character(0) > try(system("convert tmp/592c01310302060.ps tmp/592c01310302060.png",intern=TRUE)) character(0) > try(system("convert tmp/6f68w1310302060.ps tmp/6f68w1310302060.png",intern=TRUE)) character(0) > try(system("convert tmp/7ergi1310302060.ps tmp/7ergi1310302060.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.583 0.332 1.920