Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
52 16 46 56 NA 55 50 NA 60 NA NA 67 52 55 37 54 72 51 48 NA 50 63 33 67 46 54 NA 61 33 47 69 52 NA NA NA 73 NA NA NA 51 NA 56 56 NA 66 66 73 NA NA NA NA 58 NA 61 NA 50 NA 54 NA NA 80 NA 56 56 56 53 47 NA 47 NA NA NA 51 NA 35 NA NA NA NA 53 46 67 59 NA NA 50 NA NA NA NA 47 63 NA NA 51 NA 55 38 56 NA 50 54 57 NA NA NA 49 37 NA 59 46 NA NA NA 53 48 NA NA NA 62 62 NA NA NA 67 50 54 58 NA 63 31 65 NA NA 57 NA 54 47 57 NA 41 NA 63 56 NA 50 NA 41 NA 56 NA NA 42 52 NA 44 62 NA 50 NA NA 66 62 NA 47 NA NA 60 NA NA 45 58 51 NA 30 46 51 56 NA 44 NA NA 42 NA 44 NA NA 46 NA 50 54 NA NA NA 55 NA 59 54 66 55 NA 51 NA 42 NA NA NA NA NA NA 63 43 65 NA NA NA NA NA 67 52 52 69 NA 46 NA 40 70 NA 77 NA NA NA NA NA NA 41 NA NA 51 69 60 45 NA 39 51 NA NA NA NA 51 NA NA NA NA NA NA 51 65 51 NA NA 58 NA 54 52 72 50 65 NA NA NA NA NA 50 NA NA 52 NA NA 59 NA 52 45 NA 70 NA 71 NA 58 39 46 NA 67 44 NA 41 68 63 57 NA 39 NA NA 38 NA 51 59 NA NA 47 50 57 21 NA 51 37 67 43 NA NA 40 NA 58 64 NA 58 NA 59 NA NA 58 41 56 63 NA NA 58 NA 47 NA 62 60 50 46 44 58 56 43 NA 54 NA NA 66 62 58 67 25 56 53 NA 59 46 49 NA 76 33 49 53 NA 72 51 NA NA 51 54 52 59 NA 67 NA 58 NA NA NA NA 53 NA NA NA NA NA NA NA NA NA NA NA NA
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