我想知道使用fixest计算和表示标准化系数的最佳方法是什么。

library(parameters)
  library(effectsize)
  library(fixest)
  
  m <- lm(rating ~ complaints, data = attitude) 
  standardize_parameters(m, method="basic")# works

  m <- feols(rating ~ complaints, data = attitude) 
standardize_parameters(m, method="basic")# Error in stats::model.frame(model)[[1]] : subscript out of bounds

我还尝试了modelsummary方法，但它显示了未标准化的系数，没有错误。

library(parameters)
   library(effectsize)
   
   m <- lm(rating ~ complaints, data = attitude)
   
   modelsummary(m, standardize="refit") # works, coeffs are different

   m <- feols(rating ~ complaints, data = attitude)

   modelsummary(m, standardize="refit")# doesn't work, coeffs are the same

任何关于如何优雅而轻松地从fixest估计结果中提取标准化系数的见解或建议都将非常感谢。我的目标是复制在Stata中使用listcoef包的古老做法。非常感谢本文中提到的包的作者！
编辑：“”〉软件包版本（“模型摘要”）[1]“1.1.0.9000”