A Coefficient of Agreement Adjusted for Bias in Paired Ordered Categorical Data

The quality of ordered categorical recordings is determined from repeated measurements on the same subject in order to assess the level of agreement between raters, scales or occasions. The presented rating-invariant method for ordered categorical data provides means of analysing the quality of single-item rating scales, irrespective of the number of possible response values and the marginal distributions. Marginal heterogeneity implies systematic disagreement, so-called bias. An augmented ranking approach is the basis for the separation of inter-rater disagreement into systematic and random components. Correlation between pairs of augmented rank values provides a measure of agreement to the best common ordering of paired classifications, given inter-rater bias. The essential differences in interpretation and applicability of the proposed coefficient of agreement and the Spearman rank-order correlation for ordered categorical data are discussed.     

Fitting Regression Models to Ordinal Data

This paper deals with the analysis of ordinal data by means of a threshold model. Maximum likelihood estimation is discussed and two examples are used to illustrate the methods.     

Compound Regression Models for Ordered Categorical Data

A general class of sequential models for the analysis of ordered categorical variables is developed and discussed. The models apply if the ordinal response may be subdivided into two or more meaningful sets of response categories. The parametrization explicitly makes use of this subdivision. The models furnish a linear alternative to non-linear models which incorporate a scale parameter. They are shown to be special cases of multivariate generalized linear models. Applications are discussed with the use of several examples.     

A Symmetric Extended Logistic Model with Applications to Experimental Toxicity Data

The fit of the logit and probit models for quantal response data can be improved by embedding these classical models within a richer parametric family indexed by one or two shape parameters. In this paper, a symmetric extended logistic model indexed by a shape parameter λ is discussed with application to dose response curves. The usual maximum likelihood method is employed to estimate the parameters of the model. The need to include the shape parameter λ is illustrated by analyzing a set of real experimental data and comparing the fit of the extended logistic model to those obtained by the standard logit and probit models.     

Statistical Methods for Ordered Categorical Data Based on a Constrained Odds Model

The proportional odds model provides a powerful tool for analysing ordered categorical data and setting sample size, although for many clinical trials its validity is questionable. The purpose of this paper is to present a new class of constrained odds models which includes the proportional odds model. The efficient score and Fisher's information are derived from the profile likelihood for the constrained odds model. These results are new even for the special case of proportional odds where the resulting statistics define the Mann‐Whitney test. A strategy is described involving selecting one of these models in advance, requiring assumptions as strong as those underlying proportional odds, but allowing a choice of such models. The accuracy of the new procedure and its power are evaluated.     

Threshold Model for Misclassified Binary Responses with Applications to Animal Breeding

A Bayesian procedure for misclassified binary data was developed. An animal breeding simulation indicated that, when error of classification was ignored, the variance between clusters was inferred incorrectly. Data were reanalyzed assuming that the probability of misclassification was either known or unknown. In the first case, input parameter values were recovered in the analysis. When the probability was unknown, there was a slight bias; the true probability of misclassification and the true number of miscoded observations appeared within high credibility regions. An analysis of fertility in dairy cows is presented.     

Empirical Likelihood Semiparametric Regression Analysis for Longitudinal Data

A semiparametric regression model for longitudinal data is considered.The empirical likelihood method is used to estimate the regressioncoefficients and the baseline function, and to construct confidenceregions and intervals. It is proved that the maximum empiricallikelihood estimator of the regression coefficients achievesasymptotic efficiency and the estimator of the baseline functionattains asymptotic normality when a bias correction is made.Two calibrated empirical likelihood approaches to inferencefor the baseline function are developed. We propose a groupwiseempirical likelihood procedure to handle the inter-series dependencefor the longitudinal semiparametric regression model, and employbias correction to construct the empirical likelihood ratiofunctions for the parameters of interest. This leads us to provea nonparametric version of Wilks' theorem. Compared with methodsbased on normal approximations, the empirical likelihood doesnot require consistent estimators for the asymptotic varianceand bias. A simulation compares the empirical likelihood andnormal-based methods in terms of coverage accuracies and averageareas/lengths of confidence regions/intervals.     

A Dynamic Logistic Regression Model for Multiple Spell Failure Time Data

The paper deals with discrete-time regression models to analyze multistate-multiepisode failure time data. The covariate process may include fixed and external as well as internal time dependent covariates. The effects of the covariates may differ among different kinds of failures and among successive episodes. A dynamic form of the logistic regression model is investigated and maximum likelihood estimation of the regression coefficients is discussed. In the last section we give an application of the model to the analysis of survival time after breast cancer operation.     

