Journal of Probability and Statistics

Volume 2015 (2015), Article ID 612024, 7 pages

http://dx.doi.org/10.1155/2015/612024

## Posterior Analysis of State Space Model with Spherical Symmetricity

Department of Statistics, University of Delhi, Delhi 110007, India

Received 28 June 2015; Revised 15 October 2015; Accepted 18 October 2015

Academic Editor: Shein-chung Chow

Copyright © 2015 Ranjita Pandey. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The present work investigates state space model with nonnormal disturbances when the deviation from normality has been observed only with respect to kurtosis and the distribution of disturbances continues to follow a symmetric family of distributions. Spherically symmetric distribution is used to approximate behavior of symmetric nonnormal disturbances for discrete time series. The conditional posterior densities of the involved parameters are derived, which are further utilized in Gibbs sampler scheme for estimating the marginal posterior densities. The state space model with disturbances following multivariate- distribution, which is a particular case of spherically symmetric distribution, is discussed.

#### 1. Introduction

Often in the analysis of various time series dynamic or state space models, the distribution of disturbances is assumed to be normal. However, the assumption of normal disturbances may not be satisfied in many practical applications of these models. For example, deviations from normality have been observed in stock price data where behavior of disturbances appears to be leptokurtic; see, for example, Kendall [1] and Fama [2]. Praetz [3] suggested scaled- distribution to explain leptokurtic behavior of stock price data. If the tails of the distribution are flatter than the normal distribution, the multivariate- distribution may provide a more realistic model for the disturbances. Zellner [4] revealed that dependent but uncorrelated responses could be analyzed using a student- model. Fraser [5] demonstrated the robustness of the student- family as opposed to the normal distribution based on empirical studies and further suggested that a student- model holds equally good for normal model responses. Another advantage of student- family is that it handles outliers easily; see Sutradhar and Ali [6]. Student- model is also advocated by Haq and Khan [7], Khan and Haq [8], and Chib et al. [9]. Prucha and Kelejian [10] observed that normal model based analysis is heavily influenced by extreme observations and deviations from assumptions and ignores sample information beyond the first two moments. Robustness for Bayes prediction under multivariate- case has also been studied by Chib et al. [11], which is a particular member of the spherically symmetric family. Jammalamadaka et al. [12] have, however, shown that predictive inferences are completely uninfluenced by departure from the normality assumption in the direction of the spherically symmetric family. Properties of the spherical distribution have been studied by Kelker [13] and Chmielewski [14]. In carrying out the Bayesian analysis of various dynamic models involving multidimensional parameter vector, it is often difficult to derive explicit expressions for the marginal posterior distributions of various parameters. In such circumstances, Markov Chain Monte Carlo (MCMC) technique, also known as Gibbs sampler, provides a useful way for estimating the marginal posterior densities. Klock and Van Dijk [15], Naylor and Smith [16], and Geweke [17–19] applied numerical MCMC for the Bayesian analysis of various time series models. Tiwari et al. [20] considered a general state space or linear dynamic model with normal disturbances and carried out the Bayesian analysis of the model using Gibbs sampler. In the present exposition, the assumption of normality of disturbances is relaxed and it is assumed that the disturbances follow distribution belonging to the family of spherically symmetric distributions.

#### 2. Model and Assumptions

The following form of general state space (or linear dynamic model) is considered:where is the observation on the response variable, is a regression vector, is a vector of parameters generated by the first-order autoregressive model in (1), and is a known transition matrix. The initial value of parameter at , that is, , is assumed to be unknown with its uncertainty expressed by the following normal distribution , , where denotes a Gamma distribution with density function defined by and is a positive definite, symmetric matrix. The hyperparameters , , , and are assumed to be known. The errors and are assumed to be independent for all and . The joint distribution of is assumed to be spherically symmetric distribution with joint density function given bywhere is a positive measurable function and pdf and cdf of are represented by and , respectively. We further assume that the errors follow the normal distribution with mean and variance , where is a positive definite symmetric known matrix.

Thus, and . Thus, .

