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(-1.02788E+13
+ ,-1.18964E+13
+ ,-62057298249
+ ,7.24304E+12
+ ,2.41619E+12
+ ,-4.05419E+12
+ ,-8.6353E+12
+ ,-6.7979E+12
+ ,8.19387E+12
+ ,1.32084E+13
+ ,3.01757E+12
+ ,3.47692E+12
+ ,1.08919E+13
+ ,1.12154E+13
+ ,4.60923E+12
+ ,3.63868E+12
+ ,3.63868E+12
+ ,6.41451E+12
+ ,5.01673E+12
+ ,1.02511E+13
+ ,5.53413E+12
+ ,7.44932E+12
+ ,6.317E+12
+ ,2.22116E+12
+ ,-6.24218E+12
+ ,-5.78283E+12
+ ,-8.37098E+12
+ ,-7.37449E+12
+ ,-3.81579E+12
+ ,5.25789E+11
+ ,-5.40126E+12
+ ,-1.93064E+13
+ ,-1.53402E+13
+ ,-1.37226E+13
+ ,-1.13221E+13
+ ,-1.4745E+13
+ ,-1.15616E+13
+ ,-1.28557E+13
+ ,-9.10933E+12
+ ,-1.0054E+13
+ ,-1.18074E+13
+ ,-1.36966E+13
+ ,-1.49067E+13
+ ,-1.32311E+13
+ ,-1.35546E+13
+ ,-7.38183E+12
+ ,-4.60599E+12
+ ,-1.83016E+12
+ ,-2.82664E+12
+ ,-1.18312E+12
+ ,1.45689E+12
+ ,3.28809E+12
+ ,9.97208E+12
+ ,1.24244E+13
+ ,1.16996E+13
+ ,6.87892E+12
+ ,7.66179E+12
+ ,7.17651E+12
+ ,3.80549E+12
+ ,3.96725E+12
+ ,-4.79368E+12
+ ,-2.82664E+12
+ ,-3.9071E+12
+ ,-2.72293E+12
+ ,8.87629E+11
+ ,4.52412E+12
+ ,1.23707E+12
+ ,-1.45478E+12
+ ,-2.93654E+12
+ ,-9.51682E+12
+ ,-7.22627E+12
+ ,-9.30321E+12
+ ,-6.41747E+12
+ ,-4.06446E+11
+ ,8.67801E+12
+ ,1.25862E+13
+ ,8.08282E+12
+ ,6.16764E+12
+ ,1.04939E+12
+ ,-9.95432E+11
+ ,2.05207E+12
+ ,7.90133E+12
+ ,6.63319E+12
+ ,1.10007E+13
+ ,2.21362E+13
+ ,1.69599E+13
+ ,1.24306E+13
+ ,5.17735E+12
+ ,8.61702E+11
+ ,5.90217E+12
+ ,1.36148E+13
+ ,1.91208E+13
+ ,1.1628E+13
+ ,1.23269E+13
+ ,1.8875E+13
+ ,1.87133E+13
+ ,1.2217E+13
+ ,1.17317E+13
+ ,8.36068E+12
+ ,1.06772E+13
+ ,1.73612E+13
+ ,2.03247E+13
+ ,2.17805E+13
+ ,1.69018E+13
+ ,1.50644E+13
+ ,1.4631E+13
+ ,1.22564E+13
+ ,5.62431E+12
+ ,7.37773E+12
+ ,6.13551E+12
+ ,3.81283E+12
+ ,1.00634E+13
+ ,1.01412E+13
+ ,2.61627E+12
+ ,-1.39559E+12
+ ,-5.54329E+12
+ ,-1.26348E+13
+ ,-5.65319E+12
+ ,-7.0831E+12
+ ,3.45719E+12
+ ,8.14821E+12
+ ,7.33941E+12
+ ,9.2546E+12
+ ,-4.08745E+12
+ ,-9.07606E+12
+ ,-3.66022E+12
+ ,-4.76891E+11
+ ,-2.14634E+12
+ ,1.79394E+12
+ ,-3.69968E+12
+ ,5.31986E+11
+ ,-2.7681E+11
+ ,-9.90459E+12
+ ,-4.86412E+12
+ ,-5.51116E+12
+ ,7.72718E+12
+ ,1.61646E+13
+ ,2.54627E+13
+ ,1.69734E+13
+ ,1.72969E+13
+ ,1.5278E+13
+ ,8.73606E+12
+ ,9.79059E+12
+ ,4.75012E+12
+ ,8.01124E+12
+ ,7.83921E+11
+ ,1.13337E+12
+ ,4.98966E+12
+ ,1.278E+13
+ ,4.74393E+12
+ ,4.8279E+12
+ ,6.48088E+11
+ ,-1.86636E+11
+ ,1.35318E+12
+ ,-3.17609E+12
+ ,-2.5809E+12
+ ,-8.27476E+11
+ ,5.85651E+12
+ ,1.13242E+13
+ ,3.25597E+12
+ ,-5.88653E+12
+ ,-8.97921E+11
+ ,-1.0054E+13
+ ,-1.03194E+13
+ ,-7.27192E+12
+ ,-1.80061E+13
+ ,-1.16456E+13
+ ,-7.78933E+12
+ ,-8.9735E+12
+ ,-1.3988E+13
+ ,-2.2076E+13
+ ,-1.50426E+13
+ ,-1.09986E+13
+ ,-1.41301E+13
+ ,-2.11894E+13
+ ,-2.07041E+13
+ ,-2.21859E+13
+ ,-1.16975E+13
+ ,-1.33151E+13
+ ,-1.98373E+13
+ ,-1.69775E+13
+ ,-1.36065E+13
+ ,-2.44507E+12
+ ,-6.81257E+12
+ ,-8.51414E+12
+ ,6.99943E+11
+ ,7.7593E+12
+ ,9.29292E+12
+ ,1.52521E+13
+ ,1.40679E+13
+ ,8.86569E+12
+ ,5.76634E+12
+ ,4.76985E+12
+ ,3.66347E+12
+ ,4.2006E+12
+ ,4.17467E+12
+ ,1.36086E+13
+ ,1.61387E+13
+ ,1.28257E+13
+ ,1.75426E+13
+ ,1.258E+13
+ ,1.44449E+12
+ ,-9.2057E+12
+ ,-1.3651E+13
+ ,-2.1655E+13
+ ,-2.45666E+13
+ ,-3.40067E+13
+ ,-2.81056E+13
+ ,-4.32332E+13
+ ,-5.47959E+13
+ ,-6.39779E+13
+ ,-5.62394E+13
+ ,-5.46798E+13
+ ,-3.73259E+13
+ ,-9.74903E+12
+ ,-4.19116E+12
+ ,-5.2333E+12
+ ,-2.96247E+12
+ ,1.68821E+13
+ ,1.85515E+13
+ ,3.62611E+13
+ ,4.73965E+13
+ ,5.34138E+13
+ ,3.81825E+13)
> par2 = '36'
> par1 = '0'
> #'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!)
> par1 <- as.numeric(par1)
> par2 <- as.numeric(par2)
> x <- as.ts(x)
> library(lattice)
> postscript(file="/var/www/html/rcomp/tmp/1xgql1258881421.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/2r9vr1258881421.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/3s9fi1258881421.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/4rk231258881421.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> qqnorm(x)
> qqline(x)
> grid()
> dev.off()
null device 
          1 
> if (par2 > 0)
+ {
+ postscript(file="/var/www/html/rcomp/tmp/508sz1258881421.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='Lag plot (k=1), lowess, and regression line')
+ lines(lowess(z))
+ abline(lm(z))
+ dev.off()
+ if (par2 > 1) {
+ postscript(file="/var/www/html/rcomp/tmp/6che91258881421.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
+ dum <- cbind(lag(x,k=par2),x)
+ dum
+ dum1 <- dum[(par2+1):length(x),]
+ dum1
+ z <- as.data.frame(dum1)
+ z
+ mylagtitle <- 'Lag plot (k='
+ mylagtitle <- paste(mylagtitle,par2,sep='')
+ mylagtitle <- paste(mylagtitle,'), and lowess',sep='')
+ plot(z,main=mylagtitle)
+ lines(lowess(z))
+ dev.off()
+ }
+ postscript(file="/var/www/html/rcomp/tmp/72uee1258881421.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,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. 
-6.398e+13 -7.484e+12  1.185e+12  1.289e+06  8.440e+12  5.341e+13 
> 
> #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/80di51258881421.tab") 
> 
> system("convert tmp/1xgql1258881421.ps tmp/1xgql1258881421.png")
> system("convert tmp/2r9vr1258881421.ps tmp/2r9vr1258881421.png")
> system("convert tmp/3s9fi1258881421.ps tmp/3s9fi1258881421.png")
> system("convert tmp/4rk231258881421.ps tmp/4rk231258881421.png")
> system("convert tmp/508sz1258881421.ps tmp/508sz1258881421.png")
> system("convert tmp/6che91258881421.ps tmp/6che91258881421.png")
> system("convert tmp/72uee1258881421.ps tmp/72uee1258881421.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.547   1.094   2.518