Linear calibra and conditional inference

Approximate conditional inference is developed for the linear calibration problem. It is shown that this problem can be transformed so that the primary parameter is an angle, the nuisance parameter is a radial distance, and the density is rotationally symmetric. Were the nuisance parameter known, exact location confidence intervals would be available by location of structural arguments. A confidence distribution is used to average out the nuisance parameter yielding an approximate confidence interval that involves a precision indicator derived from the radial distance. Some difficulties with the ordinary solution are avoided by the conditional procedure.     

Algorithms for approximate conditional inference

This paper presents a method for listing the sample space for a conditional distribution in a discrete generalized linear model. This tabulation is used in conjunction with saddlepoint methods to approximate the associated conditional probabilities. These probabilities are used to calculate conditional p-values.     

Fiducial inference for Birnbaum-Saunders distribution

In this paper, we consider fiducial inference for the unknown parameters of the Birnbaum-Saunders distribution. Two generalized fiducial distributions of the parameters are obtained. One is based on the inverse of the structural equation, and the fiducial estimates of the parameters are obtained by a simulation method. The other is based on the method of [Hannig J. Generalized fiducial inference via discretization. Stat. Sinica. 2013;23:489–514], then we use adaptive rejection Metropolis sampling to get the fiducial estimates. We compare the fiducial estimates with the maximum likelihood estimates and Bayesian estimates by simulations. Two real data sets are analysed for illustration.     

Inverse Gaussian Distribution for Modeling Conditional Durations in Finance

The durations between market activities such as trades and quotes provide useful information on the underlying assets while analyzing financial time series. In this article, we propose a stochastic conditional duration model based on the inverse Gaussian distribution. The non-monotonic nature of the failure rate of the inverse Gaussian distribution makes it suitable for modeling the durations in financial time series. The parameters of the proposed model are estimated by an efficient importance sampling method. A simulation experiment is conducted to check the performance of the estimators. These estimates are used to compute estimated hazard functions and to compare with the empirical hazard functions. Finally, a real data analysis is provided to illustrate the practical utility of the models.     

Discrepancy measures for item fit analysis in item response theory

Item response theory (IRT) models are commonly used in educational and psychological testing to assess the (latent) ability of examinees and the effectiveness of the test items in measuring this underlying trait. The focus of this paper is on the assessment of item fit for unidimensional IRT models for dichotomous items using a Bayesian method. This paper will illustrate and compare the effectiveness of several discrepancy measures, used within the posterior predictive model check procedure, in detecting misfitted items. The effectiveness of the different discrepancy measures are illustrated in a simulation study using artificially altered simulated data. Using the best discrepancy measure among those studied, this method was applied to real data coming from a mathematics placement exam.     

Approximate Studentization with marginal and conditional inference

Many inference problems lead naturally to a marginal or conditional measure of departure that depends on a nuisance parameter. As a device for first-order elimination of the nuisance parameter, we suggest averaging with respect to an exact or approximate confidence distribution function. It is shown that for many standard problems where an exact answer is available by other methods, the averaging method reproduces the exact answer. Moreover, for the gamma-mean problem, where the exact answer is not explicitly available, the averaging method gives results that agree closely with those obtained from higher-order asymptotic methods. Examples are discussed; detailed asymptotic calculations will be examined elsewhere.     

Approximate conditional inference and the linear functional model

Approximate conditional inference is developed for the slope parameter of the linear functional model with two variables. It is shown that the model can be transformed so that the slope parameter becomes an angle and nuisance parameters are radial distances. If the nuisance parameters are known an exact confidence interval based on a location-type conditional distribution is available for the angle. More gen¬erally, confidence distributions are used to average the conditional distribution over the nuisance parameters yielding an approximate conditional confidence interval that reflects the precision indicated by the data. An example is analyzed.     

Bayesian inference for circular distributions with unknown normalising constants

Very often, the likelihoods for circular data sets are of quite complicated forms, and the functional forms of the normalising constants, which depend upon the unknown parameters, are unknown. This latter problem generally precludes rigorous, exact inference (both classical and Bayesian) for circular data.Noting the paucity of literature on Bayesian circular data analysis, and also because realistic data analysis is naturally permitted by the Bayesian paradigm, we address the above problem taking a Bayesian perspective. In particular, we propose a methodology that combines importance sampling and Markov chain Monte Carlo (MCMC) in a very effective manner to sample from the posterior distribution of the parameters, given the circular data. With simulation study and real data analysis, we demonstrate the considerable reliability and flexibility of our proposed methodology in analysing circular data.     

