Distribution-free confidence intervals for quantiles and tolerance intervals in terms of k-records

In this paper, we consider the problem of determining non-parametric confidence intervals for quantiles when available data are in the form of k-records. Distribution-free confidence intervals as well as lower and upper confidence limits are derived for fixed quantiles of an arbitrary unknown distribution based on k-records of an independent and identically distributed sequence from that distribution. The construction of tolerance intervals and limits based on k-records is also discussed. An exact expression for the confidence coefficient of these intervals are derived. Some tables are also provided to assist in choosing the appropriate k-records for the construction of these confidence intervals and tolerance intervals. Some simulation results are presented to point out some of the features and properties of these intervals. Finally, the data, representing the records of the amount of annual rainfall in inches recorded at Los Angeles Civic Center, are used to illustrate all the results developed in this paper and also to demonstrate the improvements that they provide on those based on either the usual records or the current records.     

On characterization results based on the number of observations near the k-records

In this paper, using the completeness properties of the sequence of functions {h_n(x)=(?log x)ⁿ, 0<x<1, n≥1}, some characterization results are established. The results are based on the number of observations near the k-records. It is shown that the equality of the moment of the appropriate subsequence of the number of observations near to upper and lower k-records is a characteristic property of symmetric distributions. Since ordinary record values are contained in the k-records, the results hold for usual records as a particular case.     

Simultaneous Closeness of k-Records

In this article, we study the simultaneous Pitman closeness of upper (and lower) k-records to a common parameter of the parent distribution and obtain general expressions for the corresponding probabilities. Since the usual record values are contained in the k-records, the corresponding results for the usual records can be deduced as special cases. The results are then applied to location-scale families with an emphasis on population quantiles. Exact expressions are derived for some common distributions such as exponential and uniform. In each case, the simultaneous closest k-record to a population quantile is determined.     

On Moments of k-Records from Discrete Distributions

In this article, we establish recurrence relations satisfied by first and second moments of k-records arising from discrete distributions. Then we use these relations to obtain means, variances, and correlation coefficients of geometric k-records. We consider all the three known types of k-records: strong, ordinary, and weak.     

Some results on conditional expectations of lower record values

In this paper, we propose a measure for obtaining the expectation of time between two lower k-records under the condition that the greater one is given. Several properties of the proposed measure are derived. Some characterization results and stochastic comparisons based on the new measure are also provided.     

Comparison of Records and Maxima

We establish best upper bounds on the expected differences of records and sample maxima, and kth records and kth maxima based on sequences of independent random variables with identical continuous distribution and finite variance. The bounds are expressed in terms of the standard deviation units of the parent distribution. We also provide conditions for attaining the bounds.     

A latent variable model for analyzing mixed longitudinal (k,l)-inflated count and ordinal responses

A random effects model for analyzing mixed longitudinal count and ordinal data is presented where the count response is inflated in two points (k and l) and an (k,l)-Inflated Power series distribution is used as its distribution. A full likelihood-based approach is used to obtain maximum likelihood estimates of parameters of the model. For data with non-ignorable missing values models with probit model for missing mechanism are used.The dependence between longitudinal sequences of responses and inflation parameters are investigated using a random effects approach. Also, to investigate the correlation between mixed ordinal and count responses of each individuals at each time, a shared random effect is used. In order to assess the performance of the model, a simulation study is performed for a case that the count response has (k,l)-Inflated Binomial distribution. Performance comparisons of count-ordinal random effect model, Zero-Inflated ordinal random effects model and (k,l)-Inflated ordinal random effects model are also given. The model is applied to a real social data set from the first two waves of the national longitudinal study of adolescent to adult health (Add Health study). In this data set, the joint responses are the number of days in a month that each individual smoked as the count response and the general health condition of each individual as the ordinal response. For the count response there is incidence of excess values of 0 and 30.     

A new method for fast computing unbiased estimators of cumulants

