Wessa.net

Send output to:


Data:
19064 18993 18921 18772 20246 20168 19064 18330 18401 18401 18480 18622 18843 18843 18701 18330 20246 20538 20097 19064 19506 18843 19142 19285 19434 19064 19142 18622 20246 20759 20318 19506 20389 19434 20318 20246 20467 19655 20538 20467 21792 21493 20318 19726 20538 19434 20246 20389 20688 20026 20389 20610 21422 20759 19876 18921 19805 17375 18551 19213 19876 18921 18921 18921 19434 18701 17739 16934 17518 15238 16635 17447 17596 16784 16855 16635 17375 16855 15830 15089 16342 13621 15388 16193 16193 15238 14355 14284 15089 14355 12959 11997 13030 10601 12809 13984 14355 13543 12517 13251 13543 13322 11113 10088 10821 8613 10893 11705 12367 11263 10230 10821 11113 10529 8321 7359 8242 5813 8463 10088


Sample Range: (leave blank to include all observations)
From:
To:

# periods


Chart options



R Code
par1 <- as.numeric(par1) (n <- length(x)) (np <- floor(n / par1)) arr <- array(NA,dim=c(par1,np+1)) ari <- array(0,dim=par1) j <- 0 for (i in 1:n) { j = j + 1 ari[j] = ari[j] + 1 arr[j,ari[j]] <- x[i] if (j == par1) j = 0 } ari arr arr.sd <- array(NA,dim=par1) arr.range <- array(NA,dim=par1) arr.iqr <- array(NA,dim=par1) for (j in 1:par1) { arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE)) arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE) arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE) } overall.sd <- sqrt(var(x)) overall.range <- max(x) - min(x) overall.iqr <- quantile(x,0.75) - quantile(x,0.25) bitmap(file='plot1.png') plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index') mtext(paste('# blocks = ',np)) abline(overall.sd,0) dev.off() bitmap(file='plot2.png') plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index') mtext(paste('# blocks = ',np)) abline(overall.range,0) dev.off() bitmap(file='plot3.png') plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index') mtext(paste('# blocks = ',np)) abline(overall.iqr,0) dev.off() bitmap(file='plot4.png') 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')) dev.off() bitmap(file='plot5.png') 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')) dev.off() bitmap(file='plot6.png') z <- data.frame(cbind(arr.sd,arr.range,arr.iqr)) names(z) <- list('S.D.','Range','IQR') (boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots')) dev.off()

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	Big Analytics Cloud Computing Center