A Variance Components Model for Analyzing Repeated Measures Designs with Non-Stationary Error Structure

When data are collected in the form of multiple measurements on several subjects, they are often analyzed as repeated measures data with some stationary error structure assumed for the errors. For data with non-stationary error structure, the multivariate model is often used. The multivariate model imposes restrictions that are often not met in practice by data of such type. At the same time, they ignore valuable information in the data that are related to time dependencies and time relations. In this paper, we propose a model that is a reparametrization of the multivariate model and is suitable to analyze general repeated measures designs with non-stationary error structure. The model is shown to be a variance components model whose components are estimated using the method of maximum likelihood. Several other properties of the model are derived and discussed including tests of significance. Finally, an example on neurological data is included to demonstrate its application in biological sciences.     

On Generalized Log-Logistic Model for Censored Survival Data

In the analysis of survival data with parametric models, it is well known that the Weibull model is not suitable for modeling cases where the hazard rate is non-monotonic. For such cases, log-logistic model is frequently used. However, due to the symmetric property of the log-logistic model, it may be poor for the cases where the hazard rate is skewed or heavily tailed. In this paper, we suggest a generalization of the log-logistic model by introducing a shape parameter. This generalized model is then applied to fit the lung cancer data of Prentice (1973). The results seem to improve over those obtained by using the log-logistic model.     

Analyzing Incomplete Data Subject to a Threshold using Empirical Likelihood Methods: An Application to a Pneumonia Risk Study in an ICU Setting

Summary .  The initial detection of ventilator-associated pneumonia (VAP) for inpatients at an intensive care unit needs composite symptom evaluation using clinical criteria such as the clinical pulmonary infection score (CPIS). When CPIS is above a threshold value, bronchoalveolar lavage (BAL) is performed to confirm the diagnosis by counting actual bacterial pathogens. Thus, CPIS and BAL results are closely related and both are important indicators of pneumonia whereas BAL data are incomplete. To compare the pneumonia risks among treatment groups for such incomplete data, we derive a method that combines nonparametric empirical likelihood ratio techniques with classical testing for parametric models. This technique augments the study power by enabling us to use any observed data. The asymptotic property of the proposed method is investigated theoretically. Monte Carlo simulations confirm both the asymptotic results and good power properties of the proposed method. The method is applied to the actual data obtained in clinical practice settings and compares VAP risks among treatment groups.     

A Life Table Regression Analysis for Complex Survey Data

A methodology for obtaining maximum likelihood estimates of life table regression coefficients from complex survey data is presented. Certain of the issues of writing a likelihood for survey data are presented and discussed. The proposed methodology includes consideration of the sampling design in any inference by using design based variance estimates for the parameters. An example is given using data from the 1973 United States National Survey of Family Growth.     

Analyzing Variety by Environment Data Using Multiplicative Mixed Models and Adjustments for Spatial Field Trend

The recommendation of new plant varieties for commercial use requires reliable and accurate predictions of the average yield of each variety across a range of target environments and knowledge of important interactions with the environment. This information is obtained from series of plant variety trials, also known as multi-environment trials (MET). Cullis, Gogel, Verbyla, and Thompson (1998) presented a spatial mixed model approach for the analysis of MET data. In this paper we extend the analysis to include multiplicative models for the variety effects in each environment. The multiplicative model corresponds to that used in the multivariate technique of factor analysis. It allows a separate genetic variance for each environment and provides a parsimonious and interpretable model for the genetic covariances between environments. The model can be regarded as a random effects analogue of AMMI (additive main effects and multiplicative interactions). We illustrate the method using a large set of MET data from a South Australian barley breeding program.     

Ascertainment: An Overwiew of the Classical Segregation Analysis Model for Independent Sibships

Several different methodologies for parameter estimation under various ascertainment sampling schemes have been proposed in the past. In this article, some of the methodologies that have been proposed for independent sibships under the classical segregation analysis model are synthesized, and the general likelihoods derived for single, multiple and complete ascertainment. The issue of incorporating the sibship size distribution into the analysis is addressed, and the effect of conditioning the likelihood on the observed sibship sizes is discussed. It is shown that when the number of probands in a sibship is not specified, the corresponding likelihood can be used for a broader class of ascertainment schemes than is subsumed by the classical model.     

