Estimation - Note 5

Method of Moment estimation for models on block graphs

Stefka Asenova

2023-02-11

The estimator is the same as the MME for trees described in Vignette “Estimation - Note 1”.


For application of this estimator, see Vignette “Code - Note 5”.


The idea is to find the edge weights \(\delta=(\delta^2_e, e\in E)\) which minimizes the distance between the empirical and the theoretical covariance matrices: \[\begin{equation} \hat{\delta}^{\mathrm{MM}}_{n,k} = \arg\min_{\delta\in(0,\infty)^{E}} \sum_{u\in U} \| \hat{\Sigma}_{W_u, u}-\Sigma_{W_u,u}(\delta) \|_F^2\, . \end{equation}\]

where

\[\begin{equation*} \hat{X}_{v,i} = \frac{1}{1-\hat{F}_{v,n}(\xi_{v,i})}, \qquad v \in U, \quad i = 1, \ldots, n. \end{equation*}\]

\[\begin{equation*} \hat{\Sigma}_{W_u,u} = \frac{1}{|I_u|}\sum_{i\in I_u}(\Delta_{uv,i}-\hat{\mu}_{W_u,u}, v\in W_u\setminus u) (\Delta_{uv,i}-\hat{\mu}_{W_u,u}, v\in W_u\setminus u)^\top\, . %\end{split} \end{equation*}\]

The non-parametric estimators \(\hat{\mu}\) and \(\hat{\Sigma}\) has been suggested in Engelke et al. (2015).

References

Engelke, S., A. Malinowski, Z. Kabluchko, and M. Schlather. 2015. “Estimation of Hüsler-Reiss Distributions and Brown-Resnick Processes.” Journal of the Royal Statistical Society. B 77: 239–65.