Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_density.wasp
Title produced by softwareKernel Density Estimation
Date of computationThu, 18 Dec 2014 14:12:16 +0000
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/Dec/18/t1418911948m0rvmn6424ur5hu.htm/, Retrieved Sun, 19 May 2024 20:28:42 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=270962, Retrieved Sun, 19 May 2024 20:28:42 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact72
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kernel Density Estimation] [I1: Intrinsic mot...] [2014-12-18 14:12:16] [f02c6c9412fee5ce04bac553459224aa] [Current]
Feedback Forum

Post a new message
Dataseries X:
24
22
26
11
24
28
23
19
18
23
17
15
21
20
26
19
28
21
19
20
17
20
17
21
12
23
22
22
21
20
18
21
24
22
20
17
16
19
23
22
15
21
18
23
20
21
21
22
15
19
18
20
18
22
25
23
21
19
21
16
21
22
18
4
22
17
20
18
19
20
15
24
21
19
19
27
23
23
20
17
21
23
22
20
16
20
6
15
18
21
23
20
20
18
25
16
20
14
22
20
17
22
22
20
17
22
17
22
21
25
19
24
17
22
22
17
26
19
20
19
21
24
21
19
13
27
22
21
22
22
21
19
11
19
21
19
8
23
17
25
24
22
23
17
24
22
21
19
19
16
23
23
20
24
25
20
23
21
23
11
27
22
16
18
23
24
20
20
14
23
16
19
22
21
20
24
17
19
22
16
20
14
18
13
17
19
22
18
21
18
22
21
15
18
23
23
19
20
24
26
23
25
19
23
24
20
22
24
20
21
27
26
17
26
23
13
23
18
21
21
22
15
18
25
23
21
23
19
22
23
15
24
23
19
23
23
16
23
19
19
22
26
22
24
21
18
19
27
14
26
20
13
26
19
21
26
25
20
21
19
27
13
17
18
20
22
18
6
22
15
17
10
21
23
18
20
27
13
20
20
22
20
24
23
19
22
24
21
19
20
16
17
25
16
23
20
23
22
15
16
23
25
14
18
15
27
22
26
25
19
19
24
19
25
20
17
13
20
20
24
25
20
20
20
11
23
19
26
26
20
23
19
23
27
16
26
26
20
12
20
18
28
24
24
24
12
23
16
18
25
18
7
19
17
20
19
10
25
23
23
23
20
24
17
21
22
14
19
21
24
11
21
22
22
19
15
20
22
20
18
20




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'George Udny Yule' @ yule.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 & 'George Udny Yule' @ yule.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270962&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]'George Udny Yule' @ yule.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270962&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270962&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'George Udny Yule' @ yule.wessa.net







Properties of Density Trace
Bandwidth1.02311092051039
#Observations381

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270962&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
Bandwidth1.02311092051039
#Observations381







Maximum Density Values
Kernelx-valuemax. density
Gaussian20.21705398218970.112202739923564
Epanechnikov20.98379106986060.11301713725714
Rectangular21.33766972570870.127376739501418
Triangular20.04011465426570.115097042594063
Biweight20.86583151791120.112001257878203
Cosine20.68889218998720.111724280912834
Optcosine20.98379106986060.112994614193565

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 20.2170539821897 & 0.112202739923564 \tabularnewline
Epanechnikov & 20.9837910698606 & 0.11301713725714 \tabularnewline
Rectangular & 21.3376697257087 & 0.127376739501418 \tabularnewline
Triangular & 20.0401146542657 & 0.115097042594063 \tabularnewline
Biweight & 20.8658315179112 & 0.112001257878203 \tabularnewline
Cosine & 20.6888921899872 & 0.111724280912834 \tabularnewline
Optcosine & 20.9837910698606 & 0.112994614193565 \tabularnewline
 \hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270962&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]20.2170539821897[/C][C]0.112202739923564[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]20.9837910698606[/C][C]0.11301713725714[/C][/ROW]
