将原始分布函数拟合到 R 中的数据

Question

假设我有一个随机数据集（使用帕累托分布生成）：

dataset <- rpareto(1000, shape=3, scale=100)

我有一个函数 f(x)，它基本上是一个帕累托分布，参数形状为 k 和尺度 sigma：

但我也添加了另一个参数 (m) 以获得：

如何测试我的第二个 f(x)（添加参数 m）与我的数据集的匹配程度？ 过去，我使用过 fitdistr 包函数 fistdist() 使用 mle 方法将帕累托函数拟合到我的数据中，但据我所知，这个包只能用于测试现有函数的拟合，而不是像矿。