FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_bidensity.wasp
Title produced by softwareBivariate Kernel Density Estimation
Date of computationFri, 05 Oct 2012 05:56:10 -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/2012/Oct/05/t1349431130momp9y1l21n26a8.htm/, Retrieved Fri, 03 May 2024 03:31:05 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=172173, Retrieved Fri, 03 May 2024 03:31:05 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Bivariate Kernel Density Estimation] [math scores] [2012-10-05 09:56:10] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
26
31
27
NA
NA
NA
NA
NA
NA
27
NA
23
NA
27
NA
NA
17
25
27
NA
NA
NA
26
NA
NA
23
27
NA
20
25
NA
NA
NA
NA
NA
24
NA
NA
NA
NA
NA
23
20
25
NA
26
26
20
21
NA
28
NA
26
NA
NA
NA
NA
23
NA
NA
30
12
35
25
NA
NA
NA
NA
NA
NA
NA
25
NA
18
NA
NA
NA
NA
NA
NA
20
24
NA
NA
26
NA
30
NA
23
NA
27
NA
26
13
18
27
NA
NA
NA
31
28
28
29
29
NA
NA
NA
23
NA
28
NA
NA
NA
NA
20
NA
NA
24
21
30
25
22
NA
23
NA
NA
NA
26
23
NA
NA
32
NA
26
18
NA
NA
NA
12
NA
NA
34
NA
33
NA
NA
28
26
NA
NA
24
24
25
23
NA
24
NA
NA
NA
22
NA
NA
24
30
24
NA
NA
NA
27
28
31
19
29
21
NA
NA
NA
29
25
17
NA
NA
NA
29
NA
NA
NA
NA
NA
NA
NA
27
NA
17
NA
27
NA
25
NA
NA
NA
NA
NA
14
NA
NA
NA
27
22
NA
NA
27
NA
NA
29
22
NA
26
NA
NA
25
22
23
NA
NA
NA
NA
21
NA
NA
25
31
NA
NA
27
28
23
NA
26
NA
NA
NA
NA
NA
NA
NA
23
13
NA
NA
22
25
NA
16
12
NA
NA
NA
NA
NA
NA
NA
NA
20
NA
NA
20
23
23
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
24
NA
NA
NA
NA
24
NA
NA
NA
NA
Dataseries Y:
Download CSV
Histogram 
NA
NA
NA
NA
NA
NA
NA
NA
0
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
19
NA
NA
NA
NA
NA
6
NA
NA
NA
12
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
14
19
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
8
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
8
NA
NA
NA
NA
NA
NA
NA
NA
NA
5
NA
1
16
12
18
NA
NA
10
9
14
NA
NA
4
NA
NA
1
NA
NA
NA
NA
11
NA
1
15
11
NA
NA
16
NA
NA
NA
NA
NA
NA
7
NA
NA
NA
NA
8
NA
NA
NA
NA
15
NA
NA
NA
NA
NA
9
17
NA
NA
NA
NA
NA
6
NA
5
NA
2
NA
13
NA
NA
NA
NA
NA
NA
13
NA
NA
18
8
18
15
14
9
3
2
15
8
2
15
15
6
0
20
16
17
5
16
17
13
9
10
18
12
11
2
3
2
10
9
4
NA
8
5
15
6
13
6
12
2
9
12
13
8
4
5
9
16
7
1
12
5
7
12
11
NA
7
16
15
2
11
18
NA
7
10
8
10
NA
18
3
1
15
5
13
14
11
11
8
16
NA
19
6
2
8
0
14
15
14
4
0
6
16
NA
9
17
14

Tables (Output of Computation)





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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=172173&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 time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ fisher.wessa.net






Bandwidth
x axis1.69592486315862
y axis2.23179844340726
Correlation
correlation used in KDE-0.0268514164841797
correlation(x,y)-0.0268514164841797

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

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

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

Bandwidth
x axis1.69592486315862
y axis2.23179844340726
Correlation
correlation used in KDE-0.0268514164841797
correlation(x,y)-0.0268514164841797


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 50 ; par2 = 50 ; par3 = 0 ; par4 = 0 ; par5 = 0 ; par6 = Y ; par7 = Y ; par8 = terrain.colors ;
Parameters (R input):
par1 = 50 ; par2 = 50 ; par3 = 0 ; par4 = 0 ; par5 = 0 ; par6 = Y ; par7 = Y ; par8 = terrain.colors ;
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')
x <- x[!is.na(y)]
y <- y[!is.na(y)]
y <- y[!is.na(x)]
x <- x[!is.na(x)]
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
if (par8 == 'terrain.colors') mycol <- terrain.colors(100)
if (par8 == 'rainbow') mycol <- rainbow(100)
if (par8 == 'heat.colors') mycol <- heat.colors(100)
if (par8 == 'topo.colors') mycol <- topo.colors(100)
if (par8 == 'cm.colors') mycol <- cm.colors(100)
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=mycol, 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')