R version 2.10.1 (2009-12-14)
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(589,561,640,656,727,697,640,599,568,577,553,582,600,566,653,673,742,716,660,617,583,587,565,598,628,618,688,705,770,736,678,639,604,611,594,634,658,622,709,722,782,756,702,653,615,621,602,635,677,635,736,755,811,798,735,697,661,667,645,688,713,667,762,784,837,817,767,722,681,687,660,698,717,696,775,796,858,826,783,740,701,706,677,711,734,690,785,805,871,845,801,764,725,723,690,734,750,707,807,824,886,859,819,783,740,747,711,751,804,756,860,878,942,913,869,834,790,800,763,800,826,799,890,900,961,935,894,855,809,810,766,805,821,773,883,898,957,924,881,837,784,791,760,802,828,778,889,902,969,947,908,867,815,812,773,813,834,782,892,903,966,937,896,858,817,827,797,843)
> par4 = '12'
> par3 = 'Maandelijkse melkproductie'
> par2 = 'http://www.robjhyndman.com/TSDL/index.htm'
> par1 = 'Maandelijkse melkproductie per koe'
> ylimmax = ''
> ylimmin = ''
> ylab = 'Aantal melk'
> xlab = 'Maand'
> main = 'Maandelijkse melkproductie per koe'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Dr. Ian E. Holliday
> #To cite this work: Ian E. Holliday, 2009, 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: 
> #Technical description: 
> if (par4 != 'No season') {
+ par4 <- as.numeric(par4)
+ if (par4 < 4) par4 <- 12
+ }
> summary(x)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  553.0   677.8   761.0   754.7   824.5   969.0 
> n <- length(x)
> postscript(file="/var/www/rcomp/tmp/151sv1275129335.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> if (par4=='No season') {
+ plot(x,col=2,type='b',main=main,xlab=xlab,ylab=ylab,xaxt='n')
+ axis(1,at=seq(1,n,10))
+ }
> if (par4!='No season') {
+ plot(x,col=2,type='b',main=main,xlab=xlab,ylab=ylab,xaxt='n')
+ axis(1,at=seq(1,n,par4))
+ grid(nx=0,ny=NULL,col='black')
+ abline(v=seq(1,n,par4),col='black',lty='dotted')
+ }
> dev.off()
null device 
          1 
> 
> #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,'Univariate Dataseries',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Name of dataseries',header=TRUE)
> a<-table.element(a,par1)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Source',header=TRUE)
> a<-table.element(a,par2)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Description',header=TRUE)
> a<-table.element(a,par3)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Number of observations',header=TRUE)
> a<-table.element(a,length(x))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/2ck771275129335.tab") 
> 
> try(system("convert tmp/151sv1275129335.ps tmp/151sv1275129335.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.450   0.310   0.708