Bayesian inference and prediction analysis of the power law process based on a natural conjugate prior

Under a natural conjugate prior with four hyperparameters, the importance sampling (IS) technique is applied to the Bayesian analysis of the power law process (PLP). Samples of the parameters of the PLP are obtained from IS. Based on these samples, not only the posterior analysis of parameters and some parameter functions in the PLP are performed conveniently, but also single-sample and two-sample prediction procedures are constructed easily. Furthermore, the sensitivity of the posterior mean of the parameter functions in the PLP is studied with respect to the hyperparameters of the natural conjugate prior and it can guide the selections of the hyperparameters directly. Coupled this sensitivity with the relations between the prior moments and the hyperparameters in the natural conjugate prior, it is possible to give directions about the selections of the prior moments to a certain degree. After some numerical experiments illustrate the rationality and feasibility of the proposed methods, an engineering example demonstrates its application.     

Conditional inference in finite population sampling under a calibration model

Several estimators, including the classical and the regression estimators of finite population mean, are compared, both theoretically and empirically, under a calibration model, where the dependent variable(y), and not the independent variable(x), can be observed for all units of the finite population. It is shown asymptotically that when conditioned on x, the bias of the classical estimator may be much smaller than that of the regression estimators; whereas when conditioned on y, the regression estimator may have much smaller conditional bias than the classical estimator. Since all the y's(not x's) can be observed, it seems appropriate to make comparison under the conditional distribution of each estimator with y fixed. In this case, the regression estimator has smaller variance, smaller conditional bias, and the conditional coverage probability closer to its nominal level     

A nonparametric Bayesian model for inference in related longitudinal studies

Summary.  We discuss a method for combining different but related longitudinal studies to improve predictive precision. The motivation is to borrow strength across clinical studies in which the same measurements are collected at different frequencies. Key features of the data are heterogeneous populations and an unbalanced design across three studies of interest. The first two studies are phase I studies with very detailed observations on a relatively small number of patients. The third study is a large phase III study with over 1500 enrolled patients, but with relatively few measurements on each patient. Patients receive different doses of several drugs in the studies, with the phase III study containing significantly less toxic treatments. Thus, the main challenges for the analysis are to accommodate heterogeneous population distributions and to formalize borrowing strength across the studies and across the various treatment levels. We describe a hierarchical extension over suitable semiparametric longitudinal data models to achieve the inferential goal. A nonparametric random-effects model accommodates the heterogeneity of the population of patients. A hierarchical extension allows borrowing strength across different studies and different levels of treatment by introducing dependence across these nonparametric random-effects distributions. Dependence is introduced by building an analysis of variance (ANOVA) like structure over the random-effects distributions for different studies and treatment combinations. Model structure and parameter interpretation are similar to standard ANOVA models. Instead of the unknown normal means as in standard ANOVA models, however, the basic objects of inference are random distributions, namely the unknown population distributions under each study. The analysis is based on a mixture of Dirichlet processes model as the underlying semiparametric model.     

Bayesian inference for the randomly censored Weibull distribution

In this paper, we consider the Bayesian inference of the unknown parameters of the randomly censored Weibull distribution. A joint conjugate prior on the model parameters does not exist; we assume that the parameters have independent gamma priors. Since closed-form expressions for the Bayes estimators cannot be obtained, we use Lindley's approximation, importance sampling and Gibbs sampling techniques to obtain the approximate Bayes estimates and the corresponding credible intervals. A simulation study is performed to observe the behaviour of the proposed estimators. A real data analysis is presented for illustrative purposes.     

Bayesian inference for bisexual galton-watson processes

In this work, an approach to the Bayesian estimation in a bisexual Galton-Watson process is considered. First we study an important parametric case assuming offspring distribution belonging to the bivariate series power family of distributions and then, we continue to investigate the nonparametric case. In both situations, Bayes estimators under weighted squared error loss function, for means, variances and covariance of the off spring distribution are obtained. For the superadditive case, the Bayes estimation of the asymptotic growth rate is also considered. Illustrative examples are given.     

