Spatial regression models have garnered significant attention across several disciplines, including functional magnetic resonance imaging analysis, econometrics, home price analysis, and many other domains. The phenomenon of sparsity is often found in nature, when a limited number of factors contribute significantly to the overall variation. Spatial regression models frequently use sparsity to indicate less complex computational and more superficial covariance structures. The spatial error model is a significant spatial regression model that focuses on the geographical dependence present in the error terms rather than the response variable. This study proposes an effective approach for estimating the vector of regression coefficients in the spatial error model, taking into consideration of the prior knowledge that some coefficients are insignificant and there is multicollinearity among the regressors. It also introduces pretest and shrinkage ridge estimators for spatial error regression models, evaluating their performance compared to traditional maximum likelihood estimators. It also assesses their efficacy using real-world data and bootstrapping techniques for comparison purposes.

This article deals with a finite-source stock-dependent stochastic inventory system with heterogeneous servers and retrial facility. The system can store up to S number of items, and the lifetime of each item follows an exponential distribution. The primary customers arrive in a finite waiting hall from a finite source and get service from the multi-channel (c-heterogeneous servers). The customer’s arrival rate depends on the present stock level. Those customers who arrive while the waiting hall is full enter into the finite orbit. Also, customers in the waiting hall can become impatient and enter the orbit. The (s,Q) ordering policy is adopted for the replenishment of stocks. In the steady state, the joint probability distribution of a number of items in the inventory, busy servers, customers in the waiting hall, and orbit is computed. The comparative numerical analysis of heterogeneous and homogeneous servers on the expected total cost, the average number of customers in the waiting hall and orbit, the fraction of successful retrial rate, and the average impatient customer rate are also carried out.

The aim of this paper is to study a semi-functional partial linear regression model (SFPLR) for spatial data with responses missing at random. The estimators are constructed by the kernel method, and some asymptotic properties such as probability convergence rates of the nonparametric component and asymptotic distribution of the parametric and nonparametric components are established under certain conditions. Next, the performances and the superiority of these estimators are presented and examined using a study on simulated data and on real data by carrying out a comparison between our semi-functional partially linear model with MAR estimator (SFPLRM), the semi-functional partially linear model with the full-case estimator (SFPLRC) and the nonparametric functional model estimator with MAR (FNPM). The results show that the proposed estimators outperform existing estimators as the number of random missing data increases.

This paper aims to contribute to refining the e-values for testing precise hypotheses, especially when dealing with nuisance parameters, leveraging the effectiveness of asymptotic expansions of the posterior. The proposed approach offers the advantage of bypassing the need for elicitation of priors and reference functions for the nuisance parameters and the multidimensional integration step. For this purpose, starting from a Laplace approximation, a posterior distribution for the parameter of interest only is considered and then a suitable objective matching prior is introduced, ensuring that the posterior mode aligns with an equivariant frequentist estimator. Consequently, both Highest Probability Density credible sets and the e-value remain invariant. Some targeted and challenging examples are discussed.

Empirical equations representing, interpolating and smoothing groups of measurement results by regression methods are in widespread use in metrology and other fields of science and engineering. Standard equations for the propagation of measurement uncertainties to values computed from such equations are available but suffer from a lack of general acceptance and are only infrequently applied in practice. One reason for the slow uptake is the lack of clear methods to account for systematic error. In this paper, uncertainty propagation equations in terms of covariance matrices are generalized to allow for systematic errors in least-squares regression, and effects of the resulting uncertainties are investigated analytically. A stochastic ensemble model of systematic error is proposed for the computation of the non-diagonal elements of the weight matrix of the Generalized Least-Squares method (GLS) from the measurement uncertainty. A GLS projector formalism is described which separates the effects of measurement scatter on values calculated by the equation from those on the related “residuals”, i.e., the residual errors in the fitted data. The same projectors act similarly on the associated uncertainties and covariance matrices and permit the effects of systematic errors on the simulation covariance to be quantified. It is demonstrated that, in order to include systematic errors in the uncertainty estimates of GLS, covariance matrices may be substituted by novel, specifically defined dispersion matrices that are not specified yet by the GUM[1]. Systematic measurement errors may be estimated from various sources; a particularly easy way suggested here involves analyzing the structure in the regression residuals. It is demonstrated that GLS projectors exhibit inherent features of “error blindness” also with respect to systematic errors. [1] GUM: Guide to the Expression of Uncertainty in Measurement, http://www.bipm.org/en/publications/guides

We present a new penalized method for estimation in sparse logistic regression models with group structure. Group sparsity suggests us to consider the Group Lasso penalty. Being different from penalized log-likelihood estimation, our method can be viewed as a penalized weighted score function method. Under some mild conditions, we provide non-asymptotic oracle inequalities promoting group sparsity of predictors. A modified block coordinate descent algorithm based on a weighted score function is also employed. The net advantage of our algorithm over the existing Group Lasso-type procedures is that the tuning parameter can be pre-specified. The simulations show that this algorithm is considerably faster and more stable than competing methods. Finally, we illustrate our methodology with two real data sets.

Financial time series and other human-driven, non-natural processes are known to exhibit fat-tailed outcome distributions. That is, such processes demonstrate a greater tendency for extreme outcomes than the normal distribution or other natural distributional processes would predict. We examine the mathematical expectation, or simply "expectation," traditionally the probability-weighted outcome, regarded since the seventeenth century as the mathematical definition of "expectation." However, when considering the "expectation" of an individual confronted with a finite sequence of outcomes, particularly existential outcomes (e.g., a trader with a limited time to perform or lose his position in a trading operation), we find this individual "expects" the median terminal outcome over those finite trials, with the classical seventeenth-century definition being the asymptotic limit as trials increase to this. Since such finite-sequence "expectations" often differ in values from the classic one, so do the optimal allocations (e.g., growth-optimal). We examine these for fat-tailed distributions. The focus is on implementation, and the techniques described can be applied to all distributional forms.

The distribution laws of various natural and anthropogenic processes in the world around us are stochastic in nature. The development of mathematics and, in particular, of stochastic modelling allows us to study regularities in such processes. In practice, stochastic modelling finds a huge number of applications in various fields, including finance and economics. In this work, some particular applications of stochastic processes in finance are examined in the conditions of financial crisis. More specifically, Autoregressive Integrated Moving Average (ARIMA) models and Modified Ordinary Differential Equations (ODE) models, which have been previously developed by the authors to predict assets’ prices of four Bulgarian companies, are validated on a time period during the crisis. Estimated rates of return are calculated from the models for one period ahead. The errors are estimated and the models are compared. The predicted return values with each of the two approaches are used to derive optimal risk portfolios based on the Markowitz Model. The resulting portfolios are compared in terms of distribution (weights of the stocks), risk and rate of return.

Measurements of systems taken along a continuous functional dimension, such as time or space, are ubiquitous in many fields, from the physical and biological sciences to economics and engineering. Such measurements can be viewed as realisations of an underlying smooth process sampled over the continuum. However, traditional methods for independence testing and causal learning are not directly applicable to such data, as they do not take into account the dependence along the functional dimension. By using specifically designed kernels, we introduce statistical tests for bivariate, joint, and conditional independence for functional variables. Our method not only extends the applicability to functional data of the Hilbert-Schmidt independence criterion (HSIC) and its d-variate version (d-HSIC), but also allows us to introduce a test for conditional independence by defining a novel statistic for the conditional permutation test (CPT) based on the Hilbert-Schmidt conditional independence criterion (HSCIC), with optimised regularisation strength estimated through an evaluation rejection rate. Our empirical results of the size and power of these tests on synthetic functional data show good performance, and we then exemplify their application to several constraint- and regression-based causal structure learning problems, including both synthetic examples and real socio-economic data.

This paper presents a new discrete counterpart of the mixture exponential distribution by utilizing the survival discretization method. The moment-generating function and associated moment measures are discussed. The distribution’s hazard rate function can assume increasing or decreasing forms, making it adaptable for diverse fields requiring count data modelling. The paper explores four distinct parameter estimation methods and assesses their performance through Monte Carlo simulations. The applicability of this distribution extends to time series analysis, particularly within the framework of the first-order integer-valued autoregressive process. Additionally, the paper explores quality control applications, addressing serial dependence challenges in count data encountered in production and market management. The performance of two distinct control charts, the cumulative sum chart and the exponentially weighted moving average chart, is evaluated for their effectiveness in detecting shifts in the process means under various models. A bivariate Markov chain approach is used to estimate the average run lengths of these charts, offering valuable insights for implementation. Design recommendations for achieving robustness in-control chart applications are provided. The effectiveness of the proposed models and charts is illustrated using a real data, demonstrating their practical superiority.

Random Forest (RF) is a widely used data prediction and variable selection technique. However, the variable selection aspect of RF can become unreliable when there are more irrelevant variables than relevant ones. In response, we introduced the Bayesian Random Forest (BRF) method specifically designed for high-dimensional datasets with a sparse covariate structure. Our research demonstrates that BRF possesses the oracle property, which means it achieves strong selection consistency without compromising efficiency or bias.

This article deals with the queueing-inventory system, composed of c junior servers, a senior server, two finite waiting halls, and an infinite orbit. On occasion, junior servers encounter challenges during customer service. In these instances, they approach the senior server for guidance in resolving the issue. Suppose the senior server is engaged with another junior server. The approaching junior servers await their turn in a finite waiting area with a capacity of c for consultation. Concerning this, we study the performance of junior servers approaching the senior server in the retrial queueing-inventory model with the two finite waiting halls dedicated to the primary customers and the junior servers for consultation. We formulate a level-dependent QBD process and solve its steady-state probability vector using Neuts and Rao’s truncation method. The stability condition of the system is derived, and the R matrix is computed. Optimum total cost has been obtained, and sensitivity analyses, which include expected total cost, waiting time of customers in waiting hall and orbit, number of busy servers, and fraction of successful retrial rate of the model, are done numerically.

