Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_density.wasp
Title produced by softwareKernel Density Estimation
Date of computationSun, 17 Aug 2014 17:15:31 +0100
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Aug/17/t1408292237g8huk3y1yyio3zg.htm/, Retrieved Thu, 31 Oct 2024 23:17:31 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=235617, Retrieved Thu, 31 Oct 2024 23:17:31 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsVan Reusel Raphaël
Estimated Impact127
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kernel Density Estimation] [Tijdreeks 1] [2014-08-17 16:15:31] [bf566d88435d8cc6ce5d208f6f8dd684] [Current]
Feedback Forum

Post a new message
Dataseries X:
1095
1085
1075
1054
1261
1250
1095
992
1002
1002
1013
1033
1095
1075
1106
1157
1447
1447
1385
1323
1374
1436
1447
1478
1571
1509
1509
1602
1860
1881
1829
1705
1798
1798
1808
1860
1901
1922
1922
1984
2222
2284
2294
2139
2222
2191
2129
2263
2294
2242
2253
2325
2594
2728
2728
2666
2759
2666
2614
2811
2842
2769
2955
3028
3245
3389
3369
3358
3441
3431
3307
3493
3555
3493
3751
3875
4164
4278
4247
4185
4237
4299
4092
4257
4361
4319
4588
4681
5074
5146
5053
5105
5136
5167
4970
5156
5259
5156
5456
5549
5952
6014
6034
6138
6138
6179
5993
6086
6148
6034
6365
6427
6840
6913
7016
7109
7119
7130
6944
7130




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Sir Maurice George Kendall' @ kendall.wessa.net

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 3 seconds \tabularnewline
R Server & 'Sir Maurice George Kendall' @ kendall.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235617&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Maurice George Kendall' @ kendall.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235617&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235617&T=0

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Sir Maurice George Kendall' @ kendall.wessa.net







Properties of Density Trace
Bandwidth652.451752346688
#Observations120

\begin{tabular}{lllllllll}
\hline
Properties of Density Trace \tabularnewline
Bandwidth & 652.451752346688 \tabularnewline
#Observations & 120 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235617&T=1

[TABLE]
[ROW][C]Properties of Density Trace[/C][/ROW]
[ROW][C]Bandwidth[/C][C]652.451752346688[/C][/ROW]
[ROW][C]#Observations[/C][C]120[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235617&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235617&T=1

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Properties of Density Trace
Bandwidth652.451752346688
#Observations120







Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular1926.520370298060.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 1788.81200709148 & 0.000231131863960289 \tabularnewline
Epanechnikov & 1867.50250035238 & 0.000227595691871377 \tabularnewline
Rectangular & 1926.52037029806 & 0.000235971581393028 \tabularnewline
Triangular & 1887.1751236676 & 0.000227733788516041 \tabularnewline
Biweight & 1887.1751236676 & 0.000228645319075803 \tabularnewline
Cosine & 1867.50250035238 & 0.000228910261980461 \tabularnewline
Optcosine & 1926.52037029806 & 0.000227597387157637 \tabularnewline
 \hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235617&T=2

