Earthquakes occurrences estimation through a parametric semi-Markov approach

The estimation of earthquakes’ occurrences prediction in seismic areas is a challenging problem in seismology and earthquake engineering. Indeed, the prevention and the quantification of possible damage provoked by destructive earthquakes are directly linked to this kind of prevision. In our paper, we adopt a parametric semi-Markov approach. This model assumes that a sequence of earthquakes is seen as a Markov process and besides it permits to take into consideration the more realistic assumption of events’ dependence in space and time. The elapsed time between two consecutive events is modeled as a general Weibull distribution. We determine then the transition probabilities and the so-called crossing states probabilities. We conclude then with a Monte Carlo simulation and the model is validated through a large database containing real data.     

Optimal financing and dividend policy with Markovian switching regimes

This work investigates an optimal financing and dividend problem for an insurer whose surplus process is modulated by an observable continuous-time and finite-state Markov chain. We assume that the insurer should never go bankrupt by issuing new equity. The goal of the insurer is to maximize the expected present value of the dividends payout minus the discounted cost of equity issuance. We obtain the optimal policies and explicit expressions for the value functions when the risk reserve process is modeled by both upward jump model and its diffusion approximation. Numerical illustrations of the sensitivities of the model parameters are provided.     

Efficient heterogeneous sampling for stochastic simulation with an illustration in health care applications

In modeling disease transmission, contacts are assumed to have different infection rates. A proper simulation must model the heterogeneity in the transmission rates. In this article, we present a computationally efficient algorithm that can be applied to a population with heterogeneous transmission rates. We conducted a simulation study to show that the algorithm is more efficient than other algorithms for sampling the disease transmission in a subset of the heterogeneous population. We use a valid stochastic model of pandemic influenza to illustrate the algorithm and to estimate the overall infection attack rates of influenza A (H1N1) in a Canadian city.     

Hierarchical likelihood approach to non-Gaussian factor analysis

Factor models, structural equation models (SEMs) and random-effect models share the common feature that they assume latent or unobserved random variables. Factor models and SEMs allow well developed procedures for a rich class of covariance models with many parameters, while random-effect models allow well developed procedures for non-normal models including heavy-tailed distributions for responses and random effects. In this paper, we show how these two developments can be combined to result in an extremely rich class of models, which can be beneficial to both areas. A new fitting procedures for binary factor models and a robust estimation approach for continuous factor models are proposed.     

Reversible jump Markov chain Monte Carlo algorithms for Bayesian variable selection in logistic mixed models

In this article, to reduce computational load in performing Bayesian variable selection, we used a variant of reversible jump Markov chain Monte Carlo methods, and the Holmes and Held (HH) algorithm, to sample model index variables in logistic mixed models involving a large number of explanatory variables. Furthermore, we proposed a simple proposal distribution for model index variables, and used a simulation study and real example to compare the performance of the HH algorithm with our proposed and existing proposal distributions. The results show that the HH algorithm with our proposed proposal distribution is a computationally efficient and reliable selection method.     

Kumaraswamy Marshall-Olkin Exponential distribution

In this paper, a new generalization of the Kumaraswamy distribution namely, the Kumaraswamy Marshall-Olkin Exponential distribution (KwMOE) is introduced and studied. Various properties are explored. The structural analysis includes various aspects such as limiting behaviour, shape properties, moments, quantiles, mean deviation, Renyi entropy, order statistics and stochastic ordering. Some useful characterizations of the family are also obtained. The method of maximum likelihood is used to estimate the model parameters. Monte Carlo simulation study is being conducted. An application to a real data set is presented for illustrative purposes.     

Anticipated Backward Stochastic Differential Equation with Reflection

In this article, we deal with anticipated backward stochastic differential equation with reflecting boundary. The existence and uniqueness of solution is obtained for equation with Lipschitz and non-Lipschitz generator.     

基于结构方程模型的农民工可雇佣性的决定因素研究——基于中部六省的调查

农民工是中国人力资源的重要组成部分。从智能、技能、体能和心能四个维度构建农民工可雇佣性模型,根据调查问卷收集数据,运用结构方程验证了农民工可雇佣性模型的合理性,研究了农民工可雇佣性及其各维度间的关系。结果表明,各维度对农民工可雇佣性影响程度依次为:技能、心能、智能和体能。农民工因年龄、文化程度不同,在可雇佣能力、技能、心能和智能方面存在显著性差异,但是在体能上没有显著差异。     

Accelerating inference for diffusions observed with measurement error and large sample sizes using approximate Bayesian computation

In recent years, dynamical modelling has been provided with a range of breakthrough methods to perform exact Bayesian inference. However, it is often computationally unfeasible to apply exact statistical methodologies in the context of large data sets and complex models. This paper considers a nonlinear stochastic differential equation model observed with correlated measurement errors and an application to protein folding modelling. An approximate Bayesian computation (ABC)-MCMC algorithm is suggested to allow inference for model parameters within reasonable time constraints. The ABC algorithm uses simulations of ‘subsamples’ from the assumed data-generating model as well as a so-called ‘early-rejection’ strategy to speed up computations in the ABC-MCMC sampler. Using a considerate amount of subsamples does not seem to degrade the quality of the inferential results for the considered applications. A simulation study is conducted to compare our strategy with exact Bayesian inference, the latter resulting two orders of magnitude slower than ABC-MCMC for the considered set-up. Finally, the ABC algorithm is applied to a large size protein data. The suggested methodology is fairly general and not limited to the exemplified model and data.     

Estimation of error probabilities in stochastic dominance

The process of using data to infer the existence of stochastic dominance is subject to sampling error. Kroll and Levy (1980), among others, have presented simulation results for several normal and lognormal distributions which show high error probabilities for a wide range of parameter values. This paper continues this line of research and uses simulation to estimate error probabilities. Distributions considered are a pair of normals and a pair of lognormals. Analysis of these distributions is made computationally feasible through theoretical results which reduce the number of parameters of the pair of distributions from four to two.