Plot histogram and lognormal fit on the same time
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 at 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 at 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 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]]
______________________________________________ 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.