Wessa.net

Send output to:


Data:
84728 84412 84092 83430 89981 89635 84728 81466 81781 81781 82133 82764 83746 83746 83115 81466 89981 91279 89319 84728 86692 83746 85075 85710 86372 84728 85075 82764 89981 92261 90301 86692 90617 86372 90301 89981 90964 87355 91279 90964 96852 95523 90301 87670 91279 86372 89981 90617 91946 89004 90617 91599 95208 92261 88337 84092 88021 77221 82448 85390 88337 84092 84092 84092 86372 83115 78839 75261 77857 67724 73933 77541 78204 74595 74910 73933 77221 74910 70355 67062 72630 60537 68390 71968 71968 67724 63799 63484 67062 63799 57595 53319 57911 47115 56928 62150 63799 60191 55631 58893 60191 59208 49391 44835 48093 38280 48413 52022 54964 50057 45466 48093 49391 46795 36982 32706 36631 25835 37613 44835


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	0 seconds
R Server	Big Analytics Cloud Computing Center