Functional box plots satisfy two needs; visualization of functional data, and the calculation of important box plot statistics. Data visualization illuminates key characteristics of functional sets missed by statistical tests and summary statistics. The calculation of box plot statistics for functional sets permits a novel comparison more suited to functional data. The functional box plot uses a depth

(2021). Statement of Retraction: A novel alpha power transformed exponential distribution with real-life applications. Journal of Applied Statistics. Ahead of Print.

ABSTRACT Her Majesty's Revenue and Customs (HMRC) has the ambitious target of making tax digital for all its customers and collecting tax in a more efficient, effective and accurate manner for both the government and UK taxpayers. Self-assessment tax returns, the biggest key business event for HMRC, is also one of the most popular digital services with over 90% of the approximately 12 million taxpayers

This paper proposes a dynamic infectious disease model for COVID-19 daily counts data and estimate the model using the Langevinized EnKF algorithm, which is scalable for large-scale spatio-temporal data, converges to the right filtering distribution, and is thus suitable for performing statistical inference and quantifying uncertainty for the underlying dynamic system. Under the framework of the proposed

In this paper, we propose and explore a novel semiparametric approach to analyzing longitudinal count data. We address the issue of missingness in longitudinal data and propose a weighted generalized estimations equations approach to fitting marginal mean response models for count responses with dropouts. Also, we investigate a spline regression approach to approximating the curvilinear relationship

We demonstrate how a linear factor model with latent variables can be used to estimate correlations between the outcomes of clinical trials. These correlations are needed for many policy questions of drug/vaccine development (such as calculating the optimal size of financial incentives) and the literature so far has relied on expert opinions. We apply our methodology to the case of vaccines and show

In this study, we provide simulation-based exploration and characterization of the two most crucial kernel density functionals that play a central role in kernel density estimation, considering the probability density functions that are members of the location-scale family. Kernel density functional estimates are known to rely on the choice of preliminary bandwidth. Normal-scale estimators are commonly

According to the Atlas of Human Development in Brazil, the income dimension of Municipal Human Development Index (MHDI-I) is an indicator that shows the population's ability in a municipality to ensure a minimum standard of living to provide their basic needs, such as water, food and shelter. In public policy, one of the research objectives is to identify social and economic variables that are associated

The difference in age structure and aging population level was an important factor that caused the difference in COVID-19’s case fatality rate (CFR) in various regions. To eliminate the age effect on estimating the CFR of COVID-19, our study applied nonlinear logistic model and maximum likelihood method to fit the age-fatality curves of COVID-19 in different countries and regions. We further computed

Responses from the paired organs are generally highly correlated in bilateral studies, statistical procedures ignoring the correlation could lead to incorrect results. Note the intraclass correlation in the study of combined unilateral and bilateral outcomes; 11 confidence intervals (CIs) including 7 asymptotic CIs and 4 Bootstrap-resampling CIs for assessing the equivalence of 2 treatments are derived

Cox model and traditional frailty models assume that all individuals will eventually experience the event of interest. This assumption is often overlooked, and situations will arise where it is not realistic. We introduce Compound Poisson frailty model for survival analysis to deal with populations in which some of the individuals will not experience the event of interest. This model assumes that the

The current large amounts of data and advanced technologies have produced new types of complex data, such as histogram-valued data. The paper focuses on classification problems when predictors are observed as or aggregated into histograms. Because conventional classification methods take vectors as input, a natural approach converts histograms into vector-valued data using summary values, such as the

The paper proposes a joint mixture model to model non-ignorable drop-out in longitudinal cohort studies of mental health outcomes. The model combines a (non)-linear growth curve model for the time-dependent outcomes and a discrete-time survival model for the drop-out with random effects shared by the two sub-models. The mixture part of the model takes into account population heterogeneity by accounting

Marketing mix models (MMMs) are statistical models for measuring the effectiveness of various marketing activities such as promotion, media advertisement, etc. In this research, we propose a comprehensive marketing mix model that captures the hierarchical structure and the carryover, shape and scale effects of certain marketing activities, as well as sign restrictions on certain coefficients that are

In this paper, we consider a two-stage sequential estimation procedure to estimate the parameters of a cumulative exposure model under an accelerated testing scenario. In particular, we focus on a step-stress model where the stress level changes after a pre-specified number of failures occur, which is also random. This is termed as a ‘random stress change time’ in the literature. We further aim to

