Density estimate and histogram computation

gistfile1.r

# SCRIPT_STR('
x <- .arg1
xanz <- .arg2

density_estimate <- density(rep(x, xanz))
dens <- sapply(x, function(i){
  which.min(abs(i - density_estimate$x))
})

hist_density <-  hist(rep(x, xanz), freq=FALSE, 
    breaks=c(x,max(x)+median(diff(x))), right = F)

return(paste(density_estimate$y[dens],hist_density$density,sep="~"))
# ',
# MAX([x bin]), COUNT([x axis (raw)]) )