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 computationTue, 22 Nov 2016 15:40: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/2016/Nov/22/t14798256710m5xp2x2nvsadow.htm/, Retrieved Sun, 05 May 2024 10:35:32 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=296954, Retrieved Sun, 05 May 2024 10:35:32 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact49
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kernel Density Estimation] [Kernel Density E...] [2016-11-22 14:40:31] [563c2945bc7c763925d38f2fb19cdb55] [Current]
Feedback Forum

Post a new message
Dataseries X:
7695
7730
7740
7755
7755
7765
7755
7725
7740
7745
7775
7815
7825
7850
7855
7905
7920
7950
7970
7980
7995
8000
8015
8010
8080
8100
8120
8145
8150
8175
8180
8165
8170
8170
8145
8145
8170
8180
8190
8220
8225
8235
8255
8250
8255
8265
8265
8280
8285
8290
8295
8310
8315
8325
8315
8325
8345
8355
8370
8450
8445
8400
8510
8565
8595
8605
8615
8670
8695
8710
8740
8800
8880
8915
8975
9000
9050
9105
9140
9155
9215
9245
9290
9225
9315
9355
9395
9495
9450
9455
9490
9505
9510
9535
9555
9590
9650
9665
9720
9760
9780
9785
9800
9825
9815
9815
9840
9795
9880
9905
9925
9960
9975
9990
9975





Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R ServerBig Analytics Cloud Computing Center
R Framework error message
Warning: there are blank lines in the 'Data' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time3 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
R Framework error message & 
Warning: there are blank lines in the 'Data' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.
\tabularnewline \hline \end{tabular} %Source: https://freestatistics.org/blog/index.php?pk=296954&T=0

[TABLE]
[ROW]
Summary of computational transaction[/C][/ROW] [ROW]Raw Input[/C] view raw input (R code) [/C][/ROW] [ROW]Raw Output[/C]view raw output of R engine [/C][/ROW] [ROW]Computing time[/C]3 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [ROW]R Framework error message[/C][C]
Warning: there are blank lines in the 'Data' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.
[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=296954&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=296954&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 Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R ServerBig Analytics Cloud Computing Center
R Framework error message
Warning: there are blank lines in the 'Data' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.







Properties of Density Trace
Bandwidth252.739333407582
#Observations115

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=296954&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
Bandwidth252.739333407582
#Observations115







Maximum Density Values
Kernelx-valuemax. density
Gaussian8174.93929150710.00064790280991076
Epanechnikov8152.562954909570.000627642239591485
Rectangular8130.186618312050.000632752695359211
Triangular8174.93929150710.000645935283743609
Biweight8152.562954909570.000632153455171087
Cosine8152.562954909570.000634133585762994
Optcosine8152.562954909570.000629189979257561

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 8174.9392915071 & 0.00064790280991076 \tabularnewline
Epanechnikov & 8152.56295490957 & 0.000627642239591485 \tabularnewline
Rectangular & 8130.18661831205 & 0.000632752695359211 \tabularnewline
Triangular & 8174.9392915071 & 0.000645935283743609 \tabularnewline
Biweight & 8152.56295490957 & 0.000632153455171087 \tabularnewline
Cosine & 8152.56295490957 & 0.000634133585762994 \tabularnewline
Optcosine & 8152.56295490957 & 0.000629189979257561 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=296954&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]8174.9392915071[/C][C]0.00064790280991076[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]8152.56295490957[/C][C]0.000627642239591485[/C][/ROW]
[ROW][C]Rectangular[/C][C]8130.18661831205[/C][C]0.000632752695359211[/C][/ROW]
[ROW][C]Triangular[/C][C]8174.9392915071[/C][C]0.000645935283743609[/C][/ROW]
[ROW][C]Biweight[/C][C]8152.56295490957[/C][C]0.000632153455171087[/C][/ROW]
[ROW][C]Cosine[/C][C]8152.56295490957[/C][C]0.000634133585762994[/C][/ROW]
[ROW][C]Optcosine[/C][C]8152.56295490957[/C][C]0.000629189979257561[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=296954&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=296954&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
Gaussian8174.93929150710.00064790280991076
Epanechnikov8152.562954909570.000627642239591485
Rectangular8130.186618312050.000632752695359211
Triangular8174.93929150710.000645935283743609
Biweight8152.562954909570.000632153455171087
Cosine8152.562954909570.000634133585762994
Optcosine8152.562954909570.000629189979257561







Kernel Density Values
x-valueGaussianEpanechnikovRectangularTriangularBiweightCosineOptcosine
Kernel Density Values are not shown

\begin{tabular}{lllllllll}
\hline
Kernel Density Values \tabularnewline
x-value & Gaussian & Epanechnikov & Rectangular & Triangular & Biweight & Cosine & Optcosine \tabularnewline
Kernel Density Values are not shown \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=296954&T=3

[TABLE]
[ROW][C]Kernel Density Values[/C][/ROW]
[ROW][C]x-value[/C][C]Gaussian[/C][C]Epanechnikov[/C][C]Rectangular[/C][C]Triangular[/C][C]Biweight[/C][C]Cosine[/C][C]Optcosine[/C][/ROW]
[ROW][C]Kernel Density Values are not shown[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=296954&T=3

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

As an alternative you can also use a QR Code:  

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

Kernel Density Values
x-valueGaussianEpanechnikovRectangularTriangularBiweightCosineOptcosine
Kernel Density Values are not shown



Parameters (Session):
par1 = 200 ; par2 = 5 ; par3 = 0 ; par4 = P1 P5 Q1 Q3 P95 P99 ;
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')
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(ab,file='mytable123.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(b,file='mytable2a.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='mytable1.tab')