Hi, I would like to plot the histogram of data and fit it with a lognormal distribution. The ideal, would be to superimpose the fit on the histogram and write the results of the fit on the figure. Right now, I was able to plot the histogram and fit the density with a lognormal, but I can't combine all together. Here is the code I wrote : histdata <- hist(dataframe$data) library(MASS) fitdistr(histdata$density, "lognormal") Can you help me ? Best regards, -- *Eric Leroy* /Responsable de la plateforme microscopie ?lectronique/ /Correspondant S?curit? des Syst?mes Informatiques de l'ICMPE/ ICMPE - UMR 7182 - CNRS - UPEC 2/8, rue Henri Dunant 94320 Thiais T?l : 01.49.78.12.09 Fax : 01.49.78.12.03 courriel : leroy at icmpe.cnrs.fr <mailto:leroy at icmpe.cnrs.fr> Page Web : : http://www.icmpe.cnrs.fr <http://www.icmpe.cnrs.fr>
Two solutions not exactly equivalent ; data <- rlnorm(100, meanlog = 1, sdlog = 1) histdata <- hist(data, ylim=c(0, 100)) library(MASS) f <- fitdistr(data, "lognormal") f$estimate lines(x=seq(from=0, to=50, by=0.1), ? y=dlnorm(x=seq(from=0, to=50, by=0.1), meanlog = f$estimate["meanlog"], sdlog = f$estimate["sdlog"])*400 ) library(HelpersMG) m <- modeled.hist(breaks=histdata$breaks, FUN=plnorm, ????????????????????????????? meanlog = f$estimate["meanlog"], sdlog = f$estimate["sdlog"], sum = 100) points(m$x, m$y, pch=19, col="red") Marc Girondot Le 21/01/2021 ? 12:54, Eric Leroy a ?crit?:> Hi, > > I would like to plot the histogram of data and fit it with a lognormal > distribution. > > The ideal, would be to superimpose the fit on the histogram and write > the results of the fit on the figure. > > Right now, I was able to plot the histogram and fit the density with a > lognormal, but I can't combine all together. > > Here is the code I wrote : > > histdata <- hist(dataframe$data) > > library(MASS) > > fitdistr(histdata$density, "lognormal") > > Can you help me ? > > Best regards, > > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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.[[alternative HTML version deleted]]
In future, you should try to search before posting. I realize that getting good search terms can sometimes be tricky, but not in this case: 'plot density with histogram in R' or similar on rseek.org or just on your usual search platform brought up several ways to do it. As a (slightly offtopic) side issue, be careful with this: both histograms and smooth density estimates are just that: density **estimates** that depend on both the data and various parameters, e.g. the location of the bins for histograms and kernel for kernel density estimates. You may therefore get plots where the two do not appear to "match" when you use the defaults for each. Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Thu, Jan 21, 2021 at 3:54 AM Eric Leroy <leroy at icmpe.cnrs.fr> wrote:> Hi, > > I would like to plot the histogram of data and fit it with a lognormal > distribution. > > The ideal, would be to superimpose the fit on the histogram and write > the results of the fit on the figure. > > Right now, I was able to plot the histogram and fit the density with a > lognormal, but I can't combine all together. > > Here is the code I wrote : > > histdata <- hist(dataframe$data) > > library(MASS) > > fitdistr(histdata$density, "lognormal") > > Can you help me ? > > Best regards, > > -- > > *Eric Leroy* > > /Responsable de la plateforme microscopie ?lectronique/ > > /Correspondant S?curit? des Syst?mes Informatiques de l'ICMPE/ > > ICMPE - UMR 7182 - CNRS - UPEC > > 2/8, rue Henri Dunant > > 94320 Thiais > > T?l : 01.49.78.12.09 > > Fax : 01.49.78.12.03 > > courriel : leroy at icmpe.cnrs.fr <mailto:leroy at icmpe.cnrs.fr> > > Page Web : : http://www.icmpe.cnrs.fr <http://www.icmpe.cnrs.fr> > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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. >[[alternative HTML version deleted]]