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(1149822 + ,1086979 + ,1276674 + ,1522522 + ,1742117 + ,1737275 + ,1979900 + ,2061036 + ,1867943 + ,1707752 + ,1298756 + ,1281814 + ,1281151 + ,1164976 + ,1454329 + ,1645288 + ,1817743 + ,1895785 + ,2236311 + ,2295951 + ,2087315 + ,1980891 + ,1465446 + ,1445026 + ,1488120 + ,1338333 + ,1715789 + ,1806090 + ,2083316 + ,2092278 + ,2430800 + ,2424894 + ,2299016 + ,2130688 + ,1652221 + ,1608162 + ,1647074 + ,1479691 + ,1884978 + ,2007898 + ,2208954 + ,2217164 + ,2534291 + ,2560312 + ,2429069 + ,2315077 + ,1799608 + ,1772590 + ,1744799 + ,1659093 + ,2099821 + ,2135736 + ,2427894 + ,2468882 + ,2703217 + ,2766841 + ,2655236 + ,2550373 + ,2052097 + ,1998055 + ,1920748 + ,1876694 + ,2380930 + ,2467402 + ,2770771 + ,2781340 + ,3143926 + ,3172235 + ,2952540 + ,2920877 + ,2384552 + ,2248987 + ,2208616 + ,2178756 + ,2632870 + ,2706905 + ,3029745 + ,3015402 + ,3391414 + ,3507805 + ,3177852 + ,3142961 + ,2545815 + ,2414007 + ,2372578 + ,2332664 + ,2825328 + ,2901478 + ,3263955 + ,3226738 + ,3610786 + ,3709274 + ,3467185 + ,3449646 + ,2802951 + ,2462530 + ,2490645 + ,2561520 + ,3067554 + ,3226951 + ,3546493 + ,3492787 + ,3952263 + ,3932072 + ,3720284 + ,3651555 + ,2914972 + ,2713514 + ,2703997 + ,2591373 + ,3163748 + ,3355137 + ,3613702 + ,3686773 + ,4098716 + ,4063517 + ,3551489 + ,3226663 + ,2656842 + ,2597484 + ,2572399 + ,2596631 + ,3165225 + ,3303145 + ,3698247 + ,3668631 + ,4130433 + ,4131400 + ,3864358 + ,3721110 + ,2892532 + ,2843451 + ,2747502 + ,2668775 + ,3018602 + ,3013392 + ,3393657 + ,3544233 + ,4075832 + ,4032923 + ,3734509 + ,3761285 + ,2970090 + ,2847849 + ,2741680 + ,2830639 + ,3257673 + ,3480085 + ,3843271 + ,3796961 + ,4337767 + ,4243630 + ,3927202 + ,3915296 + ,3087396 + ,2963792 + ,2955792 + ,2829925 + ,3281195 + ,3548011 + ,4059648 + ,3941175 + ,4528594 + ,4433151 + ,4145737 + ,4077132 + ,3198519 + ,3078660 + ,3028202 + ,2858642 + ,3398954 + ,3808883 + ,4175961 + ,4227542 + ,4744616 + ,4608012 + ,4295049 + ,4201144 + ,3353276 + ,3286851 + ,3169889 + ,3051720 + ,3695426 + ,3905501 + ,4296458 + ,4246247 + ,4921849 + ,4821446 + ,4425064 + ,4379099 + ,3472889 + ,3359160 + ,3200944 + ,3153170 + ,3741498 + ,3918719 + ,4403449 + ,4400407 + ,4847473 + ,4716136 + ,4297440 + ,4272253 + ,3271834 + ,3168388 + ,2911748 + ,2720999 + ,3199918 + ,3672623 + ,3892013 + ,3850845 + ,4532467 + ,4484739 + ,4014972 + ,3983758 + ,3158459 + ,3100569 + ,2935404 + ,2855719 + ,3465611 + ,3006985 + ,4095110 + ,4104793 + ,4730788 + ,4642726 + ,4246919 + ,4308117) > 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] 226 > (np <- floor(n / par1)) [1] 18 > 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] 19 19 19 19 19 19 19 19 19 19 18 18 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 1149822 1281151 1488120 1647074 1744799 1920748 2208616 2372578 2490645 [2,] 1086979 1164976 1338333 1479691 1659093 1876694 2178756 2332664 2561520 [3,] 1276674 1454329 1715789 1884978 2099821 2380930 2632870 2825328 3067554 [4,] 1522522 1645288 1806090 2007898 2135736 2467402 2706905 2901478 3226951 [5,] 1742117 1817743 2083316 2208954 2427894 2770771 3029745 3263955 3546493 [6,] 1737275 1895785 2092278 2217164 2468882 2781340 3015402 3226738 3492787 [7,] 1979900 2236311 2430800 2534291 2703217 3143926 3391414 3610786 3952263 [8,] 2061036 2295951 2424894 2560312 2766841 3172235 3507805 3709274 3932072 [9,] 1867943 2087315 2299016 2429069 2655236 2952540 3177852 3467185 3720284 [10,] 1707752 1980891 2130688 2315077 2550373 2920877 3142961 3449646 3651555 [11,] 1298756 1465446 1652221 1799608 2052097 2384552 2545815 2802951 2914972 [12,] 1281814 1445026 1608162 1772590 1998055 2248987 2414007 2462530 2713514 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] 2703997 2572399 2747502 2741680 2955792 3028202 3169889 3200944 2911748 [2,] 2591373 2596631 2668775 2830639 2829925 2858642 3051720 3153170 2720999 [3,] 3163748 3165225 3018602 3257673 3281195 3398954 3695426 3741498 3199918 [4,] 3355137 3303145 3013392 3480085 3548011 3808883 3905501 3918719 3672623 [5,] 3613702 3698247 3393657 3843271 4059648 4175961 4296458 4403449 3892013 [6,] 3686773 3668631 3544233 3796961 3941175 4227542 4246247 4400407 3850845 [7,] 4098716 4130433 4075832 4337767 4528594 4744616 4921849 4847473 4532467 [8,] 4063517 4131400 4032923 4243630 4433151 4608012 4821446 4716136 4484739 [9,] 3551489 3864358 3734509 3927202 4145737 4295049 4425064 4297440 4014972 [10,] 3226663 3721110 3761285 3915296 4077132 4201144 4379099 4272253 3983758 [11,] 2656842 2892532 2970090 3087396 3198519 3353276 3472889 3271834 3158459 [12,] 2597484 2843451 2847849 2963792 3078660 3286851 3359160 3168388 3100569 [,19] [1,] 2935404 [2,] 2855719 [3,] 3465611 [4,] 3006985 [5,] 4095110 [6,] 4104793 [7,] 4730788 [8,] 4642726 [9,] 4246919 [10,] 4308117 [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -62843 -116175 -149787 -167383 -85706 -44054 -29860 -39914 70875 [2,] 189695 289353 377456 405287 440728 504236 454114 492664 506034 [3,] 245848 190959 90301 122920 35915 86472 74035 76150 159397 [4,] 219595 172455 277226 201056 292158 303369 322840 362477 319542 [5,] -4842 78042 8962 8210 40988 10569 -14343 -37217 -53706 [6,] 242625 340526 338522 317127 234335 362586 376012 384048 459476 [7,] 81136 59640 -5906 26021 63624 28309 116391 98488 -20191 [8,] -193093 -208636 -125878 -131243 -111605 -219695 -329953 -242089 -211788 [9,] -160191 -106424 -168328 -113992 -104863 -31663 -34891 -17539 -68729 [10,] -408996 -515445 -478467 -515469 -498276 -536325 -597146 -646695 -736583 [11,] -16942 -20420 -44059 -27018 -54042 -135565 -131808 -340421 -201458 [12,] -663 43094 38912 -27791 -77307 -40371 -41429 28115 -9517 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] -112624 24232 -78727 88959 -125867 -169560 -118169 -47774 -190749 [2,] 572375 568594 349827 427034 451270 540312 643706 588328 478919 [3,] 191389 137920 -5210 222412 266816 409929 210075 177221 472705 [4,] 258565 395102 380265 363186 511637 367078 390957 484730 219390 [5,] 73071 -29616 150576 -46310 -118473 51581 -50211 -3042 -41168 [6,] 411943 461802 531599 540806 587419 517074 675602 447066 681622 [7,] -35199 967 -42909 -94137 -95443 -136604 -100403 -131337 -47728 [8,] -512028 -267042 -298414 -316428 -287414 -312963 -396382 -418696 -469767 [9,] -324826 -143248 26776 -11906 -68605 -93905 -45965 -25187 -31214 [10,] -569821 -828578 -791195 -827900 -878613 -847868 -906210 -1000419 -825299 [11,] -59358 -49081 -122241 -123604 -119859 -66425 -113729 -103446 -57890 [12,] -25085 -95949 -106169 -8000 -50458 -116962 -158216 -256640 -165165 [,19] [1,] -79685 [2,] 609892 [3,] -458626 [4,] 1088125 [5,] 9683 [6,] 625995 [7,] -88062 [8,] -395807 [9,] 61198 [10,] NA [11,] NA [12,] 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/1larq1291985469.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/www/html/freestat/rcomp/tmp/2larq1291985469.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/www/html/freestat/rcomp/tmp/3e28t1291985469.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/www/html/freestat/rcomp/tmp/4e28t1291985469.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,] 1149822 1086979 1276674 1522522 1742117 1737275 1979900 2061036 1867943 [2,] 1832774 1767894 2240376 2301569 2599332 2625111 2923572 2969538 2803888 [3,] 2572399 2591373 3067554 3013392 3546493 3544233 4075832 4032923 3720284 [4,] 2923576 2830282 3269434 3514048 3975830 3896010 4530530 4458945 4080354 [5,] 3200944 3153170 3741498 3918719 4403449 4400407 4921849 4821446 4425064 [,10] [,11] [,12] [1,] 1707752 1298756 1281814 [2,] 2735625 2052097 1998055 [3,] 3651555 2847742 2655499 [4,] 4030445 3158459 3078660 [5,] 4379099 3472889 3359160 $n [1] 19 19 19 19 19 19 19 19 19 19 18 18 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 2177008 2206282 2694544 2573896 3047544 3083562 3493347 3493047 3257594 [2,] 2967790 2976464 3440564 3452888 4045442 4004904 4658317 4572799 4182974 [,10] [,11] [,12] [1,] 3182213 2435722 2253071 [2,] 4120897 3259761 3057927 $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(1149822, 1832773.5, 2572399, 2923576, 3200944, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5obpv1291985469.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,] -190749 289353.0 -5210.0 172455.0 -118473.0 234335.0 -136604.0 -512028 [2,] -122018 416160.5 81311.0 267895.5 -39192.5 351556.0 -91099.5 -362880 [3,] -79685 478919.0 159397.0 322840.0 -3042.0 447066.0 -20191.0 -287414 [4,] -41984 554453.0 216243.5 385611.0 25778.5 536202.5 43974.5 -210212 [5,] 70875 643706.0 409929.0 511637.0 78042.0 681622.0 116391.0 -111605 [,9] [,10] [,11] [,12] [1,] -168328.0 -1000419 -201458.0 -165165 [2,] -110208.0 -828578 -123604.0 -106169 [3,] -68605.0 -691639 -84935.5 -40900 [4,] -28200.5 -515469 -49081.0 -8000 [5,] 61198.0 -408996 -16942.0 43094 $n [1] 19 19 19 19 19 19 19 19 19 18 18 18 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -108695.47 428791.2 110487.1 280170.9 -26592.48 380135.9 -69152.2 [2,] -50674.53 529046.8 208306.9 365509.1 20508.48 513996.1 28770.2 [,8] [,9] [,10] [,11] [,12] [1,] -342752.6 -98330.82 -808243.8 -112688.58 -77459.075 [2,] -232075.4 -38879.18 -575034.2 -57182.42 -4340.925 $out [1] 88959 189695 472705 -458626 1088125 150576 -324826 -340421 -256640 $group [1] 1 2 3 3 4 5 9 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-190749, -122018, -79685, -41984, 70875, 289353, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6obpv1291985469.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,] 1086979 1164976 1338333 1479691 1659093 1876694 2178756 2332664 2490645 [2,] 1279244 1449678 1630192 1786099 2025076 2314958 2479911 2632740 2814243 [3,] 1615137 1731516 1944703 2108426 2281815 2619086 2861154 3064108 3359869 [4,] 1805030 2034103 2214852 2372073 2602804 2936708 3160406 3458416 3685920 [5,] 2061036 2295951 2430800 2560312 2766841 3172235 3507805 3709274 3952263 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] 2591373 2572399 2668775 2741680 2829925 2858642 3051720 3153170 2720999 [2,] 2680420 2867992 2908970 3025594 3138590 3320064 3416024 3236389 3129514 [3,] 3290900 3485888 3206130 3638523 3744593 3992422 4075874 4095486 3761734 [4,] 3650238 3792734 3747897 3921249 4111434 4261296 4402082 4401928 3999365 [5,] 4098716 4131400 4075832 4337767 4528594 4744616 4921849 4847473 4532467 [,19] [1,] 2855719 [2,] 3006985 [3,] 4099952 [4,] 4308117 [5,] 4730788 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 1375322 1464955 1678035 1841159 2018309 2335502 2550775 2687512 2962291 [2,] 1854952 1998076 2211371 2375693 2545321 2902671 3171532 3440704 3757447 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] 2848560 3064107 2823489 3230009 3300872 3563120 3626127 3563876 3364989 [2,] 3733240 3907669 3588770 4047037 4188314 4421724 4525621 4627096 4158479 [,19] [1,] 3449854 [2,] 4750049 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" NA Warning message: In bxp(list(stats = c(1086979, 1279244, 1615137, 1805030, 2061036, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7obpv1291985469.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,] 2307174 2572399 2299088 [2,] 2560254 2751620 2596659 [3,] 3099875 3305894 3084184 [4,] 3390917 3685920 3412578 [5,] 3733234 4075832 3727051 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 2721004 2879753 2712038 [2,] 3478746 3732034 3456331 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(2307173.63157895, 2560254, 3099875.13157895, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1larq1291985469.ps tmp/1larq1291985469.png",intern=TRUE)) character(0) > try(system("convert tmp/2larq1291985469.ps tmp/2larq1291985469.png",intern=TRUE)) character(0) > try(system("convert tmp/3e28t1291985469.ps tmp/3e28t1291985469.png",intern=TRUE)) character(0) > try(system("convert tmp/4e28t1291985469.ps tmp/4e28t1291985469.png",intern=TRUE)) character(0) > try(system("convert tmp/5obpv1291985469.ps tmp/5obpv1291985469.png",intern=TRUE)) character(0) > try(system("convert tmp/6obpv1291985469.ps tmp/6obpv1291985469.png",intern=TRUE)) character(0) > try(system("convert tmp/7obpv1291985469.ps tmp/7obpv1291985469.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.982 1.473 2.254