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(1384.2,1368.9,-275.1,-408.9,-37.5,171.5,671.8,-18.5,231.6,747.5,1505.7,-83.6,1173.2,1452.1,777,-52.8,861.2,735.2,1073.6,966.9,1189.8,1093.5,1782.7,-70.4,1471.6,1273.8,900.8,-910.2,299.8,460.2,677.2,937.1,1265.4,1275.6,1582.6,-154.2,1667.6,1083.1,891.7,-26.5,423.4,662.8,711.4,993.3,1133.2,343.9,1415.8,-531.8,1193.6,1201.3,805.6,-164.8,327.3,223.7,675.8,949.7,704.4,265.6,1206,-558.2,1066.8,977.8,207.1,-980.7,-586.4,-24.3,-417.5,104.7,749.5,842.3,1176,-730.3,911.6,662.1,539.1,-236,286.9,497.4,912,519.4,260.1,945.2,412.7,-54.2,592.8,179.9,-548.6,-1685.8,-2041,-1048.7,-1708.4,-1550.7,-1650.2,-911.3)
> par2 = '0'
> 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/1czat1229888757.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,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/2v7681229888757.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> hist(x)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/3jcv61229888757.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,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/4adqs1229888757.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> qqnorm(x)
> grid()
> dev.off()
null device 
          1 
> if (par2 > 0)
+ {
+ postscript(file="/var/www/html/rcomp/tmp/5qc101229888757.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,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/6sk3n1229888757.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
+ acf(x,lag.max=par2,main='Autocorrelation Function')
+ grid()
+ dev.off()
+ }
> summary(x)
    Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
-2041.00   -48.97   627.40   410.50   989.40  1783.00 
> 
> #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/7m97v1229888757.tab") 
> 
> system("convert tmp/1czat1229888757.ps tmp/1czat1229888757.png")
> system("convert tmp/2v7681229888757.ps tmp/2v7681229888757.png")
> system("convert tmp/3jcv61229888757.ps tmp/3jcv61229888757.png")
> system("convert tmp/4adqs1229888757.ps tmp/4adqs1229888757.png")
> system("convert tmp/5qc101229888757.ps tmp/5qc101229888757.png")
convert: unable to open image `tmp/5qc101229888757.ps': No such file or directory.
convert: missing an image filename `tmp/5qc101229888757.png'.
> system("convert tmp/6sk3n1229888757.ps tmp/6sk3n1229888757.png")
convert: unable to open image `tmp/6sk3n1229888757.ps': No such file or directory.
convert: missing an image filename `tmp/6sk3n1229888757.png'.
> 
> 
> proc.time()
   user  system elapsed 
  1.060   0.644   1.670