A new method for the analysis of time to ankylosis complication on a dataset of replanted teeth is proposed. In this context of left-censored, interval-censored and right-censored data, a Cox model with piecewise constant baseline hazard is introduced. Estimation is carried out with the expectation maximisation (EM) algorithm by treating the true event times as unobserved variables. This estimation

The correct evaluation of the reproductive number R for COVID-19 is central in the quantification of the potential scope of the pandemic and the selection of an appropriate course of action. In most models, R is modeled as a constant - effectively averaging out the inherent variability of the transmission process due to varying individual contact rates, population densities, or temporal factors amongst

Clusters form the basis of a number of research study designs including survey and experimental studies. Cluster-based designs can be less costly but also less efficient than individual-based designs due to correlation between individuals within the same cluster. Their design typically relies on ad hoc choices of correlation parameters, and is insensitive to variations in cluster design. This article

The increasing inequality of private income and wealth requires the redistribution of financial resources. Thus, several financial support schemes allocate budget across countries or regions. This work shows how to estimate private wealth at low regional levels by means of a modified Fay–Herriot approach that deals with (a) unit and item non-response, especially with used multiple imputation, (b) the

The extreme value distribution was developed for modeling extreme-order statistics or extreme events. In this study, we discuss the distribution of the largest extreme. The main objective of this paper is to determine the best estimators of the unknown parameters of the extreme value distribution. Thus, both classical and Bayesian methods are used. The classical estimation methods under consideration

Logistic regression is estimated by maximizing the log-likelihood objective function formulated under the assumption of maximizing the overall accuracy. That does not apply to the imbalanced data. The resulting models tend to be biased towards the majority class (i.e. non-event), which can bring great loss in practice. One strategy for mitigating such bias is to penalize the misclassification costs

Discriminative subspace clustering (DSC) can make full use of linear discriminant analysis (LDA) to reduce the dimension of data and achieve effective clustering high-dimension data by clustering low-dimension data in discriminant subspace. However, most existing DSC algorithms do not consider the noise and outliers that may be contained in data sets, and when they are applied to the data sets with

The homogeneously weighted moving average (HWMA) control chart is a new memory-type chart that allocates a specific weight to the current sample and the remaining weight is distributed equally to the previous samples. In this paper, the HWMA control chart is proposed for monitoring count data. This chart is based on the Conway–Maxwell (COM) distribution, which can be used to model under-spread and

The National Heart, Lung and Blood Institute Growth and Health Study (NGHS) is a large longitudinal study of childhood health. A main objective of the study is to estimate the joint distributions of cardiovascular risk outcomes at any two time points conditioning on a large number of covariates. Existing multivariate longitudinal methods are not suitable for outcomes at multiple time points. We present

The Canonical Correlation Analysis (CCA) estimates the correlation between two vector variables by maximizing the correlation of linear combinations of their respective components. Here, the CCA is used to find correlation patterns in the last five successive, per pairs, earthquakes ( M≥4.0) preceding 271 main shocks ( M≥5.5) that occurred in the Greek territory during 1964–2018. The vector variables

The Generalized F-test is derived based on the Generalized P-value Method to test the equality of normally distributed group means under unequal variances. There are two approaches to compute the p-value of the GF test, based on beta and chi-squared random numbers. From prior art in the literature, it appears that the two computation approaches of the Generalized tests are equivalent. In this study

Marshall–Olkin bivariate exponential distribution is used to statistically infer the adaptive type II progressive hybrid censored data under dependent competition risk model. For complex censored data with only partial failure reasons observed, maximum likelihood estimation and approximate confidence interval based on Fisher information are established. At the same time, Bayesian estimation is performed

Quality control is an outstanding area of production management. The effectiveness of applied quality control methods strongly depends on the performance of the measurement system. Many researchers aimed to analyze the effect of measurement errors on conformity or process control and proposed solutions to treat measurement uncertainty. Although both risk-based conformity control and process control

Aim of the contribution is analyzing potential events that may negatively impact individuals, assets, and/or the environment, and making judgments about the perceived personal and social riskiness of Covid-19 compared to other hazards belonging to health (AIDS, cancer, infarction), environmental (climate change), behavioral (serious car accidents), and technological (nuclear weapons) domains. The comparative

