R version 2.12.1 (2010-12-16) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- c(193.230,199.068,195.076,191.563,191.067,186.665,185.508,184.371,183.046,175.714,175.768,171.029,170.465,170.102,156.389,124.291,99.360,86.675,85.056,128.236,164.257,162.401,152.779,156.005,153.387,153.190,148.840,144.211,145.953,145.542,150.271,147.489,143.824,134.754,131.736,126.304,125.511,125.495,130.133,126.257,110.323,98.417,105.749,120.665,124.075,127.245,146.731,144.979,148.210,144.670,142.970,142.524,146.142,146.522,148.128,148.798,150.181,152.388,155.694,160.662,155.520,158.262,154.338,158.196,160.371,154.856,150.636,145.899,141.242,140.834,141.119,139.104,134.437,129.425,123.155,119.273,120.472,121.523,121.983,123.658,124.794,124.827,120.382,117.395,115.790,114.283,117.271,117.448,118.764,120.550,123.554,125.412,124.182,119.828,115.361,114.226,115.214,115.864,114.276,113.469) > par1 = '0' > ylab = 'Density' > xlab = 'Value (€)' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: > if (par1 == '0') bw <- 'nrd0' > if (par1 != '0') bw <- as.numeric(par1) > postscript(file="/var/www/html/freestat/rcomp/tmp/1pqbk1305821783.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > mydensity1<-density(x,bw=bw,kernel='gaussian',na.rm=TRUE) > plot(mydensity1,main='Gaussian Kernel',xlab=xlab,ylab=ylab) Warning messages: 1: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 2: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 3: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 4: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 5: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 6: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for > grid() > dev.off() null device 1 > mydensity1 Call: density.default(x = x, bw = bw, kernel = "gaussian", na.rm = TRUE) Data: x (100 obs.); Bandwidth 'bw' = 8.717 x y Min. : 58.91 Min. :7.428e-06 1st Qu.:100.48 1st Qu.:9.545e-04 Median :142.06 Median :4.106e-03 Mean :142.06 Mean :6.007e-03 3rd Qu.:183.64 3rd Qu.:1.273e-02 Max. :225.22 Max. :1.610e-02 > postscript(file="/var/www/html/freestat/rcomp/tmp/266pa1305821783.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > mydensity2<-density(x,bw=bw,kernel='epanechnikov',na.rm=TRUE) > plot(mydensity2,main='Epanechnikov Kernel',xlab=xlab,ylab=ylab) Warning messages: 1: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 2: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 3: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 4: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 5: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 6: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/3z9n41305821783.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > mydensity3<-density(x,bw=bw,kernel='rectangular',na.rm=TRUE) > plot(mydensity3,main='Rectangular Kernel',xlab=xlab,ylab=ylab) Warning messages: 1: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 2: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 3: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 4: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 5: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 6: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/4dm391305821783.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > mydensity4<-density(x,bw=bw,kernel='triangular',na.rm=TRUE) > plot(mydensity4,main='Triangular Kernel',xlab=xlab,ylab=ylab) Warning messages: 1: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 2: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 3: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 4: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 5: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 6: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5ou4q1305821783.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > mydensity5<-density(x,bw=bw,kernel='biweight',na.rm=TRUE) > plot(mydensity5,main='Biweight Kernel',xlab=xlab,ylab=ylab) Warning messages: 1: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 2: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 3: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 4: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 5: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 6: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6w5yp1305821783.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > mydensity6<-density(x,bw=bw,kernel='cosine',na.rm=TRUE) > plot(mydensity6,main='Cosine Kernel',xlab=xlab,ylab=ylab) Warning messages: 1: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 2: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 3: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 4: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 5: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 6: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7p9jn1305821783.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > mydensity7<-density(x,bw=bw,kernel='optcosine',na.rm=TRUE) > plot(mydensity7,main='Optcosine Kernel',xlab=xlab,ylab=ylab) Warning messages: 1: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 2: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 3: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 4: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for 5: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for <82> 6: In title(...) : conversion failure on 'Value (€)' in 'mbcsToSbcs': dot substituted for > grid() > dev.off() null device 1 > > #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/rcomp/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="/var/www/html/freestat/rcomp/tmp/82x0a1305821783.tab") > > try(system("convert tmp/1pqbk1305821783.ps tmp/1pqbk1305821783.png",intern=TRUE)) character(0) > try(system("convert tmp/266pa1305821783.ps tmp/266pa1305821783.png",intern=TRUE)) character(0) > try(system("convert tmp/3z9n41305821783.ps tmp/3z9n41305821783.png",intern=TRUE)) character(0) > try(system("convert tmp/4dm391305821783.ps tmp/4dm391305821783.png",intern=TRUE)) character(0) > try(system("convert tmp/5ou4q1305821783.ps tmp/5ou4q1305821783.png",intern=TRUE)) character(0) > try(system("convert tmp/6w5yp1305821783.ps tmp/6w5yp1305821783.png",intern=TRUE)) character(0) > try(system("convert tmp/7p9jn1305821783.ps tmp/7p9jn1305821783.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.088 0.229 2.328