Authors of this paper were educated and have spent many years in Lomonosov Moscow State University (MSU), in the Interfaculty Laboratory of Statistical Methods (ILSM). Prof. Yu. K. Belyaev was student of Kolmogorov, and for several years his deputy in the ILSM, created and led by the great A. N. Kolmogorov. Now Prof. Yuri Belyaev is Emeritus Professor at Umea University (Sweden). Asaf Hajiyev was a PhD student in Kolmogorov's ILSM. In this paper some unusual incidents about the legendary Kolmogorov, and people around him are presented. Some stories were taken from the books and the internet (see references) dedicated to the memory of Kolmogorov; some were told by his students, and some happened with authors of this paper.

Men are cruel but Man is kind.

– A. N. Kolmogorov

Motivated by the classical De Bruijn's identity for the additive Gaussian noise channel, in this paper we consider a generalized setting where the channel is modelled via stochastic differential equations driven by fractional Brownian motion with Hurst parameter H ∈ (0, 1). We derive generalized De Bruijn's identity for Shannon entropy and Kullback–Leibler divergence by means of Itô's formula, and present two applications. In the first application we demonstrate its equivalence with Stein's identity for Gaussian distributions, while in the second application, we show that for H ∈ (0, 1/2], the entropy power is concave in time while for H ∈ (1/2, 1) it is convex in time when the initial distribution is Gaussian. Compared with the classical case of H = 1/2, the time parameter plays an interesting and significant role in the analysis of these quantities.

For a system that is subject to shocks, it is assumed that the distribution of the magnitudes of shocks changes after the first shock of size at least d1, and the system fails upon the occurrence of the first shock above a critical level d2 (> d1). In this paper, the distribution of the lifetime of such a system is studied when the times between successive shocks follow matrix-exponential distribution. In particular, it is shown that the system's lifetime has matrix-exponential distribution when the intershock times follow Erlang distribution. The model is extended to the case when the system fails upon the occurrence of l consecutive critical shocks.

We study a single server queue under a processor-sharing type of scheduling policy, where the weights for determining the sharing are given by functions of each job's remaining service (processing) amount, and obtain a fluid limit for the scaled measure-valued system descriptors.

The sample squared Sharpe ratio (SSR) is a critical statistic of the risk-return tradeoff. We show that sensitive upper-tail probabilities arise when the sample SSR is employed to test the mean-variance efficiency under different test statistics. Assuming the error's normality with a nonzero mean, we integrate the sample SSR and the arbitrage regression into a noncentral chi-square (χ2) test. We find that the distribution of the sample SSR based on the regression error is to the left of the F-distribution when assuming the returns' normality. Compared to two benchmarks that use the noncentral F-distribution and the central F-statistic, the χ2-statistic is more effective, competitive, significant, and locally robust when used to reject the upper-tailed mean-variance efficiency test using the usual parameters (sample size, portfolio size, and SSR).

In this paper, we study the equilibrium valuation for currency options in a setting of the two-country Lucas-type economy. Different from the continuous model in Bakshi and Chen [1], we propose a discontinuous model with jump processes. Empirical findings reveal that the jump components in each country's money supply can be decomposed into the simultaneous co-jump component and the country-specific jump component. Each of the jump components is modeled with a Poisson process whose jump intensity follows a mean reversion stochastic process. By solving a partial integro-differential equation (PIDE), we get a closed-form solution to the PIDE for a European call currency option. The numerical results show that the derived option pricing formula is efficient for practical use. Importantly, we find that the co-jump has a significant impact on option price and implied volatility.

Log-concavity [log-convexity] and their various properties play an increasingly important role in probability, statistics, operations research and other fields. In this paper, we first establish general preservation theorems of log-concavity and log-convexity under operator $\phi \longmapsto T(\phi , \theta )=\mathbb {E}[\phi (X_\theta )]$, θ ∈ Θ, where Θ is an interval of real numbers or an interval of integers, and the random variable $X_\theta$ has a distribution function belonging to the family $\{F_\theta , \theta \in \Theta \}$ possessing the semi-group property. The proofs are based on the theory of stochastic comparisons and weighted distributions. The main results are applied to some special operators, for example, operators occurring in reliability, Bernstein-type operators and Beta-type operators. Several known results in the literature are recovered.

This paper discusses stochastic comparisons for the residual and past lifetimes of coherent systems with dependent and identically distributed (d.i.d.) components under random monitoring in terms of the hazard rate, the reversed hazard rate, and the likelihood ratio orders. Some stochastic comparisons results are also established on the residual lifetimes of coherent systems under random observation times when all of the components are alive at that time. Sufficient conditions are established in terms of the aging properties of the components and the distortion functions induced from the system structure and dependence among components lifetimes. Numerical examples are provided to illustrate the theoretical results as well.

The relative aging is an important notion which is useful to measure how a system ages relative to another one. Among the existing stochastic orders, there are two important orders describing the relative aging of two systems, namely, aging faster orders in the cumulative hazard and the cumulative reversed hazard rate functions. In this paper, we give some sufficient conditions under which one coherent system ages faster than another one with respect to the aforementioned stochastic orders. Further, we show that the proposed sufficient conditions are satisfied for k-out-of-n systems. Moreover, some numerical examples are given to illustrate the applications of proposed results.

The aim of this paper is to introduce a new stochastic order based on the residual lifetimes of two nonnegative dependent random variables and the stochastic precedence order. We develop some characterizations and preservation properties of this stochastic order. In addition, we study some of its reliability properties and its relation with other existing stochastic orders. One of the possible applications in reliability theory has also been discussed.

