我正在通过循环遍历所有因变量名称来进行多个 TukeyHSD 测试，如 Eric Lecoutre 在此线程中提供的代码中所示：
它工作得很好，所以感谢埃里克！但在输出中，变量名称写为数字：[[ 1 ]] 表示第一个变量，[[ 2 ]] 表示第二个变量，...[[ 137 ]] 表示第 137 个变量，等等。
如何将这些数字更改为输出中的变量名（以便我可以将输出捕获到文件中并在输出文件中保留有关变量名的信息）？
具体来说，在 Eric Lecoutre 回答的输出中：
[[ 1 ]]
均值的 Tukey 多重比较
95% family-wise confidence level
拟合：aov（公式 = 模型）
例如，如果循环的 2 个因变量的名称是“sepal.width”和“sepal.length”，我想将输出的 [[ 1 ]] 更改为 [[ sepal.width ]]，[[ 2 ]] 的输出转换为 [[ sepal.length ]]。
如何使脚本循环遍历数字并将数字更改为变量名称，以便输出如下所示（我省略了统计结果以简化输出）：
[[萼片宽度]]
均值的 Tukey 多重比较
95% family-wise confidence level
拟合：aov（公式 = 模型）
[[萼片长度]]
均值的 Tukey 多重比较
95% family-wise confidence level
拟合：aov（公式 = 模型）
非常感谢！