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 computationThu, 05 Feb 2015 18:46:32 +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/2015/Feb/05/t14231621082smuf4fkbvcuaqw.htm/, Retrieved Sat, 18 May 2024 18:02:41 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=276776, Retrieved Sat, 18 May 2024 18:02:41 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact79

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] [Dichtheidsgrafiek...] [2015-02-05 18:46:32] [9fb880b65d31df281cc201e5682ad8c7] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
1145
1057
1218
1146
1150
983
1013
960
925
1087
1063
1049
1153
972
1111
985
1005
820
976
982
960
842
1008
1086
1207
958
1040
800
886
906
892
908
1025
1108
1097
1074
981
920
1065
994
1021
864
864
837
920
1085
1048
1112
1080
1140
1159
1044
871
807
1110
1078
1079
1247
1136
1066
1073
976
1073
1144
1023
694
1052
1124
1104
1183
1320
1227

Tables (Output of Computation)





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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=276776&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 time3 seconds
R Server'Gwilym Jenkins' @ jenkins.wessa.net






Properties of Density Trace
Bandwidth42.4031055991634
#Observations72

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=276776&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
Bandwidth42.4031055991634
#Observations72






Maximum Density Values
Kernelx-valuemax. density
Gaussian1073.332910750310.00350814895377133
Epanechnikov1066.441179763260.00341506048371275
Rectangular1042.32012130860.00358841804003867
Triangular1071.609978003550.00351391097917905
Biweight1068.164112510020.00344577221945038
Cosine1068.164112510020.00345688934201392
Optcosine1066.441179763260.00343329741209571

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 1073.33291075031 & 0.00350814895377133 \tabularnewline
Epanechnikov & 1066.44117976326 & 0.00341506048371275 \tabularnewline
Rectangular & 1042.3201213086 & 0.00358841804003867 \tabularnewline
Triangular & 1071.60997800355 & 0.00351391097917905 \tabularnewline
Biweight & 1068.16411251002 & 0.00344577221945038 \tabularnewline
Cosine & 1068.16411251002 & 0.00345688934201392 \tabularnewline
Optcosine & 1066.44117976326 & 0.00343329741209571 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=276776&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]1073.33291075031[/C][C]0.00350814895377133[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]1066.44117976326[/C][C]0.00341506048371275[/C][/ROW]
[ROW][C]Rectangular[/C][C]1042.3201213086[/C][C]0.00358841804003867[/C][/ROW]
[ROW][C]Triangular[/C][C]1071.60997800355[/C][C]0.00351391097917905[/C][/ROW]
[ROW][C]Biweight[/C][C]1068.16411251002[/C][C]0.00344577221945038[/C][/ROW]
[ROW][C]Cosine[/C][C]1068.16411251002[/C][C]0.00345688934201392[/C][/ROW]
[ROW][C]Optcosine[/C][C]1066.44117976326[/C][C]0.00343329741209571[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=276776&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=276776&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
Gaussian1073.332910750310.00350814895377133
Epanechnikov1066.441179763260.00341506048371275
Rectangular1042.32012130860.00358841804003867
Triangular1071.609978003550.00351391097917905
Biweight1068.164112510020.00344577221945038
Cosine1068.164112510020.00345688934201392
Optcosine1066.441179763260.00343329741209571


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):
par3 <- '512'
par2 <- 'no'
par1 <- '0'
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')
}