[ROW][C]Rectangular[/C][C]21.3376697257087[/C][C]0.127376739501418[/C][/ROW]
[ROW][C]Triangular[/C][C]20.0401146542657[/C][C]0.115097042594063[/C][/ROW]
[ROW][C]Biweight[/C][C]20.8658315179112[/C][C]0.112001257878203[/C][/ROW]
[ROW][C]Cosine[/C][C]20.6888921899872[/C][C]0.111724280912834[/C][/ROW]
[ROW][C]Optcosine[/C][C]20.9837910698606[/C][C]0.112994614193565[/C][/ROW]
 [/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270962&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270962&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
Gaussian20.21705398218970.112202739923564
Epanechnikov20.98379106986060.11301713725714
Rectangular21.33766972570870.127376739501418
Triangular20.04011465426570.115097042594063
Biweight20.86583151791120.112001257878203
Cosine20.68889218998720.111724280912834
Optcosine20.98379106986060.112994614193565







Maximum Density Values
Kernelx-valuemax. density
Gaussian20.21705398218970.112202739923564
Epanechnikov20.98379106986060.11301713725714
Rectangular21.39664950168340.127376739501418
Triangular20.04011465426570.115097042594063
Biweight20.86583151791120.112001257878203
Cosine20.68889218998720.111724280912834
Optcosine20.98379106986060.112994614193565

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 20.2170539821897 & 0.112202739923564 \tabularnewline
Epanechnikov & 20.9837910698606 & 0.11301713725714 \tabularnewline
Rectangular & 21.3966495016834 & 0.127376739501418 \tabularnewline
Triangular & 20.0401146542657 & 0.115097042594063 \tabularnewline
Biweight & 20.8658315179112 & 0.112001257878203 \tabularnewline
Cosine & 20.6888921899872 & 0.111724280912834 \tabularnewline
Optcosine & 20.9837910698606 & 0.112994614193565 \tabularnewline
 \hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270962&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]20.2170539821897[/C][C]0.112202739923564[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]20.9837910698606[/C][C]0.11301713725714[/C][/ROW]
[ROW][C]Rectangular[/C][C]21.3966495016834[/C][C]0.127376739501418[/C][/ROW]
[ROW][C]Triangular[/C][C]20.0401146542657[/C][C]0.115097042594063[/C][/ROW]
[ROW][C]Biweight[/C][C]20.8658315179112[/C][C]0.112001257878203[/C][/ROW]
[ROW][C]Cosine[/C][C]20.6888921899872[/C][C]0.111724280912834[/C][/ROW]
[ROW][C]Optcosine[/C][C]20.9837910698606[/C][C]0.112994614193565[/C][/ROW]
 [/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270962&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270962&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
Gaussian20.21705398218970.112202739923564
Epanechnikov20.98379106986060.11301713725714
Rectangular21.39664950168340.127376739501418
Triangular20.04011465426570.115097042594063
Biweight20.86583151791120.112001257878203
Cosine20.68889218998720.111724280912834
Optcosine20.98379106986060.112994614193565







Maximum Density Values
Kernelx-valuemax. density
Gaussian20.21705398218970.112202739923564
Epanechnikov20.98379106986060.11301713725714
Rectangular21.45562927765810.127376739501418
Triangular20.04011465426570.115097042594063
Biweight20.86583151791120.112001257878203
Cosine20.68889218998720.111724280912834
Optcosine20.98379106986060.112994614193565

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 20.2170539821897 & 0.112202739923564 \tabularnewline
Epanechnikov & 20.9837910698606 & 0.11301713725714 \tabularnewline
Rectangular & 21.4556292776581 & 0.127376739501418 \tabularnewline
Triangular & 20.0401146542657 & 0.115097042594063 \tabularnewline
Biweight & 20.8658315179112 & 0.112001257878203 \tabularnewline
Cosine & 20.6888921899872 & 0.111724280912834 \tabularnewline
Optcosine & 20.9837910698606 & 0.112994614193565 \tabularnewline
 \hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270962&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]20.2170539821897[/C][C]0.112202739923564[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]20.9837910698606[/C][C]0.11301713725714[/C][/ROW]
