set.seed(121343) u <- rnorm(100)
#default histogram hist(u) | |
#with shading hist(u, density=20) | |
#with specific number of bins hist(u, density=20, breaks=20) | |
# proportion, instead of frequency # also specifying y-axis hist(u, density=20, breaks=-3:3, ylim=c(0,.5), prob=TRUE) | |
# overlay normal curve with x-lab and ylim # colored normal curve m<-mean(u) std<-sqrt(var(u)) hist(u, density=20, breaks=20, prob=TRUE, xlab="x-variable", ylim=c(0, 0.7), main="normal curve over histogram") curve(dnorm(x, mean=m, sd=std), col="darkblue", lwd=2, add=TRUE) | |
hist(u, density=20, breaks=20, prob=TRUE, xlab="x-variable", ylim=c(0, 0.8), main="Density curve over histogram") lines(density(u), col = "blue") | |
# hist(u) is an object # names(uh) will show all of its components uh<-hist(u) plot(uh, ylim=c(0, 40), col="lightgray", xlab="", main="Histogram of u") text(uh$mids, uh$counts+2, label=c(uh$counts)) |
No comments:
Post a Comment