我正在使用nnet运行一个多项logit模型，然后想在modelsummary表中显示结果，因子水平作为列：

library("nnet")

multi <- multinom(DV ~ var1, data = wave1)

alpha <- summary(multi)$coefficients / summary(multi)$standard.errors
abs(alpha)
# Calculate p values by hand
p_values <- (1 - pnorm(abs(alpha))) * 2
p_values

library("modelsummary")
modelsummary_wide(multi, output = "gt", stars = TRUE)

以前我使用modelsummary_wide（）实现了这一点，但现在显然不推荐使用，我应该在modelsummary（）中使用shape参数。
但是，我真的不知道怎么做。有人能给我一些提示吗？
谢谢