In this article we introduce an extension of the Akash distribution. We use the slash 1 methodology to make the kurtosis of the Akash distribution more flexible. We study the general 2 density of this new distribution, some properties, moments, coefficients of asymmetry and kurtosis. 3 Statistical inference is performed using the methods of moments and maximum likelihood via the 4 EM algorithm. A simulation study is carried out to observe the behavior of the maximum likelihood 5 estimator. An application to a real data set with high kurtosis is considered, where it is shown that 6 the new distribution fits better than other extensions of the Akash distribution.

In this work, we propose a generalized partially functional linear model, which not only models the relationship between multiple scalar and functional predictors and response, but also automatically estimates the link function. Specifically, we use the functional principal component analysis method to reduce the dimensionality of functional predictors, estimate the regression coefficients using the maximum likelihood estimation method, estimate the link function using the method of local linear regression, iteratively obtain the final estimator, and prove the asymptotic normality of the estimator. The asymptotic normality is illustrated through simulation experiments. Finally, the proposed model is applied to study the influence of environmental, economic, and medical levels on life expectancy in China. In the study functional predictors are the daily air quality index, temperature, and humidity of 58 cities in 2020 and scalar predictors are GDP and number of beds in hospitals.

In this study, we present different methods of estimating fuzzy reliability of a two-parameter Rayleigh distribution via the maximum likelihood estimator, median first-order statistics estimator, quartile estimator, L-moment estimator, and mixed Thompson-type estimator. The mean-square error MSE as a measurement for comparing the considered methods using simulation through different values for the parameters and unalike sample sizes is used. The results of simulation show that the fuzziness values are better than the real values for all sample sizes, as well as the fuzzy reliability at the estimation of the Maximum

We introduce the lower chaos grade of a real-valued function F defined on the Markov triple (E,μ,Γ), where μ is a probability measure and Γ is the carré du champ operator. As an application of this concept, we obtain the better estimate of the four moments theorem for Markov diffusion generators worked by Bourguin et al (2019). For our purpose, we need to find the largest number except zero in the set of eigenvalues corresponding to its eigenfunction in the case where the square of a random variable F, coming from a Markov triple structure, can be expressed as a sum of eigenfunctions, We give some examples of eigenfunctions of the diffusion generators such as Ornstein-Uhlenbeck, Jacobi and Romanovski-Routh. In particular, two bounds, called the four moments theorem and fourth moment theorem respectively, will be provided for the normal approximation of the case where a random variable F comes from eigenfunctions of a Jacobi generator.

Positive recurrence for a single-server queueing system is established under generalised intensity conditions of service which does not assume existence of the density distribution function of service but a certain integral type lower bound as a sufficient condition. Positive recurrence implies existence of the invariant distribution and a guaranteed slow convergence to it in the total variation metric.

In this paper, we identify several relative ordering properties of the modified proportional hazard rate and modified proportional reversed hazard rate models. For this purpose, we use two well-known relative orderings, namely the relative hazard rate ordering and the relative reversed hazard rate ordering. The investigation is to see how a relative ordering between two possible base distributions for the response distributions in these models is preserved when the parameters of the underlying models are changed. We will give some examples to illustrate the results and the conditions under which they are obtained.

Financial technologies, stemming from the application of artificial intelligence to big data in finance, are continuously expanding, across different markets and financial services. While financial technologies bring many opportunities, such as reduced costs and extended inclusion, they also bring risks, among which cyber risks, which are constantly increasing and are difficult to measure. Among the difficulties in measurement lies the existence of interdependence among different cyber risks. The study of interdependence and possible contagion channels between cyber attacks to different institutions and economic sectors is indeed increasingly important to ensure economic and financial sustainability. Against this backdrop, this paper proposes a multivariate model for count time series of cyber risk events, in which the time-varying intensity parameter determining the probability that a cyber attack occurs evolves according to general autoregressive score models, taking both time and sectorial dependence into account. The model is particularly suitable for studying how the behaviors of different markets or sectors are interconnected and it constitutes a new approach to the multivariate analysis of count time series of cyber loss events.

In this paper, we investigate the cumulative historical entropy of the lifetime of an n-component coherent system under the assumption that all system components have failed at time t. Using the system signature, we can calculate the cumulative historical entropy of the system’s lifespan, which is an important aspect in determining a system’s predictability. This research produces a number of significant results, including expressions, bounds, and ordering aspects for this measure. In addition, we offer a method for finding a preferred system based on cumulative Kullback-Leibler discriminating information that is closely related to the parallel system. These findings provide important insights into the predictability and dependability of coherent systems, emphasizing the importance of this field.

This paper studies single change-point detection in the volatility of a class of parametric conditional heteroscedastic autoregressive nonlinear (CHARN) models. The conditional least-squares (CLS) estimators of the parameters are defined and are proved to be consistent. A Kolmogorov-Smirnov type-test for change-point detection is constructed and its null distribution is provided. An estimator of the change-point location is defined. Its consistency and its limiting distribution are studied in detail. A simulation experiment is carried out to assess the performance of the results which are also applied to two sets of real data.

The impacts of COVID-19 (Novel Coronavirus) epidemic cannot have been more severe, with the globe experiencing both economic and health crises. This study examines the trends and correlations in number of COVID-19-related deaths and number of COVID-19-infected patients in all 37 regions of Tamil Nadu state, India, in the month of August 2020 on the basis of panel regression model,

Most existing data synthesis methods are designed to tackle problems such as dataset imbalance, data anonymization and insufficient sample size. There is a lack of effective synthesis methods for the limited number of datasets which contain a large of features and unknown noise to expand the size of the dataset. We propose a data synthesis method, named Adaptive Subspace Interpolation for Sample Optimization (ASISO). The idea is to divide the original feature space into several subspaces with an equal number of samples, and then perform interpolation for the samples in the adjacent subspaces. This method can adaptively adjust the size of the dataset containing unknown noise, and the expanded data typically contain minimal error with actual. Moreover, it adjusts the structure of the samples, which can significantly reduce the proportion of samples with large errors. In addition, the hyperparameters of this method have an intuitive explanation and usually require little calibration. Experimental results on artificial data and benchmark data sets demonstrate that ASISO is a robust and stable method to optimize samples.

In this paper we study a special kind of finite mixture model. The sample drawn from the model consisits of three parts. The first two parts are drawn from specified distributions f1 and f2, while the third one is drawn from the mixture. A problem of interest is whether the two distributions f1 and f2 are the same. To test this hypothesis, we first define the regular location and scale family of distributions and assume that f1 and f2 are regular denisty functions. Then the hypothesis transforms to the equalities of the location and scale parameters, respectively. To utilize the information in the sample, we use Bayes’ theorem to obtain the posterior distribution and give the sampling method. Then we propose the posterior p-value to test the hypothesis. The simulation studies show that our posterior p-value largely improves the power in both normal and logistic cases while nicely controls the Type-I error. A real halibut dataset is used to illustrate the validity of our method.

We study the robustness of the cμ-rule for the optimal allocation of a resource consisting of one unreliable server to parallel queues with two different classes of customers. The customers in queues can be served with respect to a FCFS retrial discipline, when the customers at the head of queues repeatedly tries to occupy the server in a random time. It is proved that for the scheduling problem in the system without arrivals the cμ-rule minimizes the total average cost. For the system with arrivals it is difficult directly to prove the optimality of the same policy with explicit relations. We derived for an infinite-buffer model a static control policy that also prescribes for the certain values of system parameters the service exclusively for the class-i customers if both of queues are not empty with the aim to minimize the average cost per unit of time. It is also shown that in a finite-buffer case the cμ-rule fails.

An empirical investigation was carried out to study the relationship between the soil characteristics, including the Nitrogen (N), Phosphorus (P), Potassium (K), potential of Hydrogen (pH), Electrical conductivity (EC), Sulfur (S,), Zinc (Zn), Iron (Fe), Copper (Cu), Manganese (Mn) and Boron (B) using principal component analysis (PCA), Factor Analysis (FA), and Canonical Correlation analysis (CCA) for data reduction multivariate techniques. The first five main components accounted for 24.22%, 18.52%, 16.1%, 9.97%, and 9.27% of sample variances, respectively. All five principal components (PCs) accounted for 78.07% of the sample variance. The soil characteristics viz., P and K have highly dominated the first PC; the soil characteristics N and pH have dominated the second PC; K, EC, and Cu have dominated the third PC; Fe, Mn, and pH have dominated in the fourth PC; whereas the parameter Cu dominated the fifth PC. The first factor showed a strong negative loading on S and a strong favorable loading on Cu; P, K, and EC have substantial positive loadings in the second factor. Significant positive loadings on Mn, Fe, Cu, pH, N, and K are present in the third component. The fourth factor had a significantly positive pH, Fe, N, P, K, and B loadings. Fe, P, and B have much weight in the fifth factor. The soil characteristics viz., N, Zn, pH, K, Fe, and Mn have dominated the first PC; the soil characteristics B, P, S, Zn, and Fe have dominated the second PC; B and P have dominated the third PC whereas the only character Cu dominated the fourth PC. Cu and OC have much weight in the fourth factor. The canonical redundancies for dependent and independent sets are 12% and 17%, respectively. The Stewart-Love canonical redundancy index is 70.58%, which means that the first linear combination of the X-set explains 70.58% of the total variance in the Y-set.

To estimate the degree of quantum entanglement, it is important to understand the statistical behavior of functions of spectrum of density matrices such as von Neumann entropy, quantum purity, and entanglement capacity. These entangle- ment metrics over different generic state ensembles have been studied intensively in the literature. As an alternative metric, in this work we study sum of square root spectrum of density matrices, which is relevant to negativity and fidelity in quantum information processing. In particular, we derive the exact mean and vari- ance of sum of square root spectrum over the Bures-Hall generic state ensemble extending known results obtained recently over the Hilbert-Schmidt ensemble.