#### 3. Conditional Posterior Densities

The density function of the observation vector is given by In other words, the conditional distribution of given is normal with mean and covariance matrix . Since it is difficult to derive the explicit expressions for the marginal posterior densities of the parameters of the model, Gibbs sampler scheme can be used for estimating the marginal posterior densities. The conditional posterior densities of different parameters derived in this section can be utilised in Gibbs sampler scheme. The following joint distribution of given is obtained as follows:Let us write, ; , for

Theorem 1. *The conditional posterior density of is given by*

*Proof. *For , using expression (4), we getorwhere is the normalizing constant. For obtaining the value of , we havefor , where , and using expression (4) we getwhere the normalizing constant is obtained asfor , and utilizing expression (4), we getorwhere is the normalizing constant given by

*Remark 2. *It is interesting to observe from the above theorems that the conditional posterior density of given depends only upon and the conditional posterior density of given depends only upon whereas the conditional posterior density of given depends only upon . Therefore, the conditional posterior densities for depend upon (one or two) values of adjacent to and .

Now, we define

Theorem 3. *The conditional posterior of is obtained as*

*Proof. *Following (4), we observe thatThe normalizing constant is given by

*Remark 4. *The above derived expressions for the conditional posterior densities of , **, **, and involve integration with respect to . Sometimes, it may not be possible to obtain these integrals in neat form. In such situation, we consider as a parameter with prior density function and derive the conditional posterior density functions of given , given , and given The posterior conditional density functions of and can be straightaway obtained as follows:Hence, the conditional posterior density of is normal and conditional posterior density of is gamma.

Theorem 5. *The conditional posterior density of given is given by*

*Proof. *The joint density function of is given by Hence, the conditional posterior density function of is given as

#### 4. An Illustration

If follows -distribution with degrees of freedom and , , thenwhich is the pdf of a multivariate- distribution. For , the distribution reduces to the multivariate Cauchy distribution. For multivariate- distribution,Further, the conditional posterior density of given is and conditional posterior density of is gamma. Further,Hence, the conditional posterior density of is given by

#### 5. Implementation of Gibbs Sampler

Let the generated Gibbs samples be denoted by , . Then by employing Gibbs sampler posterior density of given can be estimated asRecall that the estimated posterior density of depends on the two values adjacent to . Hence, an estimate of is the mean of the estimated density obtained as , where is the number of iterations during implementation of the Gibbs sampler and and are the values of and based on . Then the fitted value of for our model is . If one is interested in one step ahead prediction of the unknown value , the Gibbs algorithm can be employed with the following modification. We consider as an additional parameter to the Gibbs sampler. The conditional density of is given by . Then Gibbs sampler can be run to find the following estimates of the marginal densities of and predictive density of , respectively:Under the squared error loss, can be predicted by the mean of the estimate of predictive density , which is given byAlternatively, one step ahead prediction via Gibbs sampler is given by . The above procedure can be adopted for predicting ; .

#### 6. Empirical Analysis

Stock price changes over a 3-month period are studied by analyzing records from database (of yahoo finance) documenting stocks of Amazon.com—an electronic commerce company. The data consists of daily latest intraday delayed opening prices, quoted in US dollars, taken between 6 July, 2015, and 2 October, 2015. The early diagnostics reflect that the data follows the random walk plus noise model which is an elementary form of state space model involving just one state variable. Specification of and in the model and knowledge of variances and covariances of the disturbance term are a prerequisite for implementing state space models. Based on the external facts and a preliminary examination of the data, we assumed that , , , , (a restrictive assumption), (restrictive assumption), , , . Analysis for model fitting is based on the first 60 observations. Subsequent 5 observations are excluded from the analysis and are used further in the predictive analysis. We then executed the Gibbs sampler for iterations, using parallel replications per iteration. A few initial iteration values in the chain were discarded to allow for dissipation of the initialized values called the burn-in values. The Kalman filter is fitted using . The longitudinal profile of the original 60 data points superimposed with the graph of Kalman filter model (via Gibbs sampler) is plotted in Figure 1.