Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_density.wasp
Title produced by softwareKernel Density Estimation
Date of computationTue, 04 Oct 2011 11:37:58 -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/2011/Oct/04/t1317742725iog7wqnrmu1j332.htm/, Retrieved Thu, 16 May 2024 12:50:43 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=126123, Retrieved Thu, 16 May 2024 12:50:43 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKDG2011W2EC
Estimated Impact71
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Data Series] [Tong] [2011-09-27 16:12:35] [74f6ba82782ee27f10a153e0e3d780e9]
- RMPD    [Kernel Density Estimation] [Gemiddelde consum...] [2011-10-04 15:37:58] [bd8cebb9d7961275d2f6ed94788b7e5f] [Current]
Feedback Forum

Post a new message
Dataseries X:
20.98
20.1
20.61
20.27
20.08
23.58
22.31
22.89
21.78
22.19
22.58
22.78
25.06
25.16
25.47
25.34
24.2
25.32
25.57
25.76
24.79
23.14
22.66
22.06
24.26
23.15
22.92
21.43
21.56
23.48
24.35
24.83
24.19
23.58
23.58
24.35
27.18
25.69
24.81
23.26
23.49
26.86
27.12
27.66
26.26
25.51
24.63
23.57
27.63
25.85
26.09
24.47
24.19
25.09
25.26
25.58
24.76
25.02
24.24
24.14
28.69
26.74
26.48
24.45
23.88
26.58
26.23
28.63
26.81
26.56
26.64
26.8
28.37
27.13
28.44
28.62
27.28
31.32
31.26
31.41
31.76
32.72
32.15
33.62
35.97
33.78
33.77
32.75
32.55
33.22
32.88
31.56
30.27
28.65
27.89
27.07
30.8
28.38
27.5
28
28.02
29.2
27.59
27.22
27.16
26.31
25.67
26.41
28.34
25.43
23.72
23.33
23.8
27.7
26.28
27.51
27.93
28.76
28.65
29.52
31.23
27.9
27.87
27.52
27.59
31.2
30.22
30.62
31.52
30.59
31.42
31.95






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
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 Output & view raw output of R engine  \tabularnewline
Computing time & 3 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ jenkins.wessa.net \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=126123&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]
[ROW][C]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=126123&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=126123&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'Gwilym Jenkins' @ jenkins.wessa.net
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
Bandwidth1.07627310854918
#Observations132

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

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

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



Parameters (Session):
par1 = 0 ;
Parameters (R input):
par1 = 0 ;
R code (references can be found in the software module):
if (par1 == '0') bw <- 'nrd0'
if (par1 != '0') bw <- as.numeric(par1)
bitmap(file='density1.png')
mydensity1<-density(x,bw=bw,kernel='gaussian',na.rm=TRUE)
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)
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)
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)
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)
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)
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)
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')