Distributions with bounded support show considerable sparsity over those with unbounded support, despite the fact that there are a number of real-world contexts where observations take values from a bounded range (proportions, percentages, and fractions are typical examples). For proportion modelling, a flexible family of two-parameter distribution functions associated with the exponential distribution is proposed here. Mathematical and statistical properties of the novel distribution are examined, including quantiles, mode, moments, hazard rate function, and its characterization. The parameter estimation procedure using the maximum likelihood method was carried out, and applications to environmental and engineering data were also considered. To this end, various statistical tests are used, along with some other information criterion indicators to determine how well the model fits the data. The proposed model is found to be the most efficient plan in most cases for the datasets considered.

Generalized Bayes is a Bayesian approach based on the learning rate parameter η. In this study, we examine the effect of parameter η on the estimation results considering joint type-I and type-II hybrid censored samples from k exponential populations. In addition to the learning rate parameter, we consider two loss functions, the Linex and general entropy loss functions in the Bayesian approach. Monte Carlo simulations are performed to compare the performances of the estimation results under losses and different values of η. An illustrative example is performed to study the effect of the learning rate parameter and the different losses with different parameters.

In this paper a more flexible extension of the Fréchet distribution is introduced. The new distribution is defined by means of the stochastic representation as the quotient of two independent random variables, a Fréchet distribution and the power of a random variable with uniform distribution in the interval (0,1). We will call this new extension the Slash Fréchet distribution and one of its main characteristics is that its tails are heavier than the Fréchet distribution. The general density of this distribution and some basic properties are determined. Its moments, skewness coefficients and kurtosis are calculated. In addition, the estimation of the model parameters is obtained by the method of moments and maximum likelihood. Finally, two applications with real data are performed by fitting the new model and comparing it with the Fréchet distribution.

In many research fields, statistical probability models are often used to analyze real-world data. However, data from many fields, such as the environment, economics, and health care, do not necessarily fit traditional models. New empirical models need to be developed to improve their fit. In this paper, we explore a further extension of the quasi-Lindley model. Maximum likelihood, least square error, Anderson-Darling algorithm, and expectation-maximization algorithm are four techniques for estimating the parameters under study. All techniques provide accurate and reliable estimates of the parameters. However, the mean square error of the expectation maximization approach was lower. The usefulness of the proposed model was demonstrated by analyzing a dynamical systems data set, and the analysis shows that it outperforms the other models in all statistical models considered.

In estimating logistic regression models, convergence of the maximization algorithm is critical, however, this may fail. Numerous bias correction methods for maximum likelihood estimates of the parameters have been conducted for cases of complete data sets and also longitudinal models. For binary response fixed effects panel data model, the conditional logit estimator is consistent for balanced data. When faced with missing covariates problem, researchers opt for various imputation techniques to make the data complete and without loss of generality consistent estimates still suffice asymptotically. For maximum likelihood estimates of the parameters for logistic regression in cases of imputed covariates, the optimal choice of an imputation technique which yields the best estimates with minimum variance is still elusive. The main aim of this paper is to examine the behaviour of the Hessian matrix with optimal values of the imputed covariates vector which will make the Newton-Raphson algorithm to converge faster through a reduced absolute value of the product of the score function and the inverse fisher information component. We focus on a method used to modify the conditional likelihood function through partitioning of the covariate matrix. We also confirm from the moduli of the Hessian matrices that the log likelihood of a panel data logistic model has a global maximum as the parameter estimates. Simulation results reveal that model based simulation perform better than classical imputation techniques yielding estimates with smaller bias and higher precision for the conditional maximum likelihood estimation of nonlinear panel models with single fixed effects.

Analyzing big data poses a great challenge for numerous researchers to explore the data structure. Dimension reduction methods can be used to reduce data dimensionality, taking it from occupying a high-dimensional space to existing in a lower-dimensional space while retaining as much information as possible. Principal Component Analysis is one of the most popular used to reduce the dimensional space. The bootstrap sample is obtained by randomly sampling n times with replacement from the original sample, the method provides easy tool to understand the interactive component and develop the process. There are not enough researches discussed the stability of PCA method using bootstrap method. In this paper, the bootstrap method is used to analyze the stability of PCA results and to estimate the number of PCA in efficient way. The method is used to estimate the number of PCA that needed to classify the data set, and the effectiveness of the discussed techniques is demonstrated through real data sets.

This paper focuses on a joint model to analyse longitudinal proportion and survival data. After transforming the longitudinal proportional data by a logit function, we adopt a partially linear mixed-effect model for it, in which nonlinear function of time is fitted via using B-splines technique and a centered Dirichlet Process Mixture Model (CDPMM) is specified for a general distribution of random effects. The survival data is assumed to a Cox proportional hazard model, the sharing random effects joint model is developed for the two types of data. Combining the Gibbs sampler and the Metropolis-Hastings algorithm, we propose a Bayesian Lasso (BLasso) method to simultaneously estimate unknown parameters and select important covariates. Simulation studies are conducted to investigate the finite sample performance of the proposed methods. An example from the MA.5 research experiment is used to illustrate the proposed methodologies.

Nonparametric density estimation for nonnegative data is considered in a situation where a random sample is not directly available but the data are instead observed from the length-biased sampling. Due to the boundary bias problem of the location-scale kernel, the approach in this paper is an application of asymmetric kernel. Two nonparametric density estimators are proposed. The mean integrated squared error, strong consistency, and asymptotic normality of the estimators are investigated. Some simulations illustrate the finite sample performance of the estimators.

In econometrics and finance, volatility modelling has long been a specialised field for addressing a variety of issues pertaining to the risk and uncertainties of an asset. This study presents a robust framework, through a step-by-step design, that is relevant for effective Monte Carlo simulation (MCS) with empirical verifications to estimate volatility using the Generalized Autoregressive Score (GAS) model. The framework describes an organised approach to the MCS experiment that includes "background (optional), defining the aim, research questions, method of implementation, and summarised conclusion". The method of implementation is a workflow that consists of writing the code, setting the seed, setting the true parameter a priori, data generation process, and performance assessment through meta-statistics. Among the findings, it is experimentally demonstrated in the study that the GAS model with a lower unconditional shape parameter value can generate a dataset that adequately reflects the behaviour of financial time series data, relevant for volatility estimation. This dynamic framework is intended to help interested users on MCS experiments utilising the GAS model for reliable volatility calculations in finance and other areas.

This study aims to apply a Recommendation System with Content Based Filtering method with Jaccard Similarity and Cosine Similarity algorithms on the E-Learning Platform. Recommendation systems deal with how to provide personalized recommendations to users efficiently. The Content Based Filtering method with Jaccard Similarity and Cosine Similarity algorithms can be used to calculate the similarity value between E-Courses on the E-Learning Platform. Implementation for Recommendation System using Google Colaboratory with Python programming language. In the application of the Recommendation System dataset, Coursera Free Dataset consists of 975 instances. The recommendation results use the Jaccard Similarity algorithm with an average similarity value of 0.3 while the value of Cosine Similarity with the average similarity value is 0.6 where the similarity value of Cosine Similarity is higher. Based on the results of the Mean Absolute Error in the low recommendation system, the average MAE value for all iterations of Jaccard Similarity algorithm is 0.013 and for the Cosine Similarity algorithm the average MAE value for all iterations is 0.014. This shows that the Recommendation System with Jaccard Similarity and Cosine Similarity algorithms can be used on the E-Learning Platform to provide efficiency solutions for personalized recommendations.

Asymmetry in the upper and lower tails is an important feature in modeling bivariate distributions. This article focuses on the log ratio between the tail probabilities at upper and lower corners as a measure of tail asymmetry. Asymptotic behavior of this measure at extremely large and small thresholds is explored with particular emphasis on the skew-normal copula. Our numerical studies reveal that, when the correlation or skewness parameters are around at the boundary values, some asymptotic tail approximations of the skew-normal copulas proposed in the literature are not suitable to compute the measure of tail asymmetry with practically extremal thresholds.

In recent decades, there have been numerous endeavors to develop a novel category of survival distributions possessing enhanced flexibility through the extension of existing distributions. This article constructs and validates the statistical properties of a novel survival distribution in order to obtain an alternative distribution that is suitable for analyzing survival data by presenting the novel mixture of the Fréchet distribution along with statistical properties such as the probability density function (PDF), cumulative distribution function (CDF), rth ordinary moment, skewness, kurtosis, moment-generating function, mean, variance, mode, survival function, hazard function, and asymptotic behavior, as well as constructing the estimators of the unknown parameter by employing the expectation-maximization (EM) algorithm, and simulated annealing. Additionally, the performance of the proposed estimators was compared with bias, mean squared errors (MSE), and simulated variances, and given an illustrative example of the proposed distribution to the survival data set in order to show that the proposed distribution is appropriate for the right-skewed data. This will be extremely advantageous in survival analysis.

In this paper, the problem of predicting future failure times based on a jointly type-II censored sample from k exponential populations is considered. The Bayesian prediction intervals and point predictors were then obtained. Generalized Bayes is a Bayesian study based on a learning-rate parameter. This study investigated the effects of the learning rate parameters on the prediction results. The loss functions of squared error, Linex, and general entropy were used as point predictors. Monte Carlo simulations were performed to show the effectiveness of the learning rate parameter in improving the results of prediction intervals and point predictors.

High-dimensional measurement error data are becoming more prevalent across various fields. Research on measurement error regression models has gained increasing attention due to the risk of drawing inaccurate conclusions if measurement errors are ignored. When the dimension p is larger than the sample size n, it is challenging to develop statistical inference methods for high-dimensional measurement error regression models due to the existence of bias, nonconvexity of objective function, high computational cost and many other difficulties. Over the past few years, some works have overcome the aforementioned difficulties and proposed several statistical inference methods. This paper mainly reviews the current development on estimation, hypothesis testing and variable screening methods for high-dimensional measurement error regression models and shows the theoretical results of these methods with some directions worthy of exploring for future research.

