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.     

Urn Sampling and the Proportional Hazard Model

This paper investigates the urn sampling analogue for the score statistic relating survival to covariates assuming a proportional hazard model. The exact permutation distribution can be calculated as well as the exact low order moments for arbitrary censoring patterns. The asymptotic distribution of the score statistic is an easy consequence. The method is naturally extended to deal with the multivariate case, time varying covariates and interval censoring. Finally the relationship between the censoring process, the survival times and covariates are studied considering different reference sets for the distribution of the score statistic. Some assumptions about the censoring process are investigated and as a consequence the effect of censoring is clarified.     

Regularized proportional odds models

The proportional odds model (POM) is commonly used in regression analysis to predict the outcome for an ordinal response variable. The maximum likelihood estimation (MLE) approach is typically used to obtain the parameter estimates. The likelihood estimates do not exist when the number of parameters, p, is greater than the number of observations n. The MLE also does not exist if there are no overlapping observations in the data. In a situation where the number of parameters is less than the sample size but p is approaching to n, the likelihood estimates may not exist, and if they exist they may have quite large standard errors. An estimation method is proposed to address the last two issues, i.e. complete separation and the case when p approaches n, but not the case when p>n. The proposed method does not use any penalty term but uses pseudo-observations to regularize the observed responses by downgrading their effect so that they become close to the underlying probabilities. The estimates can be computed easily with all commonly used statistical packages supporting the fitting of POMs with weights. Estimates are compared with MLE in a simulation study and an application to the real data.     

Generalizing the standardized hazard ratio to multivariate proportional hazards regression,with an application to clinical~genomic studies

The standardized hazard ratio for univariate proportional hazards regression is generalized as a scalar to multivariate proportional hazards regression. Estimators of the standardized log hazard ratio are developed, with corrections for bias and for regression to the mean in high-dimensional analyses. Tests of point and interval null hypotheses and confidence intervals are constructed. Cohort sampling study designs, commonly used in prospective–retrospective clinical genomic studies, are accommodated.     

A New Estimator for Cox Proportional Hazard Regression Model in Presence of Collinearity

We propose a new approach to estimate the parameters of the Cox proportional hazards model in the presence of collinearity. Generally, a maximum partial likelihood estimator is used to estimate parameters for the Cox proportional hazards model. However, the maximum partial likelihood estimators can be seriously affected by the presence of collinearity since the parameter estimates result in large variances.

In this study, we develop a Liu-type estimator for Cox proportional hazards model parameters and compare it with a ridge regression estimator based on the scalar mean squared error (MSE). Finally, we evaluate its performance through a simulation study.     

Development of predictive signatures for treatment selection in precision medicine with survival outcomes

For survival endpoints in subgroup selection, a score conversion model is often used to convert the set of biomarkers for each patient into a univariate score and using the median of the univariate scores to divide the patients into biomarker‐positive and biomarker‐negative subgroups. However, this may lead to bias in patient subgroup identification regarding the 2 issues: (1) treatment is equally effective for all patients and/or there is no subgroup difference; (2) the median value of the univariate scores as a cutoff may be inappropriate if the sizes of the 2 subgroups are differ substantially. We utilize a univariate composite score method to convert the set of patient's candidate biomarkers to a univariate response score. We propose applying the likelihood ratio test (LRT) to assess homogeneity of the sampled patients to address the first issue. In the context of identification of the subgroup of responders in adaptive design to demonstrate improvement of treatment efficacy (adaptive power), we suggest that subgroup selection is carried out if the LRT is significant. For the second issue, we utilize a likelihood‐based change‐point algorithm to find an optimal cutoff. Our simulation study shows that type I error generally is controlled, while the overall adaptive power to detect treatment effects sacrifices approximately 4.5% for the simulation designs considered by performing the LRT; furthermore, the change‐point algorithm outperforms the median cutoff considerably when the subgroup sizes differ substantially.     

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.     

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.     

Migration pattern in Bangladesh: a covariate-dependent Markov model

Internal migration is one of the major components of rapid and unplanned growth of towns and cities especially in the developing countries. This paper describes the transition pattern of internal out migration in Bangladesh and some sociodemographic factors influencing such migration in the country using a covariate-dependent Markov model. Four types of migration behavior namely, rural to rural, rural to urban, urban to rural and urban to urban are under consideration of this paper. Defining two discrete states, urban and rural, each of such transition can be characterized by a stochastic process; hence we use a two-state Markov chain for this purpose. We find that age, sex, division and reason of migration are significantly associated with internal migration in Bangladesh. The major findings include that any type of migration, rural to rural, rural to urban, urban to rural and urban to urban, mostly take place at the ages of 15–30 as well as at the ages of 0–15; females have higher odds than males to make a migration; Dhaka, Rajshahi and Chittagong divisions have remarkably higher migration rate as compared to Barisal and Sylhet division; and the professional reason is the main reason for rural to urban migration.