Nonparametric product partition models for multiple change-points analysis

ABSTRACT

We propose an extension of parametric product partition models. We name our proposal nonparametric product partition models because we associate a random measure instead of a parametric kernel to each set within a random partition. Our methodology does not impose any specific form on the marginal distribution of the observations, allowing us to detect shifts of behaviour even when dealing with heavy-tailed or skewed distributions. We propose a suitable loss function and find the partition of the data having minimum expected loss. We then apply our nonparametric procedure to multiple change-point analysis and compare it with PPMs and with other methodologies that have recently appeared in the literature. Also, in the context of missing data, we exploit the product partition structure in order to estimate the distribution function of each missing value, allowing us to detect change points using the loss function mentioned above. Finally, we present applications to financial as well as genetic data.     

A Whirlwind of Change Data: How Libraries and Vendors Translate Data Into Discovery and Access

Abstract

Libraries and vendors share change data to ensure successful discovery and access for users. Change data have increased dramatically and affected data communication, causing barriers to access. Serials data used to be simple: There were new, ceased, and title changes. Data were efficiently exchanged. Now data are complex and difficult to communicate in the digital age. This presentation will discuss the challenges of change data and the ways libraries and vendors are teaming up to calm the whirlwind of change.     

A Bayesian random effects model for analyzing mixed negative binomial and continuous longitudinal responses

ABSTRACT

A general Bayesian random effects model for analyzing longitudinal mixed correlated continuous and negative binomial responses with and without missing data is presented. This Bayesian model, given some random effects, uses a normal distribution for the continuous response and a negative binomial distribution for the count response. A Markov Chain Monte Carlo sampling algorithm is described for estimating the posterior distribution of the parameters. This Bayesian model is illustrated by a simulation study. For sensitivity analysis to investigate the change of parameter estimates with respect to the perturbation from missing at random to not missing at random assumption, the use of posterior curvature is proposed. The model is applied to a medical data, obtained from an observational study on women, where the correlated responses are the negative binomial response of joint damage and continuous response of body mass index. The simultaneous effects of some covariates on both responses are also investigated.     

The Extremal Dependence Measure and Asymptotic Independence

Abstract

Extremal dependence analysis assesses the tendency of large values of components of a random vector to occur simultaneously. This kind of dependence information can be qualitatively different than what is given by correlation which averages over the total body of the joint distribution. Also, correlation may be completely inappropriate for heavy tailed data. We study the extremal dependence measure (EDM), a measure of the tendency of large values of components of a random vector to occur simultaneously and show consistency of an estimator of the EDM. We also show asymptotic normality of an idealized estimator in a restricted case of multivariate regular variation where scaling functions do not have to be estimated.     

The distribution of the product of two independent generalized trapezoidal random variables

ABSTRACT

In this article, we derive the probability density function (pdf) of the product of two independent generalized trapezoidal random variables having different supports, in closed form, by considering all possible cases. We also show that the results for the product of two triangular and uniform random variables follow as special cases of our main result. As an illustration, we obtain pdf of product for a suitably constrained set of parameters and plot some graphs using MATLAB, which express variation in pdf with change in different parameters of the generalized trapezoidal distribution.     

Asymptotic normality of a conditional hazard function estimate in the single index for quasi-associated data

Abstract

The main goal of this paper is to study the estimation of the conditional hazard function of a scalar response variable Y given a hilbertian random variable X in functional single-index model. We construct an estimator of this nonparametric function and we study its asymptotic properties, under quasi-associated structure. Precisely, we establish the asymptotic normality of the constructed estimator. We carried out simulation experiments to examine the behavior of this asymptotic property over finite sample data.     

The LAD estimation of the change-point linear model with randomly censored data

Abstract

In this paper, a change-point linear model with randomly censored data is investigated. We propose the least absolute deviation estimation procedure for regression and change-point parameters simultaneously. The asymptotic properties of the change-point and regression parameter estimators are obtained. We show that the resulting regression parameter estimator is asymptotically normal, and the change-point estimator converges weakly to the minimizer of a given random process. The extensive simulation studies and the analysis of an acute myocardial infarction data set are conducted to illustrate the finite sample performance of the proposed method.     

Some Estimators of Finite Population Variance of Stratified Sample Mean

We consider a continuous-time branching random walk on Z ^d , where the particles are born and die at a single lattice point (the source of branching). The underlying random walk is assumed to be symmetric. Moreover, corresponding transition rates of the random walk have heavy tails. As a result, the variance of the jumps is infinite, and a random walk may be transient even on low-dimensional lattices (d = 1, 2). Conditions of transience for a random walk on Z ^d and limit theorems for the numbers of particles both at an arbitrary point of the lattice and on the entire lattice are obtained.     

Regression model for surrogate data in high dimensional statistics

Abstract

This paper deals with the problem of estimating the regression of a surrogated scalar response variable given a functional random one. We construct an estimator of the regression operator by using, in addition to the available (true) response data, a surrogate data. We then establish some asymptotic properties of the constructed estimator in terms of the almost-complete and the quadratic mean convergences. Notice that the obtained results generalize a part of the results obtained in the finite dimensional framework. Finally, an illustration on the applicability of our results on both simulated data and a real dataset was realized. We have thus shown the superiority of our estimator on classical estimators when we are lacking complete data.     

Testing change in the variance with epidemic alternatives

ABSTRACT

In this paper we consider the dyadic increments statistics (of type DI) based on independent not identically distributed or α-mixing random variables. We obtain their limit distributions under the null hypothesis and we present application for testing epidemic change in the variance in each case. Finally, numerical simulations are done to illustrate these results.     

Empirical likelihood inference for a semiparametric hazards regression model

ABSTRACT

