Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
154 170 170 156 178 174 160 166 175 171 161 188 158 170 165 163 180 170 166 178 189 160 164 183 160 171 174 174 183 176 154 158 199 156 175 181 153 161 175 175 183 181 159 158 194 154 173 186 148 149 175 173 189 175 161 150 200 153 173 181 158 141 174 173 196 165 151 149 198 144 166 178 158 130 181 170 188 155 158 153 210 151 169 185 159 130 181 164 189 145 161 154 210 149 164 185 165 131 173 165 185 144 156 158 210 144 164 184
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation