R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-redhat-linux-gnu (64-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(7992,6114,5965,8460,8323,6333,5675,10090,9035,6976,6459,10896,9978,7466,7199,10977,9412,6341,7784,11911,10079,7721,8197,12038,11963,8033,8618,13625,11734,8895,8727,13974,12583,9525,9662,15490,13839,10047,9788,14978,13045,9489,8741,13149,14106,9998,10034,15081,13266,9997,9027,14324,13149,11209,10332,15354,13800,11786,10550,16114,13255,11403,10269,14009,15847,12967,11328,15814,18626,13219,13818,18062,15722,12111,11702,15589,14852,13612,12380,15501,16322,12157,11124,14621,14035,11159,10944,15824,14378,11816,12233,17344,16812,12181,13275,18458,17375,14609,13323,18327,16053,15070,13806,18245,17461,14999,16022,20564,16372,15854,15115,18207,19488,16644,18631,21093,22212,19762,19403,21227,23176,20823,20647,21336,23458,22003,21647,26416,25226,24723,19945,24040,25034,24885,21168,23541,26019,24657,20599,24534,28717,26138,22968,26577,28660,30430,27356,25454,30194)
> par1 = '0'
> ylab = 'Density'
> xlab = 'Ton'
> if (par1 == '0') bw <- 'nrd0'
> if (par1 != '0') bw <- as.numeric(par1)
> postscript(file="/var/www/wessaorg/rcomp/tmp/1xju61298369036.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> mydensity1<-density(x,bw=bw,kernel='gaussian',na.rm=TRUE)
> plot(mydensity1,main='Gaussian Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> mydensity1

Call:
	density.default(x = x, bw = bw, kernel = "gaussian", na.rm = TRUE)

Data: x (149 obs.);	Bandwidth 'bw' = 1908

       x                  y            
 Min.   :  -50.34   Min.   :2.800e-08  
 1st Qu.: 9001.08   1st Qu.:5.024e-06  
 Median :18052.50   Median :2.420e-05  
 Mean   :18052.50   Mean   :2.759e-05  
 3rd Qu.:27103.92   3rd Qu.:4.591e-05  
 Max.   :36155.34   Max.   :7.181e-05  
> postscript(file="/var/www/wessaorg/rcomp/tmp/2tf4w1298369036.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> mydensity2<-density(x,bw=bw,kernel='epanechnikov',na.rm=TRUE)
> plot(mydensity2,main='Epanechnikov Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/32z531298369036.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> mydensity3<-density(x,bw=bw,kernel='rectangular',na.rm=TRUE)
> plot(mydensity3,main='Rectangular Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/4q8bw1298369036.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> mydensity4<-density(x,bw=bw,kernel='triangular',na.rm=TRUE)
> plot(mydensity4,main='Triangular Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/50qs91298369036.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> mydensity5<-density(x,bw=bw,kernel='biweight',na.rm=TRUE)
> plot(mydensity5,main='Biweight Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/62nzr1298369036.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> mydensity6<-density(x,bw=bw,kernel='cosine',na.rm=TRUE)
> plot(mydensity6,main='Cosine Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/7hozi1298369036.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> mydensity7<-density(x,bw=bw,kernel='optcosine',na.rm=TRUE)
> plot(mydensity7,main='Optcosine Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/wessaorg/rcomp/createtable")
> 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Properties of Density Trace',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Bandwidth',header=TRUE)
> a<-table.element(a,mydensity1$bw)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'#Observations',header=TRUE)
> a<-table.element(a,mydensity1$n)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/wessaorg/rcomp/tmp/8de9g1298369036.tab") 
> 
> try(system("convert tmp/1xju61298369036.ps tmp/1xju61298369036.png",intern=TRUE))
character(0)
> try(system("convert tmp/2tf4w1298369036.ps tmp/2tf4w1298369036.png",intern=TRUE))
character(0)
> try(system("convert tmp/32z531298369036.ps tmp/32z531298369036.png",intern=TRUE))
character(0)
> try(system("convert tmp/4q8bw1298369036.ps tmp/4q8bw1298369036.png",intern=TRUE))
character(0)
> try(system("convert tmp/50qs91298369036.ps tmp/50qs91298369036.png",intern=TRUE))
character(0)
> try(system("convert tmp/62nzr1298369036.ps tmp/62nzr1298369036.png",intern=TRUE))
character(0)
> try(system("convert tmp/7hozi1298369036.ps tmp/7hozi1298369036.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.680   0.310   2.092