如何使用R函数基于trajectory变量统计face类别的频次分布？

我的数据集摘要如下：

trajectory      context_face                     face    
 Min.   :  1   Pass     :404   Length:1204        Frown           :230  
 1st Qu.: 76   Undecided:400   Class :character   Scrunch         : 62  
 Median :151   Allow    :400   Mode  :character   Gasp            : 87  
 Mean   :151                                      Smile           :394  
 3rd Qu.:226                                      Neutral         :258  
 Max.   :301                                      Pout            : 25  
                                                  Wide-opened eyes:148 

我需要按trajectory的Pass、Undecided、Allow三个类别分组，统计各类face（如Frown、Scrunch等）在每个分组下的出现频次，最终得到行是face类型、列是trajectory分组的交叉统计表。

我试过用describeBy(df, df$trajectory)但没得到想要的结果，我的数据集已经是长格式。期望得到如下格式的输出（数值为示例）：

Pass     Undecided      Allow
Frown                   10         2             18
Scrunch                 19         20            4
Gasp                    23         18            14
Smile                   19         11            6
Neutral
Pout
Wide-opened eyes

请问在R中应该用哪个函数实现这个需求？感谢解答，这是我的第一个问题，若表述有不周之处还请谅解。

内容的提问来源于stack exchange，提问作者Ola