Dopasowałem rozkład normalny za pomocą funkcji fitdist
z pakietu . Używając , i , możemy wykreślić odpowiednio, jak podano poniżej.Tworzenie działek dopasowanych z ggplot2
set.seed(12345)
df <- rnorm(n=10, mean = 0, sd =1)
library(fitdistrplus)
fm1 <-fitdist(data = df, distr = "norm")
summary(fm1)
denscomp(ft = fm1, legendtext = "Normal")
qqcomp(ft = fm1, legendtext = "Normal")
cdfcomp(ft = fm1, legendtext = "Normal")
ppcomp(ft = fm1, legendtext = "Normal")
Jestem bardzo zainteresowany, aby te działki fitdist
z ggplot2
. MWE jest poniżej:
qplot(df, geom = 'blank') +
geom_line(aes(y = ..density.., colour = 'Empirical'), stat = 'density') +
geom_histogram(aes(y = ..density..), fill = 'gray90', colour = 'gray40') +
geom_line(stat = 'function', fun = dnorm,
args = as.list(fm1$estimate), aes(colour = 'Normal')) +
scale_colour_manual(name = 'Density', values = c('red', 'blue'))
ggplot(data=df, aes(sample = df)) + stat_qq(dist = "norm", dparam = fm1$estimate)
Goraco wdzięczny, jeśli ktoś daje mi wskazówki, aby te fitdist
działek z ggplot2
. Dzięki
Gdyby to nie oznaczało nagrody, głosowałbym, by zamknąć jako zbyt szeroki. Każdy wykres powinien być innym pytaniem (chociaż nie musisz pytać o każdy wykres, jeśli masz odpowiedź na jedną lub dwie z nich). – Roland