R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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(2560
+ ,2491
+ ,2380
+ ,2291
+ ,2079
+ ,1929
+ ,1851
+ ,1607
+ ,1661
+ ,2259
+ ,1668
+ ,2011
+ ,1944
+ ,1958
+ ,1844
+ ,1868
+ ,1701
+ ,2338
+ ,2018
+ ,1302
+ ,2168
+ ,2139
+ ,1560
+ ,2093
+ ,1973
+ ,2090
+ ,2811
+ ,1984
+ ,1849
+ ,2433
+ ,2071
+ ,1855
+ ,1756
+ ,1898
+ ,1770
+ ,1969
+ ,1769
+ ,2139
+ ,3013
+ ,2061
+ ,2132
+ ,2973
+ ,2081
+ ,2257
+ ,2075
+ ,2084
+ ,1747
+ ,2092
+ ,1919
+ ,2551
+ ,2643
+ ,2153
+ ,2496
+ ,2645
+ ,2035
+ ,2294
+ ,2205
+ ,2044
+ ,1762
+ ,1897
+ ,1821
+ ,1905
+ ,2111
+ ,1643
+ ,1956
+ ,1977
+ ,1685
+ ,1393
+ ,1574
+ ,1793
+ ,1562
+ ,1510
+ ,1675
+ ,1965
+ ,2173
+ ,2395
+ ,2197
+ ,2257
+ ,2885
+ ,1594
+ ,1950
+ ,1772
+ ,1280
+ ,1724
+ ,1473
+ ,1461
+ ,1576
+ ,1900
+ ,1618
+ ,2303
+ ,1994
+ ,1575
+ ,1893
+ ,1788
+ ,1817
+ ,3233
+ ,727
+ ,1121
+ ,1665
+ ,1401
+ ,1415
+ ,2058
+ ,1544
+ ,1379
+ ,1402
+ ,1313
+ ,1296
+ ,1398
+ ,1288
+ ,1563
+ ,1972
+ ,1496
+ ,1481
+ ,1819
+ ,1479
+ ,1635
+ ,1511
+ ,1547
+ ,1388
+ ,1958
+ ,1390
+ ,1597
+ ,1842
+ ,1396
+ ,1671
+ ,1385
+ ,1632
+ ,1313
+ ,1300
+ ,1431
+ ,1398
+ ,1198
+ ,1292
+ ,1434
+ ,1660
+ ,1837
+ ,1455
+ ,1315
+ ,1642
+ ,1069
+ ,1209
+ ,1586
+ ,1122
+ ,1063
+ ,1125
+ ,1414
+ ,1347
+ ,1403
+ ,1299
+ ,1547
+ ,1515
+ ,1247
+ ,1639
+ ,1296
+ ,1063
+ ,1282
+ ,1365
+ ,1268
+ ,1532
+ ,1455
+ ,1393
+ ,1515
+ ,1510
+ ,1225
+ ,1577
+ ,1417
+ ,1224
+ ,1693
+ ,1633
+ ,1639
+ ,1914
+ ,1586
+ ,1552
+ ,2081
+ ,1500
+ ,1437
+ ,1470
+ ,1849
+ ,1387
+ ,1592
+ ,1589
+ ,1798
+ ,1935
+ ,1887
+ ,2027
+ ,2080
+ ,1556
+ ,1682
+ ,1785
+ ,1869
+ ,1781
+ ,2082
+ ,2570
+ ,1862
+ ,1936
+ ,1504
+ ,1765
+ ,1607
+ ,1577
+ ,1493
+ ,1615
+ ,1700
+ ,1335
+ ,1523
+ ,1623
+ ,1540
+ ,1637
+ ,1524
+ ,1419
+ ,1821
+ ,1593
+ ,1357
+ ,1263
+ ,1750
+ ,1405
+ ,1393
+ ,1639
+ ,1679
+ ,1551
+ ,1744
+ ,1429
+ ,1784)
> par2 = '12'
> par1 = '0'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: Wessa P., (2007), Univariate Explorative Data Analysis (v1.0.5) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_edauni.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description
> par1 <- as.numeric(par1)
> par2 <- as.numeric(par2)
> x <- as.ts(x)
> library(lattice)
> postscript(file="/var/www/html/rcomp/tmp/13r8x1229179744.ps",horizontal=F,pagecentre=F,paper="special",width=13.888888888889,height=5.5555555555556) 
> plot(x,type='l',main='Run Sequence Plot',xlab='time or index',ylab='value')
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/2bas11229179744.ps",horizontal=F,pagecentre=F,paper="special",width=13.888888888889,height=5.5555555555556) 
> hist(x)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/3au591229179744.ps",horizontal=F,pagecentre=F,paper="special",width=13.888888888889,height=5.5555555555556) 
> if (par1 > 0)
+ {
+ densityplot(~x,col='black',main=paste('Density Plot   bw = ',par1),bw=par1)
+ } else {
+ densityplot(~x,col='black',main='Density Plot')
+ }
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/481hv1229179744.ps",horizontal=F,pagecentre=F,paper="special",width=13.888888888889,height=5.5555555555556) 
> qqnorm(x)
> grid()
> dev.off()
null device 
          1 
> if (par2 > 0)
+ {
+ postscript(file="/var/www/html/rcomp/tmp/50flv1229179744.ps",horizontal=F,pagecentre=F,paper="special",width=13.888888888889,height=5.5555555555556) 
+ dum <- cbind(lag(x,k=1),x)
+ dum
+ dum1 <- dum[2:length(x),]
+ dum1
+ z <- as.data.frame(dum1)
+ z
+ plot(z,main=paste('Lag plot, lowess, and regression line'))
+ lines(lowess(z))
+ abline(lm(z))
+ dev.off()
+ postscript(file="/var/www/html/rcomp/tmp/67usf1229179744.ps",horizontal=F,pagecentre=F,paper="special",width=13.888888888889,height=5.5555555555556) 
+ acf(x,lag.max=par2,main='Autocorrelation Function')
+ grid()
+ dev.off()
+ }
null device 
          1 
> summary(x)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
    727    1463    1666    1741    1963    3233 
> 
> #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,'Descriptive Statistics',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'# observations',header=TRUE)
> a<-table.element(a,length(x))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'minimum',header=TRUE)
> a<-table.element(a,min(x))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Q1',header=TRUE)
> a<-table.element(a,quantile(x,0.25))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'median',header=TRUE)
> a<-table.element(a,median(x))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'mean',header=TRUE)
> a<-table.element(a,mean(x))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Q3',header=TRUE)
> a<-table.element(a,quantile(x,0.75))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'maximum',header=TRUE)
> a<-table.element(a,max(x))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/71ihz1229179744.tab") 
> 
> system("convert tmp/13r8x1229179744.ps tmp/13r8x1229179744.png")
> system("convert tmp/2bas11229179744.ps tmp/2bas11229179744.png")
> system("convert tmp/3au591229179744.ps tmp/3au591229179744.png")
> system("convert tmp/481hv1229179744.ps tmp/481hv1229179744.png")
> system("convert tmp/50flv1229179744.ps tmp/50flv1229179744.png")
> system("convert tmp/67usf1229179744.ps tmp/67usf1229179744.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.601   1.133   1.876