R version 3.2.2 (2015-08-14) -- "Fire Safety"
Copyright (C) 2015 The R Foundation for Statistical Computing
Platform: x86_64-pc-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(350840
+ ,354950
+ ,342750
+ ,318610
+ ,303890
+ ,293480
+ ,286820
+ ,283650
+ ,276190
+ ,268670
+ ,275960
+ ,287040
+ ,291500
+ ,287540
+ ,272340
+ ,256570
+ ,246190
+ ,237340
+ ,235890
+ ,232790
+ ,226100
+ ,221110
+ ,221950
+ ,222990
+ ,232840
+ ,231380
+ ,221690
+ ,210190
+ ,202210
+ ,195210
+ ,194050
+ ,192360
+ ,187170
+ ,185930
+ ,188840
+ ,189350
+ ,194320
+ ,195850
+ ,188840
+ ,182380
+ ,177200
+ ,169140
+ ,171900
+ ,167320
+ ,164130
+ ,162550
+ ,166260
+ ,166010
+ ,177070
+ ,180800
+ ,179570
+ ,178130
+ ,178180
+ ,179430
+ ,179130
+ ,182220
+ ,183650
+ ,187110
+ ,197780
+ ,204000
+ ,221330
+ ,232500
+ ,237570
+ ,241100
+ ,244420
+ ,247240
+ ,254300
+ ,258970
+ ,262230
+ ,268880
+ ,277950
+ ,286630
+ ,299830
+ ,308090
+ ,311400
+ ,313520
+ ,310980
+ ,314910
+ ,310950
+ ,310720
+ ,310380
+ ,310570
+ ,314790
+ ,323400
+ ,335010
+ ,341630
+ ,340820
+ ,336280
+ ,325490
+ ,323750
+ ,317510
+ ,313890
+ ,308610
+ ,303720
+ ,303090
+ ,305140
+ ,304040
+ ,307100
+ ,304330
+ ,294710
+ ,286890
+ ,279050
+ ,271860
+ ,266710
+ ,259590
+ ,253830
+ ,250640
+ ,249140
+ ,250840
+ ,247590
+ ,237830
+ ,226380
+ ,217230
+ ,211420
+ ,207620
+ ,204310
+ ,197490
+ ,193580
+ ,192330
+ ,191970
+ ,196070
+ ,191940
+ ,185620
+ ,179410
+ ,173920
+ ,169190
+ ,166840
+ ,165170
+ ,161450
+ ,160830
+ ,163670
+ ,170830
+ ,182690
+ ,190940
+ ,197770
+ ,205090
+ ,210720
+ ,220210
+ ,229730
+ ,237070
+ ,241620
+ ,250370
+ ,258570
+ ,269860
+ ,283220
+ ,289610
+ ,281770
+ ,274700
+ ,267650
+ ,261380
+ ,260500
+ ,260730
+ ,254200
+ ,250450
+ ,253380
+ ,263740
+ ,276240
+ ,273820
+ ,265890
+ ,258400
+ ,253520
+ ,250710
+ ,252850
+ ,255260
+ ,251170
+ ,252500
+ ,257780
+ ,269900
+ ,291590
+ ,298870
+ ,295570
+ ,292100
+ ,290870
+ ,290580
+ ,297970
+ ,304010
+ ,304340
+ ,309850
+ ,322320
+ ,340170
+ ,369280
+ ,376690
+ ,379700
+ ,379520
+ ,377770
+ ,381560
+ ,394580
+ ,399320
+ ,400370
+ ,408200
+ ,419070
+ ,437730
+ ,460490
+ ,460070
+ ,453830
+ ,443300
+ ,435550)
> par3 = '512'
> par2 = 'no'
> par1 = '0'
> ylab = 'Density'
> xlab = 'Maanden (01-1998 t/m 12-2014)'
> par3 <- '512'
> par2 <- 'no'
> par1 <- '0'
> #'GNU S' R Code compiled by R2WASP v. 1.2.327 (Sat, 26 Sep 2015 10:34:04 +0100)
> #Author: root
> #To cite this work: Wessa, P. (2015), Kernel Density Estimation (v1.0.12) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_density.wasp/
> #Source of accompanying publication: http://www.xycoon.com/density_trace.htm
> #
> if (par1 == '0') bw <- 'nrd0'
> if (par1 != '0') bw <- as.numeric(par1)
> par3 <- as.numeric(par3)
> mydensity <- array(NA, dim=c(par3,8))
> postscript(file="/var/wessaorg/rcomp/tmp/12q1q1443897171.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)
> mydensity[,8] = signif(mydensity1$x,3)
> mydensity[,1] = signif(mydensity1$y,3)
> 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 (197 obs.);	Bandwidth 'bw' = 2.13e+04

       x                y            
 Min.   : 96926   Min.   :2.588e-09  
 1st Qu.:203793   1st Qu.:5.465e-07  
 Median :310660   Median :1.305e-06  
 Mean   :310660   Mean   :2.337e-06  
 3rd Qu.:417527   3rd Qu.:4.552e-06  
 Max.   :524394   Max.   :5.697e-06  
> postscript(file="/var/wessaorg/rcomp/tmp/2dof01443897171.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)
> mydensity[,2] = signif(mydensity2$y,3)
> plot(mydensity2,main='Epanechnikov Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/3oohc1443897171.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)
> mydensity[,3] = signif(mydensity3$y,3)
> plot(mydensity3,main='Rectangular Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/4pve81443897171.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)
> mydensity[,4] = signif(mydensity4$y,3)
> plot(mydensity4,main='Triangular Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5v6zs1443897171.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)
> mydensity[,5] = signif(mydensity5$y,3)
> plot(mydensity5,main='Biweight Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6b4jd1443897171.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)
> mydensity[,6] = signif(mydensity6$y,3)
> plot(mydensity6,main='Cosine Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/7d8fb1443897171.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)
> mydensity[,7] = signif(mydensity7$y,3)
> plot(mydensity7,main='Optcosine Kernel',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> 
> #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/rcomp/createtable")
> 
> ab<-table.start()
> ab<-table.row.start(ab)
> ab<-table.element(ab,'Properties of Density Trace',2,TRUE)
> ab<-table.row.end(ab)
> ab<-table.row.start(ab)
> ab<-table.element(ab,'Bandwidth',header=TRUE)
> ab<-table.element(ab,mydensity1$bw)
> ab<-table.row.end(ab)
> ab<-table.row.start(ab)
> ab<-table.element(ab,'#Observations',header=TRUE)
> ab<-table.element(ab,mydensity1$n)
> ab<-table.row.end(ab)
> ab<-table.end(ab)
> a <- ab
> table.save(a,file="/var/wessaorg/rcomp/tmp/8u6si1443897171.tab") 
> b<-table.start()
> b<-table.row.start(b)
> b<-table.element(b,'Maximum Density Values',3,TRUE)
> b<-table.row.end(b)
> b<-table.row.start(b)
> b<-table.element(b,'Kernel',1,TRUE)
> b<-table.element(b,'x-value',1,TRUE)
> b<-table.element(b,'max. density',1,TRUE)
> b<-table.row.end(b)
> b<-table.row.start(b)
> b<-table.element(b,'Gaussian',1,TRUE)
> b<-table.element(b,mydensity1$x[mydensity1$y==max(mydensity1$y)],1)
> b<-table.element(b,mydensity1$y[mydensity1$y==max(mydensity1$y)],1)
> b<-table.row.end(b)
> b<-table.row.start(b)
> b<-table.element(b,'Epanechnikov',1,TRUE)
> b<-table.element(b,mydensity2$x[mydensity2$y==max(mydensity2$y)],1)
> b<-table.element(b,mydensity2$y[mydensity2$y==max(mydensity2$y)],1)
> b<-table.row.end(b)
> b<-table.row.start(b)
> b<-table.element(b,'Rectangular',1,TRUE)
> b<-table.element(b,mydensity3$x[mydensity3$y==max(mydensity3$y)],1)
> b<-table.element(b,mydensity3$y[mydensity3$y==max(mydensity3$y)],1)
> b<-table.row.end(b)
> b<-table.row.start(b)
> b<-table.element(b,'Triangular',1,TRUE)
> b<-table.element(b,mydensity4$x[mydensity4$y==max(mydensity4$y)],1)
> b<-table.element(b,mydensity4$y[mydensity4$y==max(mydensity4$y)],1)
> b<-table.row.end(b)
> b<-table.row.start(b)
> b<-table.element(b,'Biweight',1,TRUE)
> b<-table.element(b,mydensity5$x[mydensity5$y==max(mydensity5$y)],1)
> b<-table.element(b,mydensity5$y[mydensity5$y==max(mydensity5$y)],1)
> b<-table.row.end(b)
> b<-table.row.start(b)
> b<-table.element(b,'Cosine',1,TRUE)
> b<-table.element(b,mydensity6$x[mydensity6$y==max(mydensity6$y)],1)
> b<-table.element(b,mydensity6$y[mydensity6$y==max(mydensity6$y)],1)
> b<-table.row.end(b)
> b<-table.row.start(b)
> b<-table.element(b,'Optcosine',1,TRUE)
> b<-table.element(b,mydensity7$x[mydensity7$y==max(mydensity7$y)],1)
> b<-table.element(b,mydensity7$y[mydensity7$y==max(mydensity7$y)],1)
> b<-table.row.end(b)
> b<-table.end(b)
> a <- b[1]
> table.save(a,file="/var/wessaorg/rcomp/tmp/9k8rw1443897171.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Kernel Density Values',8,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'x-value',1,TRUE)
> a<-table.element(a,'Gaussian',1,TRUE)
> a<-table.element(a,'Epanechnikov',1,TRUE)
> a<-table.element(a,'Rectangular',1,TRUE)
> a<-table.element(a,'Triangular',1,TRUE)
> a<-table.element(a,'Biweight',1,TRUE)
> a<-table.element(a,'Cosine',1,TRUE)
> a<-table.element(a,'Optcosine',1,TRUE)
> a<-table.row.end(a)
> if (par2=='yes') {
+ for(i in 1:par3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mydensity[i,8],1,TRUE)
+ for(j in 1:7) {
+ a<-table.element(a,mydensity[i,j],1)
+ }
+ a<-table.row.end(a)
+ }
+ } else {
+ a<-table.row.start(a)
+ a<-table.element(a,'Kernel Density Values are not shown',8)
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/10241s1443897171.tab") 
> 
> try(system("convert tmp/12q1q1443897171.ps tmp/12q1q1443897171.png",intern=TRUE))
character(0)
> try(system("convert tmp/2dof01443897171.ps tmp/2dof01443897171.png",intern=TRUE))
character(0)
> try(system("convert tmp/3oohc1443897171.ps tmp/3oohc1443897171.png",intern=TRUE))
character(0)
> try(system("convert tmp/4pve81443897171.ps tmp/4pve81443897171.png",intern=TRUE))
character(0)
> try(system("convert tmp/5v6zs1443897171.ps tmp/5v6zs1443897171.png",intern=TRUE))
character(0)
> try(system("convert tmp/6b4jd1443897171.ps tmp/6b4jd1443897171.png",intern=TRUE))
character(0)
> try(system("convert tmp/7d8fb1443897171.ps tmp/7d8fb1443897171.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.615   0.495   3.137