R语言怎么分组计算协方差阵？

数据框df，按照某一因子的不同取值，可以分成几个子数据框，然后分别计算这些子数据框的协方差（计算前将因子变量剔除），请问不用循环可以做吗？比如我尝试了下dplyr包中group_by函数等分组函数，但是没办法做到

重要提示：提问者不能发表回复，可以通过评论与回答者沟通，沟通后可以通过编辑功能完善问题描述，以便后续其他人能够更容易理解问题.