The information matrix of the bivariate extended skew-normal distribution

For the extended skew-normal distribution, which represents an extension of the normal (or Gaussian) distribution, we focus on the properties of the log-likelihood function and derived quantities in the bivariate case. Specifically, we derive explicit expressions for the score function and the information matrix, in the observed and the expected form; these do not appear to have been examined before in the literature. Ensuing numerical work illustrates some relevant aspects of the expected information matrix.