R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"
Copyright (C) 2013 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(-33.39
+ ,-31.80
+ ,-28.19
+ ,-26.42
+ ,-24.66
+ ,-22.54
+ ,-22.23
+ ,-21.80
+ ,-21.67
+ ,-21.61
+ ,-20.41
+ ,-20.15
+ ,-19.83
+ ,-19.22
+ ,-19.12
+ ,-18.85
+ ,-18.12
+ ,-18.07
+ ,-17.93
+ ,-17.44
+ ,-17.42
+ ,-17.25
+ ,-16.99
+ ,-16.85
+ ,-16.31
+ ,-15.80
+ ,-15.80
+ ,-15.80
+ ,-15.66
+ ,-15.60
+ ,-14.77
+ ,-14.75
+ ,-14.17
+ ,-13.71
+ ,-13.65
+ ,-13.54
+ ,-13.30
+ ,-13.25
+ ,-12.57
+ ,-12.54
+ ,-12.13
+ ,-11.82
+ ,-11.68
+ ,-11.52
+ ,-11.44
+ ,-11.35
+ ,-11.24
+ ,-11.21
+ ,-11.10
+ ,-10.72
+ ,-10.67
+ ,-10.63
+ ,-10.51
+ ,-10.43
+ ,-10.42
+ ,-10.41
+ ,-9.96
+ ,-9.93
+ ,-9.89
+ ,-9.82
+ ,-9.71
+ ,-9.59
+ ,-9.51
+ ,-9.45
+ ,-9.33
+ ,-9.26
+ ,-8.81
+ ,-8.81
+ ,-8.74
+ ,-8.72
+ ,-8.66
+ ,-8.64
+ ,-8.63
+ ,-8.62
+ ,-8.59
+ ,-8.57
+ ,-8.47
+ ,-8.39
+ ,-8.31
+ ,-8.22
+ ,-8.07
+ ,-8.06
+ ,-8.01
+ ,-7.89
+ ,-7.77
+ ,-7.68
+ ,-7.65
+ ,-7.45
+ ,-7.38
+ ,-7.37
+ ,-7.33
+ ,-7.16
+ ,-7.06
+ ,-7.05
+ ,-7.04
+ ,-7.02
+ ,-6.83
+ ,-6.81
+ ,-6.80
+ ,-6.76
+ ,-6.75
+ ,-6.72
+ ,-6.71
+ ,-6.71
+ ,-6.70
+ ,-6.66
+ ,-6.65
+ ,-6.55
+ ,-6.54
+ ,-6.53
+ ,-6.53
+ ,-6.49
+ ,-6.42
+ ,-6.40
+ ,-6.35
+ ,-6.32
+ ,-6.23
+ ,-6.20
+ ,-6.13
+ ,-6.07
+ ,-6.06
+ ,-6.05
+ ,-6.01
+ ,-5.99
+ ,-5.93
+ ,-5.68
+ ,-5.63
+ ,-5.60
+ ,-5.58
+ ,-5.52
+ ,-5.42
+ ,-5.39
+ ,-5.38
+ ,-5.35
+ ,-5.30
+ ,-5.27
+ ,-5.17
+ ,-5.17
+ ,-5.15
+ ,-5.07
+ ,-5.05
+ ,-5.04
+ ,-5.00
+ ,-4.99
+ ,-4.99
+ ,-4.97
+ ,-4.86
+ ,-4.82
+ ,-4.81
+ ,-4.65
+ ,-4.64
+ ,-4.51
+ ,-4.50
+ ,-4.49
+ ,-4.41
+ ,-4.41
+ ,-4.40
+ ,-4.40
+ ,-4.40
+ ,-4.39
+ ,-4.32
+ ,-4.22
+ ,-4.14
+ ,-4.09
+ ,-4.09
+ ,-4.05
+ ,-4.02
+ ,-3.92
+ ,-3.92
+ ,-3.91
+ ,-3.86
+ ,-3.83
+ ,-3.79
+ ,-3.71
+ ,-3.68
+ ,-3.67
+ ,-3.67
+ ,-3.62
+ ,-3.56
+ ,-3.54
+ ,-3.53
+ ,-3.52
+ ,-3.51
+ ,-3.42
+ ,-3.40
+ ,-3.36
+ ,-3.33
+ ,-3.31
+ ,-3.28
+ ,-3.26
+ ,-3.25
+ ,-3.23
+ ,-3.23
+ ,-3.23
+ ,-3.16
+ ,-3.15
+ ,-3.13
+ ,-3.10
+ ,-3.09
+ ,-3.06
+ ,-3.05
+ ,-3.04
+ ,-3.04
+ ,-3.00
+ ,-2.93
+ ,-2.92
+ ,-2.84
+ ,-2.81
+ ,-2.80
+ ,-2.80
+ ,-2.78
+ ,-2.78
+ ,-2.77
+ ,-2.76
+ ,-2.74
+ ,-2.74
+ ,-2.72
+ ,-2.72
+ ,-2.67
+ ,-2.65
+ ,-2.65
+ ,-2.65
+ ,-2.64
+ ,-2.63
+ ,-2.62
+ ,-2.61
+ ,-2.61
+ ,-2.60
+ ,-2.59
+ ,-2.58
+ ,-2.57
+ ,-2.55
+ ,-2.52
+ ,-2.51
+ ,-2.51
+ ,-2.48
+ ,-2.44
+ ,-2.44
+ ,-2.43
+ ,-2.43
+ ,-2.38
+ ,-2.37
+ ,-2.33
+ ,-2.33
+ ,-2.33
+ ,-2.32
+ ,-2.32
+ ,-2.31
+ ,-2.28
+ ,-2.28
+ ,-2.26
+ ,-2.23
+ ,-2.18
+ ,-2.14
+ ,-2.13
+ ,-2.13
+ ,-2.11
+ ,-2.05
+ ,-2.04
+ ,-2.03
+ ,-1.96
+ ,-1.90
+ ,-1.89
+ ,-1.89
+ ,-1.86
+ ,-1.83
+ ,-1.83
+ ,-1.83
+ ,-1.79
+ ,-1.76
+ ,-1.74
+ ,-1.68
+ ,-1.67
+ ,-1.66
+ ,-1.66
+ ,-1.66
+ ,-1.64
+ ,-1.62
+ ,-1.61
+ ,-1.61
+ ,-1.60
+ ,-1.57
+ ,-1.56
+ ,-1.53
+ ,-1.51
+ ,-1.49
+ ,-1.46
+ ,-1.43
+ ,-1.43
+ ,-1.43
+ ,-1.42
+ ,-1.39
+ ,-1.32
+ ,-1.32
+ ,-1.32
+ ,-1.31
+ ,-1.30
+ ,-1.30
+ ,-1.29
+ ,-1.29
+ ,-1.27
+ ,-1.27
+ ,-1.26
+ ,-1.22
+ ,-1.21
+ ,-1.19
+ ,-1.18
+ ,-1.16
+ ,-1.14
+ ,-1.13
+ ,-1.13
+ ,-1.12
+ ,-1.10
+ ,-1.09
+ ,-1.09
+ ,-1.08
+ ,-1.04
+ ,-1.00
+ ,-0.97
+ ,-0.96
+ ,-0.95
+ ,-0.94
+ ,-0.94
+ ,-0.93
+ ,-0.91
+ ,-0.88
+ ,-0.86
+ ,-0.82
+ ,-0.77
+ ,-0.76
+ ,-0.75
+ ,-0.74
+ ,-0.72
+ ,-0.72
+ ,-0.71
+ ,-0.71
+ ,-0.66
+ ,-0.66
+ ,-0.65
+ ,-0.65
+ ,-0.64
+ ,-0.64
+ ,-0.62
+ ,-0.62
+ ,-0.60
+ ,-0.58
+ ,-0.57
+ ,-0.57
+ ,-0.57
+ ,-0.56
+ ,-0.54
+ ,-0.53
+ ,-0.48
+ ,-0.47
+ ,-0.47
+ ,-0.45
+ ,-0.36
+ ,-0.36
+ ,-0.35
+ ,-0.34
+ ,-0.33
+ ,-0.32
+ ,-0.31
+ ,-0.30
+ ,-0.30
+ ,-0.29
+ ,-0.28
+ ,-0.25
+ ,-0.24
+ ,-0.22
+ ,-0.22
+ ,-0.21
+ ,-0.21
+ ,-0.19
+ ,-0.17
+ ,-0.15
+ ,-0.15
+ ,-0.15
+ ,-0.09
+ ,-0.08
+ ,-0.07
+ ,-0.07
+ ,-0.06
+ ,-0.05
+ ,-0.01
+ ,0.02
+ ,0.02
+ ,0.03
+ ,0.05
+ ,0.05
+ ,0.08
+ ,0.09
+ ,0.09
+ ,0.10
+ ,0.14
+ ,0.17
+ ,0.18
+ ,0.20
+ ,0.22
+ ,0.24
+ ,0.25
+ ,0.25
+ ,0.26
+ ,0.26
+ ,0.29
+ ,0.31
+ ,0.31
+ ,0.33
+ ,0.33
+ ,0.33
+ ,0.36
+ ,0.37
+ ,0.37
+ ,0.38
+ ,0.39
+ ,0.40
+ ,0.41
+ ,0.43
+ ,0.43
+ ,0.44
+ ,0.45
+ ,0.45
+ ,0.48
+ ,0.48
+ ,0.48
+ ,0.48
+ ,0.50
+ ,0.50
+ ,0.51
+ ,0.52
+ ,0.52
+ ,0.53
+ ,0.54
+ ,0.56
+ ,0.57
+ ,0.57
+ ,0.60
+ ,0.61
+ ,0.61
+ ,0.66
+ ,0.66
+ ,0.67
+ ,0.69
+ ,0.71
+ ,0.71
+ ,0.72
+ ,0.73
+ ,0.73
+ ,0.74
+ ,0.74
+ ,0.74
+ ,0.75
+ ,0.76
+ ,0.76
+ ,0.77
+ ,0.78
+ ,0.78
+ ,0.79
+ ,0.80
+ ,0.80
+ ,0.82
+ ,0.88
+ ,0.88
+ ,0.94
+ ,0.97
+ ,0.99
+ ,1.03
+ ,1.04
+ ,1.06
+ ,1.09
+ ,1.13
+ ,1.14
+ ,1.15
+ ,1.15
+ ,1.23
+ ,1.25
+ ,1.25
+ ,1.28
+ ,1.29
+ ,1.30
+ ,1.30
+ ,1.32
+ ,1.32
+ ,1.32
+ ,1.33
+ ,1.34
+ ,1.36
+ ,1.36
+ ,1.37
+ ,1.37
+ ,1.39
+ ,1.39
+ ,1.42
+ ,1.43
+ ,1.45
+ ,1.45
+ ,1.46
+ ,1.46
+ ,1.47
+ ,1.47
+ ,1.47
+ ,1.49
+ ,1.49
+ ,1.50
+ ,1.52
+ ,1.53
+ ,1.54
+ ,1.55
+ ,1.55
+ ,1.58
+ ,1.59
+ ,1.59
+ ,1.61
+ ,1.63
+ ,1.65
+ ,1.65
+ ,1.65
+ ,1.70
+ ,1.72
+ ,1.73
+ ,1.73
+ ,1.74
+ ,1.75
+ ,1.78
+ ,1.79
+ ,1.81
+ ,1.81
+ ,1.83
+ ,1.84
+ ,1.84
+ ,1.84
+ ,1.85
+ ,1.86
+ ,1.87
+ ,1.87
+ ,1.88
+ ,1.88
+ ,1.88
+ ,1.89
+ ,1.91
+ ,1.91
+ ,1.92
+ ,1.93
+ ,1.94
+ ,1.95
+ ,1.95
+ ,1.97
+ ,1.98
+ ,1.98
+ ,2.03
+ ,2.05
+ ,2.09
+ ,2.14
+ ,2.15
+ ,2.15
+ ,2.16
+ ,2.16
+ ,2.16
+ ,2.17
+ ,2.20
+ ,2.20
+ ,2.21
+ ,2.22
+ ,2.24
+ ,2.27
+ ,2.28
+ ,2.28
+ ,2.29
+ ,2.38
+ ,2.43
+ ,2.44
+ ,2.45
+ ,2.45
+ ,2.46
+ ,2.48
+ ,2.48
+ ,2.50
+ ,2.55
+ ,2.57
+ ,2.59
+ ,2.60
+ ,2.60
+ ,2.62
+ ,2.62
+ ,2.63
+ ,2.64
+ ,2.66
+ ,2.67
+ ,2.67
+ ,2.68
+ ,2.68
+ ,2.69
+ ,2.69
+ ,2.69
+ ,2.70
+ ,2.73
+ ,2.75
+ ,2.77
+ ,2.78
+ ,2.78
+ ,2.80
+ ,2.81
+ ,2.81
+ ,2.82
+ ,2.85
+ ,2.85
+ ,2.88
+ ,2.89
+ ,2.89
+ ,2.90
+ ,2.91
+ ,2.93
+ ,2.93
+ ,2.94
+ ,2.96
+ ,2.97
+ ,3.00
+ ,3.04
+ ,3.09
+ ,3.11
+ ,3.13
+ ,3.13
+ ,3.13
+ ,3.15
+ ,3.15
+ ,3.15
+ ,3.15
+ ,3.16
+ ,3.17
+ ,3.18
+ ,3.24
+ ,3.27
+ ,3.30
+ ,3.31
+ ,3.32
+ ,3.33
+ ,3.37
+ ,3.37
+ ,3.38
+ ,3.38
+ ,3.39
+ ,3.40
+ ,3.40
+ ,3.41
+ ,3.42
+ ,3.43
+ ,3.45
+ ,3.46
+ ,3.46
+ ,3.54
+ ,3.55
+ ,3.58
+ ,3.61
+ ,3.65
+ ,3.67
+ ,3.68
+ ,3.70
+ ,3.71
+ ,3.72
+ ,3.74
+ ,3.75
+ ,3.75
+ ,3.80
+ ,3.81
+ ,3.81
+ ,3.85
+ ,3.85
+ ,3.86
+ ,3.88
+ ,3.89
+ ,3.91
+ ,3.91
+ ,3.96
+ ,3.96
+ ,3.98
+ ,4.00
+ ,4.03
+ ,4.03
+ ,4.03
+ ,4.03
+ ,4.05
+ ,4.05
+ ,4.05
+ ,4.06
+ ,4.09
+ ,4.11
+ ,4.13
+ ,4.14
+ ,4.16
+ ,4.17
+ ,4.17
+ ,4.22
+ ,4.22
+ ,4.23
+ ,4.23
+ ,4.26
+ ,4.27
+ ,4.28
+ ,4.28
+ ,4.31
+ ,4.32
+ ,4.34
+ ,4.40
+ ,4.40
+ ,4.41
+ ,4.43
+ ,4.44
+ ,4.44
+ ,4.50
+ ,4.51
+ ,4.56
+ ,4.58
+ ,4.59
+ ,4.61
+ ,4.62
+ ,4.62
+ ,4.63
+ ,4.63
+ ,4.64
+ ,4.65
+ ,4.67
+ ,4.69
+ ,4.73
+ ,4.74
+ ,4.75
+ ,4.79
+ ,4.79
+ ,4.80
+ ,4.81
+ ,4.85
+ ,4.85
+ ,4.87
+ ,4.89
+ ,4.96
+ ,4.96
+ ,4.96
+ ,4.98
+ ,5.02
+ ,5.02
+ ,5.03
+ ,5.04
+ ,5.06
+ ,5.09
+ ,5.11
+ ,5.11
+ ,5.14
+ ,5.15
+ ,5.16
+ ,5.22
+ ,5.23
+ ,5.28
+ ,5.28
+ ,5.33
+ ,5.37
+ ,5.37
+ ,5.38
+ ,5.40
+ ,5.43
+ ,5.49
+ ,5.49
+ ,5.53
+ ,5.54
+ ,5.56
+ ,5.60
+ ,5.63
+ ,5.71
+ ,5.72
+ ,5.74
+ ,5.75
+ ,5.76
+ ,5.78
+ ,5.79
+ ,5.80
+ ,5.84
+ ,5.86
+ ,5.87
+ ,5.89
+ ,5.91
+ ,5.92
+ ,5.92
+ ,5.95
+ ,5.96
+ ,6.01
+ ,6.03
+ ,6.04
+ ,6.07
+ ,6.09
+ ,6.10
+ ,6.10
+ ,6.11
+ ,6.16
+ ,6.17
+ ,6.20
+ ,6.21
+ ,6.23
+ ,6.24
+ ,6.29
+ ,6.30
+ ,6.34
+ ,6.36
+ ,6.38
+ ,6.38
+ ,6.38
+ ,6.44
+ ,6.47
+ ,6.49
+ ,6.50
+ ,6.56
+ ,6.72
+ ,6.77
+ ,6.83
+ ,6.83
+ ,6.85
+ ,6.87
+ ,6.90
+ ,6.91
+ ,6.92
+ ,6.95
+ ,6.97
+ ,7.01
+ ,7.09
+ ,7.13
+ ,7.14
+ ,7.15
+ ,7.17
+ ,7.19
+ ,7.20
+ ,7.23
+ ,7.23
+ ,7.28
+ ,7.32
+ ,7.34
+ ,7.38
+ ,7.40
+ ,7.40
+ ,7.43
+ ,7.44
+ ,7.45
+ ,7.50
+ ,7.52
+ ,7.52
+ ,7.57
+ ,7.61
+ ,7.63
+ ,7.74
+ ,7.78
+ ,7.79
+ ,7.84
+ ,7.86
+ ,7.99
+ ,8.02
+ ,8.02
+ ,8.02
+ ,8.06
+ ,8.11
+ ,8.18
+ ,8.20
+ ,8.23
+ ,8.31
+ ,8.37
+ ,8.44
+ ,8.47
+ ,8.47
+ ,8.47
+ ,8.51
+ ,8.54
+ ,8.56
+ ,8.57
+ ,8.65
+ ,8.73
+ ,8.73
+ ,8.74
+ ,8.75
+ ,8.78
+ ,8.80
+ ,8.80
+ ,8.82
+ ,8.88
+ ,9.01
+ ,9.11
+ ,9.13
+ ,9.14
+ ,9.22
+ ,9.34
+ ,9.41
+ ,9.46
+ ,9.50
+ ,9.56
+ ,9.60
+ ,9.64
+ ,9.68
+ ,9.70
+ ,9.75
+ ,9.76
+ ,9.77
+ ,9.83
+ ,9.94
+ ,9.95
+ ,9.96
+ ,10.01
+ ,10.06
+ ,10.12
+ ,10.27
+ ,10.31
+ ,10.39
+ ,10.39
+ ,10.52
+ ,10.71
+ ,10.92
+ ,11.23
+ ,11.25
+ ,11.26
+ ,11.46
+ ,11.48
+ ,11.54
+ ,11.55
+ ,11.62
+ ,11.67
+ ,11.70
+ ,11.81
+ ,11.82
+ ,11.90
+ ,11.92
+ ,11.99
+ ,12.09
+ ,12.15
+ ,12.19
+ ,12.24
+ ,12.28
+ ,12.34
+ ,12.35
+ ,12.41
+ ,12.42
+ ,12.50
+ ,12.53
+ ,12.57
+ ,12.67
+ ,12.73
+ ,12.99
+ ,13.16
+ ,13.19
+ ,13.23
+ ,13.23
+ ,13.24
+ ,13.42
+ ,13.70
+ ,13.71
+ ,13.75
+ ,13.80
+ ,13.83
+ ,13.88
+ ,13.95
+ ,13.96
+ ,14.11
+ ,14.12
+ ,14.17
+ ,14.47
+ ,14.50
+ ,15.11
+ ,15.41
+ ,15.48
+ ,15.53
+ ,16.54
+ ,16.69
+ ,16.88
+ ,17.04
+ ,17.71
+ ,17.75
+ ,18.03
+ ,18.19
+ ,18.22
+ ,18.26
+ ,19.51
+ ,20.12
+ ,21.54
+ ,22.34
+ ,22.71
+ ,22.87
+ ,22.98
+ ,23.16
+ ,23.22
+ ,24.91
+ ,25.67
+ ,25.96
+ ,25.97
+ ,26.08
+ ,26.17
+ ,29.02
+ ,29.08
+ ,32.84
+ ,38.11
+ ,41.16
+ ,42.41
+ ,43.86
+ ,43.98
+ ,52.50
+ ,58.24
+ ,59.63
+ ,64.91
+ ,107.05
+ ,112.31
+ ,112.32)
> par3 = '512'
> par2 = 'yes'
> par1 = '0'
> ylab = 'Density'
> xlab = 'Value'
> par3 <- '512'
> par2 <- 'yes'
> par1 <- '0'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa, P. (2012), Kernel Density Estimation (v1.0.11) 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/1ds2k1416399142.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 (1002 obs.);	Bandwidth 'bw' = 1.27

       x                 y            
 Min.   :-37.201   Min.   :0.000e+00  
 1st Qu.:  1.132   1st Qu.:1.345e-05  
 Median : 39.465   Median :3.184e-04  
 Mean   : 39.465   Mean   :6.515e-03  
 3rd Qu.: 77.798   3rd Qu.:2.626e-03  
 Max.   :116.131   Max.   :7.389e-02  