This study examines the nonlinear relationship between BMI and earnings for workers in China using Bayesian semiparametric methods. Markov chain Monte Carlo (MCMC) methods are used to obtain the posterior distribution. We stratify the whole sample into four subsamples based on gender and type of residence area. Using longitudinal data from the China Health and Nutrition Survey (CHNS) from 1989 to 2011

Modeling cyber risks has been an important but challenging task in the domain of cyber security, which is mainly caused by the high dimensionality and heavy tails of risk patterns. Those obstacles have hindered the development of statistical modeling of the multivariate cyber risks. In this work, we propose a novel approach for modeling the multivariate cyber risks which relies on the deep learning

In this paper, we propose a Susceptible–Infected–Removal (SIR) model with time fused coefficients. In particular, our proposed model discovers the underlying time homogeneity pattern for the SIR model's transmission rate and removal rate via Bayesian shrinkage priors. MCMC sampling for the proposed method is facilitated by the nimble package in R. Extensive simulation studies are carried out to examine

Methodological development and application of joint models for longitudinal and time-to-event data have mostly coupled a single longitudinal outcome-based linear mixed-effects model with normal distribution and Cox proportional hazards model. In practice, however, (i) profile of subject's longitudinal response may follow a `broken-stick nonlinear' (piecewise) trajectory. Such multiple phases are an

Model-based learning of organism dynamics is challenging, particularly when modeling count correlated data. In this paper, we adapt the multivariate Poisson distribution to model nematode dynamics. This distribution relaxes the mean-equal-variance property of the univariate Poisson distribution and allows recovery of the correlation among nematode genera. An observational dataset with 68 soil samples

This paper aims to discuss the Bayesian estimation approach for the zero-inflated cure class of models, which extends the standard cure model by accommodating zero-inflated data in the survival analysis context. A comprehensive simulation study is carried out to assess the performance of the estimation procedure. A new estimation methodology is illustrated using a real dataset related to women diagnosed

This paper deals with the statistical inference of the unknown parameters of three-parameter exponentiated power Lindley distribution under adaptive progressive type-II censored samples. The maximum likelihood estimator (MLE) cannot be expressed explicitly, hence approximate MLEs are conducted using the Newton–Raphson method. Bayesian estimation is studied and the Markov Chain Monte Carlo method is

In multivariate longitudinal studies, several outcomes are repeatedly measured for each subject over time. The data structure of these studies creates two types of associations which should take into account by the model: association of outcomes at a given time point and association among repeated measurements over time for a specific outcome. In our approach, because of some advantageous arisen from

Aiming at the problem of fault detection in data collection in wireless sensor networks, this paper combines evolutionary computing and machine learning to propose a productive technical solution. We choose the classical particle swarm optimization (PSO) and improve it, including the introduction of a biological population model to control the population size, and the addition of a parallel mechanism

Inflated data and over-dispersion are two common problems when modeling count data with traditional Poisson regression models. In this study, we propose a latent class inflated Poisson (LCIP) regression model to solve the unobserved heterogeneity that leads to inflations and over-dispersion. The performance of the model estimation is evaluated through simulation studies. We illustrate the usefulness

In the censored data exploration, the classical linear regression model which assumes normally distributed random errors is perhaps one of the commonly used frameworks. However, practical studies have often criticized the classical linear regression model because of its sensitivity to departure from the normality and partial nonlinearity. This paper proposes to solve these potential issues simultaneously

In the statistical literature, several discrete distributions have been developed so far. However, in this progressive technological era, the data generated from different fields is getting complicated day by day, making it difficult to analyze this real data through the various discrete distributions available in the existing literature. In this context, we have proposed a new flexible family of discrete

Modelling is challenging topic and using parametric models is important stage to reach flexible function for modelling. Weibull distribution has shape and scale parameters which play the main role for modelling. Bimodality parameter is added and so bimodal Weibull distribution can capture real data set with bimodality which can be actually combination of two populations. The properties of the proposed

Over the past few months, the outbreak of Coronavirus disease (COVID-19) has been expanding over the world. A reliable and accurate dataset of the cases is vital for scientists to conduct related research and policy-makers to make better decisions. We collect the United States COVID-19 daily reported data from four open sources: the New York Times, the COVID-19 Data Repository by Johns Hopkins University

(2021). Editorial to special issue Frontiers of Data Analysis. Journal of Applied Statistics: Vol. 48, Special Issue: Frontiers of Data Analysis, pp. 1349-1351.

