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