r - 我怎样才能找到一组数据的分布，然后进一步传播这个分布？

Question

我的问题是我有一组数据，我想对其进行分布拟合，然后在找到分布后，对其运行蒙特卡洛模拟以传播找到的分布。

我的第一段代码是:

require(fitdistrplus)
example1<-c(29,23,29,25,26,29,29,27,25,25,25,26,28,25,29,28,28,26,28,25,29,26,30) 
f1<-fitdist(example1,rgamma,method="mle")

如果我然后使用命令

print(f1)

它告诉我形状是 204.00 Gamma 分布的比率为 7.568

(请注意，目前我拟合分布的数字是任意的，我通常会有数百个观察值来拟合分布)。

我现在需要帮助的地方是当我使用包 mc2d 中的代码来传播此分发时，如下所示:

require(mc2d)
ndunc(1000)
fitted<-mcstoc(rgamma, type="U", shape=204.00, rate=7.569)

目前，我必须手动输入形状并根据“fitdist”命令的先前“打印”输入上述函数。

我的问题是，有没有办法让 mcstoc 命令自动从 fitdist 命令中获取形状和速率，这样我就不必中断代码来手动执行此操作？或者，如果 fitdistrplus 包和 mc2d 包不可能，那么是否有另一个包可以为我做这件事？

非常感谢!

Answer 1

f1$estimate[1]
#   shape 
#204.0008 

f1$estimate[2]
#    rate 
#7.567762

fitted<-mcstoc(rgamma, type="U", shape=f1$estimate[1], rate=f1$estimate[2])