The granting process is based on the probability that the applicant will refund his/her loan given his/her characteristics. This probability, also called score, is learnt based on a dataset in which rejected applicants are excluded. Thus, the population on which the score is used is different from the learning population. Many ``reject inference'' methods try to exploit the data available from the

Scholars have been interested in the politicization of humans rights within the United Nations for some time. However, previous research typically looks at simple associations between voting coalitions and observable variables, such as geographic location or membership in international organizations. Our study is the first attempt at estimating the latent coalition structure based on the voting data

The relationship between the response variable and one or more independent variables refers to the quality characteristic in some statistical quality control applications, which is called profile. Most research dealt with the monitoring of profiles in single-stage processes considering a basic assumption of normality. However, some processes are made up of several sub-processes; thus, the effect of

In this article, we first propose the modified Hannan–Rissanen Method for estimating the parameters of autoregressive moving average (ARMA) process with symmetric stable noise and symmetric stable generalized autoregressive conditional heteroskedastic (GARCH) noise. Next, we propose the modified empirical characteristic function method for the estimation of GARCH parameters with symmetric stable noise

This paper presents methods of estimation of the parameters and acceleration factor for Nadarajah–Haghighi distribution based on constant-stress partially accelerated life tests. Based on progressive Type-II censoring, Maximum likelihood and Bayes estimates of the model parameters and acceleration factor are established, respectively. In addition, approximate confidence interval are constructed via

It will be recalled that the classical bivariate normal distributions have normal marginals and normal conditionals. It is natural to ask whether a similar phenomenon can be encountered involving Poisson marginals and conditionals. However, it is known, from research on conditionally specified models, that Poisson marginals will be encountered, together with both conditionals being of the Poisson form

Sporting careers observed over a preset time interval can be partitioned into two distinct subsamples. These samples consist of individuals whose careers had already commenced at the start of the time interval (prevalent subsample) and individuals whose careers began during the time interval (incident subsample) as well the respective individual-level covariate data such as salary, height, weight,

In parametric estimates, the maximum likelihood estimation method is the most popular method widely used in the social sciences and psychology, although it is biased in situations where sample sizes are small or the data are heavily censored. Therefore, the main objective of this research is to improve this estimation method using the Runge–Kutta technique. The improved method was applied to derive

This paper revitalizes the investigation of the classical cusp catastrophe model in catastrophe theory and tackles the unsolved statistical inference problem concerning stochastic cusp differential equation. This model is challenging because its associated transition density hence the likelihood function is analytically intractable. We propose a novel Bayesian approach combining Hamiltonian Monte Carlo

When the observed proportion of zeros in a data set consisting of binary outcome data is larger than expected under a regular logistic regression model, it is frequently suggested to use a zero-inflated Bernoulli (ZIB) regression model. A spline-based ZIB regression model is proposed to describe the potentially nonlinear effect of a continuous covariate. A spline is used to approximate the unknown

Parameters of a linear regression model can be estimated with the help of traditional methods like generalized least squares and mixed estimator. However, recent developments increased the importance of big data sets, which have much more predictors than observations where some predictors have no impact on the dependent variable. The estimation and model selection problem of big datasets can be solved

To correctly measure the effect of mortality rates on the stability of insurance and pension provider's financial risk, longevity risk should be considered. This paper aims to investigate the future mortality and longevity risk with different age structures for different countries. Lee–Carter mortality model is used on the historical census data to forecast future mortality rates. Turkey, Germany,

The widely used Fellegi–Sunter model for probabilistic record linkage does not leverage information contained in field values and consequently leads to identical classification of match status regardless of whether records agree on rare or common values. Since agreement on rare values is less likely to occur by chance than agreement on common values, records agreeing on rare values are more likely

The work proposes a new family of survival models called the Odd log-logistic generalized Neyman type A long-term. We consider different activation schemes in which the number of factors M has the Neyman type A distribution and the time of occurrence of an event follows the odd log-logistic generalized family. The parameters are estimated by the classical and Bayesian methods. We investigate the mean

The kth ( 1

Unsupervised document classification for imbalanced data sets poses a major challenge. To obtain accurate classification results, training data sets are often created manually by humans which requires expert knowledge, time and money. Depending on the imbalance of the data set, this approach also either requires human labelling of all of the data or it fails to adequately recognize underrepresented