Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
12.9 7.4 12.2 12.8 7.4 6.7 12.6 14.8 13.3 11.1 8.2 11.4 6.4 10.6 12 6.3 11.3 11.9 9.3 9.6 10 6.4 13.8 10.8 13.8 11.7 10.9 16.1 13.4 9.9 11.5 8.3 11.7 6.1 9 9.7 10.8 10.3 10.4 12.7 9.3 11.8 5.9 11.4 13 10.8 12.3 11.3 11.8 7.9 12.7 12.3 11.6 6.7 10.9 12.1 13.3 10.1 5.7 14.3 8 13.3 9.3 12.5 7.6 15.9 9.2 9.1 11.1 13 14.5 12.2 12.3 11.4 8.8 14.6 7.3 12.6 13 12.6 13.2 9.9 7.7 10.5 13.4 10.9 4.3 10.3 11.8 11.2 11.4 8.6 13.2 12.6 5.6 9.9 8.8 7.7 9 7.3 11.4 13.6 7.9 10.7 10.3 8.3 9.6 14.2 8.5 13.5 4.9 6.4 9.6 11.6 11.1 4.35 12.7 18.1 17.85 16.6 12.6 17.1 19.1 16.1 13.35 18.4 14.7 10.6 12.6 16.2 13.6 18.9 14.1 14.5 16.15 14.75 14.8 12.45 12.65 17.35 8.6 18.4 16.1 11.6 17.75 15.25 17.65 15.6 16.35 17.65 13.6 11.7 14.35 14.75 18.25 9.9 16 18.25 16.85 14.6 13.85 18.95 15.6 14.85 11.75 18.45 15.9 17.1 16.1 19.9 10.95 18.45 15.1 15 11.35 15.95 18.1 14.6 15.4 15.4 17.6 13.35 19.1 15.35 7.6 13.4 13.9 19.1 15.25 12.9 16.1 17.35 13.15 12.15 12.6 10.35 15.4 9.6 18.2 13.6 14.85 14.75 14.1 14.9 16.25 19.25 13.6 13.6 15.65 12.75 14.6 9.85 12.65 11.9 19.2 16.6 11.2 15.25 11.9 13.2 16.35 12.4 15.85 14.35 18.15 11.15 15.65 17.75 7.65 12.35 15.6 19.3 15.2 17.1 15.6 18.4 19.05 18.55 19.1 13.1 12.85 9.5 4.5 11.85 13.6 11.7 12.4 13.35 11.4 14.9 19.9 17.75 11.2 14.6 17.6 14.05 16.1 13.35 11.85 11.95 14.75 15.15 13.2 16.85 7.85 7.7 12.6 7.85 10.95 12.35 9.95 14.9 16.65 13.4 13.95 15.7 16.85 10.95 15.35 12.2 15.1 17.75 15.2 14.6 16.65 8.1
Sample Range:
(leave blank to include all observations)
From:
To:
Smoothing Bandwidth (leave zero for default)
Display values
(?)
no
no
yes
Number of density points
(?)
Chart options
Label y-axis:
Label x-axis:
R Code
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') 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') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Maximum Density Values',3,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Kernel',1,TRUE) a<-table.element(a,'x-value',1,TRUE) a<-table.element(a,'max. density',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Gaussian',1,TRUE) a<-table.element(a,mydensity1$x[mydensity1$y==max(mydensity1$y)],1) a<-table.element(a,mydensity1$y[mydensity1$y==max(mydensity1$y)],1) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Epanechnikov',1,TRUE) a<-table.element(a,mydensity2$x[mydensity2$y==max(mydensity2$y)],1) a<-table.element(a,mydensity2$y[mydensity2$y==max(mydensity2$y)],1) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Rectangular',1,TRUE) a<-table.element(a,mydensity3$x[mydensity3$y==max(mydensity3$y)],1) a<-table.element(a,mydensity3$y[mydensity3$y==max(mydensity3$y)],1) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Triangular',1,TRUE) a<-table.element(a,mydensity4$x[mydensity4$y==max(mydensity4$y)],1) a<-table.element(a,mydensity4$y[mydensity4$y==max(mydensity4$y)],1) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Biweight',1,TRUE) a<-table.element(a,mydensity5$x[mydensity5$y==max(mydensity5$y)],1) a<-table.element(a,mydensity5$y[mydensity5$y==max(mydensity5$y)],1) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Cosine',1,TRUE) a<-table.element(a,mydensity6$x[mydensity6$y==max(mydensity6$y)],1) a<-table.element(a,mydensity6$y[mydensity6$y==max(mydensity6$y)],1) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Optcosine',1,TRUE) a<-table.element(a,mydensity7$x[mydensity7$y==max(mydensity7$y)],1) a<-table.element(a,mydensity7$y[mydensity7$y==max(mydensity7$y)],1) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable2.tab') if (par2=='yes') { 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) 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) } a<-table.end(a) table.save(a,file='mytable1.tab') }
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation