Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
84.11 85.86 86.27 87.26 88.89 89.68 90.12 90.18 90.32 91.15 91.2 92.69 95.46 97.02 97.24 97.28 97.57 97.64 98.07 98.44 98.46 98.67 98.86 99.09 99.21 99.26 99.5 99.51 99.56 99.57 99.63 99.79 99.81 99.88 99.97 99.98 100.08 100.14 100.24 100.34 100.39 100.4 100.42 100.42 100.45 100.46 100.47 100.48 100.51 100.52 100.55 100.59 100.67 100.72 100.96 101.05 101.29 101.37 101.42 101.46 101.53 101.59 101.64 101.65 101.66 101.72 101.77 101.99 102.07 102.17 102.34 102.34 102.37 102.4 102.4 102.4 102.42 102.42 102.44 102.44 102.45 102.47 102.48 102.51 102.53 102.57 102.59 102.61 102.67 102.69 102.73 102.75 102.78 102.83 102.84 102.87 102.89 102.91 102.92 102.93 102.93 102.94 102.96 102.97 102.99 103.02 103.04 103.04 103.05 103.06 103.08 103.09 103.1 103.11 103.13 103.14 103.16 103.17 103.18 103.19 103.2 103.2 103.22 103.22 103.23 103.25 103.26 103.26 103.29 103.31 103.32 103.34 103.37 103.38 103.41 103.41 103.42 103.45 103.47 103.47 103.47 103.53 103.54 103.55 103.61 103.64 103.68 103.69 103.71 103.71 103.74 103.77 103.81 103.81 103.83 103.86 103.89 103.9 103.92 103.93 103.93 103.94 103.97 103.97 103.97 103.98 103.98 103.99 104 104.03 104.05 104.05 104.06 104.09 104.1 104.13 104.15 104.16 104.2 104.21 104.24 104.25 104.28 104.28 104.28 104.31 104.31 104.32 104.35 104.39 104.46 104.46 104.59 104.61 104.64 104.64 104.65 104.67 104.7 104.73 104.74 104.75 104.76 104.8 104.8 104.8 104.9 104.95 104.95 104.95 104.97 104.97 105 105 105 105.02 105.03 105.03 105.03 105.07 105.07 105.09 105.1 105.16 105.17 105.18 105.2 105.21 105.25 105.27 105.28 105.29 105.3 105.3 105.36 105.37 105.37 105.4 105.41 105.41 105.43 105.45 105.51 105.66 105.67 105.71 105.78 106.01 106.03 106.07 106.89 106.94 107 107.17 107.23 107.68 107.69 107.7 107.72 107.8 107.82 107.85 108.03 108.04 108.09 108.11 108.15 108.24 108.91 108.94 108.96 109.4 109.68 113.65 114.31 114.32 114.72 115.07 116.45 116.47 116.47 116.48 116.76 116.84 116.92 117.34 132.09 132.09 132.09 132.09 132.09 132.09 132.09 132.09 132.09 132.09 132.09 132.09 132.09
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