You're absolutely right. Was thinking exponential for some reason.... The rootogram is quite nice; I've never seen one before. Thanks!
Michael On Thu, Nov 3, 2011 at 5:39 PM, David L Carlson <[email protected]> wrote: > The lines() command doesn't work and histogram combines categories unless > you specify the number. How about a barplot > > Lam <- 3 > X <- table(rpois(500, Lam)) > Max <- length(X)-1 > barplot(rbind(X, 500*dpois(0:Max, Lam)), beside=TRUE, > legend.text=c("Observed", "Expected")) > > or a rootogram > library(vcd) > rootogram(X, dpois(0:Max, Lam)*500) > > > -----Original Message----- > From: [email protected] [mailto:[email protected]] On > Behalf Of R. Michael Weylandt <[email protected]> > Sent: Thursday, November 03, 2011 2:55 PM > To: kerry1912 > Cc: [email protected] > Subject: Re: [R] Histograms in R > > Try something like this > > Lam <- 3 > X <- rpois(500, Lam) > hist(X, freq = F) > x <- seq(min(X), max(X), length = 500) > lines(x, dpois(x, Lam), col=2) > > Adapt as necessary > > Michael > > On Nov 3, 2011, at 8:03 AM, kerry1912 <[email protected]> wrote: > >> We have a histogram of our observed response and we want to overlay the >> corresponding poisson distribution with respect to our poisson model. >> >> >> -- >> View this message in context: > http://r.789695.n4.nabble.com/Histograms-in-R-tp3985397p3985397.html >> Sent from the R help mailing list archive at Nabble.com. >> >> ______________________________________________ >> [email protected] mailing list >> https://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. > > ______________________________________________ > [email protected] mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
