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, 30 Oct 2009 08:12:02 -0600
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/Oct/30/t1256911990x1nizup98swgnnj.htm/, Retrieved Mon, 29 Apr 2024 03:19:34 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=52106, Retrieved Mon, 29 Apr 2024 03:19:34 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact109
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Bivariate Kernel Density Estimation] [workshop 4 review] [2009-10-30 14:12:02] [6198946fb53eb5eb18db46bb758f7fde] [Current]
Feedback Forum

Post a new message
Dataseries X:
12.28411643
12.27179952
12.25231194
12.22609584
12.20725694
12.22423668
12.39452537
12.46961698
12.40214917
12.40680698
12.3612585
12.37673061
12.3683743
12.34737254
12.33351554
12.30898116
12.2883068
12.29888881
12.4546734
12.49066871
12.42369699
12.39302513
12.34254658
12.35365178
12.29809105
12.29825521
12.26721478
12.27732762
12.26028617
12.26177695
12.38859418
12.3910407
12.34835315
12.24910819
12.19562829
12.17868089
12.17860381
12.13181235
12.09013416
12.05877841
12.02755584
12.03203468
12.21750702
12.21794193
12.11459312
12.06379454
12.0224142
12.05115417
12.04179334
12.00873085
11.99482061
11.98050115
11.92504174
11.97107157
12.13634736
12.13923281
12.06870533
12.03979371
12.02471263
12.08974667
12.13571456
12.15069735
12.16297202
12.14900482
12.13352883
12.1823382
12.31730986
12.33226026
12.27431991
Dataseries Y:
13.23983525
13.2372159
13.22732141
13.20596712
13.19312375
13.2043525
13.29384304
13.32246427
13.32548869
13.32338238
13.2949157
13.29707942
13.28934284
13.28682472
13.27838887
13.2588155
13.24891849
13.25168492
13.33865954
13.3517021
13.3506648
13.32467872
13.29699545
13.29990855
13.29363747
13.28799984
13.27043677
13.26074175
13.25824825
13.25853977
13.33677613
13.34680049
13.33733926
13.28384427
13.24589342
13.23081232
13.23672011
13.21558775
13.18378392
13.17290249
13.14419923
13.11968379
13.22737543
13.24385685
13.20238824
13.17508865
13.141864
13.15048036
13.15564727
13.13737218
13.10749742
13.10265647
13.05911895
13.07648696
13.17756911
13.18738303
13.15762434
13.13463576
13.12608843
13.15413526
13.17727193
13.1855973
13.19249001
13.19288708
13.16847758
13.19228861
13.28204729
13.29899209
13.2716847




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'George Udny Yule' @ 72.249.76.132

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=52106&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 time3 seconds
R Server'George Udny Yule' @ 72.249.76.132







Bandwidth
x axis0.0626501896862142
y axis0.0332127254092337
Correlation
correlation used in KDE0.96802082210868
correlation(x,y)0.96802082210868

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=52106&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 axis0.0626501896862142
y axis0.0332127254092337
Correlation
correlation used in KDE0.96802082210868
correlation(x,y)0.96802082210868



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