Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
615 680 680 625 710 695 640 665 700 685 645 750 630 680 660 650 720 680 665 710 755 640 655 730 640 685 695 695 730 705 615 630 795 625 700 725 610 645 700 700 730 725 635 630 775 615 690 745 590 595 700 690 755 700 645 600 800 610 690 725 630 565 695 690 785 660 605 595 790 575 665 710 630 520 725 680 750 620 630 610 840 605 675 740 635 520 725 655 755 580 645 615 840 595 655 740 660 525 690 660 740 575 625 630 840 575 655 735
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') 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')
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