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. 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(1345.9,1629.3,1629.0,1549.2,1499.4,1440.1,1124.5,1279.5,1266.4,1420.0,1547.0,1392.9,1167.6,1402.9,1377.6,1548.2,1348.0,1397.1,1249.1,1257.6,1319.1,1344.3,1439.3,1151.9,1174.5,1125.9,1230.0,1325.4,1080.4,1426.6,1200.9,1183.6,1130.4,1196.8,1126.3,1079.8,1109.7,1055.2,1232.9,1144.5,1156.0,1268.1,1153.4,1238.2,996.6,1237.2,1154.0,1037.0,1035.9,1076.7,1192.2,1107.0,1106.9,1332.0,974.0,954.2,913.4,1049.1,1088.0,995.1) > 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] 60 > (np <- floor(n / par1)) [1] 5 > 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] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1345.9 1167.6 1174.5 1109.7 1035.9 NA [2,] 1629.3 1402.9 1125.9 1055.2 1076.7 NA [3,] 1629.0 1377.6 1230.0 1232.9 1192.2 NA [4,] 1549.2 1548.2 1325.4 1144.5 1107.0 NA [5,] 1499.4 1348.0 1080.4 1156.0 1106.9 NA [6,] 1440.1 1397.1 1426.6 1268.1 1332.0 NA [7,] 1124.5 1249.1 1200.9 1153.4 974.0 NA [8,] 1279.5 1257.6 1183.6 1238.2 954.2 NA [9,] 1266.4 1319.1 1130.4 996.6 913.4 NA [10,] 1420.0 1344.3 1196.8 1237.2 1049.1 NA [11,] 1547.0 1439.3 1126.3 1154.0 1088.0 NA [12,] 1392.9 1151.9 1079.8 1037.0 995.1 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 283.4 235.3 -48.6 -54.5 40.8 NA [2,] -0.3 -25.3 104.1 177.7 115.5 NA [3,] -79.8 170.6 95.4 -88.4 -85.2 NA [4,] -49.8 -200.2 -245.0 11.5 -0.1 NA [5,] -59.3 49.1 346.2 112.1 225.1 NA [6,] -315.6 -148.0 -225.7 -114.7 -358.0 NA [7,] 155.0 8.5 -17.3 84.8 -19.8 NA [8,] -13.1 61.5 -53.2 -241.6 -40.8 NA [9,] 153.6 25.2 66.4 240.6 135.7 NA [10,] 127.0 95.0 -70.5 -83.2 38.9 NA [11,] -154.1 -287.4 -46.5 -117.0 -92.9 NA [12,] -225.3 22.6 29.9 -1.1 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/178891229599071.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/23qyg1229599071.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/3c3ps1229599071.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/4db8d1229599071.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] [1,] 1035.9 1055.2 1192.2 1107.0 1080.4 1268.1 1124.5 1183.6 913.4 1049.1 [2,] 1109.7 1076.7 1230.0 1144.5 1106.9 1332.0 1124.5 1183.6 996.6 1196.8 [3,] 1167.6 1125.9 1232.9 1325.4 1156.0 1397.1 1153.4 1238.2 1130.4 1237.2 [4,] 1174.5 1402.9 1377.6 1548.2 1348.0 1426.6 1200.9 1257.6 1266.4 1344.3 [5,] 1174.5 1629.3 1377.6 1549.2 1499.4 1440.1 1249.1 1279.5 1319.1 1420.0 [,11] [,12] [1,] 1088.0 995.1 [2,] 1126.3 1037.0 [3,] 1154.0 1079.8 [4,] 1439.3 1151.9 [5,] 1547.0 1151.9 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1121.812 895.4079 1128.606 1040.147 985.6393 1330.256 1099.416 1185.912 [2,] 1213.388 1356.3921 1337.194 1610.653 1326.3607 1463.944 1207.384 1290.488 [,9] [,10] [,11] [,12] [1,] 939.76 1132.977 932.835 998.612 [2,] 1321.04 1341.423 1375.165 1160.988 $out [1] 1345.9 1629.0 974.0 954.2 1392.9 $group [1] 1 3 7 8 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(1035.9, 1109.7, 1167.6, 1174.5, 1174.5, 1055.2, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5ezax1229599071.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] [1,] -54.5 -25.3 -88.4 -245.0 -59.3 -358.0 -19.8 -53.2 25.2 -83.2 -154.1 [2,] -48.6 -0.3 -85.2 -200.2 49.1 -315.6 -17.3 -53.2 66.4 -70.5 -154.1 [3,] 40.8 104.1 -79.8 -49.8 112.1 -225.7 8.5 -40.8 135.7 38.9 -117.0 [4,] 235.3 115.5 95.4 -0.1 225.1 -148.0 84.8 -13.1 153.6 95.0 -92.9 [5,] 283.4 177.7 170.6 11.5 346.2 -114.7 155.0 -13.1 240.6 127.0 -46.5 [,12] [1,] -225.30 [2,] -113.20 [3,] 10.75 [4,] 26.25 [5,] 29.90 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -159.8030 22.27601 -207.41151 -191.19016 -12.26116 -344.1257 -63.6436 [2,] 241.4030 185.92399 47.81151 91.59016 236.46116 -107.2743 80.6436 [,8] [,9] [,10] [,11] [,12] [1,] -69.13456 74.0847 -78.04188 -160.24377 -99.4155 [2,] -12.46544 197.3153 155.84188 -73.75623 120.9155 $out [1] 61.5 -241.6 -287.4 $group [1] 8 8 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-54.5, -48.5999999999999, 40.8000000000000, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6t72b1229599071.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] [1,] 1124.50 1151.90 1079.80 996.60 913.40 NA [2,] 1312.70 1253.35 1126.10 1082.45 984.55 NA [3,] 1430.05 1346.15 1179.05 1153.70 1062.90 NA [4,] 1548.10 1400.00 1215.45 1235.05 1106.95 NA [5,] 1629.30 1548.20 1325.40 1268.10 1192.20 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1322.682 1279.262 1138.297 1084.098 1007.073 NA [2,] 1537.418 1413.038 1219.803 1223.302 1118.727 NA $out [1] 1426.6 1332.0 $group [1] 3 5 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(1124.5, 1312.7, 1430.05, 1548.1, 1629.3, 1151.9, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/79h701229599071.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,] 1125.18 1079.8 1094.450 [2,] 1153.55 1141.9 1152.400 [3,] 1243.81 1161.8 1233.625 [4,] 1301.63 1237.7 1293.300 [5,] 1372.78 1325.4 1379.300 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1176.270 1118.105 1169.360 [2,] 1311.350 1205.495 1297.890 $out [1] 1397.1 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1125.18, 1153.55, 1243.81, 1301.63, 1372.78, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/178891229599071.ps tmp/178891229599071.png") > system("convert tmp/23qyg1229599071.ps tmp/23qyg1229599071.png") > system("convert tmp/3c3ps1229599071.ps tmp/3c3ps1229599071.png") > system("convert tmp/4db8d1229599071.ps tmp/4db8d1229599071.png") > system("convert tmp/5ezax1229599071.ps tmp/5ezax1229599071.png") > system("convert tmp/6t72b1229599071.ps tmp/6t72b1229599071.png") > system("convert tmp/79h701229599071.ps tmp/79h701229599071.png") > > > proc.time() user system elapsed 1.863 1.465 2.196