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(1213.8,1245.6,1306.3,1255.8,1257.6,1287.8,1300.4,1320.9,1370.8,1327.3,1320,1345.3,1346.7,1395.4,1462,1491.6,1461.8,1477.9,1490.3,1521.1,1561.9,1552.6,1523.6,1548.3,1552.4,1587,1621.3,1648.7,1641.8,1650.6,1688.6,1670.7,1682.2,1678.9,1650.6,1662.4,1664.5,1683.2,1736.2,1747.6,1749,1759.7,1793.6,1817.4,1858.4,1839.9,1809.1,1877.7,1880.3,1930.9,2039.3,1992.7,1987.8,1984.4,2016.5,2016.7,2064.1,2031.5,2000.3,2057.8,2041.2,2093.2,2158.3,2128.8,2131.9,2170.3,2190.8,2217.7,2254.4,2223.3,2210.5,2250.8,2249.1,2288.6,2329.2,2313.8,2309.8,2345.9,2361.3,2372,2410.4,2398.5,2362.3,2419.1,2421.6,2465,2480.5,2506.1,2506.6,2525.8,2550,2578.3,2807.8,2815.3,2767.7,2815.4,2838.8,2864,2948.6,2922.8,2917.2,2936.8,2993.4,3007.8,3046.3,3011.5,2958.6,3019.8,2998.5,3040.4,3166,3110,3099.2,3150.3,3163.6,3182.6,3244.4,3223.2,3143.6,3217,3182.3,3217.2,3262.5,3227.9,3171.6,3219,3195.4,3221.6,3262.1,3179.5,3133.6,3219.2,3245,3265.3,3312.5,3383.6,3386.3,3411.1,3467.2,3487.7,3575.5,3571.5,3582.3,3637.1,3685)
> par1 = '0'
> ylab = 'Density'
> xlab = 'Value'
> #'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: 
> if (par1 == '0') bw <- 'nrd0'
> if (par1 != '0') bw <- as.numeric(par1)
> postscript(file="/var/www/html/rcomp/tmp/1gahd1260695416.ps",horizontal=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 (145 obs.);	Bandwidth 'bw' = 236.1

       x                y            
 Min.   : 505.5   Min.   :3.100e-07  
 1st Qu.:1477.4   1st Qu.:5.549e-05  
 Median :2449.4   Median :3.051e-04  
 Mean   :2449.4   Mean   :2.569e-04  
 3rd Qu.:3421.4   3rd Qu.:4.273e-04  
 Max.   :4393.3   Max.   :4.768e-04  
> postscript(file="/var/www/html/rcomp/tmp/2k7tt1260695416.ps",horizontal=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/html/rcomp/tmp/3awfm1260695416.ps",horizontal=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/html/rcomp/tmp/4e4cs1260695416.ps",horizontal=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/html/rcomp/tmp/53n3v1260695416.ps",horizontal=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/html/rcomp/tmp/62uo11260695416.ps",horizontal=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/html/rcomp/tmp/7ltqi1260695416.ps",horizontal=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/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,'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/html/rcomp/tmp/865je1260695416.tab") 
> 
> try(system("convert tmp/1gahd1260695416.ps tmp/1gahd1260695416.png",intern=TRUE))
character(0)
> try(system("convert tmp/2k7tt1260695416.ps tmp/2k7tt1260695416.png",intern=TRUE))
character(0)
> try(system("convert tmp/3awfm1260695416.ps tmp/3awfm1260695416.png",intern=TRUE))
character(0)
> try(system("convert tmp/4e4cs1260695416.ps tmp/4e4cs1260695416.png",intern=TRUE))
character(0)
> try(system("convert tmp/53n3v1260695416.ps tmp/53n3v1260695416.png",intern=TRUE))
character(0)
> try(system("convert tmp/62uo11260695416.ps tmp/62uo11260695416.png",intern=TRUE))
character(0)
> try(system("convert tmp/7ltqi1260695416.ps tmp/7ltqi1260695416.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.304   1.074   2.678