Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
-9.689129094 -15.9156891 -20.03555923 -12.44264761 -2.468740667 -23.97386436 -3.02331487 -8.793616067 3.890809296 -1.836563963 -1.709285854 4.14432084 -12.03301028 -4.733406668 -5.29923341 -1.665398243 2.528592702 -2.076860346 2.08284686 1.362751533 17.9705551 -13.16742762 9.297539752 -0.631604907 -9.595566233 -4.093437696 -8.465289457 15.56005027 8.440539387 5.496633362 -1.542161774 -4.424920628 9.358249829 5.491823987 21.87928431 6.722498409 7.093000066 -17.05020325 0.883582789 15.92450566 3.185445132 -11.26933787 -14.58358037 -29.97174154 0.914837609 6.714446914 8.903871839 23.13841004 26.84332477 8.351341468 39.34532234 16.99485377 13.45289299 29.686769 -9.403389414 -4.311467414 2.792229304 -21.98399447 -28.17586045 -14.08502665
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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation