Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
100 101 102.01 103.0301 104.060401 105.101005 106.1520151 107.2135352 108.2856706 109.3685273 110.4622125 111.5668347 112.682503 113.809328 114.9474213 116.0968955 117.2578645 118.4304431 119.6147476 120.810895 122.019004 123.239194 124.471586 125.7163018 126.9734649 128.2431995 129.5256315 130.8208878 132.1290967 133.4503877 134.7848915 136.1327404 137.4940679 138.8690085 140.2576986 141.6602756 143.0768784 144.5076471 145.9527236 147.4122509 148.8863734 150.3752371 151.8789895 153.3977794 154.9317572 156.4810747 158.0458855 159.6263443 161.2226078 162.8348338 164.4631822 166.107814 167.7688921 169.4465811 171.1410469 172.8524573 174.5809819 176.3267917 178.0900597 179.8709603 181.6696699 183.4863666 185.3212302 187.1744425 189.0461869 190.9366488 192.8460153 194.7744755 196.7222202 198.6894424 200.6763368 202.6831002 204.7099312 800 206.7570305 208.8246008 210.9128468 213.0219753 215.1521951 217.303717 219.4767542 221.6715217 223.8882369 226.1271193 228.3883905 230.6722744 232.9789971 235.3087871 237.661875 240.0384937 242.4388787 244.8632675 247.3119001 249.7850191 252.2828693 254.805698
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