R version 2.8.1 (2008-12-22)
Copyright (C) 2008 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(2585,3368,3210,3111,3756,4216,5225,4426,3932,3816,3661,3795,2285,2934,2985,3646,4198,4935,5618,5454,3624,2898,3802,3629,2369,2511,3079,3728,4151,4326,5054,5138,3310,3508,3790,3446,2127,2523,3017,3265,3822,4027,4420,5255,4009,3074,3465,3718,1954,2604,3626,2836,4042,3584,4225,4523,2892,2876,3420,3159,2101,2181,2724,2954,4092,3470,3990,4239,2855,2897,3433,3307,1914,2214,2320,2714,3633,3295,4377,4442,2774,2840,2828,3758,1610,1968,2248,3262,3164,2972,4041,3402,2898,2555,3056,3717,1755,2193,2198,2777,3076,3389,4231,3118,2524,2280,2862,3502,1558,1940,2226,2676,3145,3224,4117,3446,2482,2349,2986,3163,1651,1725,2622,2316,2976,3263,3951,2917,2380,2458,2883,2579,1330,1686,2457,2514,2834,2757,3425,3006,2369,2017,2507,3168,1545,1643,2112,2415,2862,2822,3260,2606,2264,2250,2545,2856,1208,1412,1964,2018,2329,2660,2923,2626,2132,1772,2526,2755,1154,1568,1965,2659,2354,2592,2714,2294,2416,2016,2799,2467,1153,1482,1818,2262,2612,2967,3179)
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = ''
> par2 = ''
> par1 = '0'
> ylab = ''
> xlab = ''
> main = ''
> #'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/1usng1234974099.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 (187 obs.);	Bandwidth 'bw' = 265.3

       x                y            
 Min.   : 357.1   Min.   :1.033e-07  
 1st Qu.:1871.3   1st Qu.:2.933e-05  
 Median :3385.5   Median :1.059e-04  
 Mean   :3385.5   Mean   :1.649e-04  
 3rd Qu.:4899.7   3rd Qu.:2.900e-04  
 Max.   :6413.9   Max.   :4.762e-04  
> postscript(file="/var/www/html/rcomp/tmp/2olhb1234974099.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/3wlo21234974099.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/47ug11234974099.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/5gyib1234974099.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/60eyc1234974099.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/7tg531234974099.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/8haly1234974099.tab") 
> 
> system("convert tmp/1usng1234974099.ps tmp/1usng1234974099.png")
> system("convert tmp/2olhb1234974099.ps tmp/2olhb1234974099.png")
> system("convert tmp/3wlo21234974099.ps tmp/3wlo21234974099.png")
> system("convert tmp/47ug11234974099.ps tmp/47ug11234974099.png")
> system("convert tmp/5gyib1234974099.ps tmp/5gyib1234974099.png")
> system("convert tmp/60eyc1234974099.ps tmp/60eyc1234974099.png")
> system("convert tmp/7tg531234974099.ps tmp/7tg531234974099.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.341   1.105   4.719