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, 23 Dec 2009 08:59:17 -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/Dec/23/t1261584043yrrmwebrh4yohtd.htm/, Retrieved Mon, 29 Apr 2024 15:03:28 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=70547, Retrieved Mon, 29 Apr 2024 15:03:28 +0000
QR Codes:

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

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] [3/11/2009] [2009-11-02 21:54:51] [b98453cac15ba1066b407e146608df68]
-   PD    [Bivariate Kernel Density Estimation] [Workshop 6: Bivar...] [2009-12-23 15:59:17] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
106.5789
101.3158
98.68421
100
102.6316
102.6316
102.6316
98.68421
98.68421
93.42105
98.68421
98.68421
100
101.3158
101.3158
103.9474
106.5789
107.8947
107.8947
107.8947
103.9474
96.05263
90.78947
86.84211
88.15789
90.78947
92.10526
93.42105
94.73684
93.42105
90.78947
92.10526
89.47368
84.21053
88.15789
86.84211
84.21053
82.89474
81.57895
85.52632
89.47368
89.47368
84.21053
80.26316
76.31579
80.26316
94.73684
96.05263
90.78947
80.26316
76.31579
81.57895
93.42105
101.3158
103.9474
101.3158
97.36842
98.68421
105.2632
106.5789
Dataseries Y:
Download CSV
Histogram 
93.76979937
89.96832101
86.16684266
87.43400211
91.23548046
93.76979937
92.50263992
87.43400211
87.43400211
86.16684266
89.96832101
91.23548046
89.96832101
88.70116156
87.43400211
89.96832101
92.50263992
95.03695882
95.03695882
95.03695882
92.50263992
88.70116156
84.89968321
82.36536431
82.36536431
82.36536431
83.63252376
86.16684266
87.43400211
87.43400211
86.16684266
86.16684266
82.36536431
77.2967265
77.2967265
74.7624076
72.2280887
74.7624076
74.7624076
77.2967265
79.83104541
78.56388596
74.7624076
72.2280887
68.42661035
70.96092925
78.56388596
79.83104541
76.02956705
70.96092925
69.6937698
74.7624076
82.36536431
86.16684266
86.16684266
82.36536431
78.56388596
78.56388596
82.36536431
84.89968321

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=70547&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 axis3.87900697481608
y axis3.10460982700658
Correlation
correlation used in KDE0.875094091504901
correlation(x,y)0.875094091504901

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=70547&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 axis3.87900697481608
y axis3.10460982700658
Correlation
correlation used in KDE0.875094091504901
correlation(x,y)0.875094091504901


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