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 computationTue, 13 Dec 2011 10:22:14 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2011/Dec/13/t1323789757u91ewojre0z9uug.htm/, Retrieved Thu, 02 May 2024 20:23:01 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=154415, Retrieved Thu, 02 May 2024 20:23:01 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact149
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Bivariate Kernel Density Estimation] [] [2011-12-13 15:22:14] [d76b387543b13b5e3afd8ff9e5fdc89f] [Current]
Feedback Forum

Post a new message
Dataseries X:
188
140
188
209
234
272
400
408
314
236
188
184
173
126
165
185
224
268
375
429
348
285
214
178
212
145
158
237
275
281
356
467
366
277
200
212
248
178
227
276
304
362
488
539
469
385
277
249
262
217
248
335
321
418
573
632
547
449
329
282
322
244
334
375
459
497
615
838
611
465
363
319
341
271
325
372
433
473
645
793
657
461
426
357
412
266
337
415
476
606
745
925
700
585
397
395
401
284
362
451
445
567
735
846
649
606
404
435
410
303
470
456
506
605
873
909
755
635
469
554
486
398
482
536
616
702
940
1018
858
670
526
516
418
289
439
492
604
646
865
1013
812
665
611
532
507
377
496
604
657
667
908
1063
833
656
536
588
544
496
618
682
745
766
1172
1300
1044
816
652
652
569
468
651
668
817
858
1166
1413
1070
822
552
731
715
536
779
874
919
974
1396
1531
1235
1012
836
955
690
657
828
869
989
1163
1403
1552
1271
1099
839
931
754
593
818
864
961
1154
1436
1532
1261
1054
874
942
769
673
808
955
1075
1170
1424
1516
1272
1048
894
1067
814
668
931
848
1009
1207
1445
1516
1187
1064
893
917
818
637
860
1057
1095
1145
1533
1510
1241
1080
885
945
822
675
909
861
845
1109
1335
1389
1095
776
752
737
696
578
751
861
918
1106
1411
1399
1192
978
844
1003
762
616
773
763
933
1049
1331
1496
1214
1025
866
1028
784
686
759
983
1057
1192
1466
1672
1252
1157
983
1056
927
804
912
996
1222
1344
1525
1689
1260
1156
1032
1380
959
772
1026
1291
1314
1365
1925
1839
1637
1357
1140
1376
970
971
1059
1239
1266
1354
1820
1794
1620
1337
1076
1453
1049
996
1085
1332
1453
1447
1796
1906
1493
1458
1173
1138
1158
852
1211
1279
1274
1544
1687
2097
1634
1384
1172
1300
1041
910
980
1089
1649
1627
1976
1966
1440
1639
1283
1299
1008
941
1186
1293
1505
1715
1979
2022
1544
Dataseries Y:
117
116
166
180
202
290
298
441
388
260
175
105
137
142
176
231
240
316
363
537
487
324
185
133
169
157
206
244
243
393
405
579
525
373
198
148
201
177
222
275
290
402
534
614
578
419
203
173
229
192
294
310
365
509
537
655
643
444
259
229
276
245
324
323
349
480
530
676
670
476
281
240
259
237
400
367
497
593
696
969
878
581
373
232
358
318
410
480
604
713
844
1134
1013
755
371
280
417
417
514
548
583
839
924
1179
1109
896
452
337
484
524
575
622
664
926
1028
1361
1304
937
505
427
580
483
625
695
729
1099
1090
1393
1261
988
525
416
516
454
629
755
706
951
1099
1444
1316
1066
585
430
669
598
714
835
912
1031
1210
1581
1416
1120
652
505
741
675
782
956
996
1259
1389
1868
1609
1385
735
577
815
798
940
1007
1094
1413
1552
2038
1762
1411
805
729
912
753
989
1137
1256
1554
1629
2024
1900
1563
905
766
952
915
1197
1242
1197
1522
1591
2128
1962
1653
987
877
990
880
1258
1240
1312
1713
1683
2220
1996
1628
1119
890
1118
1164
1364
1412
1721
1752
1794
2434
2390
1929
1352
1060
1435
1196
1478
1648
1812
2118
2211
2826
2534
2290
1367
1105
1463
1299
1576
1850
1929
2367
2508
3073
2922
2377
1627
1259
1547
1436
1905
2079
1994
2501
2569
3467
2885
2211
1597
1141
1533
1546
1967
2171
2021
2753
2626
3532
3096
2639
1653
1425
1802
1674
1970
2092
2280
2715
2971
3937
3110
2662
1728
1609
1922
1863
1945
2365
2275
2962
2930
4062
3445
2943
1879
1694
2147
1999
2266
2562
2583
2965
3142
4115
3654
2992
2031
1699
2313
1970
2382
2830
2614
3321
3418
4468
3657
3250
2174
2014
2118
2227
2563
2817
2680
3337
3559
4608
3930
3133
2042
1999
2679
2425
2693
2760
2941
3611
3779
4945
4034
2906
2132
1932
2268
2178
2317
2552
2582
2886
3283
4125
3536
2568
1802
1598
2013
1872
2227
2497
2530
3119
3411
4511
3528
2833
1760
1517
1968
1809
2104
2391
2691
3023
3188
4057
3476




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'AstonUniversity' @ aston.wessa.net

\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 & 'AstonUniversity' @ aston.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=154415&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]'AstonUniversity' @ aston.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=154415&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=154415&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'AstonUniversity' @ aston.wessa.net







Bandwidth
x axis116.246431589274
y axis193.863574023945
Correlation
correlation used in KDE0.917000033865851
correlation(x,y)0.917000033865851

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=154415&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 axis116.246431589274
y axis193.863574023945
Correlation
correlation used in KDE0.917000033865851
correlation(x,y)0.917000033865851



Parameters (Session):
par1 = 50 ; par2 = 50 ; par3 = 0 ; par4 = 0 ; par5 = 0 ; par6 = Y ; par7 = Y ; par8 = rainbow ;
Parameters (R input):
par1 = 50 ; par2 = 50 ; par3 = 0 ; par4 = 0 ; par5 = 0 ; par6 = Y ; par7 = Y ; par8 = rainbow ;
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
if (par8 == 'terrain.colors') mycol <- terrain.colors(100)
if (par8 == 'rainbow') mycol <- rainbow(100)
if (par8 == 'heat.colors') mycol <- heat.colors(100)
if (par8 == 'topo.colors') mycol <- topo.colors(100)
if (par8 == 'cm.colors') mycol <- cm.colors(100)
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=mycol, 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')