R: адаптация getDescriptionStatsBy для предоставления многослойных средств по нескольким переменным

Question

Я хотел бы создать красиво отформатированную таблицу для бумаги в том же стиле, что и getDescriptionStatsBy (либо с использованием этой функции, либо с использованием чего-то похожего). Я успешно использовал эту функцию для создания моей Таблицы 1 (выборочные характеристики), а затем для создания таблицы средних баллов моей переменной интереса (которая измеряется с помощью 6 различных инструментов), стратифицированной по одной переменной (настройка):

getTable2Stats <- function(x, digits = 3, ...){
        getDescriptionStatsBy(x = x, 
                              by = as.factor(mydata$setting),    
                              digits = digits,
                              header_count = TRUE,
                              NEJMstyle = TRUE,
                             useNA = "always",    
                              ...)
      }



      t2 <- list()
      t2[["Instrument 1"]] <- getTable2Stats(mydata$instrument1_score)
      t2[["Instrument 2"]] <- getTable2Stats(mydata$instrument1_score)
      t2[["Instrument 3"]] <- getTable2Stats(mydata$instrument1_score)
      t2[["Instrument 4"]] <- getTable2Stats(mydata$instrument1_score)
      t2[["Instrument 5"]] <- getTable2Stats(mydata$instrument1_score)
      t2[["Instrument 6"]] <- getTable2Stats(mydata$instrument1_score)


      mergeDesc(t2,
                htmlTable_args = list(css.rgroup = "",
                                      caption  = "Table 2: Mean scores by setting"))

Это создает таблицу, в которой строки - это 6 инструментов, а столбцы - настройки.

Я хотел бы составить таблицу, в которой представлены средние +/- баллы по SD для каждого из шести инструментов, стратифицированные по множеству категориальных переменных, все в одной таблице - настройка, возрастной диапазон, пол, раса и др. c. как на этом изображении (заполненном поддельными примерами):

Это означает, что по сравнению с моим кодом выше строки и столбцы должны быть переключены; и что еще более важно, стратификация будет по нескольким переменным в дополнение к настройке. Я часами пробовал разные вещи, в том числе преобразовывал свои данные в длинный формат, чтобы в каждой строке был один счет с отдельным столбцом для инструмента. Затем я установил вход «by» функции в переменную инструмента, «x» - в переменную Score и заполнил свой список (t2) переменными, по которым я хотел разбить данные. Но это не создает новую строку для каждого факторного уровня переменных, которые я хочу разделить ... Я также пытался использовать data.table и другие пакеты, но я просто не могу понять это! Любая помощь приветствуется!