Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_density.wasp
Title produced by softwareKernel Density Estimation
Date of computationThu, 27 Feb 2020 14:29:17 +0100
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2020/Feb/27/t1582810209p2iw5jscptckg18.htm/, Retrieved Wed, 21 Apr 2021 10:56:33 +0200
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=, Retrieved Wed, 21 Apr 2021 10:56:33 +0200
QR Codes:

Original text written by user:
IsPrivate?This computation is private
User-defined keywords
Estimated Impact0
Dataseries X:
Latitude	Longitude
46.98321	17.97522
46.97892	17.96649
46.97825	17.96777
46.98387	17.97988
46.98588	17.98400
46.97958	17.96915
46.98295	17.98080
46.98081	17.96923
46.97616	17.96609
46.98616	17.98382
46.98354	17.97956
46.97862	17.96679
46.97854	17.96676
46.98555	17.98307
46.98774	17.98656
46.98303	17.97294
46.98357	17.97247
46.98552	17.98348
46.98274	17.97411
46.97820	17.96686
47.00554	18.00576
46.98356	17.97569
46.97813	17.96697
46.98613	17.98380
46.99067	17.99670
46.97609	17.96685
46.98254	17.97796
46.98656	17.98382
46.98420	17.97996
46.98764	17.98693
46.99811	17.99815
46.98419	17.97988
46.98423	17.98012
47.00915	18.00625
47.01085	18.00608
46.97887	17.96663
46.98509	17.98604
46.98551	17.98783
46.99852	17.99824
47.01072	18.00714
46.98432	17.98019
46.98798	17.98674
46.99740	17.99708
46.98421	17.97992
46.98472	17.98110
47.00867	18.00816
46.98432	17.98007
46.98428	17.98003
46.97790	17.96715
46.97801	17.96719
46.98414	17.97988
46.98412	17.97986
47.00948	18.00565
46.99739	17.99753
46.98399	17.97987
46.98404	17.98001
46.98060	17.97005
46.97792	17.96905
46.98506	17.98321
46.98782	17.98653
46.99820	17.99809
46.98416	17.97994
46.98361	17.97992
46.97810	17.96784
46.98413	17.98003
46.98463	17.98090
46.98167	17.97035
46.97775	17.96790
46.98601	17.98408
46.98697	17.98744
46.99255	17.99230
46.97614	17.96630
46.98416	17.97988
46.98342	17.98012
46.97825	17.96808
46.97611	17.96636
46.98407	17.97993
46.98293	17.98006
46.97917	17.96748
46.97798	17.96738
46.98414	17.98005
46.98358	17.97986
46.98712	17.98722
46.98417	17.97993
46.98244	17.97053
46.98234	17.97072
46.98177	17.97098
46.98276	17.97051
46.98239	17.97089
46.98191	17.97258
46.98429	17.98016
46.98262	17.97062
46.98245	17.97076
46.98180	17.97074
46.98265	17.97065
46.98256	17.97050
46.98157	17.97093
46.98163	17.97157
46.97640	17.96638
46.98409	17.98018
46.98258	17.97051
46.98764	17.98650
46.98430	17.98003
46.98415	17.98018
46.98248	17.97085
46.98180	17.97081
46.98412	17.97997
46.98250	17.97086
46.98190	17.97091
46.98391	17.97889
46.98411	17.98041
46.98240	17.97073
46.98239	17.97074
46.98426	17.98016
46.98234	17.97083
46.98165	17.97063
46.98173	17.97091
46.98346	17.97621
46.98272	17.97029
46.97749	17.96677
46.97787	17.96625
46.98418	17.97994
46.98263	17.97058
46.98243	17.97076
46.98180	17.97127
46.97888	17.96668
46.98247	17.97075
46.98225	17.97074
46.98413	17.98003
46.98312	17.97973
46.98245	17.97082
46.97405	17.96243
46.98415	17.98004
46.98256	17.97092
46.98116	17.96941
46.98420	17.98006
46.98164	17.97448
46.98245	17.97106
46.98262	17.97080
46.98387	17.98002
46.98246	17.97079
46.98254	17.97054
46.97913	17.96647
46.98213	17.97079
46.98754	17.98631
46.98764	17.98665
46.98512	17.98433
46.97931	17.96686
46.98253	17.97036
46.98193	17.97115
46.98186	17.97094
46.98760	17.98650
46.98767	17.98665
46.98799	17.98663
46.98353	17.97844
46.98402	17.98048
46.98186	17.97102
46.98202	17.97080
46.98440	17.98004
46.98767	17.98683
46.98747	17.98661
46.98773	17.98656
46.98248	17.97483
46.98758	17.98647
46.98359	17.97985
46.97911	17.96644
46.98043	17.97027
46.98224	17.97034
46.98409	17.97994
46.98790	17.98685
46.98861	17.98706
46.98868	17.98746
46.98422	17.97979
46.98853	17.98699
46.98870	17.98688
46.98574	17.98436
46.98431	17.98008
46.98874	17.98701
46.97889	17.96716
46.98266	17.97112
46.98875	17.98703
46.98878	17.98705
46.98903	17.98674
46.98431	17.98007
46.98766	17.98665
46.98899	17.98710
46.98926	17.98668
46.98905	17.98662
46.98282	17.97446
46.97893	17.96713
46.98272	17.97086
46.98246	17.97077
46.98232	17.97091
46.98247	17.97084
46.98271	17.97068
46.98427	17.98008
46.98260	17.97112
46.98267	17.97110
46.98244	17.97059
46.98188	17.97488
46.98256	17.97114
46.97880	17.96702
46.98248	17.97013
46.98419	17.98016
46.98366	17.97906
46.98417	17.97970
46.98459	17.98190
46.98505	17.98197
46.98360	17.97983
46.98543	17.98357
46.98725	17.98644
46.98432	17.98035
46.98580	17.98420
46.98605	17.98382
46.98426	17.98037
46.98541	17.98407
46.98728	17.98624
46.98710	17.98633
46.98546	17.98488
46.98807	17.98708
46.98729	17.98626
46.98691	17.98607
46.98426	17.98029
46.98178	17.97086
46.98661	17.98578
46.98703	17.98602
46.98269	17.97061
46.98247	17.97040
46.98171	17.97121
46.98853	17.98676
46.98860	17.98695
46.98250	17.97062
46.98299	17.97038
46.98877	17.98653
46.98858	17.98700
46.98893	17.98668
46.98922	17.98691
46.98830	17.98658
46.98762	17.98662
46.98757	17.98684
46.98835	17.98657
46.98818	17.98650
46.98755	17.98679
46.98762	17.98673
46.98757	17.98663
46.98451	17.98078
46.98534	17.98276
46.98711	17.98609
46.98702	17.98610
46.98842	17.98680
46.98617	17.98378
46.98576	17.98415
46.98464	17.98291
46.98884	17.98674
46.98428	17.98020
46.97992	17.96998
46.98183	17.97085
46.98282	17.97029
46.98428	17.98011
46.98773	17.98672
46.98763	17.98687
46.98692	17.98597
46.98271	17.97084
46.98200	17.97090
46.97872	17.96675
46.98304	17.96993
46.98274	17.97115
46.98859	17.98824
46.98512	17.98288
46.98220	17.97065
46.98259	17.97099
46.98866	17.98802
46.98846	17.98696
46.98871	17.98635
46.98766	17.98677
46.98763	17.98675
46.98734	17.98630
46.98701	17.98589
46.98777	17.98683
46.98715	17.98620
46.98890	17.98674
46.99114	17.98900
46.99122	17.98900
46.98253	17.97137
46.97773	17.96690
46.98310	17.97476
46.98350	17.97588
46.97863	17.96697
46.98774	17.98676
46.98769	17.98677
46.98848	17.98702
46.98799	17.98672
46.98810	17.98698
46.99825	17.99710
46.98787	17.98676
46.98189	17.97077
46.98208	17.97236
46.98197	17.97084
46.98281	17.97423
46.98362	17.97994
46.97970	17.96821
46.98183	17.97099
46.98609	17.98390
46.98311	17.97679
46.98443	17.98103
46.97589	17.94757
46.98220	17.97095
46.98176	17.97080
46.98615	17.98391
46.98308	17.97679
46.98423	17.98004
46.98333	17.97946
46.98458	17.98132
46.97894	17.96730
46.98365	17.98007
46.98081	17.96953
46.97904	17.96735
46.97886	17.96658
46.98327	17.97974
46.98374	17.98035
46.98560	17.98391
46.98604	17.98409
46.98411	17.97995
46.98290	17.97653
46.98426	17.98062
46.98376	17.98016
46.98596	17.98390
46.98462	17.98154
46.98303	17.97660
46.97900	17.96898
46.98417	17.98004
46.98578	17.98394
46.98502	17.98268
46.98145	17.97027
46.98420	17.98011
46.98381	17.98005
46.98469	17.98111
46.97800	17.96769
46.98103	17.97480
46.98452	17.98150
46.98015	17.96941
46.98427	17.98022
46.98360	17.98012
46.98382	17.98000
46.98348	17.97997
46.98422	17.98008
46.98420	17.97992
46.97870	17.96678
46.98408	17.98006
46.98420	17.98015
46.97616	17.96609
46.97621	17.96615
46.98457	17.98085
46.98353	17.97990
46.97492	17.96527
46.98428	17.98015
46.98071	17.97302
46.97619	17.96610
46.97713	17.96810
46.98569	17.98352
46.98415	17.98013
46.97584	17.96725
46.97507	17.96669
46.98431	17.97978
46.97469	17.96519
46.98419	17.97997
46.97612	17.96716
46.97466	17.96519
46.98412	17.97998
46.97601	17.96632
46.98329	17.97573
46.98417	17.98006
46.98202	17.97232
46.98322	17.97471
46.98317	17.97476
46.98345	17.97978
46.98529	17.98343
46.98308	17.97698
46.98357	17.97876
46.98534	17.98379
46.98556	17.98706
46.98346	17.98746
46.98429	17.98021
46.98319	17.98027
46.98332	17.98029




Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R ServerBig Analytics Cloud Computing Center

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time1 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=&T=0

