Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
350840 354950 342750 318610 303890 293480 286820 283650 276190 268670 275960 287040 291500 287540 272340 256570 246190 237340 235890 232790 226100 221110 221950 222990 232840 231380 221690 210190 202210 195210 194050 192360 187170 185930 188840 189350 194320 195850 188840 182380 177200 169140 171900 167320 164130 162550 166260 166010 177070 180800 179570 178130 178180 179430 179130 182220 183650 187110 197780 204000 221330 232500 237570 241100 244420 247240 254300 258970 262230 268880 277950 286630 299830 308090 311400 313520 310980 314910 310950 310720 310380 310570 314790 323400 335010 341630 340820 336280 325490 323750 317510 313890 308610 303720 303090 305140 304040 307100 304330 294710 286890 279050 271860 266710 259590 253830 250640 249140 250840 247590 237830 226380 217230 211420 207620 204310 197490 193580 192330 191970 196070 191940 185620 179410 173920 169190 166840 165170 161450 160830 163670 170830 182690 190940 197770 205090 210720 220210 229730 237070 241620 250370 258570 269860 283220 289610 281770 274700 267650 261380 260500 260730 254200 250450 253380 263740 276240 273820 265890 258400 253520 250710 252850 255260 251170 252500 257780 269900 291590 298870 295570 292100 290870 290580 297970 304010 304340 309850 322320 340170 369280 376690 379700 379520 377770 381560 394580 399320 400370 408200 419070 437730 460490 460070 453830 443300 435550
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