[TABLE]
[ROW][C]Maximum Density Values[/C][/ROW]
[ROW][C]Kernel[/C][C]x-value[/C][C]max. density[/C][/ROW]
[ROW][C]Gaussian[/C][C]1788.81200709148[/C][C]0.000231131863960289[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]1867.50250035238[/C][C]0.000227595691871377[/C][/ROW]
[ROW][C]Rectangular[/C][C]1926.52037029806[/C][C]0.000235971581393028[/C][/ROW]
[ROW][C]Triangular[/C][C]1887.1751236676[/C][C]0.000227733788516041[/C][/ROW]
[ROW][C]Biweight[/C][C]1887.1751236676[/C][C]0.000228645319075803[/C][/ROW]
[ROW][C]Cosine[/C][C]1867.50250035238[/C][C]0.000228910261980461[/C][/ROW]
[ROW][C]Optcosine[/C][C]1926.52037029806[/C][C]0.000227597387157637[/C][/ROW]
 [/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235617&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235617&T=2

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular1926.520370298060.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637







Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular1985.538240243730.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 1788.81200709148 & 0.000231131863960289 \tabularnewline
Epanechnikov & 1867.50250035238 & 0.000227595691871377 \tabularnewline
Rectangular & 1985.53824024373 & 0.000235971581393028 \tabularnewline
Triangular & 1887.1751236676 & 0.000227733788516041 \tabularnewline
Biweight & 1887.1751236676 & 0.000228645319075803 \tabularnewline
Cosine & 1867.50250035238 & 0.000228910261980461 \tabularnewline
Optcosine & 1926.52037029806 & 0.000227597387157637 \tabularnewline
 \hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235617&T=3

[TABLE]
[ROW][C]Maximum Density Values[/C][/ROW]
[ROW][C]Kernel[/C][C]x-value[/C][C]max. density[/C][/ROW]
[ROW][C]Gaussian[/C][C]1788.81200709148[/C][C]0.000231131863960289[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]1867.50250035238[/C][C]0.000227595691871377[/C][/ROW]
[ROW][C]Rectangular[/C][C]1985.53824024373[/C][C]0.000235971581393028[/C][/ROW]
[ROW][C]Triangular[/C][C]1887.1751236676[/C][C]0.000227733788516041[/C][/ROW]
[ROW][C]Biweight[/C][C]1887.1751236676[/C][C]0.000228645319075803[/C][/ROW]
[ROW][C]Cosine[/C][C]1867.50250035238[/C][C]0.000228910261980461[/C][/ROW]
[ROW][C]Optcosine[/C][C]1926.52037029806[/C][C]0.000227597387157637[/C][/ROW]
 [/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235617&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235617&T=3

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular1985.538240243730.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637







Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular2005.210863558960.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 1788.81200709148 & 0.000231131863960289 \tabularnewline
Epanechnikov & 1867.50250035238 & 0.000227595691871377 \tabularnewline
Rectangular & 2005.21086355896 & 0.000235971581393028 \tabularnewline
Triangular & 1887.1751236676 & 0.000227733788516041 \tabularnewline
Biweight & 1887.1751236676 & 0.000228645319075803 \tabularnewline
Cosine & 1867.50250035238 & 0.000228910261980461 \tabularnewline
Optcosine & 1926.52037029806 & 0.000227597387157637 \tabularnewline
 \hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235617&T=4

[TABLE]
[ROW][C]Maximum Density Values[/C][/ROW]
[ROW][C]Kernel[/C][C]x-value[/C][C]max. density[/C][/ROW]
[ROW][C]Gaussian[/C][C]1788.81200709148[/C][C]0.000231131863960289[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]1867.50250035238[/C][C]0.000227595691871377[/C][/ROW]
[ROW][C]Rectangular[/C][C]2005.21086355896[/C][C]0.000235971581393028[/C][/ROW]
[ROW][C]Triangular[/C][C]1887.1751236676[/C][C]0.000227733788516041[/C][/ROW]
[ROW][C]Biweight[/C][C]1887.1751236676[/C][C]0.000228645319075803[/C][/ROW]
[ROW][C]Cosine[/C][C]1867.50250035238[/C][C]0.000228910261980461[/C][/ROW]
[ROW][C]Optcosine[/C][C]1926.52037029806[/C][C]0.000227597387157637[/C][/ROW]
 [/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235617&T=4

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235617&T=4

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular2005.210863558960.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637







Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular2024.883486874180.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 1788.81200709148 & 0.000231131863960289 \tabularnewline
Epanechnikov & 1867.50250035238 & 0.000227595691871377 \tabularnewline
Rectangular & 2024.88348687418 & 0.000235971581393028 \tabularnewline
Triangular & 1887.1751236676 & 0.000227733788516041 \tabularnewline
Biweight & 1887.1751236676 & 0.000228645319075803 \tabularnewline
Cosine & 1867.50250035238 & 0.000228910261980461 \tabularnewline
Optcosine & 1926.52037029806 & 0.000227597387157637 \tabularnewline
 \hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235617&T=5

[TABLE]
[ROW][C]Maximum Density Values[/C][/ROW]
[ROW][C]Kernel[/C][C]x-value[/C][C]max. density[/C][/ROW]
[ROW][C]Gaussian[/C][C]1788.81200709148[/C][C]0.000231131863960289[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]1867.50250035238[/C][C]0.000227595691871377[/C][/ROW]
[ROW][C]Rectangular[/C][C]2024.88348687418[/C][C]0.000235971581393028[/C][/ROW]
[ROW][C]Triangular[/C][C]1887.1751236676[/C][C]0.000227733788516041[/C][/ROW]
[ROW][C]Biweight[/C][C]1887.1751236676[/C][C]0.000228645319075803[/C][/ROW]
[ROW][C]Cosine[/C][C]1867.50250035238[/C][C]0.000228910261980461[/C][/ROW]
[ROW][C]Optcosine[/C][C]1926.52037029806[/C][C]0.000227597387157637[/C][/ROW]
 [/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235617&T=5

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235617&T=5

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular2024.883486874180.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637







Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular2044.556110189410.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 1788.81200709148 & 0.000231131863960289 \tabularnewline
Epanechnikov & 1867.50250035238 & 0.000227595691871377 \tabularnewline
Rectangular & 2044.55611018941 & 0.000235971581393028 \tabularnewline
Triangular & 1887.1751236676 & 0.000227733788516041 \tabularnewline
Biweight & 1887.1751236676 & 0.000228645319075803 \tabularnewline
Cosine & 1867.50250035238 & 0.000228910261980461 \tabularnewline
Optcosine & 1926.52037029806 & 0.000227597387157637 \tabularnewline
 \hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235617&T=6

[TABLE]
[ROW][C]Maximum Density Values[/C][/ROW]
[ROW][C]Kernel[/C][C]x-value[/C][C]max. density[/C][/ROW]
[ROW][C]Gaussian[/C][C]1788.81200709148[/C][C]0.000231131863960289[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]1867.50250035238[/C][C]0.000227595691871377[/C][/ROW]
[ROW][C]Rectangular[/C][C]2044.55611018941[/C][C]0.000235971581393028[/C][/ROW]
[ROW][C]Triangular[/C][C]1887.1751236676[/C][C]0.000227733788516041[/C][/ROW]
[ROW][C]Biweight[/C][C]1887.1751236676[/C][C]0.000228645319075803[/C][/ROW]
[ROW][C]Cosine[/C][C]1867.50250035238[/C][C]0.000228910261980461[/C][/ROW]
[ROW][C]Optcosine[/C][C]1926.52037029806[/C][C]0.000227597387157637[/C][/ROW]
 [/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235617&T=6

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235617&T=6

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular2044.556110189410.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637







Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular2064.228733504630.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 1788.81200709148 & 0.000231131863960289 \tabularnewline
Epanechnikov & 1867.50250035238 & 0.000227595691871377 \tabularnewline
Rectangular & 2064.22873350463 & 0.000235971581393028 \tabularnewline
Triangular & 1887.1751236676 & 0.000227733788516041 \tabularnewline
Biweight & 1887.1751236676 & 0.000228645319075803 \tabularnewline
Cosine & 1867.50250035238 & 0.000228910261980461 \tabularnewline
Optcosine & 1926.52037029806 & 0.000227597387157637 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235617&T=7

[TABLE]
[ROW][C]Maximum Density Values[/C][/ROW]
[ROW][C]Kernel[/C][C]x-value[/C][C]max. density[/C][/ROW]
[ROW][C]Gaussian[/C][C]1788.81200709148[/C][C]0.000231131863960289[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]1867.50250035238[/C][C]0.000227595691871377[/C][/ROW]
[ROW][C]Rectangular[/C][C]2064.22873350463[/C][C]0.000235971581393028[/C][/ROW]
[ROW][C]Triangular[/C][C]1887.1751236676[/C][C]0.000227733788516041[/C][/ROW]
[ROW][C]Biweight[/C][C]1887.1751236676[/C][C]0.000228645319075803[/C][/ROW]
[ROW][C]Cosine[/C][C]1867.50250035238[/C][C]0.000228910261980461[/C][/ROW]
[ROW][C]Optcosine[/C][C]1926.52037029806[/C][C]0.000227597387157637[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235617&T=7

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235617&T=7

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Maximum Density Values
Kernelx-valuemax. density
Gaussian1788.812007091480.000231131863960289
Epanechnikov1867.502500352380.000227595691871377
Rectangular2064.228733504630.000235971581393028
Triangular1887.17512366760.000227733788516041
Biweight1887.17512366760.000228645319075803
Cosine1867.502500352380.000228910261980461
Optcosine1926.520370298060.000227597387157637



Parameters (Session):
par1 = 0 ; par2 = no ; par3 = 512 ;
Parameters (R input):
par1 = 0 ; par2 = no ; par3 = 512 ;
R code (references can be found in the software module):
if (par1 == '0') bw <- 'nrd0'
if (par1 != '0') bw <- as.numeric(par1)
par3 <- as.numeric(par3)
mydensity <- array(NA, dim=c(par3,8))
bitmap(file='density1.png')
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()
mydensity1
bitmap(file='density2.png')
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()
bitmap(file='density3.png')
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()
bitmap(file='density4.png')
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()
bitmap(file='density5.png')
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()
bitmap(file='density6.png')
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()
bitmap(file='density7.png')
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()
load(file='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='mytable.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='mytable2.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='mytable1.tab')
}