A Weighting Approach for GEE Analysis with Missing Data

We propose a new weighting (WT) method to handle missing categorical outcomes in longitudinal data analysis using generalized estimating equations (GEE). The proposed WT provides a valid GEE estimator when the data are missing at random (MAR), and has more stable weights and shows advantage in efficiency compared to the inverse probability weighing method in the presence of small observation probabilities. The WT estimator is similar to the stabilized weighting (SWT) estimator under mild conditions, but it is more stable and efficient than SWT when the associations of the outcome with the observation probabilities and the covariate are strong.     

Relations for Moments of Progressively Type-II Censored Order Statistics from Log-Logistic Distribution with Applications to Inference

In this article, we establish several recurrence relations for the single and product moments of progressively Type-II right censored order statistics from a log-logistic distribution. The use of these relations in a systematic recursive manner would enable the computation of all the means, variances and covariances of progressively Type-II right censored order statistics from the log-logistic distribution for all sample sizes n, effective sample sizes m, and all progressive censoring schemes (R ₁,…, R _m ). The results established here generalize the corresponding results for the usual order statistics due to Balakrishnan and Malik (1987 Balakrishnan , N. , Malik , H. J. ( 1987 ). Moments of order statistics from truncated log-logistic distribution . J. Statist. Plann. Infer. 17 : 251 – 267 .[Crossref], [Web of Science ®] , [Google Scholar]) and Balakrishnan et al. (1987 Balakrishnan , N. , Malik , H. J. , Puthenpura , S. ( 1987 ). Best linear unbiased estimation of location and scale parameters of the log-logistic distribution . Commun. Statist. Theor. Meth. 16 : 3477 – 3495 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The moments so determined are then utilized to derive best linear unbiased estimators for the scale- and location-scale log-logistic distributions. A comparison of these estimates with the maximum likelihood estimates is made through Monte Carlo simulation. The best linear unbiased predictors of progressively censored failure times is then discussed briefly. Finally, a numerical example is presented to illustrate all the methods of inference developed here.     

Use of Quasi-Ranges in the Estimation of Scale Parameter of a Symmetric Distribution

In this article, we consider the problem of best linear unbiased estimation and best linear invariant estimation of the scale parameter of a symmetric distribution using quasi-ranges is considered. We also prove a sufficient condition for the non negativity of the scale estimator obtained by the above method. Further, we obtain necessary and sufficient conditions for the derived estimators to be constant multiple of the sample range.     

The Recursive Kernel Distribution Function Estimator Based on Negatively and Positively Associated Sequences

Let {X _j , j ≥ 1} be a strictly stationary negatively or positively associated sequence of real valued random variables with unknown distribution function F(x). On the basis of the random variables {X _j , j ≥ 1}, we propose a smooth recursive kernel-type estimate of F(x), and study asymptotic bias, quadratic-mean consistency and asymptotic normality of the recursive kernel-type estimator under suitable conditions.     

Bernoulli Processes with Non-Positive Correlations

We deal with sequences of dependent Bernoulli variables with non positive correlations. Some special models for 1-dependent and 2-dependent 0–1 valued variables are analyzed, namely Bernoulli variables obtained by clipping a linear combination of iid variables. Formulas describing dependence of their correlation function on the clipping parameters are derived. The lower bound for the sum of autocorrelations of Bernoulli variables obtained by clipping a Gaussian process is provided.     

Normal Deviation and Poisson Approximation of a Security Market by the Geometric Markov Renewal Processes

We consider the geometric Markov renewal processes (GMRP) as a model for a security market. Normal deviations of the geometric Markov renewal processes for ergodic averaging and double averaging schemes are derived. We introduce Poisson averaging scheme for the geometric Markov renewal processes. European call option pricing formulas for GMRP are presented.     

General Multivariate Weibull Processes

Two multivariate stationary processes with general multivariate Weibull marginals are developed and studied. The joint distribution of the two adjacent events in the processes and the distributions of the finite sample minima as well as the geometric minima are derived. The characterization properties of these two processes are also proved.     

The Scale Invariant Wigner Spectrum Estimation of Gaussian Locally Self-Similar Processes

We study locally self-similar processes (LSSPs) in Silverman’s sense. By deriving the minimum mean-square optimal kernel within Cohen’s class counterpart of time–frequency representations, we obtain an optimal estimation for the scale invariant Wigner spectrum (SIWS) of Gaussian LSSPs. The class of estimators is completely characterized in terms of kernels, so the optimal kernel minimizes the mean-square error of the estimation. We obtain the SIWS estimation for two cases: global and local, where in the local case, the kernel is allowed to vary with time and frequency. We also introduce two generalizations of LSSPs: the locally self-similar chirp process and the multicomponent LSSP, and obtain their optimal kernels. Finally, the performance and accuracy of the estimation is studied via simulation.     

A Note on Wald's Type Chi-squared Tests of Fit

The Wald's method for constructing chi-squared tests of fit has been formulated more accurately. It is shown that Wald's type statistics will follow the central chi-squared distribution if and only if the limit covariance matrix of standardized frequencies will not depend on unknown parameters. Several examples that illustrate this important fact are presented. In particular, it is shown that the goodness-of-fit statistic developed by Moore and Stubblebine does not follow the chi-squared limit distribution, and, hence, cannot be used for testing multivariate normality.     

Modeling of covariance structures of random effects and random errors in linear mixed models

Abstract

In this paper, we discuss how to model the mean and covariancestructures in linear mixed models (LMMs) simultaneously. We propose a data-driven method to modelcovariance structures of the random effects and random errors in the LMMs. Parameter estimation in the mean and covariances is considered by using EM algorithm, and standard errors of the parameter estimates are calculated through Louis’ (1982 Louis, T.A. (1982). Finding observed information using the EM algorithm. J. Royal Stat. Soc. B 44:98–130. [Google Scholar]) information principle. Kenward’s (1987 Kenward, M.G. (1987). A method for comparing profiles of repeated measurements. Appl. Stat. 36:296–308.[Crossref], [Web of Science ®] , [Google Scholar]) cattle data sets are analyzed for illustration,and comparison to the literature work is made through simulation studies. Our numerical analysis confirms the superiority of the proposed method to existing approaches in terms of Akaike information criterion.