R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> x <- c(24
+ ,25
+ ,17
+ ,18
+ ,18
+ ,16
+ ,20
+ ,16
+ ,18
+ ,17
+ ,23
+ ,30
+ ,23
+ ,18
+ ,15
+ ,12
+ ,21
+ ,15
+ ,20
+ ,31
+ ,27
+ ,34
+ ,21
+ ,31
+ ,19
+ ,16
+ ,20
+ ,21
+ ,22
+ ,17
+ ,24
+ ,25
+ ,26
+ ,25
+ ,17
+ ,32
+ ,33
+ ,13
+ ,32
+ ,25
+ ,29
+ ,22
+ ,18
+ ,17
+ ,20
+ ,15
+ ,20
+ ,33
+ ,29
+ ,23
+ ,26
+ ,18
+ ,20
+ ,11
+ ,28
+ ,26
+ ,22
+ ,17
+ ,12
+ ,14
+ ,17
+ ,21
+ ,19
+ ,18
+ ,10
+ ,29
+ ,31
+ ,19
+ ,9
+ ,20
+ ,28
+ ,19
+ ,30
+ ,29
+ ,26
+ ,23
+ ,13
+ ,21
+ ,19
+ ,28
+ ,23
+ ,18
+ ,21
+ ,20
+ ,23
+ ,21
+ ,21
+ ,15
+ ,28
+ ,19
+ ,26
+ ,10
+ ,16
+ ,22
+ ,19
+ ,31
+ ,31
+ ,29
+ ,19
+ ,22
+ ,23
+ ,15
+ ,20
+ ,18
+ ,23
+ ,25
+ ,21
+ ,24
+ ,25
+ ,17
+ ,13
+ ,28
+ ,21
+ ,25
+ ,9
+ ,16
+ ,19
+ ,17
+ ,25
+ ,20
+ ,29
+ ,14
+ ,22
+ ,15
+ ,19
+ ,20
+ ,15
+ ,20
+ ,18
+ ,33
+ ,22
+ ,16
+ ,17
+ ,16
+ ,21
+ ,26
+ ,18
+ ,18
+ ,17
+ ,22
+ ,30
+ ,30
+ ,24
+ ,21
+ ,21
+ ,29
+ ,31
+ ,20
+ ,16
+ ,22
+ ,20
+ ,28
+ ,38
+ ,22
+ ,20
+ ,17
+ ,28
+ ,22
+ ,31
+ ,14
+ ,11
+ ,6
+ ,12
+ ,8
+ ,10
+ ,10
+ ,11
+ ,16
+ ,11
+ ,13
+ ,12
+ ,8
+ ,12
+ ,11
+ ,4
+ ,9
+ ,8
+ ,8
+ ,14
+ ,15
+ ,16
+ ,9
+ ,14
+ ,11
+ ,8
+ ,9
+ ,9
+ ,9
+ ,9
+ ,10
+ ,16
+ ,11
+ ,8
+ ,9
+ ,16
+ ,11
+ ,16
+ ,12
+ ,12
+ ,14
+ ,9
+ ,10
+ ,9
+ ,10
+ ,12
+ ,14
+ ,14
+ ,10
+ ,14
+ ,16
+ ,9
+ ,10
+ ,6
+ ,8
+ ,13
+ ,10
+ ,8
+ ,7
+ ,15
+ ,9
+ ,10
+ ,12
+ ,13
+ ,10
+ ,11
+ ,8
+ ,9
+ ,13
+ ,11
+ ,8
+ ,9
+ ,9
+ ,15
+ ,9
+ ,10
+ ,14
+ ,12
+ ,12
+ ,11
+ ,14
+ ,6
+ ,12
+ ,8
+ ,14
+ ,11
+ ,10
+ ,14
+ ,12
+ ,10
+ ,14
+ ,5
+ ,11
+ ,10
+ ,9
+ ,10
+ ,16
+ ,13
+ ,9
+ ,10
+ ,10
+ ,7
+ ,9
+ ,8
+ ,14
+ ,14
+ ,8
+ ,9
+ ,14
+ ,14
+ ,8
+ ,8
+ ,8
+ ,7
+ ,6
+ ,8
+ ,6
+ ,11
+ ,14
+ ,11
+ ,11
+ ,11
+ ,14
+ ,8
+ ,20
+ ,11
+ ,8
+ ,11
+ ,10
+ ,14
+ ,11
+ ,9
+ ,9
+ ,8
+ ,10
+ ,13
+ ,13
+ ,12
+ ,8
+ ,13
+ ,14
+ ,12
+ ,14
+ ,15
+ ,13
+ ,16
+ ,9
+ ,9
+ ,9
+ ,8
+ ,7
+ ,16
+ ,11
+ ,9
+ ,11
+ ,9
+ ,14
+ ,13
+ ,16
+ ,11
+ ,7
+ ,17
+ ,10
+ ,12
+ ,12
+ ,11
+ ,11
+ ,12
+ ,13
+ ,14
+ ,16
+ ,11
+ ,10
+ ,11
+ ,15
+ ,9
+ ,11
+ ,17
+ ,17
+ ,11
+ ,18
+ ,14
+ ,10
+ ,11
+ ,15
+ ,15
+ ,13
+ ,16
+ ,13
+ ,9
+ ,18
+ ,18
+ ,12
+ ,17
+ ,9
+ ,9
+ ,12
+ ,18
+ ,12
+ ,18
+ ,14
+ ,15
+ ,16
+ ,10
+ ,11
+ ,14
+ ,9
+ ,12
+ ,17
+ ,5
+ ,12
+ ,12
+ ,6
+ ,24
+ ,12
+ ,12
+ ,14
+ ,7
+ ,13
+ ,12
+ ,13
+ ,14
+ ,8
+ ,11
+ ,9
+ ,11
+ ,13
+ ,10
+ ,11
+ ,12
+ ,9
+ ,15
+ ,18
+ ,15
+ ,12
+ ,13
+ ,14
+ ,10
+ ,13
+ ,13
+ ,11
+ ,13
+ ,16
+ ,8
+ ,16
+ ,11
+ ,9
+ ,16
+ ,12
+ ,14
+ ,8
+ ,9
+ ,15
+ ,11
+ ,21
+ ,14
+ ,18
+ ,12
+ ,13
+ ,15
+ ,12
+ ,19
+ ,15
+ ,11
+ ,11
+ ,10
+ ,13
+ ,15
+ ,12
+ ,12
+ ,16
+ ,9
+ ,18
+ ,8
+ ,13
+ ,17
+ ,9
+ ,15
+ ,8
+ ,7
+ ,12
+ ,14
+ ,6
+ ,8
+ ,17
+ ,10
+ ,11
+ ,14
+ ,11
+ ,13
+ ,12
+ ,11
+ ,9
+ ,12
+ ,20
+ ,12
+ ,13
+ ,12
+ ,12
+ ,9
+ ,15
+ ,24
+ ,7
+ ,17
+ ,11
+ ,17
+ ,11
+ ,12
+ ,14
+ ,11
+ ,16
+ ,21
+ ,14
+ ,20
+ ,13
+ ,11
+ ,15
+ ,19
+ ,12
+ ,8
+ ,8
+ ,8
+ ,9
+ ,7
+ ,4
+ ,11
+ ,7
+ ,7
+ ,12
+ ,10
+ ,10
+ ,8
+ ,8
+ ,4
+ ,9
+ ,8
+ ,7
+ ,11
+ ,9
+ ,11
+ ,13
+ ,8
+ ,8
+ ,9
+ ,6
+ ,9
+ ,9
+ ,6
+ ,6
+ ,16
+ ,5
+ ,7
+ ,9
+ ,6
+ ,6
+ ,5
+ ,12
+ ,7
+ ,10
+ ,9
+ ,8
+ ,5
+ ,8
+ ,8
+ ,10
+ ,6
+ ,8
+ ,7
+ ,4
+ ,8
+ ,8
+ ,4
+ ,20
+ ,8
+ ,8
+ ,6
+ ,4
+ ,8
+ ,9
+ ,6
+ ,7
+ ,9
+ ,5
+ ,5
+ ,8
+ ,8
+ ,6
+ ,8
+ ,7
+ ,7
+ ,9
+ ,11
+ ,6
+ ,8
+ ,6
+ ,9
+ ,8
+ ,6
+ ,10
+ ,8
+ ,8
+ ,10
+ ,5
+ ,7
+ ,5
+ ,8
+ ,14
+ ,7
+ ,8
+ ,6
+ ,5
+ ,6
+ ,10
+ ,12
+ ,9
+ ,12
+ ,7
+ ,8
+ ,10
+ ,6
+ ,10
+ ,10
+ ,10
+ ,5
+ ,7
+ ,10
+ ,11
+ ,6
+ ,7
+ ,12
+ ,11
+ ,11
+ ,11
+ ,5
+ ,8
+ ,6
+ ,9
+ ,4
+ ,4
+ ,7
+ ,11
+ ,6
+ ,7
+ ,8
+ ,4
+ ,8
+ ,9
+ ,8
+ ,11
+ ,8
+ ,5
+ ,4
+ ,8
+ ,10
+ ,6
+ ,9
+ ,9
+ ,13
+ ,9
+ ,10
+ ,20
+ ,5
+ ,11
+ ,6
+ ,9
+ ,7
+ ,9
+ ,10
+ ,9
+ ,8
+ ,7
+ ,6
+ ,13
+ ,6
+ ,8
+ ,10
+ ,16
+ ,24
+ ,25
+ ,30
+ ,19
+ ,22
+ ,22
+ ,25
+ ,23
+ ,17
+ ,21
+ ,19
+ ,19
+ ,15
+ ,16
+ ,23
+ ,27
+ ,22
+ ,14
+ ,22
+ ,23
+ ,23
+ ,21
+ ,19
+ ,18
+ ,20
+ ,23
+ ,25
+ ,19
+ ,24
+ ,22
+ ,25
+ ,26
+ ,29
+ ,32
+ ,25
+ ,29
+ ,28
+ ,17
+ ,28
+ ,29
+ ,26
+ ,25
+ ,14
+ ,25
+ ,26
+ ,20
+ ,18
+ ,32
+ ,25
+ ,25
+ ,23
+ ,21
+ ,20
+ ,15
+ ,30
+ ,24
+ ,26
+ ,24
+ ,22
+ ,14
+ ,24
+ ,24
+ ,24
+ ,24
+ ,19
+ ,31
+ ,22
+ ,27
+ ,19
+ ,25
+ ,20
+ ,21
+ ,27
+ ,23
+ ,25
+ ,20
+ ,21
+ ,22
+ ,23
+ ,25
+ ,25
+ ,17
+ ,19
+ ,25
+ ,19
+ ,20
+ ,26
+ ,23
+ ,27
+ ,17
+ ,17
+ ,19
+ ,17
+ ,22
+ ,21
+ ,32
+ ,21
+ ,21
+ ,18
+ ,18
+ ,23
+ ,19
+ ,20
+ ,21
+ ,20
+ ,17
+ ,18
+ ,19
+ ,22
+ ,15
+ ,14
+ ,18
+ ,24
+ ,35
+ ,29
+ ,21
+ ,25
+ ,20
+ ,22
+ ,13
+ ,26
+ ,17
+ ,25
+ ,20
+ ,19
+ ,21
+ ,22
+ ,24
+ ,21
+ ,26
+ ,24
+ ,16
+ ,23
+ ,18
+ ,16
+ ,26
+ ,19
+ ,21
+ ,21
+ ,22
+ ,23
+ ,29
+ ,21
+ ,21
+ ,23
+ ,27
+ ,25
+ ,21
+ ,10
+ ,20
+ ,26
+ ,24
+ ,29
+ ,19
+ ,24
+ ,19
+ ,24
+ ,22
+ ,17
+ ,26
+ ,23
+ ,25
+ ,23
+ ,19
+ ,29
+ ,25
+ ,21
+ ,22
+ ,25
+ ,24
+ ,18
+ ,22
+ ,15
+ ,22
+ ,28
+ ,20
+ ,12
+ ,24
+ ,20
+ ,21
+ ,20
+ ,21
+ ,23
+ ,28
+ ,24
+ ,24
+ ,24
+ ,23
+ ,23
+ ,29
+ ,24
+ ,18
+ ,25
+ ,21
+ ,26
+ ,22
+ ,22
+ ,22
+ ,23
+ ,30
+ ,23
+ ,17
+ ,23
+ ,23
+ ,25
+ ,24
+ ,24
+ ,23
+ ,21
+ ,24
+ ,24
+ ,28
+ ,16
+ ,20
+ ,29
+ ,27
+ ,22
+ ,28
+ ,16
+ ,25
+ ,24
+ ,28
+ ,24
+ ,23
+ ,30
+ ,24
+ ,21
+ ,25
+ ,25
+ ,22
+ ,23
+ ,26
+ ,23
+ ,25
+ ,21
+ ,25
+ ,24
+ ,29
+ ,22
+ ,27
+ ,26
+ ,22
+ ,24
+ ,27
+ ,24
+ ,24
+ ,29
+ ,22
+ ,21
+ ,24
+ ,24
+ ,23
+ ,20
+ ,27
+ ,26
+ ,25
+ ,21
+ ,21
+ ,19
+ ,21
+ ,21
+ ,16
+ ,22
+ ,29
+ ,15
+ ,17
+ ,15
+ ,21
+ ,21
+ ,19
+ ,24
+ ,20
+ ,17
+ ,23
+ ,24
+ ,14
+ ,19
+ ,24
+ ,13
+ ,22
+ ,16
+ ,19
+ ,25
+ ,25
+ ,23
+ ,24
+ ,26
+ ,26
+ ,25
+ ,18
+ ,21
+ ,26
+ ,23
+ ,23
+ ,22
+ ,20
+ ,13
+ ,24
+ ,15
+ ,14
+ ,22
+ ,10
+ ,24
+ ,22
+ ,24
+ ,19
+ ,20
+ ,13
+ ,20
+ ,22
+ ,24
+ ,29
+ ,12
+ ,20
+ ,21
+ ,24
+ ,22
+ ,20
+ ,14
+ ,18
+ ,11
+ ,12
+ ,16
+ ,18
+ ,14
+ ,14
+ ,15
+ ,15
+ ,17
+ ,19
+ ,10
+ ,16
+ ,18
+ ,14
+ ,14
+ ,17
+ ,14
+ ,16
+ ,18
+ ,11
+ ,14
+ ,12
+ ,17
+ ,9
+ ,16
+ ,14
+ ,15
+ ,11
+ ,16
+ ,13
+ ,17
+ ,15
+ ,14
+ ,16
+ ,9
+ ,15
+ ,17
+ ,13
+ ,15
+ ,16
+ ,16
+ ,12
+ ,12
+ ,11
+ ,15
+ ,15
+ ,17
+ ,13
+ ,16
+ ,14
+ ,11
+ ,12
+ ,12
+ ,15
+ ,16
+ ,15
+ ,12
+ ,12
+ ,8
+ ,13
+ ,11
+ ,14
+ ,15
+ ,10
+ ,11
+ ,12
+ ,15
+ ,15
+ ,14
+ ,16
+ ,15
+ ,15
+ ,13
+ ,12
+ ,17
+ ,13
+ ,15
+ ,13
+ ,15
+ ,16
+ ,15
+ ,16
+ ,15
+ ,14
+ ,15
+ ,14
+ ,13
+ ,7
+ ,17
+ ,13
+ ,15
+ ,14
+ ,13
+ ,16
+ ,12
+ ,14
+ ,17
+ ,15
+ ,17
+ ,12
+ ,16
+ ,11
+ ,15
+ ,9
+ ,16
+ ,15
+ ,10
+ ,10
+ ,15
+ ,11
+ ,13
+ ,14
+ ,18
+ ,16
+ ,14
+ ,14
+ ,14
+ ,14
+ ,12
+ ,14
+ ,15
+ ,15
+ ,15
+ ,13
+ ,17
+ ,17
+ ,19
+ ,15
+ ,13
+ ,9
+ ,15
+ ,15
+ ,15
+ ,16
+ ,11
+ ,14
+ ,11
+ ,15
+ ,13
+ ,15
+ ,16
+ ,14
+ ,15
+ ,16
+ ,16
+ ,11
+ ,12
+ ,9
+ ,16
+ ,13
+ ,16
+ ,12
+ ,9
+ ,13
+ ,13
+ ,14
+ ,19
+ ,13
+ ,12
+ ,13)
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> 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)
+ }
> postscript(file="/var/www/html/rcomp/tmp/10qi11291552986.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> myqqnorm <- qqnorm(x,col=2)
> qqline(x)
> grid()
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/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="/var/www/html/rcomp/tmp/24qg71291552986.tab") 
> 
> try(system("convert tmp/10qi11291552986.ps tmp/10qi11291552986.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.158   0.312   4.259