The aim of this study is to obtain the Bayes estimators and the maximum likelihood estimators (MLEs) for the unknown parameters of the Rayleigh Weibull (RW) distribution based on progres-sive type-II censored samples. The approximate Bayes estimators are calculated using the idea of Lindley and Tierney-Kadane's approximation method under the squared-error loss function when the Bayes estimators are not handed in explicit forms. In this study, the approximate Bayes esti-mates are compared with the maximum likelihood estimates in the aspect of the estimated risks (ERs) using Monte Carlo simulation. In addition, the coverage probabilities of the parametric bootstrap estimates are calculated. Real lifetime data sets belonging to the cancer types as bladder cancer, head and neck cancer, and leukemia are used to illustrate the emprical results belonging to the approximate Bayes estimates, the maximum likelihood estimates, and the parametric bootstrap intervals.

This study rolls out a robust framework relevant for simulation studies through the Generalised Autoregressive Conditional Heteroscedasticity (GARCH) model using the rugarch package. The package is thoroughly investigated, and novel findings are identified for improved and effective simulations. The focus of the study is to provide necessary simulation steps for volatility estimation that involve "background (optional), defining the aim, research questions, method of implementation, and summarised conclusion". The method of implementation is a workflow that includes writing the code, setting the seed, setting the true parameters a priori, data generation process and performance assessment through meta-statistics. This novel, easy-to-understand steps are demonstrated on financial returns using illustrative Monte Carlo simulation with empirical verification. Among the findings, the study shows that regardless of the arrangement of the seed values, the efficiency and consistency of an estimator generally remain the same as the sample size increases. The study also derived a new and flexible true-parameter-recovery measure which can be used by researchers to determine the level of recovery of the true parameter by the MCS estimator. It is anticipated that the outcomes of this study will be broadly applicable in finance, with intuitive appeal in other areas, for volatility estimation.

Control charts have been widely used for monitoring process quality in manufacturing and play an important role in triggering a signal in time when detecting a change in process quality. Many control charts in literature assume that the in-control distribution of the univariate or multivariate process data is continuous and not categorical. This research develops two exponentially weighted moving average (EWMA) proportion control charts for monitoring a process with multinomial proportions when considering both large and small sample sizes. For a large sample size , the charting statistic depends on the well-known Pearson χ2 statistic, and the control limit of the EWMA proportion chart is determined by an asymptotical chi-square distribution. For a small sample size, we derive the exact mean and variance of the Pearson χ2 statistic. Hence, the exact EWMA proportion chart is determined. The proportion chart can also be applied to monitor the distribution-free continuous multivariate process as long as each categorical proportion associated with specification limits of each quality variable is known or estimated. Lastly, we investigate the detection performance of the proposed EWMA proportion chart by numerical analyses. Real data analysis demonstrates the beneficial application of the proposed EWMA proportion charts.

In actuarial practice, the modeling of total losses tied to a certain policy is a non-trivial task. Traditional parametric models to predict total losses have limitations due to complex distributional features such as extreme skewness, zero inflation, multi-modality, etc., and the lack of explicit solutions for log-normal convolution. In the recent literature, the application of the Dirichlet process mixture for insurance loss has been proposed to eliminate the risk of model misspecification biases; however, the effect of covariates as well as missing covariates in the modeling framework is rarely studied. In this article, we propose novel connections among covariate-dependent Dirichlet process mixture, log-normal convolution, and missing covariate imputation. Assuming an individual loss is log-normally distributed, we develop a log skew-normal Dirichlet process to approximate the log-normal sum. As a generative approach, our framework models the joint of outcome and covariates, which allows to impute missing covariates under the assumption of missingness at random. The performance is assessed by applying our model to several insurance datasets, and the empirical results demonstrate the benefit of our model compared to the existing actuarial models such as the Tweedie-based generalized linear model, generalized additive model, or multivariate adaptive regression spline.

For order-of-addition experiments, the response is affected by the addition order of the experimental materials. Consequently, the main interest focuses on creating a predictive model and an optimal design for optimizing the response. Van Nostrand (1995) proposed the pairwise-order (PWO) model for detecting PWO effects. Under the PWO model, the full PWO design is optimal under various criteria but is often unaffordable because of the large run size. In this paper, we consider the D-, A- and M.S.-optimal fractional PWO designs. We first present some results on information matrices. Then, a flexible and efficient algorithm is given for generating these optimal PWO designs. Numerical simulation shows that the generated design has an appealing efficiency in comparison with the full PWO design, though with only a small fraction of runs. Several comparisons with existing designs illustrate that the generated designs achieve better efficiencies, and the best PWO designs and some selected 100% efficient PWO designs generated by the new algorithm are reported.

Polynomial recurrence bounds for a class of stochastic differential equations with a gradient type drift and an additive Wiener process are studied without Lyapunov functions.

Let {ξ1, ξ2. . . .} be a sequence of independent possibly differently distributed random variables defined on a probability space (Ω, F, P) with distribution functions {F_ξ1., F_ξ2 , . . .}. Let η be a counting random variable independent of sequence {ξ1, ξ2. . . .}. In this paper, we find conditions under which distribution function of randomly stopped sum Sη = ξ1 + ξ2 + . . . + ξη belongs to the class of generalized subexponential distributions.

We give two examples of the appearance of heavy-tailed distributions in applications to social sciences. Among these distributions are the laws of Pareto, Lotka, and some new ones. The examples are illustrated by constructing suitable toy models.

The multivariate Gaussian random fields with matrix-based scaling laws are widely used for inference in statistics and many applied areas. In such contexts interests are often in symmetry and in the rates of change of spatial surfaces in any given direction. This article analyzes the almost sure sample function behavior for operator fractional Brownian motion, including multivariate fractional Brownian motion. We obtain the estimations of small ball probability and the strongly locally nondeterministic for operator fractional Brownian motion in any given direction. Applying these estimates we obtain Chung's laws of the iterated logarithm for spatial surfaces of operator fractional Brownian motion. Our results show that the precise rates of change of spatial surfaces are completely determined by the self-similarity exponent.

In this paper, under the symmetric entropy and the scale squared error loss functions, we consider the maximum likelihood (ML) estimation and Bayesian estimation of the Shannon entropy and Rényi entropy of the two-parameter inverse Weibull distribution. In the ML estimation, the dichotomy is used to solve likelihood equation. In addition, the approximation confidence interval is given by the Delta method. Because the form of estimation results is more complex in the Bayesian estimation, the Lindley approximation method is used to achieve the numerical calculation. Finally, Monte Carlo simulations and a real data set are used to illustrate the results derived. By comparing the mean square error between the estimated value and the real value, it can be found that the performance of ML estimation of Shannon entropy is better than that of Bayesian estimation, and there is no significant difference between the performance of ML estimation of Rényi entropy and that of Bayesian estimation.

In the first part of this article, we discuss and generalize the complete convergence introduced by Hsu and Robbins (1947) to the r-complete convergence introduced by Tartakovsky (1998). We also establish its relation to the r-quick convergence first introduced by Strassen (1967) and extensively studied by Lai (1976). Our work is motivated by various statistical problems, mostly in sequential analysis. As we show in the second part, generalizing and studying these convergence modes is important not only in probability theory but also to solve challenging statistical problems in hypothesis testing and changepoint detection for general stochastic non-i.i.d. models.

The article proposes a method for determining the p-value of correlations in compositional data, i.e., those data that arise as a result of dividing the original values by their sum. Data organized in this way are typical for many fields of knowledge, but there is still no consensus on methods for interpreting correlations in such data. In a space closed by normalizing quantity, correlation coefficients behave differently than under normal conditions: their probabilities of occurrence do not coincide with those inherent in the standard scale of estimates. In the tens of the new millennium, almost all newly emerging methods for estimating correlation in compositional data began to require mandatory log-transformation of the variable values. In the method proposed here there are no log-transformations. We return to the early stages of attempting to solve the problem and rely on negative shifts in correlations in the multinomial distribution. In modeling the data, we use a hybrid method that combines the hypergeometric distribution with the distribution of any other law. During our work on the calculation method, we found that the number of degrees of freedom in compositional data measures discretely only when all the normalizing sums are equal and that it decreases when the sums are not equal, becoming a continuously varying quantity. Estimation of the number of degrees of freedom and the strength of its influence on the magnitude of the shift in the distribution of correlation coefficients is the basis of the proposed method.

L\'evy walks represent important modeling tools for variety of real life processes. Their natural scaling limits are known to be described by the so-called material fractional derivatives. So far these scaling limits were derived for spatially homogeneous walks, where Fourier and Laplace transforms represent natural tools of analysis. Here we derive the corresponding limiting equations in the case of position depending times and velocities of walks, where Fourier transforms cannot be effectively applied. In fact, we derive three different limits (specified by the way the process is stopped at an attempt to cross the boundary), leading to three different multi-dimensional versions of Caputo-Dzherbashian derivatives, which correspond to different boundary conditions for the generators of the related Feller semigroups and processes. Some other extensions and generalisations are analysed.

In this paper, we will extend the expected value of the function w.r.t the uniform probability measure on sets measurable in the Caratheodory sense to be finite for a larger class of functions, since the set of all measurable functions with infinite or undefined expected values forms a prevalent subset of the set of all measurable functions, which means "almost all" measurable functions have infinite or undefined expected values. Before we define the specific problem in section 2, we will outline some preliminary definitions. We'll then define the specific problem (along with a partial solution in section 3) to visualize the complete solution. Along the way, we will ask a series of questions to clarify our understanding of the paper.

Quantile regression model is widely used in variable relationship research of general size data, due to strong robustness and more comprehensive description of the response variables' characteristics. With the increase of data size and data dimension, there have been some studies on high-dimensional quantile regression under the classical statistical framework, including higher-efficient frequency perspective, which is however at cost of randomness quantification, or lower-efficient Bayesian method based on MCMC sampling. To overcome these problems, we propose the high-dimensional quantile regression with Spike-and-Slab Lasso penalty based on variational Bayesian (VBSSLQR), which can not only improve the computational efficiency but also measure the randomness via variational distributions. The simulation studies and real data analysis illustrate that the proposed VBSSLQR method is superior to or equivalent to other quantile and non-quantile regression methods (including Bayesian and non-Bayesian methods), and its efficiency is higher than any other method.

