Simulate from parametric distribution

rfitdistr(x, n = length(x), p = 1)