In light of recent empirical research on jump activity, this article study the calibration of a new class of stochastic volatility models that include both jumps in return and volatility. Specifically, we consider correlated jump sizes and both contemporaneous and independent arrival of jumps in return and volatility. Based on the specifications of this model, we derive a closed-form relationship between the VIX index and latent volatility. Also, we propose a closed-form logarithmic likelihood formula by using the link to the VIX index. By estimating alternative models, we find that the general counting processes setting lead to better capturing of return jump behaviors. That is, the part where the return and volatility jump simultaneously and the part that jump independently can both be captured. In addition, the size of the jumps in volatility is, on average, positive for both contemporaneous and independent arrivals. However, contemporaneous jumps in the return are negative, but independent return jumps are positive. The sub-period analysis further supports above insight, and we find that the jumps in return and volatility increased significantly during the two recent economic crises.

We show that the kth order statistic from a heterogeneous sample of n ≥ k exponential random variables is larger than that from a homogeneous exponential sample in the sense of star ordering, as conjectured by Xu and Balakrishnan [14]. As a consequence, we establish hazard rate ordering for order statistics between heterogeneous and homogeneous exponential samples, resolving an open problem of Pǎltǎnea [11]. Extensions to general spacings are also presented.

This paper probabilistically explores a class of stationary count time series models built by superpositioning (or otherwise combining) independent copies of a binary stationary sequence of zeroes and ones. Superpositioning methods have proven useful in devising stationary count time series having prespecified marginal distributions. Here, basic properties of this model class are established and the idea is further developed. Specifically, stationary series with binomial, Poisson, negative binomial, discrete uniform, and multinomial marginal distributions are constructed; other marginal distributions are possible. Our primary goal is to derive the autocovariance function of the resulting series.

We introduce a Markov queueing system with Poisson arrivals, exponential services, and jockeying between two parallel and equivalent servers. An arriving customer admits to the shortest line. Every transition, of only the last customer in line, from the longer line to the shorter line may accompanied by a certain fixed cost. Thus, a transition from the longer queue to the shorter queue occurs whenever the difference between the lines reaches a certain discrete threshold (d = 2, 3, …). In this study, we focus on the stochastic analysis of the number of transitions of an arbitrary customer.

Two main methods are used to solve continuous-time quasi birth-and-death processes: matrix geometric (MG) and probability generating functions (PGFs). MG requires a numerical solution (via successive substitutions) of a matrix quadratic equation A0 + RA1 + R2A2 = 0. PGFs involve a row vector $\vec{G}(z)$ of unknown generating functions satisfying $H(z)\vec{G}{(z)^\textrm{T}} = \vec{b}{(z)^\textrm{T}},$ where the row vector $\vec{b}(z)$ contains unknown “boundary” probabilities calculated as functions of roots of the matrix H(z). We show that: (a) H(z) and $\vec{b}(z)$ can be explicitly expressed in terms of the triple A0, A1, and A2; (b) when each matrix of the triple is lower (or upper) triangular, then (i) R can be explicitly expressed in terms of roots of $\det [H(z)]$; and (ii) the stability condition is readily extracted.

This article considers the interval availability and instant availability of the k-system. A certain relationship between the two types of availability is established. Some lower and upper bounds to interval availability are derived. It also provides a couple of conditions under which the availability of two systems can be compared. Several examples are given to show the complexity of comparisons of availability.

Dielectric breakdown in a thin oxide is presented in terms of an interacting particle system on a two-dimensional lattice. All edges in the system are initially assumed to be closed. An edge between two adjacent vertices will open according to an exponentially distributed random variable. Breakdown occurs at the time an open path connects the top layer of the lattice to the bottom layer. Using the extreme value theory, we show that the time until breakdown is asymptotically Weibull distributed.

This paper deals with stochastic comparisons of the largest order statistics arising from two sets of independent and heterogeneous gamma samples. It is shown that the weak supermajorization order between the vectors of scale parameters together with the weak submajorization order between the vectors of shape parameters imply the reversed hazard rate ordering between the corresponding maximum order statistics. We also establish sufficient conditions of the usual stochastic ordering in terms of the p-larger order between the vectors of scale parameters and the weak submajorization order between the vectors of shape parameters. Numerical examples and applications in auction theory and reliability engineering are provided to illustrate these results.

In this paper, the problem of goodness of fit test for the inverse Rayleigh distribution based on progressively Type-II censored samples is studied. We develop two test statistics via entropy and propose one new non-entropy test statistic via a pivotal method. We also study the properties of these test statistics. Critical values are obtained by simulations. Then, we do power analysis of these test statistics against various alternatives under different censoring schemes. We conclude that the tests we proposed perform well against various alternatives, especially for non-monotone hazard alternatives. Finally, one real data set is analyzed.

This paper presents the distribution of the number of customers served during a busy period for special cases of the Geo/G/1 queue when initiated with m customers. We analyze the system under the assumptions of a late arrival system with delayed access and early arrival system policies. It is not easy to invert the functional equation for the number of customers served during a busy period except for the simple case Geo/Geo/1 queue, as stated by several researchers. Using the Lagrange inversion theorem, we give an elegant solution to this equation. We find the distribution of the number of customers served during a busy period for various service-time distributions such as geometric, deterministic, binomial, negative binomial, uniform, Delaporte, discrete phase-type and interrupted Bernoulli process. We compute the mean and variance of these distributions and also give numerical results. Due to the clarity of the expressions, the computations are very fast and robust. We also show that in the limiting case, the results tend to the analogous continuous-time counterparts.