This study considers the change-point test problem for time series based on the self-normalization ratio statistic test, which is constructed using residuals obtained from a support vector regression (SVR)-autoregressive moving average (ARMA) model. Under the null hypothesis, the series is a stationary process, and our test statistic converges to a non-degenerate distribution. Under the alternative hypothesis, there are change-points in the time series, and the self-normalization test statistic diverges to infinity. The simulations show that our proposed new test has better finite sample performance than other SVR-based tests in the literature. Finally, we illustrate its usefulness by analyzing two actual data sets.

This article presents six probability distributions from the Gamma family with three parameters, for the flood frequency analysis in hydrology. The choice of the Gamma family of statistical dis-tributions was driven by its frequent use in hydrology. In the Faculty of Hydrotechnics, the im-provement of the estimation of maximum flows and including the methodological bases for the realization of a regionalization study with the linear moments method with the corrected pa-rameters was researched, being an element of novelty. The linear moments method is better than MOM because it avoids the choice of skewness depending on the origin of the flows, practiced in Romania. The L-moments method conforms to the current trend for estimating the parameters of statistical distributions. Observed data from hydrometric stations are of relatively short length, so the statistical parameters that characterize them are of a sample that requires correction. The correction of the statistical parameters is proposed, using the method of least squares based on the inverse functions of the statistical distributions expressed with the frequency factor for L-moments. All the necessary elements for their use are presented like, quantile functions, the exact and ap-proximate relations for estimating parameters and frequency factors. A flood frequency analysis case study was carried out for the Ialomita river, to verify the proposed methodology. The per-formance of this distributions is evaluated using Kiling-Gupta and Nash-Sutcliff coefficients.

This article presents the results of the consideration of the Bertrand paradox in terms of geometric probability theory. Although the method and conclusions of this article are straightforward, no equivalent studies were found when reviewing the relevant literature. This could be conditioned by the fact that the hypotheses presented in this research have low intuitive obviousness and, in contrast, could be due to the historically established agreement regarding the issue. This article shows that out of three classical solutions to the problem described by Bertrand, two methods are inconsistent with the claimed relative objectivity. Although the remaining solution (1/4) seems to be the most correct, we cannot claim that we have exhaustively ruled out all aspects that could reduce its adequacy to solve the problem.

Multilevel functional principal components analysis (mFPCA) is used to treat dynamical changes in shape in this article. Results of standard (single-level) FPCA are presented here also as a comparison. Monte Carlo (MC) simulation is used to create univariate data (i.e., a single “outcome” variable) that contains two distinct classes of trajectory with time. MC simulation is also used to create multivariate data of sixteen 2D points that represent (broadly) an eye; this data also has two distinct classes of trajectory (an eye blinking and an eye widening in surprise). This is followed-up by an application of mFPCA and single-level FPCA to “real” data consisting of twelve 3D landmark outlining the mouth that are tracked over all phases of a smile. By consideration of eigenvalues, results for the MC datasets find correctly that variation due to differences in groups between the two classes of trajectories variation are larger than variation within each group. In both cases, differences in standardized component scores between the two groups are observed, as expected. Modes of variation are shown to model the univariate MC data correctly, and good model fits are found for both the “blinking” and “surprised” trajectories for the MC “eye” data. Results for the “smile” data show that the smile trajectory is modelled correctly, namely: the corners of the mouth are drawn backwards and wider during a smile. Furthermore, the first mode of variation at level 1 of the mFPCA model show only subtle and minor changes in mouth shape due to sex, whereas the first mode of variation at level 2 of the mFPCA model govern whether the mouth is upturned or downturned. These results are all an excellent test of mFPCA, which show that mFPCA presents a viable method of modelling dynamical changes in shape.

Building fires are preventable incidents that have proven to be both deadly and costly. Addressing their root causes will lead to safer neighbourhoods for families and businesses to live and operate in. Multiple studies have established the effect of residents’ socioeconomic compositions on an area’s building fire rates; however, the existing model based on the classical stepwise regression procedure has several limitations. This paper aims to construct a more accurate predictive model of building fire rates based on a set of explanatory socioeconomic variables. In building the socioeconomic model, a backward elimination by Robust Final Predictor Error (RFPE) criterion is proposed to enhance the forecasting capability of the model. The proposed method has been implemented on the census data and the fire incident data of the South East Queensland region in Australia. A cross-validation was then conducted to assess the model’s accuracy. In addition, comparative analyses of other elimination criteria, such as p-value, adjusted R-squared, Akaike’s Information Criterion (AIC), Bayesian Information Criterion (BIC) and predicted residual error sum of squared (PRESS), were conducted. The cross-validation analyses demonstrate that the proposed criterion is a more accurate predictive model based on a couple of goodness-of-fit measures. All in all, the RFPE equation was found to be a suitable criterion for the backward elimination procedure in the socioeconomic modelling of building fires.

We consider a new model of a branching random walk on a multidimensional lattice with continuous time and one source of particle reproduction and death, as well as an infinite number of sources in which, in addition to the walk, only absorption of particles can occur. The asymptotic behavior of the integer moments of both the total number of particles and the number of particles at a lattice point is studied depending on the relationship between the model parameters. In the case of the existence of an isolated positive eigenvalue of the evolution operator of the average number of particles, a limit theorem is obtained on the exponential growth of both the total number of particles and the number of particles at a lattice point.

Bootstrapping is a flexible, powerful and well-established statistical approach to quantify the uncertainty of virtually any point estimate. While multiple versions of bootstrap confidence intervals are already available in Stata, dbs implements the double (iterated) bootstrap. Instead of relying on parametric assumptions such as the non-parametric resampling bootstrap confidence interval does, it is more flexible and derives critical values directly from that data. To do so, multiple methods are available (analytic approach, double resampling, jackknife estimation). In a comparative simulation study it is empirically demonstrated that the strengths of the double bootstrap are particularly evident for small samples (n < 100) when heteroscedasticity is present. While all other approaches result in undercoverage, only the double bootstrap reaches the target coverage level and hence avoids incorrect statistical conclusions. The computational burden is not even necessarily larger than for other bootstrap approaches.

In the era of big data with increasingly complex data structures and ever-larger data scales, matrix-type data are becoming highly valued and their applications in the fields of medicine, industry, education, geography, and astronomy are growing in extent. In recent years, significant progress has been made in the practical use of matrix variable t-distribution finite mixture models for handling data in order to address the issues of multi-subgroup structures and long data tails. In this paper, the expectation-maximization (EM) algorithm with penalized maximum likelihood is proposed to resolve the problem of the unbounded nature of the likelihood function applied to the model by considering the degeneracy of the variance-covariance matrix of this model. Our data were analyzed through simulations and real data, and the results demonstrate that our model is effective in both preventing the likelihood function from being unbounded and in ensuring the accuracy of the estimated parameters of the EM algorithm.

Many papers in recent years have been devoted to estimating the per pixel proportions of three broad classes of materials (e.g. photosynthetic vegetation, non-photosynthetic vegetation and bare soil) using data from multispectral sensors. Many of these papers use estimation methods based on the linear mixture model. Very few of these papers assess the accuracy of their estimators. I show how to produce confidence intervals (CIs) and joint confidence regions (JCRs) for the proportions associated with various linear mixture models. There are two main models, both of which assume that the coefficients in the model are non-negative. The first model assumes that the coefficients sum to 1. The second does not, but uses rescaling of the estimated coefficients to produce estimated proportions. Three variants of these two models are also analysed. JCRs are shown to be particularly informative, because they are typically better at localising the information than CIs are. The methodology is illustrated using examples from Landsat Thematic Mapper data at 1169 locations across Australia, each of which has associated field observations. There is also discussion about the extent to which the methodology can be extended to hyperspectral data.

The opioid crisis in the United States has had devastating effects on communities across the country, leading many states to pass legislation that limits the prescription of opioid medications in an effort to reduce the number of overdose deaths. This study evaluates the impact of two categories of PDMP and Pill Mill regulations on the supply of opioid prescriptions at the level of dispensers and distributors (excluding manufacturers) using ARCOS data. The study uses a difference-in-difference method with a two-way fixed design to analyze the data. The study finds that both of the regulations are associated with reductions in the volume of opioid distribution. However, the study reveals that these regulations may have unintended consequences, such as shifting the distribution of controlled substances to neighboring states. For example, in Tennessee, the implementation of Operational PDMP regulations reduces in-state distribution of opioid drugs by 3.36% (95% CI, 2.37 to 4.3), while the out-of-state distribution to Georgia, which did not have effective PDMP regulations in place, increases by 16.93% (95% CI, 16.42 to 17.44). Our studies emphasize that policymakers should consider the potential for unintended distribution shifts of opioid drugs to neighboring states with laxer regulations as well as varying impacts on different dispenser types.

For each $t\in(-1,1)$, exact values of the least upper bounds $$ H(t)=\sup_{\E X=t,\,\E X^2=1} \frac{\E\abs{X}^3}{\E \abs{X-t}^3},\quad \sup_{\E X=t,\,\E X^2=1} \frac{L_1(X)}{L_1(X-t)} $$ are obtained, where $L_1(X)=\E|X|^3/(\E X^2)^{3/2}$ is the non-cental Lyapunov ratio. It is demonstrated that these values are attained only at two-point distributions. As a corollary, S.\,Shorgin's conjecture is proved that states that the exact value is $$ \sup\frac{L_1(X)}{L_1(X-\E X)}= \frac{\sqrt{17 + 7\sqrt7}}{4} = 1.4899\ldots, $$ where the supremum is taken over all non-degenerate distributions of the random variable $X$ with the finite third moment. Also, in terms of the central Lyapunov ratio $L_1(X-\E X)$, an analog of the Berry--Esseen inequality is proved for Poisson random sums of independent identically distributed random variables with the constant $$ 0.3031\cdot H\left(\frac{\E X}{\sqrt{\E X^2}}\right) \left(1-\frac{(\E X)^2}{\E X^2}\right)^{3/2} \le 0.3031\cdot \frac{\sqrt{17 + 7\sqrt7}}{4}<0.4517.$$ where $\Law(X)$ is the common distribution of the summands.

