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, 12 Dec 2008 06:52:26 -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/2008/Dec/12/t1229090002l662a5xz1uuy2sb.htm/, Retrieved Sun, 19 May 2024 07:08:13 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=32735, Retrieved Sun, 19 May 2024 07:08:13 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact152
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Bivariate Kernel Density Estimation] [Bivariate KDP Tot...] [2008-12-12 13:52:26] [02e7fb326979b65614900650d62c19a6] [Current]
Feedback Forum

Post a new message
Dataseries X:
13807,9
14101,7
16010,3
14633,1
14478,5
15327,3
14179,5
11398,2
16111,5
15887,4
14529,3
13923,1
13960,2
14807,8
17511,5
15845,9
14594,2
17252,2
14832,8
13132,1
17665,9
16913
17318,8
16224,2
15469,6
16557,5
19414,8
17335
16525,2
18160,4
15553,8
15262,2
18581
17564,1
18948,6
17187,8
17564,8
17668,4
20811,7
17257,8
18984,2
20532,6
17082,3
16894,9
20274,9
20078,6
19900,9
17012,2
19642,9
19024
21691
18835,9
19873,4
21468,2
19406,8
18385,3
20739,3
22268,3
21569
17514,8
Dataseries Y:
2236
2084,9
2409,5
2199,3
2203,5
2254,1
1975,8
1742,2
2520,6
2438,1
2126,3
2267,5
2201,1
2128,5
2596
2458,2
2210,5
2621,2
2231,4
2103,6
2685,8
2539,3
2462,4
2693,3
2307,7
2385,9
2737,6
2653,9
2545,4
2848,8
2359,5
2488,3
2861,1
2717,9
2844
2749
2652,9
2660,2
3187,1
2774,1
3158,2
3244,6
2665,5
2820,8
2983,4
3077,4
3024,8
2731,8
3046,2
2834,8
3292,8
2946,1
3196,9
3284,2
3003
2979
3137,4
3630,2
3270,7
2942,3




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 2 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=32735&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]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=32735&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=32735&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 time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24







Bandwidth
x axis1256.65422399394
y axis203.379085726035
Correlation
correlation used in KDE0.961196026675653
correlation(x,y)0.961196026675653

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=32735&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 axis1256.65422399394
y axis203.379085726035
Correlation
correlation used in KDE0.961196026675653
correlation(x,y)0.961196026675653



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