R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-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(49.164,50.112,52.329,52.872,53.008,53.804,57.150,57.863,58.334,60.621,62.084,62.395,61.064,65.989,67.170,68.813,69.255,70.643,73.657,72.723,72.966,75.159,74.414,75.573,70.977,72.291,72.384,74.385,70.977,68.889,71.950,70.463,70.615,72.567,75.962,75.851,75.005,78.635,82.937,80.337,82.336,81.858,83.732,84.158,84.806,86.610,91.300,92.038,99.274,99.292,100.169,103.661,100.427,106.692,108.459,109.815,108.914,111.181,114.572,115.019,149.365,154.780,163.305,168.139,170.969,178.788,184.302,182.706,189.415,192.260,200.070,200.555,203.220,207.325,209.362,205.709,205.752,204.527,208.881,212.113,216.893,219.171,226.970,222.226)
> par1 = '0'
> ylab = 'kp v/e woonhuis - in duizendtallen'
> xlab = 'per kwartaal - 1990 -2010'
> #'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/wessaorg/rcomp/tmp/1a7351317807124.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 (84 obs.);	Bandwidth 'bw' = 20.92

       x                y            
 Min.   :-13.60   Min.   :6.044e-06  
 1st Qu.: 62.23   1st Qu.:8.420e-04  
 Median :138.07   Median :2.482e-03  
 Mean   :138.07   Mean   :3.293e-03  
 3rd Qu.:213.90   3rd Qu.:4.080e-03  
 Max.   :289.74   Max.   :1.046e-02  
> postscript(file="/var/wessaorg/rcomp/tmp/2mbms1317807124.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/wessaorg/rcomp/tmp/38bvd1317807124.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/wessaorg/rcomp/tmp/49cfn1317807124.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/wessaorg/rcomp/tmp/547b71317807124.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/wessaorg/rcomp/tmp/6hlcn1317807124.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/wessaorg/rcomp/tmp/79lnl1317807124.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/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/wessaorg/rcomp/tmp/8n0cx1317807124.tab") 
> 
> try(system("convert tmp/1a7351317807124.ps tmp/1a7351317807124.png",intern=TRUE))
character(0)
> try(system("convert tmp/2mbms1317807124.ps tmp/2mbms1317807124.png",intern=TRUE))
character(0)
> try(system("convert tmp/38bvd1317807124.ps tmp/38bvd1317807124.png",intern=TRUE))
character(0)
> try(system("convert tmp/49cfn1317807124.ps tmp/49cfn1317807124.png",intern=TRUE))
character(0)
> try(system("convert tmp/547b71317807124.ps tmp/547b71317807124.png",intern=TRUE))
character(0)
> try(system("convert tmp/6hlcn1317807124.ps tmp/6hlcn1317807124.png",intern=TRUE))
character(0)
> try(system("convert tmp/79lnl1317807124.ps tmp/79lnl1317807124.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.725   0.335   2.068