在R语言中在分面直方图上绘制von Bertalanffy生长曲线

我需要将通过模态进展分析得到的von Bertalanffy生长曲线（描述鱼类随时间的生长规律），叠加到按月份分面的体长频率直方图上。以下是相关信息：

当前生成的无曲线图形 #

按月份分面的体长频率直方图，每个面板对应一个月的体长分布数据，暂未叠加生长曲线。

数据样例 #

> head(L_data)
  Length Month
1   18.5    01
2   61.5    01
3   61.5    01
4   62.5    01
5   62.5    01
6   63.5    01

现有绘图代码（无曲线版本） #

geom_histogram(bins = 60, color = "black", fill = "white") +
  facet_wrap(~Month, ncol = 1) +
  theme_bw() + theme(axis.text.y = element_blank(),
                     axis.ticks = element_blank(),
                     strip.text.x = element_blank(),
                     panel.grid = element_blank(),
                     panel.border = element_blank()) +
  ylab("") + xlim(5,100)

目标效果 #

每个分面的直方图上叠加一条非线性的von Bertalanffy生长曲线，曲线参数由贝叶斯推断得出。也可接受使用R基础绘图函数实现该需求。

内容的提问来源于stack exchange，提问作者Matheus de Barros