A computer algebra package for approximate conditional inference

This paper presents a set of REDUCE procedures that make a number of existing higher-order asymptotic results available for both theoretical and practical research. Attention has been restricted to the context of exact and approximate inference for a parameter of interest conditionally either on an ancillary statistic or on a statistic partially sufficient for the nuisance parameter. In particular, the procedures apply to regression-scale models and multiparameter exponential families. Most of them support algebraic computation as well as numerical calculation for a given data set. Examples illustrate the code.     

Bayesian inference for random coefficient dynamic panel data models

We develop a hierarchical Bayesian approach for inference in random coefficient dynamic panel data models. Our approach allows for the initial values of each unit's process to be correlated with the unit-specific coefficients. We impose a stationarity assumption for each unit's process by assuming that the unit-specific autoregressive coefficient is drawn from a logitnormal distribution. Our method is shown to have favorable properties compared to the mean group estimator in a Monte Carlo study. We apply our approach to analyze energy and protein intakes among individuals from the Philippines.     

Bayesian inference for multinomial populations under stochastic ordering

In many parametric problems the use of order restrictions among the parameters can lead to improved precision. Our interest is in the study of several multinomial populations under the stochastic order restriction (SOR) for univariate situations. We use Bayesian methods to show that the SOR can lead to larger gains in precision than the method without the SOR when the SOR is reasonable. Unlike frequentist order restricted inference, our methodology permits analysis even when there is uncertainty about the SOR. Our method is sampling based, and we use simple and efficient rejection sampling. The Bayes factor in favor of the SOR is computed in a simple manner, and samples from the requisite posterior distributions are easily obtained. We use real data to illustrate the procedure, and we show that there is likely to be larger gains in precision under the SOR.     

Exact likelihood inference for two exponential populations based on a joint generalized Type-I hybrid censored sample

Following the work of Chen and Bhattacharyya [Exact confidence bounds for an exponential parameter under hybrid censoring. Comm Statist Theory Methods. 1988;17:1857–1870], several results have been developed regarding the exact likelihood inference of exponential parameters based on different forms of censored samples. In this paper, the conditional maximum likelihood estimators (MLEs) of two exponential mean parameters are derived under joint generalized Type-I hybrid censoring on the two samples. The moment generating functions (MGFs) and the exact densities of the conditional MLEs are obtained, using which exact confidence intervals are then developed for the model parameters. We also derive the means, variances, and mean squared errors of these estimates. An efficient computational method is developed based on the joint MGF. Finally, an example is presented to illustrate the methods of inference developed here.     

Sums of squares and products matrices for a non-full ranks hypothesis in the model of potthoff and roy

Formulae for sums of squares and products matrices, useful in testing a general linear hypothesis in the model of POTTHOFF and ROY, are given, Contrary to the customary approach, these formulae are expressed in original terms of the design matrices and the matrices formulating the hypothesis. They are applicable regardless of the ranks of the matrices involved, which allows to avoid a transformation of the hypothesis and a repara-metrization of the model.     

A general method of inference for two-parameter continuous distributions

Abstract

This article presents a general method of inference of the parameters of a continuous distribution with two unknown parameters. Except in a few distributions such as the normal distribution, the classical approach fails in this context to provide accurate inferences with small samples.Therefore, by taking the generalized approach to inference (cf. Weerahandi, 1995 Weerahandi, S. (1995). Exact Statistical Methods for Data Analysis. New York: Springer Verlag. [Google Scholar]), in this article we present a general method of inference to tackle practically useful two-parameter distributions such as the gamma distribution as well as distributions of theoretical interest such as the two-parameter uniform distribution. The proposed methods are exact in the sense that they are based on exact probability statements and exact expected values. The advantage of taking the generalized approach over the classical approximate inferences is shown via simulation studies.

This article has the potential to motivate much needed further research in non normal regressions, multiparameter problems, and multivariate problems for which basically there are only large sample inferences available. The approach that we take should pave the way for researchers to solve a variety of non normal problems, including ANOVA and MANOVA problems, where even the Bayesian approach fails. In the context of testing of hypotheses, the proposed method provides a superior alternative to the classical generalized likelihood ratio method.