R version 2.12.0 (2010-10-15)
Copyright (C) 2010 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(160984
+ ,147727
+ ,165384
+ ,151202
+ ,135237
+ ,143684
+ ,150327
+ ,164909
+ ,154582
+ ,162221
+ ,165019
+ ,169673
+ ,150297
+ ,146944
+ ,156522
+ ,148536
+ ,126854
+ ,99485
+ ,135265
+ ,157406
+ ,153170
+ ,147381
+ ,136753
+ ,158421
+ ,173332
+ ,121795
+ ,123438
+ ,103478
+ ,115420
+ ,120867
+ ,120728
+ ,155286
+ ,172705
+ ,156418
+ ,148870
+ ,171422
+ ,173574
+ ,134051
+ ,123753
+ ,117862
+ ,94628
+ ,133387
+ ,147477
+ ,165077
+ ,150981
+ ,167026
+ ,169155
+ ,160271
+ ,170293
+ ,151312
+ ,160860
+ ,137556
+ ,106249
+ ,108464
+ ,121223
+ ,169114
+ ,155352
+ ,159722
+ ,166638
+ ,165920
+ ,212488
+ ,161262
+ ,168545
+ ,150959
+ ,114509
+ ,128341
+ ,180244
+ ,159268
+ ,149027
+ ,173723
+ ,166053
+ ,162952
+ ,151625
+ ,126219
+ ,152071
+ ,171245
+ ,144467
+ ,137722
+ ,198451
+ ,203411
+ ,174618
+ ,192387
+ ,171786
+ ,167966
+ ,181481
+ ,141108
+ ,164991
+ ,122032
+ ,144562
+ ,163153
+ ,183695
+ ,201470
+ ,166062
+ ,203959
+ ,194742
+ ,165450
+ ,190427
+ ,147520
+ ,166109
+ ,164804
+ ,176121
+ ,149173
+ ,170451
+ ,175412
+ ,166340
+ ,209292
+ ,176079
+ ,165432
+ ,193754
+ ,139788
+ ,158378
+ ,157995
+ ,165045
+ ,113321
+ ,162527
+ ,168683
+ ,164104
+ ,177656
+ ,150257
+ ,177723
+ ,185473
+ ,142068
+ ,130617
+ ,161522
+ ,195235
+ ,151377
+ ,177385
+ ,172286
+ ,158314
+ ,142479
+ ,167690
+ ,149838
+ ,135503
+ ,158414
+ ,167871
+ ,161318
+ ,157863
+ ,141916
+ ,163330
+ ,164925
+ ,152883
+ ,175959
+ ,171660
+ ,176426
+ ,151048
+ ,133931
+ ,144500
+ ,141476
+ ,136659
+ ,138192
+ ,143229
+ ,147819
+ ,139456
+ ,150545
+ ,144622
+ ,147211
+ ,145788
+ ,127517
+ ,134908
+ ,135425
+ ,144149
+ ,134913
+ ,150898
+ ,138370
+ ,146320
+ ,158698
+ ,147410
+ ,155281
+ ,144658
+ ,131474
+ ,143098
+ ,140664
+ ,142011
+ ,139980
+ ,149512
+ ,146613
+ ,141335
+ ,146087
+ ,143822
+ ,156251
+ ,160500
+ ,146866
+ ,152560
+ ,154946
+ ,148239
+ ,150344
+ ,157722
+ ,155789
+ ,149591
+ ,157821
+ ,152981
+ ,165561
+ ,162723
+ ,153944
+ ,157191
+ ,151201
+ ,154338
+ ,151773
+ ,159417
+ ,158863
+ ,153009
+ ,157433
+ ,152753
+ ,160532
+ ,171240
+ ,149951
+ ,158530
+ ,152473
+ ,156337
+ ,155593
+ ,164848
+ ,163900
+ ,159612
+ ,165514
+ ,162519
+ ,169272
+ ,166859
+ ,152068
+ ,164349
+ ,161185
+ ,162728
+ ,158524
+ ,165934
+ ,163616
+ ,159216
+ ,167323
+ ,166526
+ ,176076
+ ,187502
+ ,164936
+ ,172274
+ ,170532
+ ,174122
+ ,165161
+ ,172175
+ ,172293
+ ,164661
+ ,173427
+ ,170402
+ ,165221
+ ,163411
+ ,150358
+ ,162993
+ ,156878
+ ,163197
+ ,162195
+ ,177127
+ ,181363
+ ,170555
+ ,179649
+ ,178388
+ ,183210)
> 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/rcomp/tmp/1mwn21292447090.ps",horizontal=F,onefile=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/rcomp/tmp/2mwn21292447090.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> hist(x)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/3mwn21292447090.ps",horizontal=F,onefile=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/rcomp/tmp/4f5nn1292447090.ps",horizontal=F,onefile=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/rcomp/tmp/5pe4q1292447090.ps",horizontal=F,onefile=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/rcomp/tmp/6pe4q1292447090.ps",horizontal=F,onefile=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/rcomp/tmp/7teke1292447090.ps",horizontal=F,onefile=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. 
  94630  146900  158200  156500  166700  212500 
> 
> #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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/rcomp/tmp/8poin1292447090.tab") 
> 
> try(system("convert tmp/1mwn21292447090.ps tmp/1mwn21292447090.png",intern=TRUE))
character(0)
> try(system("convert tmp/2mwn21292447090.ps tmp/2mwn21292447090.png",intern=TRUE))
character(0)
> try(system("convert tmp/3mwn21292447090.ps tmp/3mwn21292447090.png",intern=TRUE))
character(0)
> try(system("convert tmp/4f5nn1292447090.ps tmp/4f5nn1292447090.png",intern=TRUE))
character(0)
> try(system("convert tmp/5pe4q1292447090.ps tmp/5pe4q1292447090.png",intern=TRUE))
character(0)
> try(system("convert tmp/6pe4q1292447090.ps tmp/6pe4q1292447090.png",intern=TRUE))
character(0)
> try(system("convert tmp/7teke1292447090.ps tmp/7teke1292447090.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.000   1.190   3.164