Maximum Likelihood Methods for Fitting the Burr Type XII Distribution to Life Test Data

This paper describes mathematical and computational methodology for estimating the parameters of the Burr Type XII distribution by the method of maximum likelihood. Expressions for the asymptotic variances and covariances of the parameter estimates are given, and the modality of the log-likelihood and conditional log-likelihood functions is analyzed. As a result of this analysis for various a priori known and unknown parameter combinations, conditions are given which guarantee that the parameter estimates obtained will, indeed, be maximum likelihood estimates. An efficient numerical method for maximizing the conditional log-likelihood function is described, and mathematical expressions are given for the various numerical approximations needed to evaluate the expressions given for the asymptotic variances and covariances of the parameter estimates. The methodology discussed is applied in a numerical example to life test data arising in a clinical setting.     

Zero-Altered and other Regression Models for Count Data with Added Zeros

On occasion, generalized linear models for counts based on Poisson or overdispersed count distributions may encounter lack of fit due to disproportionately large frequencies of zeros. Three alternative types of regression models that utilize all the information and explicitly account for excess zeros are examined and given general formulations. A simple mechanism for added zeros is assumed that directly motivates one type of model, here called the added-zero type, particular forms of which have been proposed independently by D. LAMBERT (1992) and in unpublished work by the author. An original regression formulation (the zero-altered model) is presented as a reduced form of the two-part model for count data, which is also discussed. It is suggested that two-part models be used to aid in development of an added-zero model when the latter is thought to be appropriate.     

A bivariate pseudolikelihood for incomplete longitudinal binary data with nonignorable nonmonotone missingness

For analyzing longitudinal binary data with nonignorable and nonmonotone missing responses, a full likelihood method is complicated algebraically, and often requires intensive computation, especially when there are many follow-up times. As an alternative, a pseudolikelihood approach has been proposed in the literature under minimal parametric assumptions. This formulation only requires specification of the marginal distributions of the responses and missing data mechanism, and uses an independence working assumption. However, this estimator can be inefficient for estimating both time-varying and time-stationary effects under moderate to strong within-subject associations among repeated responses. In this article, we propose an alternative estimator, based on a bivariate pseudolikelihood, and demonstrate in simulations that the proposed method can be much more efficient than the previous pseudolikelihood obtained under the assumption of independence. We illustrate the method using longitudinal data on CD4 counts from two clinical trials of HIV-infected patients.     

An empirical Bayes method for estimating epistatic effects of quantitative trait loci

Summary .   The genetic variance of a quantitative trait is often controlled by the segregation of multiple interacting loci. Linear model regression analysis is usually applied to estimating and testing effects of these quantitative trait loci (QTL). Including all the main effects and the effects of interaction (epistatic effects), the dimension of the linear model can be extremely high. Variable selection via stepwise regression or stochastic search variable selection (SSVS) is the common procedure for epistatic effect QTL analysis. These methods are computationally intensive, yet they may not be optimal. The LASSO (least absolute shrinkage and selection operator) method is computationally more efficient than the above methods. As a result, it has been widely used in regression analysis for large models. However, LASSO has never been applied to genetic mapping for epistatic QTL, where the number of model effects is typically many times larger than the sample size. In this study, we developed an empirical Bayes method (E-BAYES) to map epistatic QTL under the mixed model framework. We also tested the feasibility of using LASSO to estimate epistatic effects, examined the fully Bayesian SSVS, and reevaluated the penalized likelihood (PENAL) methods in mapping epistatic QTL. Simulation studies showed that all the above methods performed satisfactorily well. However, E-BAYES appears to outperform all other methods in terms of minimizing the mean-squared error (MSE) with relatively short computing time. Application of the new method to real data was demonstrated using a barley dataset.     

An evolutionary model for maximum likelihood alignment of DNA sequences

Summary Most algorithms for the alignment of biological sequences are not derived from an evolutionary model. Consequently, these alignment algorithms lack a strong statistical basis. A maximum likelihood method for the alignment of two DNA sequences is presented. This method is based upon a statistical model of DNA sequence evolution for which we have obtained explicit transition probabilities. The evolutionary model can also be used as the basis of procedures that estimate the evolutionary parameters relevant to a pair of unaligned DNA sequences. A parameter-estimation approach which takes into account all possible alignments between two sequences is introduced; the danger of estimating evolutionary parameters from a single alignment is discussed.