Goodness-of-fit indice for `fitdistr` and `fitdist` in `R`

Question

The following code:

library(ggplot2)
library(fitdistrplus)
library(MASS)
library(extraDistr)  # for the Rayleigh distribution functions

estimate <- c(19.00000