Построение сложных функций распределения в RStudio

Прошу Вас порекомендовать пакеты в RStudio, с помощью которых можно осуществить подгонку данных по 1 - минутным логарифмическим доходностям под сложные одномерные функции распределения вероятностей. В частности, необходимы:

• обобщенное распределение ошибки (GED);
• обратное нормальное распределение (NIG);
• Variance - Gamma - распределение (VG);
• обобщенное гиперболическое распределение (GHYP);
• распределение Бурра (BURR)

Массив данных, по которым необходимо оценить параметры вышеуказанных функций распределения, достаточно велик и содержит порядка 30000 значений.

Также хотелось бы узнать, существуют ли какие-то функции или пакеты, позволяющие определить, насколько точно выбранное теоретическое распределение описывает эмпирическое распределение данных, а также сравнить несколько теоретических функций распределения между собой и выбрать из них наиболее подходящее под имеющиеся данные.

Заранее спасибо за помощь.