FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_density.wasp
Title produced by softwareKernel Density Estimation
Date of computationSun, 18 Aug 2013 07:08:46 -0400
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2013/Aug/18/t1376824157q82t72j5m40lc39.htm/, Retrieved Mon, 06 May 2024 06:20:15 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=211164, Retrieved Mon, 06 May 2024 06:20:15 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsDe Laere Dieter
Estimated Impact120

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kernel Density Estimation] [Tijdreeks 2 - Stap 4] [2013-08-18 11:08:46] [bc2cf5f41ec5ca561b7a550898b8dd0d] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
4640
4880
4400
4120
4440
4640
4680
4360
4640
4840
5000
4800
4720
4840
3800
4280
4480
4880
4680
4480
4720
5000
4960
4920
4480
5320
3960
4440
4360
4840
4880
4880
4400
4800
5280
4720
4440
5200
4240
4520
4640
5040
4840
4760
4520
4680
5480
4680
4160
5360
4200
4520
4600
4880
4840
4600
4520
4600
5760
4640
4520
5400
4200
4600
4480
4680
4400
4480
4840
4680
5480
4680
4440
5280
4240
4600
4640
4920
4560
4400
5080
4640
5520
4600
4720
5480
4320
4640
4920
4840
4520
4440
5000
4840
5480
4320
4880
5440
4480
4600
4720
5000
4160
4720
5000
4480
5720
4600

Tables (Output of Computation)





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

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

As an alternative you can also use a QR Code:  
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 time7 seconds
R Server'Gwilym Jenkins' @ jenkins.wessa.net






Properties of Density Trace
Bandwidth105.320509479822
#Observations108

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

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

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

As an alternative you can also use a QR Code:  
QR Code


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

Properties of Density Trace
Bandwidth105.320509479822
#Observations108






Maximum Density Values
Kernelx-valuemax. density
Gaussian4620.22392115130.00129620091491391
Epanechnikov4574.573612908810.00127435393042797
Rectangular4544.140074080490.00134998523354853
Triangular4640.512947036850.00130025893009424
Biweight4599.934895265750.00127726333670997
Cosine4610.079408208520.00128060590514053
Optcosine4579.64586938020.00127204758328764

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 4620.2239211513 & 0.00129620091491391 \tabularnewline
Epanechnikov & 4574.57361290881 & 0.00127435393042797 \tabularnewline
Rectangular & 4544.14007408049 & 0.00134998523354853 \tabularnewline
Triangular & 4640.51294703685 & 0.00130025893009424 \tabularnewline
Biweight & 4599.93489526575 & 0.00127726333670997 \tabularnewline
Cosine & 4610.07940820852 & 0.00128060590514053 \tabularnewline
Optcosine & 4579.6458693802 & 0.00127204758328764 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=211164&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]4620.2239211513[/C][C]0.00129620091491391[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]4574.57361290881[/C][C]0.00127435393042797[/C][/ROW]
[ROW][C]Rectangular[/C][C]4544.14007408049[/C][C]0.00134998523354853[/C][/ROW]
[ROW][C]Triangular[/C][C]4640.51294703685[/C][C]0.00130025893009424[/C][/ROW]
[ROW][C]Biweight[/C][C]4599.93489526575[/C][C]0.00127726333670997[/C][/ROW]
[ROW][C]Cosine[/C][C]4610.07940820852[/C][C]0.00128060590514053[/C][/ROW]
[ROW][C]Optcosine[/C][C]4579.6458693802[/C][C]0.00127204758328764[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=211164&T=2

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

As an alternative you can also use a QR Code:  
QR Code


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

Maximum Density Values
Kernelx-valuemax. density
Gaussian4620.22392115130.00129620091491391
Epanechnikov4574.573612908810.00127435393042797
Rectangular4544.140074080490.00134998523354853
Triangular4640.512947036850.00130025893009424
Biweight4599.934895265750.00127726333670997
Cosine4610.079408208520.00128060590514053
Optcosine4579.64586938020.00127204758328764


Figures (Output of Computation)

Input Parameters & R Code

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