In recent times, there has been a growing interest in the use of control charts as a feedback process monitoring technique. Particularly, CUSUM schemes have been proven to be efficient for monitoring processes where the magnitude of the sustained mean shifts is small or moderate. However, because of the constant control limits, the detection ability of CUSUM schemes becomes slow at the initial set-up of the process. The fast initial response FIR feature has often been used to enhance it at the process startup. Meanwhile, the dynamism of real-life problems has always encouraged a more sensitive CUSUM scheme capable of detecting process shifts more rapidly. In this paper, an improved CUSUM scheme for monitoring process mean is proposed. This scheme will be substantial for monitoring processes whose observations are obtained at a distant time interval for example hourly, daily or weekly and where the sustained shifts are assumed to be small or moderate. To demonstrate the practical applications of the proposed scheme, we present its real-life application using datasets from a bottling company and petroleum refinery laboratory.

Statistical modeling of lifetime data plays an essential role in a wide range of practical fields, such as health and engineering. There have been a lot of studies done to develop statistical models that can better describe health data than traditional models. For the first time, we pioneer a novel family of continuous probability distributions called the generalized odd beta prime generalized (GOBP-G) family of distributions. The cumulative distribution and probability density functions of the new family are presented. A new generalization of the Weibull distribution called "generalized odd beta prime-Weibull" (GOBPW) is proposed using the pioneered GOBP-G family. The mixture representations of the new distribution are defined and derived. Some formal statistical properties of the GOBPW distribution, such as the moments, moment generating function, incomplete moments, information generating function, entropies, stress-strength function, quantile function, and order statistics, are derived. The estimation of the parameters of the proposed distribution is evaluated using the maximum likelihood estimation approach. Different cancer disease data sets, such as the bladder, head and neck, acute bone, and blood cancers, are used to illustrate the applicability and usefulness of the new model and were compared using several statistical accuracy measures with that of well-established extended Weibull distributions, which are the beta modified Weibull distribution, Kumaraswamy modified Weibull distribution, gamma generalized modified Weibull distribution, gamma log-logistic Weibull distribution, and beta log-logistic Weibull distribution. The results show that the proposed model gives better results than the competitive models. This study could guide the relevant stakeholders in choosing a suitable statistical model for the health data instead of relying on traditional models to enhance decision-making.

Background. COVID-19 efforts were often ineffective in controlling the spread of the pandemic. Identifying ineffective controls during a pandemic is thus vital. Method. Utilizing publicly available data on COVID deaths in the counties of US states, we create an index to capture and interpret ineffectiveness in the efforts to reduce the spread of the pandemic in US counties. This index is based on the Intervened Poisson Distribution (IPD) introduced originally by Shanmugam. Motivation for the research idea occurred while we noticed the data dispersion of the COVID deaths is smaller than the average only in some counties. Under-dispersed data is common in statistical modeling. A novel approach we adapted in this article includes the estimation of an intervention parameter estimated through iterative non-linear optimization. Results. Twenty-five counties in California, Idaho, Minnesota, Mississippi, Montana, Nebraska, North Carolina, North Dakota, Texas, and Utah were found to be ineffective in controlling for fatalities based on the expected probability distribution. A review of the policies enacted in these areas would provide insight into ineffective prevention efforts, and some of these issues are documented in current literature. Conclusion. The IPD index an innovate way to document efficacy of interventions during pandemics. The IPD may identify ineffective efforts prior to statistical models intended to evaluate efficacy of efforts.

Sword-duels are interesting, from a game-theoretical perspective, as they correspond to sudden-death interactions between players (i.e., defeat means removal from the game), and where victory/defeat depend on not just skill, but also luck. Analyzing probabilities of victory streaks, given a certain level of self and others' skill, is thus relevant for any application domain that can be modeled in the same way. This paper takes inspiration from Miyamoto Musashi, famously undefeated for 61 duels, and implements a Markov-chain Monte-Carlo simulation approach to evaluate this scenario. Results suggest that a 61 victory streak can be probabilistically observed when skill level is roughly 6.5 times that of the average duelist, and with 95\% confidence when skill level is roughly 1000 times that of the average duelist. More generally, this paper provides a method for determining chances of victory streaks in game-theoretical sudden-death encounters, when both "skill" and "luck" contribute to the outcome of the encounter. Specific scenarios can be modeled by modifying the utilized Markov chain and adjusting sampled distributions as required.

Water resource management and flood forecasting are crucial societal and financial stakes requiring reliable predictions of flow parameters (depth, velocity), the accuracy of which is often limited by uncertainties in hydrodynamic numerical models. In this study, we assess the effect of two uncertainty sources, namely breach characteristics induced by overtopping and the roughness coefficient, on water elevations and inundation extent. A two-dimensional (2D) hydraulic solver was applied in a Monte Carlo integration framework to a reach of the Loire river (France) including about 300 physical parameters. Inundation hazard maps for different flood scenarios allowed highlighting the impact of the breach development chronology. A special attention was paid to proposing a relevant sensitivity analysis to exam the factors influencing the depth and extent of flooding. The spatial analysis of vulnerability area induced by a levee breach width exhibits that, with increasing the flood discharge, the rise of the parameter influence is accompanied by a more localized spatial effect. This argues for a local analysis to allow a clear understanding of the flood hazard. The physical interpretation, highlighted by a global sensitivity analysis, showed precisely the dependence of the flood simulation on the main factors studied, i.e. the roughness coefficients and the characteristics of the breaches.

Guessing effects frequently occur in testing data in educational or psychological applications. Different item response models have been proposed to handle guessing effects in dichotomous test items. However, it has been pointed out in the literature that the oftenemployed three-parameter logistic model poses implausible assumptions regarding the guessing process. The four-parameter guessing model has been proposed as an alternative to circumvent these conceptual issues. In this article, the four-parameter guessing model is compared with alternative item response models for handling guessing effects through a simulation study and an empirical example. It turned out that model selection for item response models should be rather based on the AIC than the BIC. However, the RMSD item fit statistic used with typical cutoff values was found to be ineffective in detecting misspecified item response models. Furthermore, sufficiently large sample sizes are required for sufficiently precise item parameter estimation. Moreover, it is argued that the criterion of statistical model fit should not be the sole criterion of model choice. The item response model used in operational practice should be valid with respect to the meaning of the ability variable and the underlying model assumptions. In this sense, the four-parameter guessing model could be the model of choice in educational large-scale assessment studies.

Noisy data is always encountered in real applications, such as bioinformatics, neuroimage and remote sensing. Existing methods mainly consider linear or generalized linear errors-in-variables regression, while relatively little attention is paid for the multivariate response case, and how to evaluate the estimation performance under perturbed covariates is still an open question. In this paper, we consider the information-theoretic limitations of estimating a low-rank matrix in the multi-response errors-in-variables regression model. By application of the information theory and statistical techniques on concentration inequalities, the minimax lower bound is provided in terms of the squared Frobenius loss, which recaptures the rate provided under the clean covariate assumption in previous literatures. Hence our result further indicates that though under the more realistic errors-in-variables situation, no more samples are required so as to achieve a rate-optimal estimation.

Directed acyclic graphs (DAGs) are nonparametric causal path diagrams that have substantial utility as principled representations of disease and healthcare pathways, and of the underlying ‘data generating mechanisms’ these pathways involve. As such, DAGs provide a valuable bridge between: the aetiological knowledge, operational insight and professional experience on which clinical training and practice depend; and the more abstract epistemological and analytical considerations required to extract robust statistical insight from health and healthcare data. DAGs are nonetheless vulnerable to imperfect biomedical paradigms, partial clinical knowledge and limited empirical data. DAGs drawn under such circumstances offer limited scope for statistical insight free from cognitive, analytical or inferential bias if: they misrepresent the data generating mechanisms involved; or ignore the important role that omitted variables (whether measured, unmeasured or unacknowledged) might play therein. To address these weaknesses and broaden the appeal and application of DAGs, this chapter provides ten simple steps that educators can use to improve the analytical competence and statistical confidence of the healthcare students, qualified practitioners and experienced researchers they support. These steps use temporal logic to draw DAGs so as to: reduce reliance on uncertain knowledge, incomplete information, flawed assumptions or guesswork; and avoid, mitigate or acknowledge the errors and biases that each of these incur. The chapter comprises an accessible, non-technical overview of the perspective and thoughtfulness required to generate temporally coherent DAGs as objective representations of the probabilistic causal paths involved in context-specific data generating mechanisms. It encourages a focus on those variables operating as potential sources of analytical or inferential bias when estimating the plausible, probabilistic causal relationship between two pre-specified variables; and specifically addresses the challenges posed by: omitted; time-variant; non-asynchronous; and temporally obscure variables. The chapter includes a worked example based on a published clinical study to demonstrate how each of the steps required to generate temporally-informed DAGs can be applied to: critically appraise the analytical decisions made during applied healthcare research; and inform the decisions required when designing, undertaking and analysing primary and secondary, prospective and retrospective research. The appendices include a summary of ten recommendations for improving the reporting and interrogability of DAGs and DAG-informed analyses.

For more than half a century, Manfred Deistler has been contributing to the construction of the rigorous theoretical foundations of the statistical analysis of time series and more general stochastic processes. Half a century of unremitting activity is not easily summarized in a few pages. In this short note, we chose to concentrate on a relatively little-known aspect of Manfred's contribution which nevertheless had quite an impact on the development of one of the most powerful tools of contemporary time series and econometrics: dynamic factor models.

