R Under development (unstable) (2026-01-22 r89323 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2026 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> require("fitdistrplus")
Loading required package: fitdistrplus
Loading required package: MASS
Loading required package: survival
> nbboot <- 100
> nbboot <- 10
> set.seed(123) # here just to make random sampling reproducible
> 
> nsample <- 10
> 
> # (1) Description of a sample from a normal distribution
> # with and without uncertainty on skewness and kurtosis estimated by bootstrap 
> #
> x1 <- rnorm(nsample)
> descdist(x1)
summary statistics
------
min:  -1.265061   max:  1.715065 
median:  -0.07983455 
mean:  0.07462564 
estimated sd:  0.9537841 
estimated skewness:  0.715477 
estimated kurtosis:  2.858524 
> descdist(x1,boot=nbboot)
summary statistics
------
min:  -1.265061   max:  1.715065 
median:  -0.07983455 
mean:  0.07462564 
estimated sd:  0.9537841 
estimated skewness:  0.715477 
estimated kurtosis:  2.858524 
> 
> # (2) Description of a sample from a beta distribution
> # with uncertainty on skewness and kurtosis estimated by bootstrap
> # with changing of default colors 
> #
> descdist(rbeta(nsample,shape1=0.05,shape2=1),boot=nbboot,
+ obs.col="blue",boot.col="orange")
summary statistics
------
min:  4.15903e-23   max:  0.4220559 
median:  5.825247e-08 
mean:  0.06196849 
estimated sd:  0.1387746 
estimated skewness:  2.425894 
estimated kurtosis:  8.742888 
> 
> # (3) Description of a sample from a gamma distribution
> # with uncertainty on skewness and kurtosis estimated by bootstrap
> # without plotting 
> #
> descdist(rgamma(nsample,shape=2,rate=1),boot=nbboot,graph=FALSE)
summary statistics
------
min:  0.8772438   max:  4.943292 
median:  2.409935 
mean:  2.296641 
estimated sd:  1.279247 
estimated skewness:  0.9159561 
estimated kurtosis:  3.67572 
> 
> # (3) Description of a sample from a Poisson distribution
> # with uncertainty on skewness and kurtosis estimated by bootstrap 
> #
> descdist(rpois(nsample,lambda=2),discrete=TRUE,boot=nbboot)
summary statistics
------
min:  0   max:  6 
median:  1 
mean:  2.1 
estimated sd:  2.378141 
estimated skewness:  0.9678028 
estimated kurtosis:  2.334913 
> 
> # (4) Description of serving size data
> # with uncertainty on skewness and kurtosis estimated by bootstrap 
> #
> data(groundbeef)
> serving <- groundbeef$serving
> descdist(serving, boot=nbboot)
summary statistics
------
min:  10   max:  200 
median:  79 
mean:  73.64567 
estimated sd:  35.88487 
estimated skewness:  0.7352745 
estimated kurtosis:  3.551384 
> 
> proc.time()
   user  system elapsed 
   1.62    0.28    1.89