Free Statistics

of Irreproducible Research!

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 computationFri, 13 Nov 2009 08:00:59 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2009/Nov/13/t1258124531lfpm65llrbf8nxi.htm/, Retrieved Sun, 05 May 2024 10:38:29 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=56772, Retrieved Sun, 05 May 2024 10:38:29 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact108
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Bivariate Kernel Density Estimation] [workshop 6] [2009-11-13 15:00:59] [6c94b261890ba36343a04d1029691995] [Current]
Feedback Forum

Post a new message
Dataseries X:
123.560
122.117
121.782
121.789
122.273
121.683
119.869
118.873
117.607
122.783
124.454
127.064
125.166
124.554
124.272
128.836
127.408
126.420
124.465
124.526
124.379
130.189
132.196
134.893
132.709
129.955
127.947
130.369
129.852
124.278
126.141
121.743
110.898
117.707
120.738
121.445
120.439
116.313
117.173
119.773
119.639
113.006
113.776
107.866
106.924
114.562
115.367
116.602
114.393
115.140
117.623
119.361
120.527
121.660
122.852
119.325
119.151
126.494
127.832
128.780
Dataseries Y:
283.042
276.687
277.915
277.128
277.103
275.037
270.150
267.140
264.993
287.259
291.186
292.300
288.186
281.477
282.656
280.190
280.408
276.836
275.216
274.352
271.311
289.802
290.726
292.300
278.506
269.826
265.861
269.034
264.176
255.198
253.353
246.057
235.372
258.556
260.993
254.663
250.643
243.422
247.105
248.541
245.039
237.080
237.085
225.554
226.839
247.934
248.333
246.969
245.098
246.263
255.765
264.319
268.347
273.046
273.963
267.430
271.993
292.710
295.881
293.299




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\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 & 0 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=56772&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]0 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=56772&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=56772&T=0

As an alternative you can also use a 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 time0 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135







Bandwidth
x axis2.90805122715653
y axis7.42998839027608
Correlation
correlation used in KDE0.800558872376785
correlation(x,y)0.800558872376785

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=56772&T=1

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Bandwidth
x axis2.90805122715653
y axis7.42998839027608
Correlation
correlation used in KDE0.800558872376785
correlation(x,y)0.800558872376785



Parameters (Session):
par1 = 50 ; par2 = 50 ; par3 = 0 ; par4 = 0 ; par5 = 0 ; par6 = Y ; par7 = Y ;
Parameters (R input):
par1 = 50 ; par2 = 50 ; par3 = 0 ; par4 = 0 ; par5 = 0 ; par6 = Y ; par7 = Y ;
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')
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
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=terrain.colors(100), 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')