Location-and-scale transformation of a random variable underpins normal distribution, but it is however fundamentally incorrect for scale estimation such as relative dispersion. In this paper, a parametrized alternative to a normal distribution, called scaloc-normal distribution, is proposed that efficiently works and is fundamentally correct with absolute and relative dispersion estimators. The Monte Carlos simulation experiment was used to generate a total of 600,000 artificial datasets in 600 different simulation scenarios from loc-normal (normal) and scaloc-normal distributions. The absolute and relative dispersion were estimated and compared from the two distributions. The results show that scaloc-normal distribution is a good parametrized alternative to loc-normal distribution, fundamentally correct and efficient with both standard deviation and coefficient of variation. The key statistical advancement from loc-normal to scaloc-normal distribution is its fundamental correctness (i.e., scale-invariant property) with an efficient relative estimator of dispersion (i.e., coefficient of variation). Parametrically, the loc-normal and scaloc-normal distributions are very different, but both have linear transformations.

Numerous approaches address over-fitting in neural networks: by imposing a penalty on the parameters of the network (L1, L2, etc); by changing the network stochastically (drop-out, Gaussian noise, etc.); or by transforming the input data (batch normalization, etc.). In contrast, we aim to ensure that a minimum amount of supporting evidence is present when fitting the model parameters to the training data. This, at the single neuron level, is equivalent to ensuring that both sides of the separating hyperplane (for a standard artificial neuron) have a minimum number of data points — noting that these points need not belong to the same class for the inner layers. We firstly benchmark the results of this approach on the standard Fashion-MINST dataset, comparing it to various regularization techniques. Interestingly, we note that by nudging each neuron to divide, at least in part, its input data, the resulting networks make use of each neuron, avoiding a hyperplane completely on one side of its input data (which is equivalent to a constant into the next layers). To illustrate this point, we study the prevalence of saturated nodes throughout training, showing that neurons are activated more frequently and earlier in training when using this regularization approach. A direct consequence of the improved neuron activation is that deep networks are now easier to train. This is crucially important when the network topology is not known a priori and fitting often remains stuck in a suboptimal local minima. We demonstrate this property by training a network of increasing depth (and constant width): most regularization approaches will result in increasingly frequent training failures (over different random seeds) whilst the proposed evidence-based regularization significantly outperforms in its ability to train deep networks.

In many Bayesian computations, first, we obtain the expression of the joint distribution of all the unknown variables given the observed data. In general, this expression is not separable in those variables. Thus, obtaining their marginals for each variable and computing the expectations are difficult and costly. This problem becomes even more difficult in high dimensional quandaries, which is an important issue in inverse problems. We may then try to propose a surrogate expression with which we can do approximate computations. Often a separable expression approximation can be useful enough. The Variational Bayesian Approximation (VBA) is a technique that approximates the joint distribution $p$ with an easier, for example separable, one $q$ by minimizing Kullback–Leibler Divergence $KL(q|p)$. When $q$ is separable in all the variables, the approximation is also called Mean Field Approximation (MFA) and so $q$ is the product of the approximated marginals. A first standard and general algorithm is alternate optimization of $KL(q|p)$ with respect to $q_i$. A second general approach is its optimization in the Riemannian manifold. However, in this paper, for practical reasons, we consider the case where $p$ is in the exponential family and so is $q$. For this case, $KL(q|p)$ becomes a function of the parameters $\thetab$ of the exponential family. Then, we can use any other optimization algorithm to obtain those parameters. In this paper, we compare three optimization algorithms: standard alternate optimization, a gradient-based algorithm and a natural gradient algorithm and study their relative performances on three examples.

Measurement uncertainty is typically expressed in terms of a symmetric interval , where denotes the measurement result and the expanded uncertainty. However, in the case of heteroscedasticity, symmetric uncertainty intervals can be misleading. In this paper, a different approach for the calculation of uncertainty intervals is introduced. This approach is applicable when a validation study has been conducted with samples with known concentrations. It will be shown how, under certain circumstances, asymmetric uncertainty intervals arise quite naturally and lead to more reliable uncertainty intervals.

This paper focuses on the so-called proportional intensity-based software reliability models (PI-SRMs), which are extensions of the common non homogeneous Poisson process (NHPP)-based SRMs, and describe the probabilistic behavior of software fault-detection process by incorporating the time-dependent software metrics data observed in the development process. Especially we generalize the seminal PI-SRM in Rinsaka, Shibata and Dohi (2006) by introducing eleven well-known fault-detection time distributions, and investigate their goodness-of-fit and predictive performances. In numerical illustrations with four data sets collected in real software development projects, we utilize the maximum likelihood estimation to estimate model parameters with three time-dependent covariates; test execution time, failure identification work and computer time-failure identification, and examine the performances of our PI SRMs in comparison with the existing NHPP-based SRMs without covariates. It is shown that our PI-STMs could give better goodness-of-fit and predictive performances in many cases.

Developing new ways to estimate probabilities can be valuable for science, statistics, and engineering. By considering the information content of different output patterns, recent work invoking algorithmic information theory has shown that a priori probability predictions based on pattern complexities can be made in a broad class of input-output maps. These algorithmic probability predictions do not depend on a detailed knowledge of how output patterns were produced, or historical statistical data. Although quantitatively fairly accurate, a main weakness of these predictions is that they are given as an upper bound on the probability of a pattern, but many low complexity, low probability patterns occur, for which the upper bound has little predictive value. Here we study this low complexity, low probability phenomenon by looking at example maps, namely a finite state transducer, natural time series data, RNA molecule structures, and polynomial curves. Some mechanisms causing low complexity, low probability behaviour are identified, and we argue this behaviour should be assumed as a default in the real world algorithmic probability studies. Additionally, we examine some applications of algorithmic probability and discuss some implications of low complexity, low probability patterns for several research areas including simplicity in physics and biology, a priori probability predictions, Solomonoff induction and Occam's razor, machine learning, and password guessing.

The twisted Wang transform distribution family, defined as the composition of parameter shifted inverse CDF function with an original CDF function, is found to be most suitable for matching low shape factor distributions, characterizing hard to fit or to simulate reinsurance portfolio losses for some perils from our previous study. Among them, the best form for matching a hard-to-fit empirical loss distribution for a specific peril, is the Exponential Fractional Extra Power 0 Distribution in (0,1) with CDF:.The simplest yet still a good form of this family is the Transformed Hyperbolic Tangent Distribution with CDF:,which has analytical formulas for the moments. The twisted Wang transform distribution family is compared and confirmed to be superior to all other well-known distribution families through extensive numerical optimization practice, distribution forms guesses, and computer-aided exploration experiments.

For the symmetrical mixture model and mixture test area, the lattice point set is used to partition, and then the corresponding test statistics can be constructed. In this paper, we first proposes the partition methods under the lattice point sets and obtains several sub-simplexes without common interior points. Furthermore, we present the method for constructing a uniform design on the simplex using the center points of these sub-simplexes. The designs satisfy the uniformity of maximum distance deviation and provide good results for the mean square error deviation. Finally, the uniformity test on the mixture region is considered and illustrated by examples.

Two extensions of the generalized double Lomax distribution are introduced - the spliced-scale generalized double Lomax distribution and the exponentiated generalized double Lomax distribution. Their properties are studied. Usefulness of the new distributions is shown by fitting them to stock returns datasets.

Considering that a fair coin has two sides and a cylindrical edge, the probability that it would fall on its edge is calculated, yielding the probability of heads or tails less than 50%. Theoretical models for a static case and for five dynamic cases, without and with rebounds, show that there is a small probability that the coin does not fall on its head or tail, depending on initial toss conditions, the coin geometry and conditions of the coin and landing surfaces. For the dynamic model with rebounds, it is found that the probability that a 50 Eurocent coin thrown from a normal height with common initial velocity conditions and appropriate surface conditions will end up on its edge is in the order of one against several thousand.

In this note, we provide some effective treatments of a general linear model with adding-up restrictions via algebraic operations of given vectors and matrices in the model, including analytic expressions of the well-known ordinary least-squares estimator (OLSE) and the best linear unbiased estimator (BLUE) of the unknown parameters in the model.

As its capital, Jakarta plays a critical role in boosting Indonesia’s economic growth and setting the precedent for broader change outside of the city. One crucial avenue of inquiry to better understand, and prepare for, the future of a country so heavily impacted by disastrous weather events is understanding the effects of climate change through data. This study investigates meteorological data collected from 1996 to 2021 and compares the application of the SARIMA and the Holt-Winters methods to predict the future influence of climatic parameters on Jakarta’s weather. The performance of the SARIMA method is proven to provide better results than the Holt-Winter models and both methods showed the best performances when forecasting the humidity data. The results of the forecast are able to demonstrate the characteristic of the climate in Jakarta, with dry season ranging from May to October and wet season ranging from November to April.

Automated software bots infiltrate online surveys and corrupt data integrity, not to mention waste researcher time and budgets. Although resources exist to help keep bots out and identify bots when they do evade survey barriers, bot attacks may be a persistent problem for online surveys for a long time to come. Bots are evolving -- even as survey designers try ever more sophisticated methods to fend them off and weed their answers out. Vigilance needs to be high and the bot generators should not be under-estimated. We recount here some bot features we encountered after our own survey was attacked that helped to identify them, and that have not been detailed elsewhere. We also discuss reasons why commonly recommended strategies for how to keep bots out may not be feasible for many scientific researchers.

