[R] Plot histogram and lognormal fit on the same time

[Bert Gunter]

OK, but that's not the point of my comment.

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 Fri, Jan 22, 2021 at 5:03 PM Abby Spurdle <spurdle.a using gmail.com> wrote:

> Sorry, Bert.
> The fitdistr function estimates parameters via maximum likelihood.
> (i.e. The "lognormal" part of this, is not a kernel).
>
>
> On Fri, Jan 22, 2021 at 5:14 AM Bert Gunter <bgunter.4567 using gmail.com>
> wrote:
> >
> > 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 using 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 using icmpe.cnrs.fr <mailto:leroy using icmpe.cnrs.fr>
> > >
> > > Page Web : : http://www.icmpe.cnrs.fr <http://www.icmpe.cnrs.fr>
> > >
> > > ______________________________________________
> > > R-help using 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]]
> >
> > ______________________________________________
> > R-help using 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]]