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 computationWed, 20 Jun 2012 08:04:15 -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/Jun/20/t1340193970d3oo5j5jkpub88u.htm/, Retrieved Sun, 28 Apr 2024 21:36:58 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=168750, Retrieved Sun, 28 Apr 2024 21:36:58 +0000
QR Codes:

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

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] [] [2012-06-20 12:04:15] [d76b387543b13b5e3afd8ff9e5fdc89f] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
182
110
162
155
141
165
101
152
127
128
121
112
142
164
156
161
176
152
155
85
164
155
88
145
122
161
150
132
171
118
96
114
140
139
184
121
139
140
146
148
112
146
126
169
75
83
163
180
181
168
94
114
107
152
186
118
133
90
127
87
121
103
50
134
89
84
163
50
98
96
123
104
122
124
128
123
76
85
121
88
116
137
66
136
159
102
110
104
107
158
126
83
48
97
63
131
93
97
105
88
89
95
70
81
107
84
129
77
134
84
58
69
93
166
40
102
97
75
67
81
146
103
48
63
147
70
64
77
146
103
63
126
41
65
30
35
56
37
30
20
49
8
22
21
23
12
13
16
18
1
12
8
4
4
0
7
0
0
0
0
0
0
0
0
Dataseries Y:
Download CSV
Histogram 
198
153
150
148
147
138
137
136
136
135
130
129
127
126
126
122
121
121
121
120
118
118
116
116
115
112
112
112
111
111
110
108
107
107
107
106
106
104
104
104
103
103
103
102
102
101
101
100
99
99
97
97
96
95
95
94
94
91
91
90
90
90
89
89
89
87
87
87
87
86
86
86
86
86
85
85
84
84
84
83
82
82
82
80
79
79
79
79
79
78
77
76
75
74
74
73
72
72
72
72
71
71
71
70
69
69
69
68
67
66
66
65
65
65
65
64
63
63
63
61
60
60
60
58
58
56
56
55
55
54
51
50
50
42
41
40
38
36
28
25
22
20
20
18
17
13
12
11
9
9
7
7
6
6
5
2
0
0
0
0
0
0
0
0

Tables (Output of Computation)





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

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






Bandwidth
x axis15.5395488336291
y axis9.79881379049444
Correlation
correlation used in KDE0.824414075319505
correlation(x,y)0.824414075319505

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=168750&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 axis15.5395488336291
y axis9.79881379049444
Correlation
correlation used in KDE0.824414075319505
correlation(x,y)0.824414075319505


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