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 14:48:33 +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/t1423147841qz63zwkv4cy3j8d.htm/, Retrieved Fri, 17 May 2024 23:55:38 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=276755, Retrieved Fri, 17 May 2024 23:55:38 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Data Series] [gemiddeld aantal ...] [2015-02-05 13:07:03] [55e6f8ecf0292dd7e776661106ae2b3f]
- RMP     [Kernel Density Estimation] [gemiddeld aantal ...] [2015-02-05 14:48:33] [d4d789bd6635ea6f10ab86fa6020c964] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232

Tables (Output of Computation)





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

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






Properties of Density Trace
Bandwidth14.5820867002872
#Observations153

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=276755&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
Bandwidth14.5820867002872
#Observations153






Maximum Density Values
Kernelx-valuemax. density
Gaussian156.2343371039160.0065423418457738
Epanechnikov167.4825180918630.00654251913498104
Rectangular115.9283552304360.00659872811808285
Triangular124.8331651792280.00654366864047559
Biweight143.1114592846430.00654236172087743
Cosine168.8885407153570.00654235300303239
Optcosine115.9283552304360.00654256107445054

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 156.234337103916 & 0.0065423418457738 \tabularnewline
Epanechnikov & 167.482518091863 & 0.00654251913498104 \tabularnewline
Rectangular & 115.928355230436 & 0.00659872811808285 \tabularnewline
Triangular & 124.833165179228 & 0.00654366864047559 \tabularnewline
Biweight & 143.111459284643 & 0.00654236172087743 \tabularnewline
Cosine & 168.888540715357 & 0.00654235300303239 \tabularnewline
Optcosine & 115.928355230436 & 0.00654256107445054 \tabularnewline
 \hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=276755&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]156.234337103916[/C][C]0.0065423418457738[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]167.482518091863[/C][C]0.00654251913498104[/C][/ROW]
[ROW][C]Rectangular[/C][C]115.928355230436[/C][C]0.00659872811808285[/C][/ROW]
[ROW][C]Triangular[/C][C]124.833165179228[/C][C]0.00654366864047559[/C][/ROW]
[ROW][C]Biweight[/C][C]143.111459284643[/C][C]0.00654236172087743[/C][/ROW]
[ROW][C]Cosine[/C][C]168.888540715357[/C][C]0.00654235300303239[/C][/ROW]
[ROW][C]Optcosine[/C][C]115.928355230436[/C][C]0.00654256107445054[/C][/ROW]
 [/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=276755&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=276755&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
Gaussian156.2343371039160.0065423418457738
Epanechnikov167.4825180918630.00654251913498104
Rectangular115.9283552304360.00659872811808285
Triangular124.8331651792280.00654366864047559
Biweight143.1114592846430.00654236172087743
Cosine168.8885407153570.00654235300303239
Optcosine115.9283552304360.00654256107445054






Maximum Density Values
Kernelx-valuemax. density
Gaussian156.2343371039160.0065423418457738
Epanechnikov167.4825180918630.00654251913498104
Rectangular115.9283552304360.00659872811808285
Triangular124.8331651792280.00654366864047559
Biweight168.8885407153570.00654236172087743
Cosine168.8885407153570.00654235300303239
Optcosine196.0716447695640.00654256107445054

\begin{tabular}{lllllllll}
\hline
Maximum Density Values \tabularnewline
Kernel & x-value & max. density \tabularnewline
Gaussian & 156.234337103916 & 0.0065423418457738 \tabularnewline
Epanechnikov & 167.482518091863 & 0.00654251913498104 \tabularnewline
Rectangular & 115.928355230436 & 0.00659872811808285 \tabularnewline
Triangular & 124.833165179228 & 0.00654366864047559 \tabularnewline
Biweight & 168.888540715357 & 0.00654236172087743 \tabularnewline
Cosine & 168.888540715357 & 0.00654235300303239 \tabularnewline
Optcosine & 196.071644769564 & 0.00654256107445054 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=276755&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]156.234337103916[/C][C]0.0065423418457738[/C][/ROW]
[ROW][C]Epanechnikov[/C][C]167.482518091863[/C][C]0.00654251913498104[/C][/ROW]
[ROW][C]Rectangular[/C][C]115.928355230436[/C][C]0.00659872811808285[/C][/ROW]
[ROW][C]Triangular[/C][C]124.833165179228[/C][C]0.00654366864047559[/C][/ROW]
[ROW][C]Biweight[/C][C]168.888540715357[/C][C]0.00654236172087743[/C][/ROW]
[ROW][C]Cosine[/C][C]168.888540715357[/C][C]0.00654235300303239[/C][/ROW]
[ROW][C]Optcosine[/C][C]196.071644769564[/C][C]0.00654256107445054[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=276755&T=3

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

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
Gaussian156.2343371039160.0065423418457738
Epanechnikov167.4825180918630.00654251913498104
Rectangular115.9283552304360.00659872811808285
Triangular124.8331651792280.00654366864047559
Biweight168.8885407153570.00654236172087743
Cosine168.8885407153570.00654235300303239
Optcosine196.0716447695640.00654256107445054


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