Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

rwasp_bidensity.wasp

Title produced by software

Bivariate Kernel Density Estimation

Date of computation

Fri, 06 Nov 2009 13:08:14 -0700

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2009/Nov/06/t1257538243ujomkoez56bq6pd.htm/, Retrieved Mon, 12 May 2025 20:13:13 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=54355, Retrieved Mon, 12 May 2025 20:13:13 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

277

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-     [Bivariate Kernel Density Estimation] [3/11/2009] [2009-11-02 21:54:51] [b98453cac15ba1066b407e146608df68]
-   PD    [Bivariate Kernel Density Estimation] [ws 6] [2009-11-06 20:08:14] [f7d3e79b917995ba1c8c80042fc22ef9] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

100
114,242116
115,8697864
108,0366226
100,3051882
98,1688708
97,5584944
105,4933876
104,8830112
105,4933876
114,0386572
88,40284842
96,64292981
113,733469
110,681587
111,190234
103,1536114
102,2380468
102,4415056
115,4628688
107,9348932
113,5300102
116,8870804
90,13224822
101,2207528
117,0905392
120,0406918
113,3265514
109,155646
107,121058
107,121058
119,9389624
112,1057986
114,3438454
119,5320448
94,60834181
105,2899288
118,311292
122,0752798
116,276704
106,5106816
111,698881
114,5473042
116,3784334
117,7009156
116,6836216
120,4476094
96,54120041
105,595117
117,0905392
115,6663276
105,7985758
95,93082401
94,09969481
94,81180061
106,5106816
95,62563581
99,7965412
104,4760936
83,82502543

Dataseries Y:

Download CSV

Histogram

100
100
93,5483871
88,17204301
89,24731183
91,39784946
92,47311828
91,39784946
88,17204301
87,09677419
84,94623656
92,47311828
93,5483871
93,5483871
91,39784946
90,32258065
91,39784946
93,5483871
93,5483871
92,47311828
91,39784946
89,24731183
86,02150538
88,17204301
87,09677419
87,09677419
86,02150538
84,94623656
84,94623656
86,02150538
86,02150538
84,94623656
86,02150538
82,79569892
77,41935484
80,64516129
78,49462366
75,2688172
75,2688172
75,2688172
77,41935484
78,49462366
76,34408602
73,11827957
68,8172043
65,59139785
69,89247312
82,79569892
84,94623656
80,64516129
74,19354839
70,96774194
74,19354839
82,79569892
86,02150538
86,02150538
82,79569892
78,49462366
79,56989247
87,09677419

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

\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 & 1 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=54355&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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=54355&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=54355&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 Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Bandwidth
x axis	3.58725546733050
y axis	3.00763238103525
Correlation
correlation used in KDE	-0.321072294081109
correlation(x,y)	-0.321072294081109

\begin{tabular}{lllllllll}
\hline
Bandwidth \tabularnewline
x axis & 3.58725546733050 \tabularnewline
y axis & 3.00763238103525 \tabularnewline
Correlation \tabularnewline
correlation used in KDE & -0.321072294081109 \tabularnewline
correlation(x,y) & -0.321072294081109 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=54355&T=1

[TABLE]
[ROW][C]Bandwidth[/C][/ROW]
[ROW][C]x axis[/C][C]3.58725546733050[/C][/ROW]
[ROW][C]y axis[/C][C]3.00763238103525[/C][/ROW]
[ROW][C]Correlation[/C][/ROW]
[ROW][C]correlation used in KDE[/C][C]-0.321072294081109[/C][/ROW]
[ROW][C]correlation(x,y)[/C][C]-0.321072294081109[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=54355&T=1

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

As an alternative you can also use a QR Code:

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

Bandwidth
x axis	3.58725546733050
y axis	3.00763238103525
Correlation
correlation used in KDE	-0.321072294081109
correlation(x,y)	-0.321072294081109

Figure 1

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = 50 ; par2 = 50 ; par3 = 0 ; par4 = 0 ; par5 = 0 ; par6 = Y ; par7 = Y ;

Parameters (R input):

par1 = 50 ; par2 = 50 ; par3 = 0 ; par4 = 0 ; par5 = 0 ; par6 = Y ; par7 = Y ;

R code (references can be found in the software module):

par1 <- as(par1,'numeric')
par2 <- as(par2,'numeric')
par3 <- as(par3,'numeric')
par4 <- as(par4,'numeric')
par5 <- as(par5,'numeric')
library('GenKern')
if (par3==0) par3 <- dpik(x)
if (par4==0) par4 <- dpik(y)
if (par5==0) par5 <- cor(x,y)
if (par1 > 500) par1 <- 500
if (par2 > 500) par2 <- 500
bitmap(file='bidensity.png')
op <- KernSur(x,y, xgridsize=par1, ygridsize=par2, correlation=par5, xbandwidth=par3, ybandwidth=par4)
image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE,main=main,xlab=xlab,ylab=ylab)
if (par6=='Y') contour(op$xords, op$yords, op$zden, add=TRUE)
if (par7=='Y') points(x,y)
(r<-lm(y ~ x))
abline(r)
box()
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Bandwidth',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'x axis',header=TRUE)
a<-table.element(a,par3)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'y axis',header=TRUE)
a<-table.element(a,par4)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Correlation',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'correlation used in KDE',header=TRUE)
a<-table.element(a,par5)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'correlation(x,y)',header=TRUE)
a<-table.element(a,cor(x,y))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code