> postscript(file="/var/wessaorg/rcomp/tmp/2wjv01416399142.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/3x3ze1416399142.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/4v51m1416399142.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/5o8sf1416399142.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/6mwnm1416399142.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/7xgzp1416399142.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")
> 
> 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/8xc6p1416399142.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Maximum Density Values',3,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Kernel',1,TRUE)
> a<-table.element(a,'x-value',1,TRUE)
> a<-table.element(a,'max. density',1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Gaussian',1,TRUE)
> a<-table.element(a,mydensity1$x[mydensity1$y==max(mydensity1$y)],1)
> a<-table.element(a,mydensity1$y[mydensity1$y==max(mydensity1$y)],1)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Epanechnikov',1,TRUE)
> a<-table.element(a,mydensity2$x[mydensity2$y==max(mydensity2$y)],1)
> a<-table.element(a,mydensity2$y[mydensity2$y==max(mydensity2$y)],1)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Rectangular',1,TRUE)
> a<-table.element(a,mydensity3$x[mydensity3$y==max(mydensity3$y)],1)
> a<-table.element(a,mydensity3$y[mydensity3$y==max(mydensity3$y)],1)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Triangular',1,TRUE)
> a<-table.element(a,mydensity4$x[mydensity4$y==max(mydensity4$y)],1)
> a<-table.element(a,mydensity4$y[mydensity4$y==max(mydensity4$y)],1)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Biweight',1,TRUE)
> a<-table.element(a,mydensity5$x[mydensity5$y==max(mydensity5$y)],1)
> a<-table.element(a,mydensity5$y[mydensity5$y==max(mydensity5$y)],1)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Cosine',1,TRUE)
> a<-table.element(a,mydensity6$x[mydensity6$y==max(mydensity6$y)],1)
> a<-table.element(a,mydensity6$y[mydensity6$y==max(mydensity6$y)],1)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Optcosine',1,TRUE)
> a<-table.element(a,mydensity7$x[mydensity7$y==max(mydensity7$y)],1)
> a<-table.element(a,mydensity7$y[mydensity7$y==max(mydensity7$y)],1)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/94btx1416399142.tab") 
> if (par2=='yes') {
+ 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)
+ 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)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/wessaorg/rcomp/tmp/10951f1416399142.tab") 
+ }
> 
> try(system("convert tmp/1ds2k1416399142.ps tmp/1ds2k1416399142.png",intern=TRUE))
character(0)
> try(system("convert tmp/2wjv01416399142.ps tmp/2wjv01416399142.png",intern=TRUE))
character(0)
> try(system("convert tmp/3x3ze1416399142.ps tmp/3x3ze1416399142.png",intern=TRUE))
character(0)
> try(system("convert tmp/4v51m1416399142.ps tmp/4v51m1416399142.png",intern=TRUE))
character(0)
> try(system("convert tmp/5o8sf1416399142.ps tmp/5o8sf1416399142.png",intern=TRUE))
character(0)
> try(system("convert tmp/6mwnm1416399142.ps tmp/6mwnm1416399142.png",intern=TRUE))
character(0)
> try(system("convert tmp/7xgzp1416399142.ps tmp/7xgzp1416399142.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
 15.092   0.515  15.639