lines density

Hi all,

I attach a picture, with my output plot.

I put this command to obtain it:

hist(MCT[,2],probability=TRUE,xlab=" ",ylab="Max",main="M distribution")
lines(density(MCT[,2]), lwd = 2)

The problem is that when I extract the "bell" is too high? And the higher
part doesn?t appear in the plot. There is some trick to adjust the line on
the plotted area?

You need to set ylim in your call to hist().  You can do it by trial and 
error, or compute the numbers first via

h <- hist(MCT[,2], plot=FALSE)

and set

ylim <- range(c(0, h$density, density(MCT[,2])$y))

before you plot:

hist(MCT[,2], probability=TRUE, ylim=ylim, xlab=" ", ylab="Max", main="M 
distribution")
lines(density(MCT[,2]), lwd = 2)

Duncan Murdoch