Rで複数の変数にわたってテーブル関数を実行し、関数を使用して結果を新しいデータセットにコンパイルするにはどうすればよいですか？

約100個の変数のデータセットがあり、これらの変数の約30個の要約を含むテーブルを作成したいと思います。これを行うために、これらの変数に対してテーブルと他の関数を手動で実行し、結果をバインドしました。ただし、30以上の変数に対してこれを行う必要があるため、関数を使用してプロセスを自動化したいと思います。

データセットの例を次に示します。

df <- data.frame(v1=c('a','b','c','c','b'),v2=c('d','d','e','e','e'),cat=c('1low','1low','2med','3high','2med'))

目標は、以下のようなテーブルを作成することです（NAなし）。

以下は私のコードです：

library(formattable)

# For var1 & var2, apply the table function and convert to dataframe so that the row labels are incorporated into dataset
var1.df <- as.data.frame(table(df$v1, df$cat))

# reshape to achive wide format (goal to view the count of each var1 level across low, med, high cats)
var1.df <- reshape(var1.df, idvar = "Var1", timevar = "Var2", direction = "wide")

# add col names
names(var1.df) <- c("vcat","low","med","high"); var1.df

# repeat above steps for next variable. in true dataset, I will need to repeat for 30 vars...
var2.df <- as.data.frame(table(df$v2, df$cat))
var2.df <- reshape(var2.df, idvar = "Var1", timevar = "Var2", direction = "wide")
names(var2.df) <- c("vcat","low","med","high")

# Create variable headings
var1.heading <- data.frame("variable 1",NA,NA,NA) # ideally, the NAs are blanks
names(var1.heading) <- c("vcat","low","med","high")

var2.heading <- data.frame("variable 2","","","")
names(var2.heading) <- c("vcat","low","med","high")

# Rbind the category headings and the table result data
table01 <- do.call("rbind", list(var1.heading, var1.df, 
                                 var2.heading, var2.df))

# Format the table for presentation
heading.list <- c("variable 1", "variable 2")
x <- formattable(table01, 
                 align =c("l","c","c","c","c"),
                 list(vcat = formatter("span", style = x ~ ifelse(x %in% heading.list, 
                                                                  style(font.weight = "bold"), NA))))

上記のコードを自動化するための以下の試みは、不完全（a）であるか、正しく実行されません（b）

# (a)
lapply(df, function(x) as.data.frame(table(x, df$cat)))

# (b)
myfxn <- function(x){
  y <- as.data.frame(table(x, df$cat))
  y <- reshape(y, idvar = "x", timevar = "Var2", direction = "wide")
  names(y) <- c("vcat","low","med","high")
}
lapply(df, myfxn(x))

さらにいくつかの変数についてこのプロセスを自動化する方法について何か提案はありますか？また、1行のデータフレームの挿入を手動で作成する以外に、テーブルにカテゴリ見出しを挿入する別の方法はありますか？これは最初のデータフレームであるため、var1.headingにNAを挿入したことに注意してください。空白の代わりに「」を挿入しようとすると（var2.headingなど）、後続のデータフレームは文字ではなく因子変数であるためバインドされませんでした。事前にどうもありがとうございました！