We propose new algorithms for generating k-statistics, multivariate k-statistics, polykays and multivariate polykays. The resulting computational times are very fast compared with procedures existing in the literature. Such speeding up is obtained by means of a symbolic method arising from the classical umbral calculus. The classical umbral calculus is a light syntax that involves only elementary rules to managing sequences of numbers or polynomials. The cornerstone of the procedures here introduced is the connection between cumulants of a random variable and a suitable compound Poisson random variable. Such a connection holds also for multivariate random variables.     

Prediction intervals for generalized-order statistics with random sample size

The problem of predicting future generalized-order statistics, by assuming the future sample size is a random variable, is discussed. A general expression for the coverage probability of the prediction intervals is derived. Since k-records and progressively type-II censored-order statistics are contained in the model of generalized-order statistics, the corresponding results for them can be deduced as special cases. When the future sample size has degenerate, binomial, Poisson and geometric distributions, numerical computations are given. The procedure for finding an optimal prediction interval is presented for each case. Finally, we apply our results to a real data set in life testing given in Lee and Wang [Statistical methods for survival data analysis. Hoboken, NJ: John Wiley and Sons; 2003, p. 58, Table 3.4] for illustrative the proposed procedure in this paper.     

Mean and cell-mean imputation resulting in the use of a semicontinuous distribution and a mixture of semicontinuous distributions

Zero-inflated models are commonly used for modeling count and continuous data with extra zeros. Inflations at one point or two points apart from zero for modeling continuous data have been discussed less than that of zero inflation. In this article, inflation at an arbitrary point α as a semicontinuous distribution is presented and the mean imputation for a continuous response is discussed as a cause of having semicontinuous data. Also, inflation at two points and generally at k arbitrary points and their relation to cell-mean imputation in the mixture of continuous distributions are studied. To analyze the imputed data, a mixture of semicontinuous distributions is used. The effects of covariates on the dependent variable in a mixture of k semicontinuous distributions with inflation at k points are also investigated. In order to find the parameter estimates, the method of expectation–maximization (EM) algorithm is used. In a real data of Iranian Households Income and Expenditure Survey (IHIES), it is shown how to obtain a proper estimate of the population variance when continuous missing at random responses are mean imputed.     

On k-match problems

Several waiting time random variables for a duplication within a memory window of size k in a sequence of {1,2,…,m}-valued random variables are investigated. The exact distributions of the waiting time random variables are derived by the method of conditional probability generating functions. In particular, the exact distribution of the waiting time for the first k-match is obtained when the underlying sequence is generated by higher order Markov dependent trials. Examples for numerical calculations are also given.     

Maximizing survival time in a random walk on an interval

A gambler buys N tokens that enable him to play N rounds of the following game. A symmetric random walk on a discrete interval { ? r, …, r} starts from the point 0. The gambler knows only the number of steps made so far, but is unaware of the current position of the walk. Once the walk hits one of the barriers ? r or r for the first time in the current round, the round ends with no payoff. The gambler can start a new round by inserting a new token, if there are any tokens left. The gambler can end the game at any time getting the payoff equal to the number of steps made in the current round. We find the optimal stopping strategy for this game and calculate the expected payoff once the optimal strategy is applied.     

Bayesian prediction of order statistics with fixed and random sample sizes based on k-record values from Pareto distribution

In this paper, the two-parameter Pareto distribution is considered and the problem of prediction of order statistics from a future sample and that of its geometric mean are discussed. The Bayesian approach is applied to construct predictors based on observed k-record values for the cases when the future sample size is fixed and when it is random. Several Bayesian prediction intervals are derived. Finally, the results of a simulation study and a numerical example are presented for illustrating all the inferential procedures developed here.     

DETERMINATION OF DOMAINS OF ATTRACTION BASED ON A SEQUENCE OF MAXIMA

