R语言 多因素协方差计算

R语言的协方差的输入必须是：aov(观测变量域名～协变量+控制变量）这样吗？
不可以像Python那样直接cov()吗？