Skew-t distributions derived from skew-normal distributions, as developed by Azzalini and several co-workers, are popular because of their theoretical foundation and the availability of computational methods in the R package sn. One difficulty with this skew-t family is that the elements of the expected information matrix do not have closed form analytic formulas. Thus, we developed a numerical integration method of computing the expected information matrix in the R package skewtInfo. The accuracy of our expected information matrix calculation method was confirmed by comparing the result with that obtained using an observed information matrix for a very large sample size. A Monte Carlo study to evaluate the accuracy of the standard errors obtained with our expected information matrix calculation method, for the case of three realistic skew-t parameter vectors, indicates that use of the expected information matrix results in standard errors as accurate as, and sometimes a little more accurate than, use of an observed information matrix.
Supplementary materials are available in addition to this article. It can be downloaded at RJ-2020-019.zip
Text and figures are licensed under Creative Commons Attribution CC BY 4.0. The figures that have been reused from other sources don't fall under this license and can be recognized by a note in their caption: "Figure from ...".
For attribution, please cite this work as
Martin, et al., "Skew-t Expected Information Matrix Evaluation and Use for Standard Error Calculations", The R Journal, 2020
BibTeX citation
@article{RJ-2020-019,
author = {Martin, R. Douglas and Uthaisaad, Chindhanai and Xia, Daniel Z.},
title = {Skew-t Expected Information Matrix Evaluation and Use for Standard Error Calculations},
journal = {The R Journal},
year = {2020},
note = {https://doi.org/10.32614/RJ-2020-019},
doi = {10.32614/RJ-2020-019},
volume = {12},
issue = {1},
issn = {2073-4859},
pages = {188-205}
}