[ROW][C]Rectangular[/C][C]21.4556292776581[/C][C]0.127376739501418[/C][/ROW]
[ROW][C]Triangular[/C][C]20.0401146542657[/C][C]0.115097042594063[/C][/ROW]
[ROW][C]Biweight[/C][C]20.8658315179112[/C][C]0.112001257878203[/C][/ROW]
[ROW][C]Cosine[/C][C]20.6888921899872[/C][C]0.111724280912834[/C][/ROW]
[ROW][C]Optcosine[/C][C]20.9837910698606[/C][C]0.112994614193565[/C][/ROW]
 [/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270962&T=4

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270962&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
Gaussian20.21705398218970.112202739923564
Epanechnikov20.98379106986060.11301713725714
Rectangular21.45562927765810.127376739501418
Triangular20.04011465426570.115097042594063
Biweight20.86583151791120.112001257878203
Cosine20.68889218998720.111724280912834
Optcosine20.98379106986060.112994614193565







Maximum Density Values
Kernelx-valuemax. density
Gaussian20.21705398218970.112202739923564
Epanechnikov20.98379106986060.11301713725714
Rectangular21.57358882960740.127376739501418
Triangular20.04011465426570.115097042594063
Biweight20.86583151791120.112001257878203
Cosine20.68889218998720.111724280912834
Optcosine20.98379106986060.112994614193565

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 20.2170539821897 & 0.112202739923564 \tabularnewline
Epanechnikov & 20.9837910698606 & 0.11301713725714 \tabularnewline
Rectangular & 21.5735888296074 & 0.127376739501418 \tabularnewline
Triangular & 20.0401146542657 & 0.115097042594063 \tabularnewline
Biweight & 20.8658315179112 & 0.112001257878203 \tabularnewline
Cosine & 20.6888921899872 & 0.111724280912834 \tabularnewline
Optcosine & 20.9837910698606 & 0.112994614193565 \tabularnewline
 \hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270962&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]20.2170539821897[/C][C]0.112202739923564[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]20.9837910698606[/C][C]0.11301713725714[/C][/ROW]
[ROW][C]Rectangular[/C][C]21.5735888296074[/C][C]0.127376739501418[/C][/ROW]
[ROW][C]Triangular[/C][C]20.0401146542657[/C][C]0.115097042594063[/C][/ROW]
[ROW][C]Biweight[/C][C]20.8658315179112[/C][C]0.112001257878203[/C][/ROW]
[ROW][C]Cosine[/C][C]20.6888921899872[/C][C]0.111724280912834[/C][/ROW]
[ROW][C]Optcosine[/C][C]20.9837910698606[/C][C]0.112994614193565[/C][/ROW]
 [/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270962&T=5

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270962&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
Gaussian20.21705398218970.112202739923564
Epanechnikov20.98379106986060.11301713725714
Rectangular21.57358882960740.127376739501418
Triangular20.04011465426570.115097042594063
Biweight20.86583151791120.112001257878203
Cosine20.68889218998720.111724280912834
Optcosine20.98379106986060.112994614193565







Maximum Density Values
Kernelx-valuemax. density
Gaussian20.21705398218970.112202739923564
Epanechnikov20.98379106986060.11301713725714
Rectangular21.63256860558210.127376739501418
Triangular20.04011465426570.115097042594063
Biweight20.86583151791120.112001257878203
Cosine20.68889218998720.111724280912834
Optcosine20.98379106986060.112994614193565

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 20.2170539821897 & 0.112202739923564 \tabularnewline
Epanechnikov & 20.9837910698606 & 0.11301713725714 \tabularnewline
Rectangular & 21.6325686055821 & 0.127376739501418 \tabularnewline
Triangular & 20.0401146542657 & 0.115097042594063 \tabularnewline
Biweight & 20.8658315179112 & 0.112001257878203 \tabularnewline
Cosine & 20.6888921899872 & 0.111724280912834 \tabularnewline
Optcosine & 20.9837910698606 & 0.112994614193565 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270962&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]20.2170539821897[/C][C]0.112202739923564[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]20.9837910698606[/C][C]0.11301713725714[/C][/ROW]
[ROW][C]Rectangular[/C][C]21.6325686055821[/C][C]0.127376739501418[/C][/ROW]
[ROW][C]Triangular[/C][C]20.0401146542657[/C][C]0.115097042594063[/C][/ROW]
[ROW][C]Biweight[/C][C]20.8658315179112[/C][C]0.112001257878203[/C][/ROW]
[ROW][C]Cosine[/C][C]20.6888921899872[/C][C]0.111724280912834[/C][/ROW]
[ROW][C]Optcosine[/C][C]20.9837910698606[/C][C]0.112994614193565[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270962&T=6

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270962&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
Gaussian20.21705398218970.112202739923564
Epanechnikov20.98379106986060.11301713725714
Rectangular21.63256860558210.127376739501418
Triangular20.04011465426570.115097042594063
Biweight20.86583151791120.112001257878203
Cosine20.68889218998720.111724280912834
Optcosine20.98379106986060.112994614193565



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')
}