In [1]:

    

data <- scan('rn01.txt')


    

In [2]:

    

x.lim <- round(c(min(data)-1, max(data)+1))
his <- hist(data, breaks="Scott", xlim=x.lim, prob=TRUE)


    

In [3]:

    

histx <- function(his, x) {
    "Return height of histogram at given x."
    b <- diff(his$breaks)[1]
    r <- length(his$breaks)
    y <- ifelse( (x < his$breaks[1]) | (x > his$breaks[r]), 0,
               his$density[(abs(x-his$breaks[1]) %/% b)+1])
    
    return(y)
}


    

In [4]:

    

b <- diff(his$breaks)[1]
r <- length(his$mids)
x1 <- c(his$mids[1]-b, his$mids, his$mids[r]+b)
his <- hist(data, breaks="Scott", xlim=x.lim, prob=TRUE)
lines(x1, histx(his, x1), lty=1, lwd=2)


    

In [ ]: