Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
11 15 19 16 24 15 17 19 19 28 26 15 26 16 24 25 22 15 21 22 27 26 26 22 21 22 20 21 20 22 21 8 22 18 20 24 17 20 23 20 22 19 15 20 22 17 14 24 17 23 25 16 18 20 18 23 24 23 13 20 20 19 22 22 15 17 19 20 22 21 21 16 20 21 20 23 15 18 22 16 17 24 13 19 20 22 19 21 15 21 24 22 20 21 19 14 25 11 17 22 20 22 15 23 20 22 16 25 18 19 25 21 22 21 22 23 20 6 15 18 24 22 21 23 20 20 18 25 16 20 14 22 26 20 17 22 22 20 17 22 17 22 21 25 11 19 24 17 22 22 17 26 19 20 19 21 24 21 19 13 24 28 27 22 23 19 18 23 21 22 17 15 21 20 26 19 28 21 19 22 21 20 19 11 17 19 20 17 21 21 12 23 22 22 21 20 18 21 24 22 20 17 19 16 19 23 8 22 23 15 17 21 25 18 23 20 21 21 24 22 22 23 17 15 24 22 19 18 21 20 19 19 16 18 23 22 23 20 24 25 25 20 23 21 23 23 11 21 27 19 21 16 22 21 22 16 18 23 24 20 20 18 4 14 22 17 23 20 18 19 20 15 24 21 19 19 27 23 23 20 17 21 23 22 16 20 16 21 19 27 13 17 18 20 22 18 6 22 15 19 17 22 10 21 21 23 18 20 27 13 20 20 22 20 24 23 19 22 24 21 19 20 16 17 25 16 23 20 23 22 15 16 20 23 24 17 19 25 14 18 22 15 27 22 26 16 25 20 19 19 24 14 18 13 19 25 20 17 17 13 20 20 24 25 19 20 20 22 18 21 20 11 18 22 21 15 23 18 23 19 23 26 19 26 20 20 23 24 26 23 19 25 23 19 27 23 24 20 16 22 26 26 24 20 20 12 21 27 26 17 20 18 28 24 24 24 12 26 23 13 23 16 23 18 25 18 18 21 7 19 21 17 22 15 20 19 10 18 25 23 25 23 21 23 19 22 23 15 23 23 24 20 23 24 17 21 19 23 22 14 19 21 23 16 23 19 19 22 26 22 24 24 11 21 21 22 22 19 18 19 27 14 15 20 22 26 20 13 26 19 20 18 20 21 26 25 20 21
Sample Range:
(leave blank to include all observations)
From:
To:
Chart options
R Code
x <-sort(x[!is.na(x)]) q1 <- function(data,n,p,i,f) { np <- n*p; i <<- floor(np) f <<- np - i qvalue <- (1-f)*data[i] + f*data[i+1] } q2 <- function(data,n,p,i,f) { np <- (n+1)*p i <<- floor(np) f <<- np - i qvalue <- (1-f)*data[i] + f*data[i+1] } q3 <- function(data,n,p,i,f) { np <- n*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i] } else { qvalue <- data[i+1] } } q4 <- function(data,n,p,i,f) { np <- n*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- (data[i]+data[i+1])/2 } else { qvalue <- data[i+1] } } q5 <- function(data,n,p,i,f) { np <- (n-1)*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i+1] } else { qvalue <- data[i+1] + f*(data[i+2]-data[i+1]) } } q6 <- function(data,n,p,i,f) { np <- n*p+0.5 i <<- floor(np) f <<- np - i qvalue <- data[i] } q7 <- function(data,n,p,i,f) { np <- (n+1)*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i] } else { qvalue <- f*data[i] + (1-f)*data[i+1] } } q8 <- function(data,n,p,i,f) { np <- (n+1)*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i] } else { if (f == 0.5) { qvalue <- (data[i]+data[i+1])/2 } else { if (f < 0.5) { qvalue <- data[i] } else { qvalue <- data[i+1] } } } } lx <- length(x) qval <- array(NA,dim=c(99,8)) mystep <- 25 mystart <- 25 if (lx>10){ mystep=10 mystart=10 } if (lx>20){ mystep=5 mystart=5 } if (lx>50){ mystep=2 mystart=2 } if (lx>=100){ mystep=1 mystart=1 } for (perc in seq(mystart,99,mystep)) { qval[perc,1] <- q1(x,lx,perc/100,i,f) qval[perc,2] <- q2(x,lx,perc/100,i,f) qval[perc,3] <- q3(x,lx,perc/100,i,f) qval[perc,4] <- q4(x,lx,perc/100,i,f) qval[perc,5] <- q5(x,lx,perc/100,i,f) qval[perc,6] <- q6(x,lx,perc/100,i,f) qval[perc,7] <- q7(x,lx,perc/100,i,f) qval[perc,8] <- q8(x,lx,perc/100,i,f) } bitmap(file='test1.png') myqqnorm <- qqnorm(x,col=2) qqline(x) grid() dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Percentiles - Ungrouped Data',9,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p',1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/method_1.htm', 'Weighted Average at Xnp',''),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/method_2.htm','Weighted Average at X(n+1)p',''),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/method_3.htm','Empirical Distribution Function',''),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/method_4.htm','Empirical Distribution Function - Averaging',''),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/method_5.htm','Empirical Distribution Function - Interpolation',''),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/method_6.htm','Closest Observation',''),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/method_7.htm','True Basic - Statistics Graphics Toolkit',''),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/method_8.htm','MS Excel (old versions)',''),1,TRUE) a<-table.row.end(a) for (perc in seq(mystart,99,mystep)) { a<-table.row.start(a) a<-table.element(a,round(perc/100,2),1,TRUE) for (j in 1:8) { a<-table.element(a,round(qval[perc,j],6)) } a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable.tab')
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