复合对数正态帕累托分布

Question

对于我正在NetLogo中研究的模型，我想将某个变量X分配为复合对数正态-对数，它采用对数正态密度直至未知阈值θ和此后具有缩放参数α的两参数帕累托密度。

以下是pdf文件：

f1（x）=（αθ^α/（1 +Φ（k））x ^（α+ 1））exp {（-α^ 2 / 2k ^ 2）ln ^ 2（x /θ）} ，0 f 2（x）＝αθ^α/（1 +Φ（k））x ^（α+ 1）。 ，θ

其中Φ（·）是标准正态分布的累积分布函数，k是常数。

如何在NetLogo中编写此代码？