R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-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(1418
+ ,869
+ ,1530
+ ,2172
+ ,901
+ ,463
+ ,3201
+ ,371
+ ,1192
+ ,1583
+ ,1439
+ ,1764
+ ,1495
+ ,1373
+ ,2187
+ ,1491
+ ,4041
+ ,1706
+ ,2152
+ ,1036
+ ,1882
+ ,1929
+ ,2242
+ ,1220
+ ,1289
+ ,2515
+ ,2147
+ ,2352
+ ,1638
+ ,1222
+ ,1812
+ ,1677
+ ,1579
+ ,1731
+ ,807
+ ,2452
+ ,829
+ ,1940
+ ,2662
+ ,186
+ ,1499
+ ,865
+ ,1793
+ ,2527
+ ,2747
+ ,1324
+ ,2702
+ ,1383
+ ,1179
+ ,2099
+ ,4308
+ ,918
+ ,1831
+ ,3373
+ ,1713
+ ,1438
+ ,496
+ ,2253
+ ,744
+ ,1161
+ ,2352
+ ,2144
+ ,4691
+ ,1112
+ ,2694
+ ,1973
+ ,1769
+ ,3148
+ ,2474
+ ,2084
+ ,1954
+ ,1226
+ ,1389
+ ,1496
+ ,2269
+ ,1833
+ ,1268
+ ,1943
+ ,893
+ ,1762
+ ,1403
+ ,1425
+ ,1857
+ ,1840
+ ,1502
+ ,1441
+ ,1420
+ ,1416
+ ,2970
+ ,1317
+ ,1644
+ ,870
+ ,1654
+ ,1054
+ ,937
+ ,3004
+ ,2008
+ ,2547
+ ,1885
+ ,1626
+ ,1468
+ ,2445
+ ,1964
+ ,1381
+ ,1369
+ ,1659
+ ,2888
+ ,1290
+ ,2845
+ ,1982
+ ,1904
+ ,1391
+ ,602
+ ,1743
+ ,1559
+ ,2014
+ ,2143
+ ,2146
+ ,874
+ ,1590
+ ,1590
+ ,1210
+ ,2072
+ ,1281
+ ,1401
+ ,834
+ ,1105
+ ,1272
+ ,1944
+ ,391
+ ,761
+ ,1605
+ ,530
+ ,1988
+ ,1386
+ ,2395
+ ,387
+ ,1742
+ ,620
+ ,449
+ ,800
+ ,1684
+ ,1050
+ ,2699
+ ,1606
+ ,1502
+ ,1204
+ ,1138
+ ,568
+ ,1459
+ ,2158
+ ,1111
+ ,1421
+ ,2833
+ ,1955
+ ,2922
+ ,1002
+ ,1060
+ ,956
+ ,2186
+ ,3604
+ ,1035
+ ,1417
+ ,3261
+ ,1587
+ ,1424
+ ,1701
+ ,1249
+ ,946
+ ,1926
+ ,3352
+ ,1641
+ ,2035
+ ,2312
+ ,1369
+ ,1577
+ ,2201
+ ,961
+ ,1900
+ ,1254
+ ,1335
+ ,1597
+ ,207
+ ,1645
+ ,2429
+ ,151
+ ,474
+ ,141
+ ,1639
+ ,872
+ ,1318
+ ,1018
+ ,1383
+ ,1314
+ ,1335
+ ,1403
+ ,910
+ ,616
+ ,1407
+ ,771
+ ,766
+ ,473
+ ,1376
+ ,1232
+ ,1521
+ ,572
+ ,1059
+ ,1544
+ ,1230
+ ,1206
+ ,1205
+ ,1255
+ ,613
+ ,721
+ ,1109
+ ,740
+ ,1126
+ ,728
+ ,689
+ ,592
+ ,995
+ ,1613
+ ,2048
+ ,705
+ ,301
+ ,1803
+ ,799
+ ,861
+ ,1186
+ ,1451
+ ,628
+ ,1161
+ ,1463
+ ,742
+ ,979
+ ,675
+ ,1241
+ ,676
+ ,1049
+ ,620
+ ,1081
+ ,1688
+ ,736
+ ,617
+ ,812
+ ,1051
+ ,1656
+ ,705
+ ,945
+ ,554
+ ,1597
+ ,982
+ ,222
+ ,1212
+ ,1143
+ ,435
+ ,532
+ ,882
+ ,608
+ ,459
+ ,578
+ ,826
+ ,509
+ ,717
+ ,637
+ ,857
+ ,830
+ ,652
+ ,707
+ ,954
+ ,1461
+ ,672
+ ,778
+ ,1141
+ ,680
+ ,1090
+ ,616
+ ,285
+ ,1145
+ ,733
+ ,888
+ ,849
+ ,1182
+ ,528
+ ,642
+ ,947
+ ,819
+ ,757
+ ,894)
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa P., (2012), Percentiles (v1.0.6) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_percentiles.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> 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/1n2on1355832498.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/2i6xv1355832498.tab") 
> 
> try(system("convert tmp/1n2on1355832498.ps tmp/1n2on1355832498.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.565   0.156   1.698