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 computationFri, 13 Nov 2009 07:50:04 -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/t1258123875d0vk9i34qcd24ll.htm/, Retrieved Sun, 05 May 2024 08:45:38 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=56736, Retrieved Sun, 05 May 2024 08:45:38 +0000
QR Codes:

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

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] [workshop 6] [2009-11-13 14:50:04] [6c94b261890ba36343a04d1029691995] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
111.632
106.707
108.827
108.413
106.249
104.861
102.382
100.320
100.228
117.089
121.523
114.948
112.831
107.605
108.928
101.993
102.850
99.925
101.536
99.450
98.305
110.159
109.483
106.810
96.279
91.982
90.276
90.999
86.622
83.117
80.367
77.550
77.443
92.844
92.175
84.822
81.632
78.872
81.485
80.651
78.192
76.844
76.335
71.415
73.889
86.822
86.371
83.469
82.662
82.880
89.406
95.378
97.657
100.247
99.180
97.493
101.628
114.585
115.669
111.311
Dataseries Y:
Download CSV
Histogram 
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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=56736&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'Gwilym Jenkins' @ 72.249.127.135






Bandwidth
x axis5.16486149133762
y axis2.90805122715653
Correlation
correlation used in KDE0.601167078652812
correlation(x,y)0.601167078652812

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=56736&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 axis5.16486149133762
y axis2.90805122715653
Correlation
correlation used in KDE0.601167078652812
correlation(x,y)0.601167078652812


Figures (Output of Computation)

Input Parameters & R Code

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