Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
1894.00 1757.00 3582.00 5321.00 5561.00 5907.00 4944.00 4966.00 3258.00 1964.00 1743.00 1262.00 2086.00 1793.00 3548.00 5672.00 6084.00 4914.00 4990.00 5139.00 3218.00 2179.00 2238.00 1442.00 2205.00 2025.00 3531.00 4977.00 7998.00 4880.00 5231.00 5202.00 3303.00 2683.00 2202.00 1376.00 2422.00 1997.00 3163.00 5964.00 5657.00 6415.00 6208.00 4500.00 2939.00 2702.00 2090.00 1504.00 2549.00 1931.00 3013.00 6204.00 5788.00 5611.00 5594.00 4647.00 3490.00 2487.00 1992.00 1507.00 2306.00 2002.00 3075.00 5331.00 5589.00 5813.00 4876.00 4665.00 3601.00 2192.00 2111.00 1580.00 2288.00 1993.00 3228.00 5000.00 5480.00 5770.00 4962.00 4685.00 3607.00 2222.00 2467.00 1594.00 2228.00 1910.00 3157.00 4809.00 6249.00 4607.00 4975.00 4784.00 3028.00 2461.00 2218.00 1351.00 2070.00 1887.00 3024.00 4596.00 6398.00 4459.00 5382.00 4359.00 2687.00 2249.00 2154.00 1169.00 2429.00 1762.00 2846.00 5627.00 5749.00 4502.00 5720.00 4403.00 2867.00 2635.00 2059.00 1511.00 2359.00 1741.00 2917.00 6249.00 5760.00 6250.00 5134.00 4831.00 3695.00 2462.00 2146.00 1579.00
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