Home » date » 2010 » Feb » 13 »

Max. prijs om naar show te gaan

*Unverified author*

R Software Module: /rwasp_density.wasp (opens new window with default values)

Title produced by software: Kernel Density Estimation

Date of computation: Sat, 13 Feb 2010 09:14:05 -0700

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68.htm/, Retrieved Mon, 20 May 2013 19:43:04 +0000

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

KDGP1W21

System-generated keywords (parent):

(pk = 0)

Estimated Impact

Dataseries X:

» Textfile « » CSV « » Stem and Leaf « » Histogram « » Kernel Density « » Harrell-Davis Quantiles « » Central Tendency « » Variability «

Output produced by software:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	3 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Properties of Density Trace
Bandwidth	0.861507799335009
#Observations	900

Charts produced by software:

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/13jxq1266077642.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/13jxq1266077642.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/2fsiv1266077642.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/2fsiv1266077642.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/3yhbd1266077642.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/3yhbd1266077642.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/44n121266077642.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/44n121266077642.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/5o8nu1266077642.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/5o8nu1266077642.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/6jur51266077642.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/6jur51266077642.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/7wkoy1266077642.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Feb/13/t1266077700lndmfyf58xnlf68/7wkoy1266077642.ps (opens in new window)

Click here to open pdf file.

Parameters (Session):

par1 = 0 ;

Parameters (R input):

par1 = 0 ;

R code (references can be found in the software module):

if (par1 == '0') bw <- 'nrd0'
if (par1 != '0') bw <- as.numeric(par1)
bitmap(file='density1.png')
mydensity1<-density(x,bw=bw,kernel='gaussian',na.rm=TRUE)
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)
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)
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)
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)
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)
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)
plot(mydensity7,main='Optcosine Kernel',xlab=xlab,ylab=ylab)
grid()
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Properties of Density Trace',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Bandwidth',header=TRUE)
a<-table.element(a,mydensity1$bw)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'#Observations',header=TRUE)
a<-table.element(a,mydensity1$n)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')