ggplot()+aes(X)+geom_density()+
stat_function(fun = dpois, colour = "red", args = list(lambda=parms$estimate))+
xlab( "confirmation time [s]") + ggtitle("Shape of curve");
leads to:
but if args = list(mean = mean(X),sd = sd(X)))
we get a nice normal distribution:
How can you explain this strange behavior?