- 注册时间
- 2013-10-24
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 8856
- 在线时间
- 小时
|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?欢迎注册
×
假设有已知一维分布的随机变量 `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如何解决呢? |
|