我有一个类似的数据框：
word freq1 freq2 freq3 freq4 freq5 freq6
1 a 141 97 103 158 71 101
2 b 84 41 51 42 19 22
3 c 15 24 33 87 102 59
4 d 9 33 41 19 93 75
5 e 1 0 4 15 21 16
6 f 91 31 42 16 75 9
我想计算负数二次分布的参数（r，p），使用R中的最大似然估计我的数据帧中的每个字。我尝试使用方法mle2，但无法获得任何结果。
例如word_a <- c(141,97,103,158,71,101)
如何使用mle方法估计此向量“word_a”的分布参数？