Skip to contents

For more information, check out the smoothic website.

Implementation of the SIC epsilon-telescope method, either using single or multi-parameter regression. Includes classical regression with normally distributed errors and robust regression, where the errors are from the Laplace distribution. The “smooth generalized normal distribution” is used, where the estimation of an additional shape parameter allows the user to move smoothly between both types of regression. See O’Neill and Burke (2022) “Robust Distributional Regression with Automatic Variable Selection” for more details on arXiv. This package also contains the data analyses from O’Neill and Burke (2021). “Variable Selection Using a Smooth Information Criterion for Multi-Parameter Regression Models” on arXiv.



You can install the released version of smoothic from CRAN with:



Install the development version from GitHub:

# install.packages("devtools")


The smoothic() function performs automatic variable selection with distributional regression.

fit <- smoothic(
  formula = y ~ .,
  data = dataset

A summary table that includes the estimated coefficients, estimated standard errors (SEE) and the value of the penalized likelihood function is returned with:


Further information and examples of implementation (including plotting of the coefficient paths - vignette("sgnd-boston")) are available in the function documentation and vignettes.