Considering the Concordant Observations in Likelihood Ratio Test for Paired Binary Data

In the formula of the likelihood ratio test on fourfold tables with matched pairs of binary data, only the two parts b and c, which represent changes, are considered; the retained parts a and d, which represent concordant observations, are not included. To develop the test by considering all the four parts and the mixture distribution of likelihood ratio chi-squares, a formula based on the entire sample is proposed. The revised formula is the same as the unrevised one when a + d is zero. The revised test is more valid than the revised McNemar's test in most cases.     

Inference of Genetic Networks from Time Course Expression Data Using Functional Regression with Lasso Penalty

Statistical inference of genetic regulatory networks is essential for understanding temporal interactions of regulatory elements inside the cells. In this work, we propose to infer the parameters of the ordinary differential equations using the techniques from functional data analysis (FDA) by regarding the observed time course expression data as continuous-time curves. For networks with a large number of genes, we take advantage of the sparsity of the networks by penalizing the linear coefficients with a L ₁ norm. The ability of the algorithm to infer network structure is demonstrated using the cell-cycle time course data for Saccharomyces cerevisiae.     

The Poisson Generalized Linear Failure Rate Model

We formulate a new cure rate survival model by assuming that the number of competing causes of the event of interest has the Poisson distribution, and the time to this event has the generalized linear failure rate distribution. A new distribution to analyze lifetime data is defined from the proposed cure rate model, and its quantile function as well as a general expansion for the moments is derived. We estimate the parameters of the model with cure rate in the presence of covariates for censored observations using maximum likelihood and derive the observed information matrix. We obtain the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and present some ways to perform global influence analysis. The usefulness of the proposed cure rate survival model is illustrated in an application to real data.     

An Improved Test for Continuous Local Martingales

We present a new test for the “continuous martingale hypothesis”. That is, a test for the hypothesis that observed data are from a process which is a continuous local martingale. The basis of the test is an embedded random walk at first passage times, obtained from the well-known representation of a continuous local martingale as a continuous time-change of Brownian motion. With a variety of simulated diffusion processes our new test shows higher power than existing tests using either the crossing tree or the quadratic variation, including the situation where non-negligible drift is present. The power of the test in the presence of jumps is also explored with a variety of simulated jump diffusion processes. The test is also applied to two sequences of high-frequency foreign exchange trade-by-trade data. In both cases the continuous martingale hypothesis is rejected at times less than hourly and we identify significant dependence in price movements at these small scales.     

Statistical Inference on Partially Linear Additive Models with Missing Response Variables and Error-prone Covariates

This paper considers statistical inference for the partially linear additive models, which are useful extensions of additive models and partially linear models. We focus on the case where some covariates are measured with additive errors, and the response variable is sometimes missing. We propose a profile least-squares estimator for the parametric component and show that the resulting estimator is asymptotically normal. To construct a confidence region for the parametric component, we also propose an empirical-likelihood-based statistic, which is shown to have a chi-squared distribution asymptotically. Furthermore, a simulation study is conducted to illustrate the performance of the proposed methods.     

Operating Characteristics of a Subset Selection Procedure Applied to a Hybrid Randomized Response Model

The operating characteristics (OCs) of a subset ranking and selection procedure are derived for the hybrid randomized response model developed by Jia and McDonald (2009 Jia, F., McDonald, G. (2009). Analyzing hybrid randomized response data with a binomial selection procedure. Commun. Statist. Theor. Meth. 38:784807.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The OCs include the probability of a correct P(CS), the individual selection probability γ_i, and the expected subset size E(S), under the slippage configuration or the equi-spaced configuration. An example comparing failure rates of contraceptive methods is used to illustrate the use of these new results.     

Bridging the gap: A likelihood function approach for the analysis of ranking data

ABSTRACT

In the parametric setting, the notion of a likelihood function forms the basis for the development of tests of hypotheses and estimation of parameters. Tests in connection with the analysis of variance stem entirely from considerations of the likelihood function. On the other hand, non parametric procedures have generally been derived without any formal mechanism and are often the result of clever intuition. In the present article, we propose a more formal approach for deriving tests involving the use of ranks. Specifically, we define a likelihood function motivated by characteristics of the ranks of the data and demonstrate that this leads to well-known tests of hypotheses. We also point to various areas of further exploration such as how co-variates may be incorporated.     

Stochastic EM algorithm of a finite mixture model from hurdle Poisson distribution with missing responses

ABSTRACT

In this article, a finite mixture model of hurdle Poisson distribution with missing outcomes is proposed, and a stochastic EM algorithm is developed for obtaining the maximum likelihood estimates of model parameters and mixing proportions. Specifically, missing data is assumed to be missing not at random (MNAR)/non ignorable missing (NINR) and the corresponding missingness mechanism is modeled through probit regression. To improve the algorithm efficiency, a stochastic step is incorporated into the E-step based on data augmentation, whereas the M-step is solved by the method of conditional maximization. A variation on Bayesian information criterion (BIC) is also proposed to compare models with different number of components with missing values. The considered model is a general model framework and it captures the important characteristics of count data analysis such as zero inflation/deflation, heterogeneity as well as missingness, providing us with more insight into the data feature and allowing for dispersion to be investigated more fully and correctly. Since the stochastic step only involves simulating samples from some standard distributions, the computational burden is alleviated. Once missing responses and latent variables are imputed to replace the conditional expectation, our approach works as part of a multiple imputation procedure. A simulation study and a real example illustrate the usefulness and effectiveness of our methodology.     

Consistency of information criteria for model selection with missing data

ABSTRACT

In this paper, we investigate the consistency of the Expectation Maximization (EM) algorithm-based information criteria for model selection with missing data. The criteria correspond to a penalization of the conditional expectation of the complete data log-likelihood given the observed data and with respect to the missing data conditional density. We present asymptotic properties related to maximum likelihood estimation in the presence of incomplete data and we provide sufficient conditions for the consistency of model selection by minimizing the information criteria. Their finite sample performance is illustrated through simulation and real data studies.     

Isotontic regression for interval censored survival data using an e-m algorithm

This paper develops a nonparametric model of the relationship between survival S and a dichotomous random variable X under the order constraint that P(X=1|S=s) is increasing (or decreasing) with s. The estimation procedure, called isotonic regression, has been studied in some depth for the case of uncensored data, but we give a methodology which is appropriate in the more general context of right, left, and interval censored data. An E-M Algorithm (Dempster et. al., 1977) is used for maximum likelihood estimation.