OK, ваш пример здесь не совсем воспроизводимый, но я думаю, что ответ, который вы хотите, будет выглядеть примерно так:
result <- los_data %>%
group_by(group, setting, diagnosis) %>%
do({
fit <- fitdist(.$my_column, "gamma")
data_frame(group=.$group[1], setting=.$setting[1], diagnosis=.$diagnosis[1], fit = list(fit))
}) %>%
ungroup()
Это даст вам фрейм данных всех приступов, с столбцы для группы, настройка, диагноз, а также столбец-столбец, который содержит подгонки для каждого из них. Поскольку это столбец списка, вам нужно будет использовать двойные скобки для извлечения отдельных подстановок. Пример:
# Get the fit in the first row
result$fit[[1]]