We investigated the empirical likelihood inference approach under a general class of semiparametric hazards regression models with survival data subject to right-censoring. An empirical likelihood ratio for the full 2p regression parameters involved in the model is obtained. We showed that it converged weakly to a random variable which could be written as a weighted sum of 2p independent chi-squared variables with one degree of freedom. Using this, we could construct a confidence region for parameters. We also suggested an adjusted version for the preceding statistic, whose limit followed a standard chi-squared distribution with 2p degrees of freedom.     

Guaranteed Local Maximum Likelihood Detection of a Change Point in Nonparametric Logistic Regression

We consider nonparametric logistic regression and propose a generalized likelihood test for detecting a threshold effect that indicates a relationship between some risk factor and a defined outcome above the threshold but none below it. One important field of application is occupational medicine and in particular, epidemiological studies. In epidemiological studies, segmented fully parametric logistic regression models are often threshold models, where it is assumed that the exposure has no influence on a response up to a possible unknown threshold, and has an effect beyond that threshold. Finding efficient methods for detection and estimation of a threshold is a very important task in these studies. This article proposes such methods in a context of nonparametric logistic regression. We use a local version of unknown likelihood functions and show that under rather common assumptions the asymptotic power of our test is one. We present a guaranteed non asymptotic upper bound for the significance level of the proposed test. If applying the test yields the acceptance of the conclusion that there was a change point (and hence a threshold limit value), we suggest using the local maximum likelihood estimator of the change point and consider the asymptotic properties of this estimator.     

Limit Distributions of Random Sums of Z+-Valued Random Variables

In this article, we formulate a transfer theorem in terms of probability generating functions and discuss two approaches to limit distributions of random sums of Z ₊-valued random variables. We then develop Z ₊-valued N-ID and ?-ID laws.     

A Resampling Approach for Under-estimating a Finite Population Total from a Censored Sample

Abstract

Suppose a finite population of N objects each of which has an unknown value μ _i  ≥ 0, i = 1, … , N of a nonnegative characteristic of interest. A random sample has been drawn, but only for a selected subset of the sample the μ-values have been observed. The subset selection procedure has been somewhat obscure, and thus the subsample is censorized rather than random. Despite that, a reliable lower bound for the population total (the sum of all μ _i ) is required which uses the statistical information contained in the data. We propose a resampling procedure to construct an under-estimate of the population total. We also consider the case when the objects of the population have unequal sampling probabilities, in particular when the population is divided into a few number of strata with constant probabilities within each stratum. A real data example illustrates the method.     

Proportional hazards model for discrete data: Some new developments

ABSTRACT

Many times, a product lifetime can be described through a non negative integer valued random variable. In this article, we propose a proportional hazards model for discrete data analogous to the version for continuous data. Some ageing properties of the model are discussed. Stochastic comparison of pair of random variables that follow the model are also made. A new test based on U-statistics is developed for testing that the proportionality parameter in the proposed model is 1. The asymptotic properties of the proposed test are studied. We present some numerical results to asses the performance of the test procedure.     

The asymptotic distributions of maximum likelihood ratio test and maximally selected χ2-test in binomial observations

We considered binomial distributed random variables whose parameters are unknown and some of those parameters need to be estimated. We studied the maximum likelihood ratio test and the maximally selected χ²-test to detect if there is a change in the distributions among the random variables. Their limit distributions under the null hypothesis and their asymptotic distributions under the alternative hypothesis were obtained when the number of the observations is fixed. We discussed the properties of the limit distribution and found an efficient way to calculate the probability of multivariate normal random variables. Finally, those results for both tests have been applied to examples of Lindisfarne's data, the Talipes Data. Our conclusions are consistent with other researchers' findings.     

Nonparametric Estimation of Distribution Functions of Nonstandard Mixtures

ABSTRACT

Nonstandard mixtures are those that result from a mixture of a discrete and a continuous random variable. They arise in practice, for example, in medical studies of exposure. Here, a random variable that models exposure might have a discrete mass point at no exposure, but otherwise may be continuous. In this article we explore estimating the distribution function associated with such a random variable from a nonparametric viewpoint. We assume that the locations of the discrete mass points are known so that we will be able to apply a classical nonparametric smoothing approach to the problem. The proposed estimator is a mixture of an empirical distribution function and a kernel estimate of a distribution function. A simple theoretical argument reveals that existing bandwidth selection algorithms can be applied to the smooth component of this estimator as well. The proposed approach is applied to two example sets of data.     

Generation of pseudo random numbers and estimation under Cox models with time-dependent covariates

ABSTRACT

We study the method for generating pseudo random numbers under various special cases of the Cox model with time-dependent covariates when the baseline hazard function may not be constant and the random variable may equal infinity with a positive probability. During our simulation studies in computing the partial likelihood estimates, in between 3% and 20% of the time with a moderate sample size, it happens that the partial likelihood estimate of the regression coefficient is ∞ for the data from the Cox model. We propose a semi-parametric estimator as a modification for such a case. We present simulation results on the asymptotic properties of the semi-parametric estimator.     

Optimization problems for consecutive-2-out-of-n:G system

Abstract

This paper considers the optimization problems for a consecutive-2-out-of-n:G system where n is considered to be fixed or random. When the number of components is constant, the optimal number of components and the optimal replacement time are discussed by minimizing the expected cost rates. Furthermore, we focus on the above discussions again when n is a random variable. We give an approximate value of MTTF and propose the preventive replacement policy, respectively.     

On some study of skew-t distributions

Abstract

In this note, through ratio of independent random variables, new families of univariate and bivariate skew-t distributions are introduced. Probability density function for each skew-t distribution will be given. We also derive explicit forms of moments of the univariate skew-t distribution and recurrence relations for its cumulative distribution function. Finally we illustrate the flexibility of this class of distributions with applications to a simulated data and the volcanos heights data.