R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

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(35.323
+ ,35.478
+ ,4.39
+ ,41.667
+ ,22.173
+ ,28.021
+ ,18.109
+ ,13.962
+ ,40.174
+ ,16.065
+ ,18.145
+ ,18.439
+ ,10.603
+ ,34.811
+ ,69.064
+ ,51.202
+ ,14.786
+ ,33.01
+ ,81.101
+ ,89.232
+ ,21.223
+ ,15.173
+ ,241.66
+ ,26.848
+ ,8.752
+ ,60.535
+ ,60.535
+ ,26.052
+ ,49.218
+ ,30.669
+ ,18.673
+ ,86
+ ,10.632
+ ,35.802
+ ,33.974
+ ,36.972
+ ,4.928
+ ,53.976
+ ,15.467
+ ,35.723
+ ,40.424
+ ,9.706
+ ,26.532
+ ,23.843
+ ,18.062
+ ,35.681
+ ,68.125
+ ,23.937
+ ,31.479
+ ,66.659
+ ,250.234
+ ,49.469
+ ,42.951
+ ,43.402
+ ,24.112
+ ,56.95
+ ,17.313
+ ,25.658
+ ,48.172
+ ,13.891
+ ,32.048
+ ,19.797
+ ,31.317
+ ,20.966
+ ,22.708
+ ,26.81
+ ,52.004
+ ,32.354
+ ,27.128
+ ,26.529
+ ,28.392
+ ,57.393
+ ,194.731
+ ,9.415
+ ,91.076
+ ,57.751
+ ,8.236
+ ,20.407
+ ,13.681
+ ,79.659
+ ,53.48
+ ,6.906
+ ,50.202
+ ,37.877
+ ,85.903
+ ,35.351
+ ,283.801
+ ,5.974
+ ,3.441
+ ,51.987
+ ,13.22
+ ,1.455
+ ,18.187
+ ,21.29
+ ,5.686
+ ,4.944
+ ,32.789
+ ,50.494
+ ,35.162
+ ,38.095
+ ,19.172
+ ,24.5
+ ,20.573
+ ,42.042
+ ,302.912
+ ,25.027
+ ,16.488
+ ,32.36
+ ,6.193
+ ,37.7
+ ,6.343
+ ,23.025
+ ,48.578
+ ,21.564
+ ,33.697
+ ,10.831
+ ,19.172
+ ,21.075
+ ,33.189
+ ,60.5
+ ,33.686
+ ,40.838
+ ,13.491
+ ,106.637
+ ,35.897
+ ,7.314
+ ,49.094
+ ,14.667
+ ,54.179
+ ,145.846
+ ,18.56
+ ,23.525
+ ,21.804
+ ,26.301
+ ,41.33
+ ,10.5
+ ,13.338
+ ,60.31
+ ,34.256
+ ,48.267
+ ,41.559
+ ,32.45
+ ,10.951
+ ,22.561
+ ,57.095
+ ,19.105
+ ,13.151
+ ,27.426
+ ,15.355
+ ,13.82
+ ,47.21
+ ,110.349
+ ,34.985
+ ,27.257
+ ,23.556
+ ,50.108
+ ,18.158
+ ,87.357
+ ,18.187
+ ,28.33
+ ,13.474
+ ,26.244)
> #'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/wessaorg/rcomp/tmp/1haqz1322789289.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/2fa821322789289.tab") 
> 
> try(system("convert tmp/1haqz1322789289.ps tmp/1haqz1322789289.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.103   0.102   1.630