Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
2115 2180 2341 2402 2453 2466 2473 2507 2513 2523 2537 2544 2564 2582 2618 2622 2648 2679 2716 2772 2783 2811 2851 2856 2904 2932 2935 2985 2995 3003 3036 3127 3176 3228 3236 3238 3251 3252 3287 3318 3330 3362 3426 3458 3483 3518 3526 3597 3696 3708 3745 3760 3814 3824 3902 3915 3955 3996 4018 4029 4037 4115 4136 4160 4160 4202 4202 4233 4328 4375 4379 4394 4401 4427 4472 4491 4542 4614 4616 4635 4645 4648 4658 4684 4698 4703 4774 4799 4801 4839 4873 4916 4920 4946 4951 4952 4979 4980 4995 5011 5014 5019 5032 5137 5174 5181 5184 5221 5224 5225 5249 5252 5265 5266 5294 5294 5296 5335 5365 5373 5380 5396 5397 5398 5403 5436 5437 5456 5488 5531 5572 5582 5589 5596 5612 5612 5620 5625 5678 5706 5718 5724 5740 5764 5792 5813 5816 5821 5842 5844 5845 5866 5869 5869 5880 5881 5898 5898 5928 5968 5974 5983 5994 6011 6014 6028 6062 6063 6070 6082 6086 6105 6109 6112 6120 6145 6164 6175 6175 6187 6188 6198 6208 6210 6255 6268 6279 6303 6342 6355 6365 6369 6375 6404 6406 6424 6426 6462 6487 6496 6510 6513 6529 6537 6548 6554 6555 6575 6588 6591 6593 6596 6596 6598 6640 6645 6678 6719 6727 6778 6801 6810 6814 6897 6897 7128 7133 7203
Sample Range:
(leave blank to include all observations)
From:
To:
# blockwidth
Chart options
R Code
par1 <- as.numeric(par1) (n <- length(x)) (np <- floor(n / par1)) 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 arr darr 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 bitmap(file='plot1.png') plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index') mtext(paste('#blocks = ',np)) abline(overall.mean,0) dev.off() bitmap(file='plot2.png') plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index') mtext(paste('#blocks = ',np)) abline(overall.median,0) dev.off() bitmap(file='plot3.png') plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index') mtext(paste('#blocks = ',np)) abline(overall.midrange,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='plot4b.png') 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')) 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.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')) dev.off()
Compute
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
Click here to blog (archive) this computation