Interval estimation in discriminant analysis for large dimension

This paper is concerned with the interval estimation for the log odds of the posterior probability that the observation vector belongs to one of two homoscedastic multivariate normal distributions (Π₁ and Π₂). We give the limiting distribution of the unbiased estimator for the log odds as the sample sizes and the dimension jointly tend to infinity, and approximate the confidence interval based on the asymptotic distribution. Small-scale simulations are performed to check the precision of the approximation.     

On square ordinal contingency tables: a comparison of social class and income mobility for the same individuals

Summary.  Square contingency tables with matching ordinal rows and columns arise in particular as empirical transition matrices and the paper considers these in the context of social class and income mobility tables. Such tables relate the socio-economic position of parents to the socio-economic position of their child in adulthood. The level of association between parental and child socio-economic position is taken as a measure of mobility. Several approaches to analysis are described and illustrated by UK data in which interest focuses on comparisons of social class and income mobility tables that are derived from the same individuals. Account is taken of the use of the same individuals in the two tables. Additionally comparisons over time are considered.     

Data-dependent probability matching priors of the second order

In this paper, we consider the preliminary test approach for the estimation of the regression parameter in a multiple regression model under a multicollinearity situation. The preliminary test two-parameter estimators based on the Wald (W), likelihood ratio, and Lagrangian multiplier tests are given, when it is suspected that the regression parameter may be restricted to a subspace and the regression error is distributed with multivariate Student's t distribution. The bias and mean square error of the proposed estimators are derived and compared. The conditions of superiority of the proposed estimators are obtained. Finally, we conclude that the optimum choice of the level of significance becomes the traditional choice by using the Wald test.     

Estimation of Causal Odds of Concordance using the Aalen Additive Model

A simple summary of a treatment effect is attractive, which is part of the explanation of the success of the Cox model when analysing time‐to‐event data since the relative risk measure is such a convenient summary measure. In practice, however, the Cox model may fail to give a reasonable fit, very often because of time‐changing treatment effect. The Aalen additive hazards model may be a good alternative as time‐changing effects are easily modelled within this model, but results are then evidently more complicated to communicate. In such situations, the odds of concordance measure (OC) is a convenient way of communicating results, and recently Martinussen & Pipper (2012) showed how a variant of the OC measure may be estimated based on the Aalen additive hazards model. In this study, we propose an estimator that should be preferred in observational studies as it always estimates the causal effect on the chosen scale, only assuming that there are no un‐measured confounders. The resulting estimator is shown to be consistent and asymptotically normal, and an estimator of its limiting variance is provided. Two real applications are provided.     

Jackknifing in Categorical Data Analysis

Estimation of nonlinear functions of a multinomial parameter vector is necessary in many categorical data problems. The first and second order jackknife are explored for the purpose of reduction of bias. The second order jackknife of a function g(.) of a multinomial parameter is shown to be asymptotically normal if all second order partials ?²g( p )?dp_i?p_j obey a Hölder condition with exponent α>1/2. Numerical results for the estimation of the log odds ratio in a 2times2 table demonstrate the efficiency of the jackknife method for reduction of mean-square-error and the construction of approximate confidence intervals.     

Comparing the relative efficiency of three sequential procedures

We consider the sequential procedures developed by Robbins and Siegmund (1974), Louis (1975) and Zoubeidi (1992) for comparing the means of two treatments. We let the procedures have equal power functions and compare their Bayes and minimax risks using the invariance property of their power functions. For each of several formulations of the problem we determine the most relatively efficient procedure and compute its expected total sample size.     

Acute Middle Ear Infection in Small Children: A Bayesian Analysis Using Multiple Time Scales

The study is based on a sample of 965 children living in Oulu region (Finland), who were monitored for acute middle ear infections from birth to the age of two years. We introduce a nonparametrically defined intensity model for ear infections, which involves both fixed and time dependent covariates, such as calendar time, current age, length of breast-feeding time until present, or current type of day care. Unmeasured heterogeneity, which manifests itself in frequent infections in some children and rare in others and which cannot be explained in terms of the known covariates, is modelled by using individual frailty parameters. A Bayesian approach is proposed to solve the inferential problem. The numerical work is carried out by Monte Carlo integration (Metropolis-Hastings algorithm).     

Modelling two way contingency tables with recursive logits and odds ratios

In this work a new type of logits and odds ratios, which includes as special cases the continuation and the reverse-continuation logits and odds ratios, are defined. We prove that these logits and odds ratios define a parameterization of the joint probabilities of a two way contingency table. The problem of testing equality and inequality constraints on these logits and odds ratios is examined with particular regard to two new hypotheses of monotone dependence. Work partially supported by a MIUR2004 grant. Preliminary findings have been presented at SIS (Società Italiana di Statistica) Annual Meeting, Torino, 2006.     

Estimation of Dynamic Bivariate Mixture Models

This note compares a Bayesian Markov chain Monte Carlo approach implemented by Watanabe with a maximum likelihood ML approach based on an efficient importance sampling procedure to estimate dynamic bivariate mixture models. In these models, stock price volatility and trading volume are jointly directed by the unobservable number of price-relevant information arrivals, which is specified as a serially correlated random variable. It is shown that the efficient importance sampling technique is extremely accurate and that it produces results that differ significantly from those reported by Watanabe.     

Bayesian estimation for the Pareto income distribution

The Bayes estimators of the Gini index, the mean income and the proportion of the population living below a prescribed income level are obtained in this paper on the basis of censored income data from a pareto income distribution. The said estimators are obtained under the assumptions of a two-parameter exponential prior distribution and the usual squared error loss function. This work is also extended to the case when the income data are grouped and the exact incomes for the individuals in the population are not available. The method for the assessment of the hyperparameters is also outlined. Finally, the results are generalized for the doubly truncated gamma prior distribution. Now deceased.