Generalized Laplacian approximations in Bayesian inference

This paper presents a new Laplacian approximation to the posterior density of η = g(θ). It has a simpler analytical form than that described by Leonard et al. (1989). The approximation derived by Leonard et al. requires a conditional information matrix R_η to be positive definite for every fixed η. However, in many cases, not all R_η are positive definite. In such cases, the computations of their approximations fail, since the approximation cannot be normalized. However, the new approximation may be modified so that the corresponding conditional information matrix can be made positive definite for every fixed η. In addition, a Bayesian procedure for contingency-table model checking is provided. An example of cross-classification between the educational level of a wife and fertility-planning status of couples is used for explanation. Various Laplacian approximations are computed and compared in this example and in an example of public school expenditures in the context of Bayesian analysis of the multiparameter Fisher-Behrens problem.     

The analysis of incomplete data using stochastic covariates

In the development of many diseases there are often associated random variables which continuously reflect the progress of a subject towards the final expression of the disease (failure). At any given time these processes, which we call stochastic covariates, may provide information about the current hazard and the remaining time to failure. Likewise, in situations when the specific times of key prior events are not known, such as the time of onset of an occult tumour or the time of infection with HIV-1, it may be possible to identify a stochastic covariate which reveals, indirectly, when the event of interest occurred. The analysis of carcinogenicity trials which involve occult tumours is usually based on the time of death or sacrifice and an indicator of tumour presence for each animal in the experiment. However, the size of an occult tumour observed at the endpoint represents data concerning tumour development which may convey additional information concerning both the tumour incidence rate and the rate of death to which tumour-bearing animals are subject. We develop a stochastic model for tumour growth and suggest different ways in which the effect of this growth on the hazard of failure might be modelled. Using a combined model for tumour growth and additive competing risks of death, we show that if this tumour size information is used, assumptions concerning tumour lethality, the context of observation or multiple sacrifice times are no longer necessary in order to estimate the tumour incidence rate. Parametric estimation based on the method of maximum likelihood is outlined and is applied to simulated data from the combined model. The results of this limited study confirm that use of the stochastic covariate tumour size results in more precise estimation of the incidence rate for occult tumours.     

On the asymptotic normality of the L1- and L2-errors in histogram density estimation

The L₁ and L₂-errors of the histogram estimate of a density f from a sample X₁,X₂,…,Xn using a cubic partition are shown to be asymptotically normal without any unnecessary conditions imposed on the density f. The asymptotic variances are shown to depend on f only through the corresponding norm of f. From this follows the asymptotic null distribution of a goodness-of-fit test based on the total variation distance, introduced by Györfi and van der Meulen (1991). This note uses the idea of partial inversion for obtaining characteristic functions of conditional distributions, which goes back at least to Bartlett (1938).     

Bayesian analysis of outlier problems using divergence measures

A Bayesian approach is presented for detecting influential observations using general divergence measures on the posterior distributions. A sampling-based approach using a Gibbs or Metropolis-within-Gibbs method is used to compute the posterior divergence measures. Four specific measures are proposed, which convey the effects of a single observation or covariate on the posterior. The technique is applied to a generalized linear model with binary response data, an overdispersed model and a nonlinear model. An asymptotic approximation using Laplace method to obtain the posterior divergence is also briefly discussed.     

India's demand for international reserve and monetary disequilibrium: Reserve adequacy under floating regime

The main aim of this study is to investigate India's demand for international reserve by focusing on the role of national monetary disequilibrium and to present new benchmarks for assessing the adequacy of international reserves. We assessed India's position in terms of reserve adequacy and found that India is well placed and has sufficient stock of international reserves to meet the minimum adequacy requirements. Also, the results reveal that the central bank is holding substantial excess reserves and the related opportunity cost (1.5% of GDP) appears to be quite considerable. Further, the estimates of reserve demand function suggest that scale of foreign trade, uncertainty and profitability considerations play significant role in determining India's long-term reserve demand policies. More importantly, validating the monetary approach to balance of payment, our results show that national monetary disequilibrium does play a crucial role in short-run reserve movements. An excess of money demand (supply) induces an inflow (outflow) of international reserves with an elasticity of 0.56 which also implies that Reserve Bank of India responds to correct the domestic money market disequilibrium; and did not just leave it completely on the mercy of reserve inflows.     

A multivariate version of Gini's rank association coefficient

In this paper, we introduce a multivariate generalization of the population version of Gini's rank association coefficient, giving a response to this open question posed in [4]. We also study some properties of this version, present the corresponding results for the sample statistic, and provide several examples.     

On Christensen's conjecture

In this paper we give solution to a conjecture appearing in Christensen (1997, p. 360) in relation to the definitions of standardized residuals in loglinear models. Work done when he was visiting as the Distinguished Lukacs Professor on leave from the Complutense University of Madrid. An erratum to this article is available at .     

Two-stage Huber estimation

In this paper we propose a new robust estimator in the context of two-stage estimation methods directed towards the correction of endogeneity problems in linear models. Our estimator is a combination of Huber estimators for each of the two stages, with scale corrections implemented using preliminary median absolute deviation estimators. In this way we obtain a two-stage estimation procedure that is an interesting compromise between concerns of simplicity of calculation, robustness and efficiency. This method compares well with other possible estimators such as two-stage least-squares (2SLS) and two-stage least-absolute-deviations (2SLAD), asymptotically and in finite samples. It is notably interesting to deal with contamination affecting more heavily the distribution tails than a few outliers and not losing as much efficiency as other popular estimators in that case, e.g. under normality. An additional originality resides in the fact that we deal with random regressors and asymmetric errors, which is not often the case in the literature on robust estimators.     

Evaluation of multilevel decision trees

The evaluation of decision trees under uncertainty is difficult because of the required nested operations of maximizing and averaging. Pure maximizing (for deterministic decision trees) or pure averaging (for probability trees) are both relatively simple because the maximum of a maximum is a maximum, and the average of an average is an average. But when the two operators are mixed, no simplification is possible, and one must evaluate the maximization and averaging operations in a nested fashion, following the structure of the tree. Nested evaluation requires large sample sizes (for data collection) or long computation times (for simulations).     

Minimum contrast estimation of random processes based on information of second and third orders

This paper presents the results on consistency and asymptotic normality of a class of minimum contrast estimators for random processes with short- or long-range dependence based on the second- and third-order cumulant spectra. Asymptotic properties of sample spectral functionals of second and third orders, which are of independent interest in view of their possible use for nonparametric estimation of processes with short- or long-range dependence, are also provided.