Objectives: This paper presents a statistical review of modelling simulations for frequency and sensitivity of COVID-19 testing paradigms. Methods: We performed a review of preprints and published articles on PubMed from January 1, 2020 – March 1, 2021 using the search terms “COVID screening testing”, “COVID testing frequency”, “COVID testing frequency screening” and “SARS-CoV-2 testing frequency”.Results: Several authors’ conclusions support the claim that test frequency and test sensitivity both play a role in reducing SARS-CoV-2 transmission. We highlight the interplay between frequency of testing, test sensitivity and the speed at which test results are available in our review. Conclusions: Evidence suggests that sensitivity and frequency of testing both play a part in decreasing transmission of disease. We conclude that, overall, test sensitivity plays less of a role in reducing disease transmission in a population compared to the frequency of testing and how quickly test results are available.

A paradigm shift away from null hypothesis significance testing seems in progress. Based on simulations, we illustrate some of the underlying motivations. First, P-values vary strongly from study to study, hence dichotomous inference using significance thresholds is usually unjustified. Second, statistically significant results have overestimated effect sizes, a bias declining with increasing statistical power. Third, statistically non-significant results have underestimated effect sizes, and this bias gets stronger with higher statistical power. Fourth, the tested statistical hypotheses usually lack biological justification and are often uninformative. Despite these problems, a screen of 48 papers from the 2020 volume of the Journal of Evolutionary Biology exemplifies that significance testing is still used almost universally in evolutionary biology. All screened studies tested default null hypotheses of zero effect with the default significance threshold of p = 0.05, none presented a pre-specified alternative hypothesis, pre-study power calculation, and the probability of ‘false negatives’(beta error). The results sections of the papers presented 49 significance tests on average (median 23, range 0–390). Of 41 studies that contained verbal descriptions of a ‘statistically non-significant’ result, 26 (63%) falsely claimed the absence of an effect. We conclude that studies in ecology and evolutionary biology are mostly exploratory and descriptive. We should thus shift from claiming to ‘test’ specific hypotheses statistically to describing and discussing many hypotheses (possible true effect sizes) that are most compatible with our data, given our statistical model. We already have the means for doing so, because we routinely present compatibility (‘confidence’) intervals covering these hypotheses.

This paper explores the tuning and results of two-part models on rich datasets provided through the Casualty Actuarial Society (CAS). These data sets include BI (bodily injury), PD (property damage) and COLL (collision) coverage, each documenting policy characteristics and claims across a four year period. The datasets are explored, including summaries of all variables, then the methods for modeling are set forth. Models are tuned and the tuning results are displayed, after which we train the final models and seek to explain select predictions. All of the code will be made available on GitHub. Data was provided by a private insurance carrier to the CAS after anonymizing the data set. This data is available to actuarial researchers for well-defined research projects that have universal benefit to the insurance industry and the public. Our hope is that the methods demonstrated here can be a good foundation for future ratemaking models to be developed and tested more efficiently.

Stage of cancer is a discrete ordinal response that indicates aggressiveness of disease and is often used by physicians to determine the type and intensity of treatment to be administered. For example, the FIGO stage in cervical cancer is based on the size and depth of the tumor as well as the level of spread. It may be of clinical relevance to identify molecular features from high-throughput genomic assays that are associated with stage of cervical cancer, to elucidate pathways related to tumor aggressiveness, identify improved molecular features that may be useful for staging, and identify therapeutic targets. High-throughput RNA-Seq data and corresponding clinical data (including stage) for cervical cancer patients has been made available through The Cancer Genome Atlas Project (TCGA). We recently described penalized Bayesian ordinal response models that can be used for variable selection for over-parameterized datasets such as the TCGA-CESC dataset. Herein, we describe our ordinalbayes R package, available from the Comprehensive R Archive Network (CRAN), which is capable of fitting cumulative logit models when the outcome is ordinal and the number of predictors exceeds the sample size, P>N, such as for TCGA data. We demonstrate use of this package through application to TCGA cervical cancer dataset. Our ordinalbayes package can be used to fit models to high-dimensional dataset and effectively performs variable selection.

Regression models are mostly used in all fields of sciences for modelling the relationship between a dependent variable and independent variable(s). The least square method is often used to estimate the parameters in a linear model because it is the best linear unbiased estimator. These estimates can only be reliable if the assumption of normality is satisfied. In some cases the dependent variable might be bimodal and shows a non-linear relationship with the independent variable(s). In this case, a non-linear model should be used. In non-linear model, the standard errors are often obtained by linearizing the nonlinear function around the parameter, assuming central limit theorem. After the linearization, the least square parameter estimates are obtained. It should be noted that the error of the non-linear model is different from that of the transformed linear model. Thus, there is a need to transform back to the original non-linear model. In this note, a novel non-linear function was developed into a non-linear regression model, called Normal-Power model. The least square method was used to estimate the parameter of the transformed model. Its usefulness in regression model was demonstrated using real data of Nigeria Economy-Tourism model.

In this work, we introduce a generalized measure of cumulative residual entropy and study its properties. We show that several existing measures of entropy such as cumulative residual entropy, weighted cumulative residual entropy and cumulative residual Tsallis entropy, are all special cases of the generalized cumulative residual entropy. We also propose a measure of generalized cumulative entropy, which includes cumulative entropy, weighted cumulative entropy and cumulative Tsallis entropy as special cases. We discuss generating function approach using which we derive different entropy measures. Finally, using the newly introduced entropy measures, we establish some relationships between entropy and extropy measures.

Fitting parameters on spline curves produces more parsimonious models, maintaining fit quality. Smoothing the splines reduces predictive variance. Individual splines are often fit by type of variable, e.g., in age-period-cohort models. Linear and cubic splines are most common. Several smoothing criteria have been used for parameter curves, with cubic splines fit by constraining the integral of splines’ squared-second derivatives popular recently. Constraining the sum of squared second differences for linear splines is analogous. Generally the degree of smoothing is selected using cross-validation. Known spline dummy-variable matrices allow regression estimation of splines, with smoothing done via constrained regression. Smoothing criteria based on sums of squares or absolute values of parameters, as in ridge regression or LASSO, improves predictive accuracy and produces splines similar to smoothing by second-derivative constraints. Variables with very low t-statistics represent points where curve-shapes barely change. Eliminating those variables leaves knots concentrated where spline shapes change. A Bayesian version of this puts shrinkage priors on spline parameters. This yields realistic joint parameter distributions, avoids problems associated with using cross-validation for parameter estimation, and readily expands to non-linear modeling, such as interactions among variable types. Regularized regression and Bayesian spline methods are compared for two example datasets.

The present era of the Internet of Things (IoT) having intelligent functionalities in solving problems pertaining to realtime mission-critical systems has brought an immense revolution in diverse fields including healthcare and navigation systems. However, to the best of our knowledge, the potential of IoT has not been fully exploited yet in the field of the energy sector. We argue that there is an immense need to shift the traditional mission-critical electric power system architecture to IoT-based fully orchestrated architecture in order to increase efficiency, as billions of investment is reserved for the energy sector globally. Since network orchestration deals with auomating the interaction between multiple components involved to execute a particular service, therefore, scheduling the relevant processes within strict deadlines becomes the core pillar of the architecture. The mission-critical systems with urgent task execution often suffer from issues of missing task deadlines. In this study, we present a novel IoT task orchestration architecture for efficient energy management of a nanogrid system that focuses on minimizing the use of nonrenewable energy resources and maximizing the use of renewable energy resources. Moreover, major components of IoT task orchestration such as task mapping and task scheduling are also enhanced using NLP and PSO optimization modules. The proposed task scheduling algorithm incorporates the optimized surplus time, and efficiently executes the energy management-related tasks contemplating to their types. The study utilizes sensors to obtain data from physical IoT devices, including photovoltaic (PV), Energy Storage System (ESS), and diesel generator (DG). The performance of the proposed model is evaluated using data set of nanogrid houses. The outcomes revealed that IoT-task orchestration has played a pivotal role in efficient energy management for nanogrid mission-critical system. Furthermore, the comparison with state-of-the-art scheduling algorithms showed that the task starvation rate is reduced to 16% and 12% when compared with RR and FEF algorithms, respectively.

To detect aberrant human behaviors from large volume of passive data collected by smartphones in real time, we propose an online anomaly detection method using Hotelling’s T-squared test. The test statistic is a weighted average, with more weight on the between-individual component when there are little data available for the individual and more weight on the within-individual component when the data are adequate. The algorithm takes only O(1) run time in each update and the required memory usage is fixed after a pre-specified number of updates. The performance of the proposed method, in terms of accuracy, sensitivity and specificity, is consistently better than or equal to the offline method that it builds upon depending on the sample size of the individual data.

Bayesian nonlinear mixed effects models for data in the form of continuous, repeated measurements from a population, also known as Bayesian hierarchical nonlinear models, are a popular platform for analysis when interest focuses on individual specific characteristics and relevant uncertainty quantification. Due to the limitation of computational power, this framework was relatively dormant until the late 1980s, but in recent years, the statistical research community saw vigorous development of new methodological and computational techniques for these models, the emergence of software, and wide application of the models in numerous industrial and academic fields. This article presents an overview of the formulation, interpretation, and implementation of Bayesian nonlinear mixed effects models and surveys recent advances and applications.

The purpose of this study was to identify the factors that affect the mortality among adult HIV/TB co-infected patients and to see the nutritional difference among mortality in residence level. Retrospective cohort studies of 417 patients which fulfill our criteria were included. Multilevel logistic regression models were used. MLwiN and SPSS software are used to estimate the parameter. The variance of the random factor in the empty model was significant which indicates that there were residential differences in TB-HIV co-infected mortality and it shows multilevel analysis was an appropriate approach for further analysis. The prevalence of HIV/TB co-infected patients' death was 12.9% in study time. Functional status, age of patients, WHO clinical stages, nutritional status, CD4 counts, regimen, and BMI were found to be significant determinants of HIV/TB co-infected mortality. In our study, patients with the bedridden category of functional status, the fourth stages of WHO clinical stages (stage IV), patients with higher age, patients whose treatments were second-line regimen and low CD4 cell counts were more at risk of death. The study also revealed that; poor nutritional status increased the risk of mortality among HIV/TB co-infected patients and it varies among the residence of the patients (rural area were more at risk).