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(12.182
+ ,12.242
+ ,12.368
+ ,12.552
+ ,12.377
+ ,12.350
+ ,12.333
+ ,12.427
+ ,12.281
+ ,12.381
+ ,12.569
+ ,12.684
+ ,12.689
+ ,12.580
+ ,12.264
+ ,12.259
+ ,12.205
+ ,12.255
+ ,12.039
+ ,11.893
+ ,11.741
+ ,12.149
+ ,12.096
+ ,12.146
+ ,11.946
+ ,11.976
+ ,12.392
+ ,12.102
+ ,12.362
+ ,12.547
+ ,12.532
+ ,12.422
+ ,12.304
+ ,12.216
+ ,12.267
+ ,12.652
+ ,12.605
+ ,12.626
+ ,12.613
+ ,12.603
+ ,12.575
+ ,12.527
+ ,12.580
+ ,12.325
+ ,12.304
+ ,12.372
+ ,12.617
+ ,12.627
+ ,12.851
+ ,12.825
+ ,12.721
+ ,12.765
+ ,12.848
+ ,12.891
+ ,12.870
+ ,12.831
+ ,12.818
+ ,13.013
+ ,13.057
+ ,12.969
+ ,13.011
+ ,12.815
+ ,12.861
+ ,12.746
+ ,12.872
+ ,12.825
+ ,12.891
+ ,12.993
+ ,12.985
+ ,13.026
+ ,12.825
+ ,12.598
+ ,12.621
+ ,12.480
+ ,12.543
+ ,12.594
+ ,12.647
+ ,12.638
+ ,12.503
+ ,12.392
+ ,12.389
+ ,12.603
+ ,12.210
+ ,12.278
+ ,12.286
+ ,12.090
+ ,12.145
+ ,12.307
+ ,12.270
+ ,12.159
+ ,12.023
+ ,12.062
+ ,11.844
+ ,11.842
+ ,11.805
+ ,11.809
+ ,11.453
+ ,11.346
+ ,11.345
+ ,11.382
+ ,11.216
+ ,11.289
+ ,11.225
+ ,11.382
+ ,11.148
+ ,11.226
+ ,11.104
+ ,11.051
+ ,10.962
+ ,11.238
+ ,11.437
+ ,11.495
+ ,11.458
+ ,11.603
+ ,11.630
+ ,11.341
+ ,11.369
+ ,11.133
+ ,11.398
+ ,11.578
+ ,11.380
+ ,11.326
+ ,11.286
+ ,11.604
+ ,11.655
+ ,11.432
+ ,11.730
+ ,11.782
+ ,11.633
+ ,11.532
+ ,11.611
+ ,11.660
+ ,11.478
+ ,11.346
+ ,11.415
+ ,11.427
+ ,11.626
+ ,11.384
+ ,11.412
+ ,11.500
+ ,11.713
+ ,11.546
+ ,11.506
+ ,11.532
+ ,11.186
+ ,11.225
+ ,11.515
+ ,11.234
+ ,11.264
+ ,11.429
+ ,11.416
+ ,10.906
+ ,11.057
+ ,10.609
+ ,11.028
+ ,11.394
+ ,11.016
+ ,10.850
+ ,10.827
+ ,11.019
+ ,11.140
+ ,10.372
+ ,10.847
+ ,10.826
+ ,10.484
+ ,10.323
+ ,9.955
+ ,9.437
+ ,9.262
+ ,8.569
+ ,8.462
+ ,9.389
+ ,9.302
+ ,8.577
+ ,8.975
+ ,8.852
+ ,9.263
+ ,9.028
+ ,8.520
+ ,8.683
+ ,8.376
+ ,8.179
+ ,9.062
+ ,9.005
+ ,9.179
+ ,9.326
+ ,9.324
+ ,9.617
+ ,9.134
+ ,8.696
+ ,8.947
+ ,8.864
+ ,8.685
+ ,8.822
+ ,8.495
+ ,8.273
+ ,8.421
+ ,7.996
+ ,7.552
+ ,8.048
+ ,8.445
+ ,8.464
+ ,8.725
+ ,8.827
+ ,8.154
+ ,8.409
+ ,8.587
+ ,8.376
+ ,8.638
+ ,8.934
+ ,8.693
+ ,8.750
+ ,8.563
+ ,8.629
+ ,8.566
+ ,8.922
+ ,8.824
+ ,8.607
+ ,8.573
+ ,8.519
+ ,8.417
+ ,8.469
+ ,8.516
+ ,8.488
+ ,8.666
+ ,8.772
+ ,9.027
+ ,8.955
+ ,8.997
+ ,8.770
+ ,8.739
+ ,8.599
+ ,8.475
+ ,8.446
+ ,8.196
+ ,8.216
+ ,8.280
+ ,7.949
+ ,8.224
+ ,8.109
+ ,8.078
+ ,8.117
+ ,8.176
+ ,8.373
+ ,8.149
+ ,8.001
+ ,7.937
+ ,8.070
+ ,7.955
+ ,8.056)
> par1 = '0'
> ylab = 'WAARDE VAN HET AANDEEL'
> xlab = 'TIJD IN WERKDAGEN'
> #'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/1s4wl1236064432.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 (250 obs.);	Bandwidth 'bw' = 0.5057

       x                y            
 Min.   : 6.035   Min.   :5.593e-05  
 1st Qu.: 8.170   1st Qu.:2.803e-02  
 Median :10.305   Median :9.865e-02  
 Mean   :10.305   Mean   :1.170e-01  
 3rd Qu.:12.439   3rd Qu.:1.966e-01  
 Max.   :14.574   Max.   :2.879e-01  
> postscript(file="/var/www/html/rcomp/tmp/25txx1236064432.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/32d9d1236064432.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/4t4tj1236064432.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/5h2w91236064432.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/67n781236064433.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/73a5v1236064433.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/8e8v11236064433.tab") 
> 
> system("convert tmp/1s4wl1236064432.ps tmp/1s4wl1236064432.png")
> system("convert tmp/25txx1236064432.ps tmp/25txx1236064432.png")
> system("convert tmp/32d9d1236064432.ps tmp/32d9d1236064432.png")
> system("convert tmp/4t4tj1236064432.ps tmp/4t4tj1236064432.png")
> system("convert tmp/5h2w91236064432.ps tmp/5h2w91236064432.png")
> system("convert tmp/67n781236064433.ps tmp/67n781236064433.png")
> system("convert tmp/73a5v1236064433.ps tmp/73a5v1236064433.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.341   1.125   1.856