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:58:43 -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/t1258124362wfsdav7y2rtlio5.htm/, Retrieved Sun, 05 May 2024 11:33:55 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=56755, Retrieved Sun, 05 May 2024 11:33:55 +0000
QR Codes:

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

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:58:43] [6c94b261890ba36343a04d1029691995] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
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:
Download CSV
Histogram 
93.028
92.285
91.685
94.260
93.858
92.437
92.980
92.099
92.803
88.551
98.334
98.329
96.455
97.109
97.687
98.512
98.673
96.028
98.014
95.580
97.838
97.760
99.913
97.588
93.942
93.656
92.881
93.120
91.063
90.930
91.946
94.624
65.484
95.862
95.530
94.574
94.677
93.845
91.533
91.214
90.922
89.563
89.945
91.850
92.505
92.437
93.876
93.561
94.119
95.264
96.089
97.160
98.644
96.266
97.938
99.757
101.550
102.449
102.416
102.587

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=56755&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=56755&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=56755&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 axis2.90805122715653
y axis1.55441324953223
Correlation
correlation used in KDE0.442676662005938
correlation(x,y)0.442676662005938

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=56755&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 axis2.90805122715653
y axis1.55441324953223
Correlation
correlation used in KDE0.442676662005938
correlation(x,y)0.442676662005938


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