A nonlinear mixed-effects model for multivariate longitudinal data with partially observed outcomes with application to HIV disease dynamics

The measurable multiple bio-markers for a disease are used as indicators for studying the response variable of interest in order to monitor and model disease progression. However, it is common for subjects to drop out of the studies prematurely resulting in unbalanced data and hence complicating the inferences involving such data. In this paper we consider a case where data are unbalanced among subjects and also within a subject because for some reason only a subset of the multiple outcomes of the response variable are observed at any one occasion. We propose a nonlinear mixed-effects model for the multivariate response variable data and derive a joint likelihood function that takes into account the partial dropout of the outcomes of the response variable. We further show how the methodology can be used in the estimation of the parameters that characterise HIV disease dynamics. An approximation technique of the parameters is also given and illustrated using a routine observational HIV dataset.     

Application of copulas to multivariate control charts

The most popular multivariate control chart for monitoring the mean of a distribution is probably the Hotelling T² rule. Unfortunately, this rule relies on the assumption that the distribution under control is Gaussian, which is rarely true in practice. The objective of this paper is to propose a new approach for the non-normal multivariate case. It consists in the construction of a tolerance region obtained from a density level set estimation. The method follows a “plug-in” approach in which the density of the observations is previously estimated. This estimation is conducted using copulas modeling, an increasingly popular tool in multivariate modeling.     

A multivariate two-factor skew model

It is well known that many data, such as the financial or demographic data, exhibit asymmetric distributions. In recent years, researchers have concentrated their efforts to model this asymmetry. Skew normal model is one of such models that are skew and yet possess many properties of the normal model. In this paper, a new multivariate skew model is proposed, along with its statistical properties. It includes the multivariate normal distribution and multivariate skew normal distribution as special cases. The quadratic form of this random vector follows a χ² distribution. The roles of the parameters in the model are investigated using contour plots of bivariate densities.     

Probabilities and moments of generalized discrete distributfons using finite difference operators

The probabilities and factorial moments of the univar iate and multivariate generalized (or compound) discrete di st r-Lbut Lons with probability generating functions H(t)=F(G(t)) and H(t₁,…,t_k)=F(G(t₁,…,t_k))or H(t₁,…,t_k) = F(G₁(t₁),…, G_k( t_k)) are derived using finite difference operators.     

A Multivariate Intervention Model for the Dutch Mint Circulation: Estimation and Monte Carlo Simulation

For each of the five Dutch coinage denominations, a transfer-function model is estimated. The output variables are monthly observations of coins in circulation. Two input variables represent transaction flows; all other inputs are step functions, representing the occurrence of interventions. Using the method of cross-correlating the residuals of the individual equations, a multivariate transfer-function model is constructed and estimated. Next Monte Carlo simulation is applied to derive expectations and variances of the yearly addition to the stock of coins until 1996. Our results enlighten on some aspects of a problem situation faced by the Dutch State Mint.     

New HEAVY Models for Fat-Tailed Realized Covariances and Returns

ABSTRACT

We develop a new score-driven model for the joint dynamics of fat-tailed realized covariance matrix observations and daily returns. The score dynamics for the unobserved true covariance matrix are robust to outliers and incidental large observations in both types of data by assuming a matrix-F distribution for the realized covariance measures and a multivariate Student's t distribution for the daily returns. The filter for the unknown covariance matrix has a computationally efficient matrix formulation, which proves beneficial for estimation and simulation purposes. We formulate parameter restrictions for stationarity and positive definiteness. Our simulation study shows that the new model is able to deal with high-dimensional settings (50 or more) and captures unobserved volatility dynamics even if the model is misspecified. We provide an empirical application to daily equity returns and realized covariance matrices up to 30 dimensions. The model statistically and economically outperforms competing multivariate volatility models out-of-sample. Supplementary materials for this article are available online.     

Asymptotic Distribution of the Estimated BDS Statistic from The Residuals of Location-Scale Type Processes

We introduce local asymptotic normality (LAN) for the log-likelihood ratio and the contiguity property of a class of location-scale processes in this paper. Using the contiguity property, we study the asymptotic distribution of the estimated Brock-Dechert-Scheinkman (BDS) statistic from the standardized residuals of fitting location-scale type processes.     

A Bayesian nonparametric estimator of a multivariate survival function

Using reinforced processes related to beta-Stacy process and generalized Pólya urn scheme jointly with a structure assumption about dependence, a Bayesian nonparametric prior and a predictive estimator for a multivariate survival function are provided. This estimator can be computed through an easy implementation of a Gibbs sampler algorithm. Moreover consistency of the estimator is studied.     

Optimal smoothing parameters for multivariate fized and adaptive kernel methods

Except in special cases optimum smoothing parameters of kernel methods are difficult to obtain for small samples, and large sample results are often used. Simulation is used to obtain finite sample optimum smoothing parameters and mean integrated square errors for the bivariate normal density. For this example, comparison is made of finite and asymptotic results, and of fixed and adaptive kernel methods. Further comparisons are made of fixed and adaptive methods by considering four other different types of density. Finally, some examples are given.     

A Procedure for Identification of Principal Variables by Least Generalized Dependence

Principal components are often used for reducing dimensions in multivariate data, but they frequently fail to provide useful results and their interpretation is rather difficult. In this article, the use of entropy optimization principles for dimensional reduction in multivariate data is proposed. Under the assumptions of multivariate normality, a four-step procedure is developed for selecting principal variables and hence discarding redundant variables. For comparative performance of the information theoretic procedure, we use simulated data with known dimensionality. Principal variables of cluster bean (Guar) are identified by applying this procedure to a real data set generated in a plant breeding experiment.