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 18:25:23 +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/t1418927268drkt9hyelr2p7vp.htm/, Retrieved Fri, 17 May 2024 10:35:58 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=271192, Retrieved Fri, 17 May 2024 10:35:58 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Bootstrap Plot - Central Tendency] [bootstrap plot] [2014-12-15 10:33:55] [7b9119a46b6eb1a22eecca7bc054a6e2]
- RMPD    [Bivariate Kernel Density Estimation] [] [2014-12-18 18:25:23] [88f8137dd67cdcd531568536dca46410] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
139
158
128
224
105
159
167
165
159
119
163
153
148
188
149
244
150
132
161
105
97
166
111
145
162
163
109
148
125
116
149
138
148
159
164
176
202
214
188
110
48
54
50
124
121
221
150
153
154
94
156
151
157
194
159
39
100
187
105
111
162
99
186
183
138
101
177
126
101
139
114
114
162
111
75
82
159
158
67
121
32
117
165
147
165
150
154
126
149
145
109
120
172
132
169
172
114
156
167
2
113
165
165
118
173
158
155
49
220
122
151
44
141
152
103
107
175
154
110
143
131
167
137
121
168
149
94
51
168
145
140
109
66
119
164
132
126
83
142
93
117
166
Dataseries Y:
Download CSV
Histogram 
12,2
7,4
6,7
12,6
13,3
11,1
8,2
11,4
6,4
10,6
11,9
6,4
13,8
13,8
11,7
10,9
9,9
9
9,7
10,8
10,3
9,3
5,9
11,4
13
10,8
11,8
10,9
13,3
10,1
12,9
14,3
12,8
14,8
9,3
12
15,9
13
14,5
7,7
4,35
6,3
12,7
9,3
10
10,8
18,1
11,5
17,85
8,3
11,7
10,4
11,8
17,1
16,1
14,7
10,6
11,3
12,7
12,6
12,5
7,6
9,2
11,1
16,2
13,6
12,3
11,4
14,1
12,6
14,5
14,75
13
14,8
12,45
12,65
13,2
19,1
13,35
17,35
8,6
16,1
17,75
18,4
16,15
18,4
15,25
17,65
16,35
17,65
14,35
13,6
18,25
14,75
9,9
16,85
16
18,25
18,95
16,1
15,6
15,4
15,4
13,35
17,1
19,1
19,1
7,6
19,25
12,75
14,75
9,85
13,6
15,25
12,4
11,9
18,15
16,35
12,35
17,75
15,6
19,3
17,1
18,4
12,85
19,05
9,5
4,5
18,55
13,6
19,1
13,35
11,7
14,05
17,6
13,35
16,1
11,95
11,85
7,7
14,6
13,2

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'Herman Ole Andreas Wold' @ wold.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 & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=271192&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]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=271192&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=271192&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'Herman Ole Andreas Wold' @ wold.wessa.net






Bandwidth
x axis8.99269987073071
y axis1.17946936914141
Correlation
correlation used in KDE0.291725182582846
correlation(x,y)0.291725182582846

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=271192&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 axis8.99269987073071
y axis1.17946936914141
Correlation
correlation used in KDE0.291725182582846
correlation(x,y)0.291725182582846


Figures (Output of Computation)

Input Parameters & R Code

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