[TABLE]
[ROW]
Summary of computational transaction[/C][/ROW] [ROW]Raw Input[/C] view raw input (R code) [/C][/ROW] [ROW]Raw Output[/C]view raw output of R engine [/C][/ROW] [ROW]Computing time[/C]1 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=&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 Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R ServerBig Analytics Cloud Computing Center



Parameters (Session):
par1 = 0 ; par2 = no ; par3 = 512 ;
Parameters (R input):
par1 = 0 ; par2 = no ; par3 = 512 ;
R code (references can be found in the software module):
par3 <- '512'
par2 <- 'no'
par1 <- '0'
if (par1 == '0') bw <- 'nrd0'
if (par1 != '0') bw <- as.numeric(par1)
par3 <- as.numeric(par3)
mydensity <- array(NA, dim=c(par3,8))
bitmap(file='density1.png')
mydensity1<-density(x,bw=bw,kernel='gaussian',na.rm=TRUE)
mydensity[,8] = signif(mydensity1$x,3)
mydensity[,1] = signif(mydensity1$y,3)
plot(mydensity1,main='Gaussian Kernel',xlab=xlab,ylab=ylab)
grid()
dev.off()
mydensity1
bitmap(file='density2.png')
mydensity2<-density(x,bw=bw,kernel='epanechnikov',na.rm=TRUE)
mydensity[,2] = signif(mydensity2$y,3)
plot(mydensity2,main='Epanechnikov Kernel',xlab=xlab,ylab=ylab)
grid()
dev.off()
bitmap(file='density3.png')
mydensity3<-density(x,bw=bw,kernel='rectangular',na.rm=TRUE)
mydensity[,3] = signif(mydensity3$y,3)
plot(mydensity3,main='Rectangular Kernel',xlab=xlab,ylab=ylab)
grid()
dev.off()
bitmap(file='density4.png')
mydensity4<-density(x,bw=bw,kernel='triangular',na.rm=TRUE)
mydensity[,4] = signif(mydensity4$y,3)
plot(mydensity4,main='Triangular Kernel',xlab=xlab,ylab=ylab)
grid()
dev.off()
bitmap(file='density5.png')
mydensity5<-density(x,bw=bw,kernel='biweight',na.rm=TRUE)
mydensity[,5] = signif(mydensity5$y,3)
plot(mydensity5,main='Biweight Kernel',xlab=xlab,ylab=ylab)
grid()
dev.off()
bitmap(file='density6.png')
mydensity6<-density(x,bw=bw,kernel='cosine',na.rm=TRUE)
mydensity[,6] = signif(mydensity6$y,3)
plot(mydensity6,main='Cosine Kernel',xlab=xlab,ylab=ylab)
grid()
dev.off()
bitmap(file='density7.png')
mydensity7<-density(x,bw=bw,kernel='optcosine',na.rm=TRUE)
mydensity[,7] = signif(mydensity7$y,3)
plot(mydensity7,main='Optcosine Kernel',xlab=xlab,ylab=ylab)
grid()
dev.off()
load(file='createtable')
ab<-table.start()
ab<-table.row.start(ab)
ab<-table.element(ab,'Properties of Density Trace',2,TRUE)
ab<-table.row.end(ab)
ab<-table.row.start(ab)
ab<-table.element(ab,'Bandwidth',header=TRUE)
ab<-table.element(ab,mydensity1$bw)
ab<-table.row.end(ab)
ab<-table.row.start(ab)
ab<-table.element(ab,'#Observations',header=TRUE)
ab<-table.element(ab,mydensity1$n)
ab<-table.row.end(ab)
ab<-table.end(ab)
a <- ab
table.save(ab,file='mytable123.tab')
b<-table.start()
b<-table.row.start(b)
b<-table.element(b,'Maximum Density Values',3,TRUE)
b<-table.row.end(b)
b<-table.row.start(b)
b<-table.element(b,'Kernel',1,TRUE)
b<-table.element(b,'x-value',1,TRUE)
b<-table.element(b,'max. density',1,TRUE)
b<-table.row.end(b)
b<-table.row.start(b)
b<-table.element(b,'Gaussian',1,TRUE)
b<-table.element(b,mydensity1$x[mydensity1$y==max(mydensity1$y)],1)
b<-table.element(b,mydensity1$y[mydensity1$y==max(mydensity1$y)],1)
b<-table.row.end(b)
b<-table.row.start(b)
b<-table.element(b,'Epanechnikov',1,TRUE)
b<-table.element(b,mydensity2$x[mydensity2$y==max(mydensity2$y)],1)
b<-table.element(b,mydensity2$y[mydensity2$y==max(mydensity2$y)],1)
b<-table.row.end(b)
b<-table.row.start(b)
b<-table.element(b,'Rectangular',1,TRUE)
b<-table.element(b,mydensity3$x[mydensity3$y==max(mydensity3$y)],1)
b<-table.element(b,mydensity3$y[mydensity3$y==max(mydensity3$y)],1)
b<-table.row.end(b)
b<-table.row.start(b)
b<-table.element(b,'Triangular',1,TRUE)
b<-table.element(b,mydensity4$x[mydensity4$y==max(mydensity4$y)],1)
b<-table.element(b,mydensity4$y[mydensity4$y==max(mydensity4$y)],1)
b<-table.row.end(b)
b<-table.row.start(b)
b<-table.element(b,'Biweight',1,TRUE)
b<-table.element(b,mydensity5$x[mydensity5$y==max(mydensity5$y)],1)
b<-table.element(b,mydensity5$y[mydensity5$y==max(mydensity5$y)],1)
b<-table.row.end(b)
b<-table.row.start(b)
b<-table.element(b,'Cosine',1,TRUE)
b<-table.element(b,mydensity6$x[mydensity6$y==max(mydensity6$y)],1)
b<-table.element(b,mydensity6$y[mydensity6$y==max(mydensity6$y)],1)
b<-table.row.end(b)
b<-table.row.start(b)
b<-table.element(b,'Optcosine',1,TRUE)
b<-table.element(b,mydensity7$x[mydensity7$y==max(mydensity7$y)],1)
b<-table.element(b,mydensity7$y[mydensity7$y==max(mydensity7$y)],1)
b<-table.row.end(b)
b<-table.end(b)
a <- b[1]
table.save(b,file='mytable2a.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Kernel Density Values',8,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'x-value',1,TRUE)
a<-table.element(a,'Gaussian',1,TRUE)
a<-table.element(a,'Epanechnikov',1,TRUE)
a<-table.element(a,'Rectangular',1,TRUE)
a<-table.element(a,'Triangular',1,TRUE)
a<-table.element(a,'Biweight',1,TRUE)
a<-table.element(a,'Cosine',1,TRUE)
a<-table.element(a,'Optcosine',1,TRUE)
a<-table.row.end(a)
if (par2=='yes') {
for(i in 1:par3) {
a<-table.row.start(a)
a<-table.element(a,mydensity[i,8],1,TRUE)
for(j in 1:7) {
a<-table.element(a,mydensity[i,j],1)
}
a<-table.row.end(a)
}
} else {
a<-table.row.start(a)
a<-table.element(a,'Kernel Density Values are not shown',8)
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')