找回密码
 欢迎注册
查看: 27448|回复: 2

[原创] 复合随机变量的分布及参数估计

[复制链接]
发表于 2014-9-16 10:46:37 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
假设有已知一维分布的随机变量 `X_i \sim F(x;\theta) \;(i=1,2,\cdots ,N)`,而 `N \sim G(x;\beta)`。其中 `F(\cdot)` 和 `G(\cdot)` 为累计分布函数(CDF), `\theta,\beta` 为参数。
如果一个新的随机变量 $$Y=\sum_{i=1}^N X_i\tag{1}$$问:
1. 在已知参数值的情况下,1) 若已知随机变量 `X_i` 两两之间的协方差矩阵(或相关系数矩阵);2) 若随机变量 `X_i` 是独立同分布的(i.i.d.),求 `Y` 的概率密度函数(PDF).
2. 已知随机变量 `Y` 的 `n` 个样本 `y_1,y_2,\cdots,y_n`,求参数 `\theta,\beta` 的估计.

问题背景

最近在处理一个 单变量跳扩散(Jump-Diffusion)随机微分方程模型的参数估计问题$$\dif S(t)=a(b-S(t))\dif t+\sqrt{\small S(t)}\sigma \dif W(t)+\dif Y(t) $$其中,`W(t)` 是标准的维纳随机过程;`Y(t)` 是跳跃项(去掉跳跃微分项后便是经济学中著名的 CIR 随机微分方程模型)。待估参数是 `a,b,\sigma` 以及 `Y(t)` 中涉及到的相关参数。

`Y(t)` 的形式是什么?它被定义为一个复合的随机过程。即在式 `(1)` 中,将求和上限 `N` 换成一个随机过程 `N(t)`,比如泊松随机过程:$$ p (N(t)=k)=\frac{(\lambda t)^k}{k!}\mathrm e^{\lambda t}$$`X` 仍为服从已知分布 `F(x;\theta)` 的随机变量。这样,`Y(t)` 的参数就是是 `N(t)`以及 `X` 中的参数 `\lambda,\theta`。

参数估计考虑使用极大似然估计(MLE)方法,但是这个跳跃项是离散求和的形式,它的微分 `\D \dif \sum_{i=1}^{N(t)} X_i` 难以写成 Ito 随机微分形式,所以得不到跳跃微分项的PDF,因而似然函数(即`S_{t_i}`的联合概率密度分布函数)就不容易写出表达式。

本文开头的问题便是上面的问题的简化情形(定常情形,即不考虑时间变化)。若将其实例化,我们可以考虑选择 `F(x;\theta)` 是参数为 `\theta` 的指数分布,`G(x;\beta)` 是参数为 `\beta` 的泊松分布。那么问题1和2如何解决呢?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-9-16 11:00:26 | 显示全部楼层
对问题背景内容了解不多,仅考虑初始的简化问题。
若干个独立同分布随机变量之和的分布规律,理论上是通过卷积进行计算的,而且不容易得出一个“好看”的结果。

点评

关键是“若干个”这个量,本身也是一个随机变量。  发表于 2014-9-16 11:39
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-24 08:04 , Processed in 0.025745 second(s), 18 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表