Suppose that the maximum of a random sample from a distribution F(x) may be obtained in each of k equally spaced observation periods. This paper proposes a test to determine the domain of attraction of F(x), and investigates the properties when the sample size is very large and perhaps unknown and k is fixed and small. The test statistic is a function of the spacings between the order statistics based on the sequence of maxima and is suggested by reference to one studied previously when inference was based on the largest k observations of a random sample. A Monte Carlo study shows that the proposed test is more powerful than its main competitor. The test is illustrated by two examples.     

Upper and lower bounds on the variances of linear combinations of the kth records

We describe a method of determining upper bounds on the variances of linear combinations of the kth records values from i.i.d. sequences, expressed in terms of variances of parent distributions. We also present conditions for which the bounds are sharp, and those for which the respective lower ones are equal to zero. A special attention is paid to the case of the kth record spacings, i.e. the differences of consecutive kth record values.     

A new test of discordancy in cylindrical data

ABSTRACT

Cylindrical data are bivariate data from the combination of circular and linear variables. However, up to now no work has been done on the detection of outlier in cylindrical data. We introduce a definition of outlier for cylindrical data and present a new test of discordancy to detect outlier in this type of data, based on the k-nearest neighbor’s distance. Cut-off points of the new test statistic based on the Johnson-Wehrly distribution are calculated and its performance is examined using simulation. A practical example is presented using wind speed and wind direction data obtained from the Malaysian Meteorological Department.     

Fisher information in generalized order statistics

A representation of the Fisher information in generalized order statistics in terms of the hazard rate of the underlying distribution function is derived under mild regularity conditions. This expression supplements results for complete, Type-II censored, and progressively Type-II censored data. As a byproduct, we find a hazard rate based representation for samples of k-records which apparently has not been known so far. Moreover, sufficient conditions for the validity of this representation in location and scale family settings are given. The result is illustrated by considering generalized order statistics based on logistic, Laplace, and extreme value distributions.     

Asymptotic tail behavior of a random sum with conditionally dependent subexponential summands

In this paper, we obtain some results for the asymptotic behavior of the tail probability of a random sum S_τ = ∑^τ_{k = 1}X_k, where the summands X_k, k = 1, 2, …, are conditionally dependent random variables with a common subexponential distribution F, and the random number τ is a non negative integer-valued random variable, independent of {X_k: k ? 1}.     

Joint (k,l)-hurdle random effects models for mixed longitudinal power series and normal outcomes

A random effects model for analyzing mixed longitudinal normal and count outcomes with and without the possibility of non ignorable missing outcomes is presented. The count response is inflated in two points (k and l) and the (k, l)-Hurdle power series is used as its distribution. The new distribution contains, as special submodels, several important distributions which are discussed, such as (k, l)-Hurdle Poisson and (k, l)-Hurdle negative binomial and (k, l)-Hurdle binomial distributions among others. Random effects are used to take into account the correlation between longitudinal outcomes and inflation parameters. A full likelihood-based approach is used to yield maximum likelihood estimates of the model parameters. A simulation study is performed in which for count outcome (k, l)-Hurdle Poisson, (k, l)-Hurdle negative binomial and (k, l)-Hurdle binomial distributions are considered. To illustrate the application of such modelling the longitudinal data of body mass index and the number of joint damage are analyzed.     

A stochastic restricted k–d class estimator

In this paper, we introduce a stochastic restricted k–d class estimator for the vector of parameters in a linear model when additional linear restrictions on the parameter vector are assumed to hold. The stochastic restricted k–d class estimator is a generalization of the ordinary mixed estimator and the k–d class estimator. We show that our new biased estimator is superior in the mean squared error matrix sense to the k–d class estimator [S. Sakall?o?lu and S. Kaçiranlar, A new biased estimator based on ridge estimation, Statist. Papers 49 (2008), pp. 669–689] and the stochastic restricted Liu estimator [H. Yang and J.W. Xu, An alternative stochastic restricted Liu estimator in linear regression, Statist. Papers 50 (2009), pp. 639–647]. Finally, a numerical example is given to show the theoretical results.