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 computationThu, 18 Dec 2014 23:11:52 +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/2014/Dec/18/t1418944331eu9lghhth8r7n9k.htm/, Retrieved Sun, 19 May 2024 18:45:56 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=271320, Retrieved Sun, 19 May 2024 18:45:56 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Boxplot and Trimmed Means] [Reddy Moores Boxp...] [2010-10-12 16:37:57] [98fd0e87c3eb04e0cc2efde01dbafab6]
- R P   [Boxplot and Trimmed Means] [Reddy-Moores Plac...] [2010-10-13 09:46:26] [98fd0e87c3eb04e0cc2efde01dbafab6]
- RMPD    [Notched Boxplots] [] [2010-10-15 11:13:23] [b98453cac15ba1066b407e146608df68]
- RMPD      [Testing Mean with unknown Variance - Critical Value] [I1 hypo] [2014-12-18 22:46:45] [0f2441473984dd788bd8148394ed78a8]
- R  D        [Testing Mean with unknown Variance - Critical Value] [I2 hypo] [2014-12-18 22:51:24] [0f2441473984dd788bd8148394ed78a8]
- RM D          [Bivariate Kernel Density Estimation] [Ext motivatie kernal] [2014-12-18 23:03:23] [0f2441473984dd788bd8148394ed78a8]
- R  D            [Bivariate Kernel Density Estimation] [E2 num] [2014-12-18 23:08:44] [0f2441473984dd788bd8148394ed78a8]
-    D                [Bivariate Kernel Density Estimation] [I3 hypo] [2014-12-18 23:11:52] [8aa9b0b9e9cdf95f84c1d02ac9593640] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
21
26
22
22
18
23
12
20
22
21
19
22
15
20
19
18
20
21
15
16
23
21
18
25
9
23
16
16
19
25
25
18
23
21
10
22
26
23
23
24
24
23
15
16
19
18
27
13
28
23
21
19
19
18
19
17
25
19
26
14
28
16
20
24
23
22
21
25
27
23
23
18
18
23
19
15
20
16
25
25
19
19
16
19
19
23
21
22
19
20
3
23
14
23
20
15
13
16
7
24
17
24
24
19
28
23
19
23
25
25
20
16
20
25
25
23
17
20
16
23
12
24
11
14
23
18
29
16
19
16
23
19
4
20
20
4
24
16
3
24
23
17
20
22
19
24
19
27
22
23
22
17
23
23
28
29
21
24
20
7
19
28
26
19
20
23
24
16
19
24
21
16
16
21
NA
28
16
23
26
29
18
19
19
16
16
16
18
22
14
20
15
22
16
15
11
15
20
21
16
17
15
16
18
25
20
24
28
22
20
27
17
22
23
22
13
19
15
20
24
18
19
15
20
13
23
24
23
19
20
22
25
26
24
27
16
15
25
27
23
21
14
24
16
22
13
17
23
22
23
26
14
24
21
16
11
NA
19
16
19
16
11
23
27
23
25
24
22
26
19
19
19
20
16
Dataseries Y:
Download CSV
Histogram 
15
18
15
14
19
15
15
18
17
21
9
17
15
21
18
19
20
19
21
21
17
19
20
17
21
19
18
19
17
18
19
19
17
19
17
14
19
13
16
17
19
18
20
19
18
19
21
18
21
21
21
14
17
18
13
15
15
16
18
20
12
20
14
20
12
20
19
18
16
16
16
19
14
20
18
21
21
15
20
15
16
20
17
21
9
17
19
19
13
21
18
20
13
18
11
15
15
16
16
16
21
17
18
18
13
20
17
15
14
15
13
19
15
20
18
6
18
21
17
20
18
21
16
17
19
21
16
16
20
14
20
18
19
15
21
17
16
18
5
16
15
16
18
16
16
20
19
18
19
15
17
12
13
21
20
18
18
17
21
14
14
16
17
18
18
16
20
16
18
20
20
21
15
18
18
21
20
15
15
9
16
21
18
18
18
16
21
17
5
20
20
20
17
16
21
14
21
13
16
18
13
19
20
20
14
18
19
14
19
20
16
16
18
17
16
21
21
19
18
13
20
15
19
19
21
20
17
15
18
19
20
17
18
15
16
20
21
14
13
17
16
18
19
19
19
19
15
18
17
15
19
21
18
15
20
11
17
16
15
18
21
20
17
21
18
21
14
18
20
21
15
15
21

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gertrude Mary Cox' @ cox.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 & 1 seconds \tabularnewline
R Server & 'Gertrude Mary Cox' @ cox.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=271320&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]'Gertrude Mary Cox' @ cox.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=271320&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=271320&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 time1 seconds
R Server'Gertrude Mary Cox' @ cox.wessa.net






Bandwidth
x axis1.34582978817889
y axis0.695839646402925
Correlation
correlation used in KDE0.0487291233871223
correlation(x,y)0.0487291233871223

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=271320&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.34582978817889
y axis0.695839646402925
Correlation
correlation used in KDE0.0487291233871223
correlation(x,y)0.0487291233871223


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
Parameters (R input):
par1 = 30 ; par2 = 30 ; 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')