FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_bidensity.wasp
Title produced by softwareBivariate Kernel Density Estimation
Date of computationSat, 06 Oct 2012 09:27:17 -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/06/t1349530191b5pmiz519fi6k68.htm/, Retrieved Tue, 30 Apr 2024 09:43:07 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=172649, Retrieved Tue, 30 Apr 2024 09:43:07 +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] [Bivariate Kernel ...] [2012-10-06 13:27:17] [4289cf790da1cc09a0cb8798de13fde3] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
41
39
30
31
34
35
39
34
36
37
38
36
38
39
33
32
36
38
39
32
32
31
39
37
39
41
36
33
33
34
31
27
37
34
34
32
29
36
29
35
37
34
38
35
38
37
38
33
36
38
32
32
32
34
32
37
39
29
37
35
30
38
34
31
34
35
36
30
39
35
38
31
34
38
34
39
37
34
28
37
33
37
35
37
32
33
38
33
29
33
31
36
35
32
29
39
37
35
37
32
38
37
36
32
33
40
38
41
36
43
Dataseries Y:
Download CSV
Histogram 
3
4
4
3
5
5
4
4
3
3
5
4
4
4
4
2
2
4
5
3
3
5
5
4
3
2
4
4
3
3
4
2
2
4
2
4
1
2
3
3
3
4
3
4
4
4
4
3
3
4
4
2
2
3
4
4
4
4
3
4
1
2
2
3
4
4
4
3
4
3
4
4
5
4
2
4
4
3
4
2
3
3
4
2
4
2
4
4
3
2
5
4
4
4
3
3
4
4
4
4
2
1
5
4
4
2
2
3
2
5

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 Maurice George Kendall' @ kendall.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 Maurice George Kendall' @ kendall.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=172649&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 Maurice George Kendall' @ kendall.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=172649&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=172649&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 Maurice George Kendall' @ kendall.wessa.net






Bandwidth
x axis1.23601999273031
y axis0.100117578839257
Correlation
correlation used in KDE0.111496000473474
correlation(x,y)0.111496000473474

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

[TABLE]
[ROW][C]Bandwidth[/C][/ROW]
[ROW][C]x axis[/C][C]1.23601999273031[/C][/ROW]
[ROW][C]y axis[/C][C]0.100117578839257[/C][/ROW]
[ROW][C]Correlation[/C][/ROW]
[ROW][C]correlation used in KDE[/C][C]0.111496000473474[/C][/ROW]
[ROW][C]correlation(x,y)[/C][C]0.111496000473474[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=172649&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=172649&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.23601999273031
y axis0.100117578839257
Correlation
correlation used in KDE0.111496000473474
correlation(x,y)0.111496000473474


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