Search

This paper consists of two parts. The first part of the paper is to propose an explicit robust estimation method for the regression coefficients in simple linear regression based on the power-weighted repeated medians technique that has a tuning constant for dealing with the trade-offs between efficiency and robustness. We then investigate the lower and upper bounds of the finite-sample breakdown point of the proposed method. The second part of the paper is to show that based on the linearization of the cumulative distribution function, the proposed method can be applied to obtain robust parameter estimators for the Weibull and Birnbaum-Saunders distributions that are commonly used in both reliability and survival analysis. Numerical studies demonstrate that the proposed method performs well in a manner that is approximately comparable with the ordinary least squares method, whereas it is far superior in the presence of data contamination that occurs frequently in practice. [ABSTRACT FROM AUTHOR]

In this paper we introduce a novel statistical framework based on the first two quantile conditional moments that facilitates effective goodness-of-fit testing for one-sided Lévy distributions. The scale-ratio framework introduced in this paper extends our previous results in which we have shown how to extract unique distribution features using conditional variance ratio for the generic class of α-stable distributions. We show that the conditional moment-based goodness-of-fit statistics are a good alternative to other methods introduced in the literature tailored to the one-sided Lévy distributions. The usefulness of our approach is verified using an empirical test power study. For completeness, we also derive the asymptotic distributions of the test statistics and show how to apply our framework to real data. [ABSTRACT FROM AUTHOR]

In medical diagnostic research, it is customary to collect multiple continuous biomarker measures to improve the accuracy of diagnostic tests. A prevalent practice is to combine the measurements of these biomarkers into one single composite score. However, incorporating those biomarker measurements into a single score depends on the combination of methods and may lose vital information needed to make an effective and accurate decision. Furthermore, a diagnostic cut-off is required for such a combined score, and it is difficult to interpret in actual clinical practice. The paper extends the classical biomarkers' accuracy and predictive values from univariate to bivariate markers. Also, we will develop a novel pseudo-measures system to maximize the vital information from multiple biomarkers. We specified these pseudo-and-or classifiers for the true positive rate, true negative rate, false-positive rate, and false-negative rate. We used them to redefine classical measures such as the Youden index, diagnostics odds ratio, likelihood ratios, and predictive values. We provide optimal cut-off point selection based on the modified Youden index with numerical illustrations and real data analysis for this paper's newly developed pseudo measures. [ABSTRACT FROM AUTHOR]

The histogram has all its bin widths equal to some non-random number arbitrary set by an analyst (EBWH). In the result, particular bin counts are random variables. This paper presents also a histogram that is constructed in a converse manner. Bin counts are all equal to some non-random number arbitrary set by an analyst (EBCH). In the result, particular bin widths are random variables. The first goal of the paper is a choose of constant bin width (of bin numbers k) in the EBWH, which maximize the similarity measure in the Monte Carlo simulation. The second goal is a choose of constant bin count in the EBCH, which maximize the similarity measure in the Monte Carlo simulation. The third goal is to present similarity measures between empirical and theoretical data. The fourth goal is the comparative analysis of two histogram methods by means of the frequency formula. The first additional goal is a tip how to proceed in EBCH when modulo(n,k)≠0. The second additional goal is the software in the form of a Mathcad file with the implementation of EBWH and EBCH. [ABSTRACT FROM AUTHOR]

In this paper, reparameterization and student-t are applied to Stochastic Volatility (SV) model. We aim to reduce the amount of autocorrelation of the SV parameters and to introduce heavy-tailed model via the Bayesian computation of the Markov Chain Monte Carlo (MCMC) samplers. This research paper helps support better MCMC estimation of the SV model for volatile Asian FX series during Covid-19. [ABSTRACT FROM AUTHOR]

Coronavirus disease 2019 (COVID-19) caused by the SARS-CoV-2 virus has spread seriously throughout the world. Predicting the spread, or the number of cases, in the future can facilitate preparation for, and prevention of, a worst-case scenario. To achieve these purposes, statistical modeling using past data is one feasible approach. This paper describes spatio-temporal modeling of COVID-19 case counts in 47 prefectures of Japan using a nonlinear random effects model, where random effects are introduced to capture the heterogeneity of a number of model parameters associated with the prefectures. The negative binomial distribution is frequently used with the Paul-Held random effects model to account for overdispersion in count data; however, the negative binomial distribution is known to be incapable of accommodating extreme observations such as those found in the COVID-19 case count data. We therefore propose use of the beta-negative binomial distribution with the Paul-Held model. This distribution is a generalization of the negative binomial distribution that has attracted much attention in recent years because it can model extreme observations with analytical tractability. The proposed beta-negative binomial model was applied to multivariate count time series data of COVID-19 cases in the 47 prefectures of Japan. Evaluation by one-step-ahead prediction showed that the proposed model can accommodate extreme observations without sacrificing predictive performance. [ABSTRACT FROM AUTHOR]

Orthogonality and optimality of run order are two important and worthy to be considered criteria in design of experiment. For the mixed-level designs commonly used in the case of unequal levels of factors in the experiment, both orthogonality and optimality of run order are taken into account at the same time in this paper. The construction methods of mixed-level designs with orthogonality and relatively minimum or maximum level changes are respectively proposed. Based on such designs, the orthogonal main-effect plans with relatively minimum or maximum level changes are obtained by the method of collapsing factor. The designs with relatively minimum or maximum level changes constructed in this paper guarantee the cost minimum or the time trend most robust in the actual experiment arrangement. [ABSTRACT FROM AUTHOR]

The paper considers parameter estimation in count data models using penalized likelihood methods. The motivating data consists of multiple independent count variables with a moderate sample size per variable. The data were collected during the assessment of oral reading fluency (ORF) in school-aged children. A sample of fourth-grade students were given one of ten available passages to read with these differing in length and difficulty. The observed number of words read incorrectly (WRI) is used to measure ORF. Three models are considered for WRI scores, namely the binomial, the zero-inflated binomial, and the beta-binomial. We aim to efficiently estimate passage difficulty, a quantity expressed as a function of the underlying model parameters. Two types of penalty functions are considered for penalized likelihood with respective goals of shrinking parameter estimates closer to zero or closer to one another. A simulation study evaluates the efficacy of the shrinkage estimates using Mean Square Error (MSE) as metric. Big reductions in MSE relative to unpenalized maximum likelihood are observed. The paper concludes with an analysis of the motivating ORF data. [ABSTRACT FROM AUTHOR]

In this paper, we consider the problem of classification of misaligned multivariate functional data. We propose to use a model-based approach for the joint registration and classification of such data. The observed functional inputs are modeled as a functional nonlinear mixed effects model containing a nonlinear functional fixed effect constructed upon warping functions to account for curve alignment, and a nonlinear functional random effects component to address the variability among subjects. The warping functions are also modeled to accommodate common effect within groups and the variability between subjects. Then, a functional logistic regression model defined upon the representation of the aligned curves and scalar inputs is used to account for curve classification. EM-based algorithms are developed to perform maximum likelihood inference of the proposed models. The identifiability of the registration model and the asymptotical properties of the proposed method are established. The performance of the proposed procedure is illustrated via simulation studies and an analysis of a hyoid bone movement data application. The statistical developments proposed in this paper were motivated by the hyoid bone movement study, the methodology is designed and presented generality and can be applied to numerous areas of scientific research. [ABSTRACT FROM AUTHOR]

When the data are stored in a distributed manner, direct application of traditional hypothesis testing procedures is often prohibitive due to communication costs and privacy concerns. This paper mainly develops and investigates a distributed two-node Kolmogorov–Smirnov hypothesis testing scheme, implemented by the divide-and-conquer strategy. In addition, this paper also provides a distributed fraud detection and a distribution-based classification for multi-node machines based on the proposed hypothesis testing scheme. The distributed fraud detection is to detect which node stores fraud data in multi-node machines and the distribution-based classification is to determine whether the multi-node distributions differ and classify different distributions. These methods can improve the accuracy of statistical inference in a distributed storage architecture. Furthermore, this paper verifies the feasibility of the proposed methods by simulation and real example studies. [ABSTRACT FROM AUTHOR]

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 for further tuning. The proposed RS-PPSO algorithm was successfully used to optimize the initial weights and biases of back propagation neural network (BPNN), shortening the training time and raising the prediction accuracy. Wireless sensor networks (WSN) has become the key supporting platform of Internet of Things (IoT). The correctness of the data collected by the sensor nodes has a great influence on the reliability, real-time performance and energy saving of the entire network. The optimized machine learning technology scheme given in this paper can effectively identify the fault data, so as to ensure the effective operation of WSN. [ABSTRACT FROM AUTHOR]

In many situations, it is common to have more than one observation per experimental unit, thus generating the experiments with repeated measures. In the modeling of such experiments, it is necessary to consider and model the intra-unit dependency structure. In the literature, there are several proposals to model positive continuous data with repeated measures. In this paper, we propose one more with the generalization of the beta prime regression model. We consider the possibility of dependence between observations of the same unit. Residuals and diagnostic tools also are discussed. To evaluate the finite-sample performance of the estimators, using different correlation matrices and distributions, we conducted a Monte Carlo simulation study. The methodology proposed is illustrated with an analysis of a real data set. Finally, we create an R package for easy access to publicly available the methodology described in this paper. [ABSTRACT FROM AUTHOR]

Climate change has become increasingly important in recent years. It is the outcome of the burning of fossil fuels that increased the concentration of atmospheric carbon dioxide (CO 2 ), over the last century. Mitigating the impacts of climate change requires a better understanding and assessment of the countries' economic decisions on the amount of CO 2 emissions. This paper assesses the variability between the different countries in the trends of CO 2 emissions and electricity consumption from 1975 to 2014, while identifying clusters of countries of similar trends over time. The novel methodology applied in this paper enables us to assess long-debated issues in climate literature. The temporal dynamic effects of electricity consumption and economic growth on CO 2 emissions across countries are studied using functional data analysis (FDA) methods. The latter have proven to be useful tools for visualising similarities and differences in the non-linear trends of CO 2 emissions without forcing linear trends and stationary relationships which can be unrealistic and misleading. The results indicate the possibility of identifying changes in the trends of CO 2 emissions and electricity consumption for a wide range of heterogeneous countries over the study period. The findings also reveal that economic growth puts a strain on the environment, where many high-income countries are still away from attaining economic-energy sustainability. [ABSTRACT FROM AUTHOR]

Polytomous categorical data are frequent in studies, that can be obtained with an individual or grouped structure. In both structures, the generalized logit model is commonly used to relate the covariates on the response variable. After fitting a model, one of the challenges is the definition of an appropriate residual and choosing diagnostic techniques. Since the polytomous variable is multivariate, raw, Pearson, or deviance residuals are vectors and their asymptotic distribution is generally unknown, which leads to difficulties in graphical visualization and interpretation. Therefore, the definition of appropriate residuals and the choice of the correct analysis in diagnostic tools is important, especially for nominal data, where a restriction of methods is observed. This paper proposes the use of randomized quantile residuals associated with individual and grouped nominal data, as well as Euclidean and Mahalanobis distance measures, as an alternative to reduce the dimension of the residuals. We developed simulation studies with both data structures associated. The half-normal plots with simulation envelopes were used to assess model performance. These studies demonstrated a good performance of the quantile residuals, and the distance measurements allowed a better interpretation of the graphical techniques. We illustrate the proposed procedures with two applications to real data. [ABSTRACT FROM AUTHOR]

In recent years, some interested data can be recorded only if the values fall within an interval range, and the responses are often subject to censoring. Attempting to perform effective statistical analysis with censored, especially heavy-tailed and asymmetric data, can be difficult. In this paper, we develop a novel linear regression model based on the proposed skewed generalized t distribution for censored data. The likelihood-based inference and diagnostic analysis are established using the Expectation/Conditional Maximization Either algorithm in conjunction with smoothing approximate functions. We derive relevant measures to perform global influence for this novel model and develop local influence analysis based on the conditional expectation of the complete-data log-likelihood function. Some useful perturbation schemes are discussed. We illustrate the finite sample performance and the robustness of the proposed method by simulation studies. The proposed model is compared with other procedures based on a real dataset, and a sensitivity analysis is also conducted. [ABSTRACT FROM AUTHOR]

In various scenarios where products and services are accompanied by warranties to ensure their reliability over a specified time, the two-parameter (shifted) exponential distribution serves as a fundamental model for time-to-event data. In modern production process, the products often come with warranties, and their quality can be manifested by the changes in the scale and origin parameters of a shifted exponential (SE) distribution. This paper introduces the Max-EWMA chart, employing maximum likelihood estimators and exponentially weighted moving average (EWMA) statistics, to jointly monitor SE distribution parameters. Additionally, we extend two additional charts, namely the Max-DEWMA and Max-TEWMA charts to enhance early-stage shift detection. Performance evaluations under zero-state and steady-state conditions compare these charts with the existing Max-CUSUM chart in terms of expected value and standard deviation of the run length (RL) distribution. Our findings reveal that among the Max-EWMA schemes, the Max-EWMA SE chart outperforms the others in terms of steady-state performance, while the Max-TEWMA chart surpasses the Max-EWMA and Max-DEWMA SE charts in respect to zero-state performance. Moreover, the proposed Max-EWMA schemes demonstrate advantages over Max-CUSUM, especially for small to moderate smoothing constants. We also provide an illustrative example to demonstrate the implementation of the proposed schemes. [ABSTRACT FROM AUTHOR]

Proportional data arise frequently in a wide variety of fields of study. Such data often exhibit extra variation such as over/under dispersion, sparseness and zero inflation. For example, the hepatitis data present both sparseness and zero inflation with 19 contributing non-zero denominators of 5 or less and with 36 having zero seropositive out of 83 annual age groups. The whitefly data consists of 640 observations with 339 zeros (53%), which demonstrates extra zero inflation. The catheter management data involve excessive zeros with over 60% zeros averagely for outcomes of 193 urinary tract infections, 194 outcomes of catheter blockages and 193 outcomes of catheter displacements. However, the existing models cannot always address such features appropriately. In this paper, a new two-parameter probability distribution called Lindley–binomial (LB) distribution is proposed to analyze the proportional data with such features. The probabilistic properties of the distribution such as moment, moment generating function are derived. The Fisher scoring algorithm and EM algorithm are presented for the computation of estimates of parameters in the proposed LB regression model. The issues on goodness of fit for the LB model are discussed. A limited simulation study is also performed to evaluate the performance of derived EM algorithms for the estimation of parameters in the model with/without covariates. The proposed model is illustrated through three aforementioned proportional datasets. [ABSTRACT FROM AUTHOR]

This paper presents a novel approach for analyzing bivariate positive data, taking into account a covariate vector and left-censored observations, by introducing a hierarchical Bayesian analysis. The proposed method assumes marginal Weibull distributions and employs either a usual Weibull likelihood or Weibull–Tobit likelihood approaches. A latent variable or frailty is included in the model to capture the possible correlation between the bivariate responses for the same sampling unit. The posterior summaries of interest are obtained through Markov Chain Monte Carlo methods. To demonstrate the effectiveness of the proposed methodology, we apply it to a bivariate data set from stellar astronomy that includes left-censored observations and covariates. Our results indicate that the new bivariate model approach, which incorporates the latent factor to capture the potential dependence between the two responses of interest, produces accurate inference results. We also compare the two models using the different likelihood approaches (Weibull or Weibull–Tobit likelihoods) in the application. Overall, our findings suggest that the proposed hierarchical Bayesian analysis is a promising approach for analyzing bivariate positive data with left-censored observations and covariate information. [ABSTRACT FROM AUTHOR]

It is well-known that the presence of serial correlation may result in an inefficient or even biased estimation in time series analysis. In this paper, we consider testing serial correlation in a general d-factor model when the model errors follow the GARCH process, which is frequently used in modeling financial data. Two empirical likelihood-based testing statistics are suggested as a way to deal with problems that might come up with infinite variance. Both statistics are shown to be chi-squared distributed asymptotically under mild conditions. Simulations confirm the excellent finite-sample performance of both tests. Finally, to emphasize the importance of using our tests, we explore the impact of the exchange rate on the stock return using both monthly and daily data from eight countries. [ABSTRACT FROM AUTHOR]

This paper presents an effort to investigate the estimations of the Weibull distribution using an improved adaptive Type-II progressive censoring scheme. This scheme effectively guarantees that the experimental time will not exceed a pre-fixed time. The point and interval estimations using two classical estimation methods, namely maximum likelihood and maximum product of spacing, are considered to estimate the unknown parameters as well as the reliability and hazard rate functions. The approximate confidence intervals of these quantities are obtained based on the asymptotic normality of the maximum likelihood and maximum product of spacing methods. The Bayesian estimations are also considered using MCMC techniques based on the two classical approaches. An extensive simulation study is implemented to compare the performance of the different methods. Further, we propose the use of various optimality criteria to find the optimal sampling scheme. Finally, one real data set is applied to show how the proposed estimators and the optimality criteria work in real-life scenarios. The numerical outcomes demonstrated that the Bayesian estimates using the likelihood and product of spacing functions performed better than the classical estimates. [ABSTRACT FROM AUTHOR]

This paper aims to detect anomalous changes in social network structure in real time and to offer early warnings by phase II monitoring social networks. First, the exponential random graph model is used to model social networks. Then, a test and online monitoring technique of the exponential random graph model is developed based on the split likelihood-ratio test after determining the model and its parameters for a specific data set. This proposed approach uses pseudo-maximum likelihood estimation and likelihood ratio to construct the test statistics, avoiding the several steps of discovering Monte Carlo Markov Chain maximum likelihood estimation through an iterative method. A bisection algorithm for the control limit is given. Simulations on three data sets Flobusiness, Kapferer and Faux.mesa.high are presented to study the performance of the procedure. Different change points and shift sizes are compared to see how they affect the average run length. A real application example on the MIT reality mining social proximity network is used to illustrate the proposed modelling and online monitoring methods. [ABSTRACT FROM AUTHOR]

Biclustering is widely used in different kinds of fields including gene information analysis, text mining, and recommendation system by effectively discovering the local correlation between samples and features. However, many biclustering algorithms will collapse when facing heavy-tailed data. In this paper, we propose a robust version of convex biclustering algorithm with Huber loss. Yet, the newly introduced robustification parameter brings an extra burden to selecting the optimal parameters. Therefore, we propose a tuning-free method for automatically selecting the optimal robustification parameter with high efficiency. The simulation study demonstrates the more fabulous performance of our proposed method than traditional biclustering methods when encountering heavy-tailed noise. A real-life biomedical application is also presented. The R package RcvxBiclustr is available at https://github.com/YifanChen3/RcvxBiclustr. [ABSTRACT FROM AUTHOR]

In this paper, we consider the estimation of intracluster correlation for ordinal data. We focus on pure-tone audiometry hearing threshold data, where thresholds are measured in 5 decibel increments. We estimate the intracluster correlation for tests from iPhone-based hearing assessment applications as a measure of test/retest reliability. We present a method to estimate the intracluster correlation using mixed effects cumulative logistic and probit models, which assume the outcome data are ordinal. This contrasts with using a mixed effects linear model which assumes that the outcome data are continuous. In simulation studies, we show that using a mixed effects linear model to estimate the intracluster correlation for ordinal data results in a negative finite sample bias, while using mixed effects cumulative logistic or probit models reduces this bias. The estimated intracluster correlation for the iPhone-based hearing assessment application is higher when using the mixed effects cumulative logistic and probit models compared to using a mixed effects linear model. When data are ordinal, using mixed effects cumulative logistic or probit models reduces the bias of intracluster correlation estimates relative to using a mixed effects linear model. [ABSTRACT FROM AUTHOR]

Comparative lifetime experiments are remarkable when the study is to ascertain the relative merits of two competing products regarding the duration of their service life. This paper considers the comparative lifetime experiments of two Gompertz populations under a balanced joint progressive Type-II censoring scheme. The lifetime distributions of the units are assumed to follow the Gompertz distribution with a common shape but different scale parameters. The maximum likelihood estimates of the unknown parameters are derived. The existence of the maximum likelihood estimates is proved. Expectation-maximization and stochastic expectation-maximization algorithms are provided to calculate the estimates. The bootstrap-p, bootstrap-t, and approximate confidence intervals are established. To obtain the Bayesian estimates, it is assumed that the prior of scale parameters is a Beta-Gamma distribution and the prior of the common shape parameter is an independent Gamma distribution. Under squared error loss and LINEX loss functions, the Metropolis-Hastings algorithm is provided to compute the Bayes estimates and the credible intervals. Further, the statistical inferences with order restriction are studied when it is known a priori that the expectation of the lifespan of one population is shorter than that of the other population. A wide range of simulation experiments is conducted to evaluate the performance of the proposed methods. Finally, the lifetimes of white organic light-emitting diodes and the breaking strengths of jute fiber of gauge lengths are analyzed to illustrate the practical application of the proposed model and methods. [ABSTRACT FROM AUTHOR]

The modality is an important topic for modelling. Using parametric models is an efficient way when real data set shows trimodality. In this paper, we propose a new class of trimodal probability distributions, that is, probability distributions that have up to three modes. Trimodality itself is achieved by applying a proper transformation to density function of certain continuous probability distributions. At first, we obtain preliminary results for an arbitrary density function $ g(x) $ g (x) and, next, we focus on the Gaussian case, studying trimodal Gaussian model more deeply. The Gaussian distribution is applied to produce the trimodal form of Gaussian known as normal distribution. The tractability of analytical expression of normal distribution and properties of the trimodal normal distribution are important reasons why we choose normal distribution. Furthermore, the existing distributions should be improved to be capable of modelling efficiently when there exists a trimodal form in a data set. After new density function is proposed, estimating its parameters is important. Since Mathematica 12.0 software has optimization tools and important modelling techniques, computational steps are performed using this software. The bootstrapped form of real data sets are applied to show the modelling ability of the proposed distribution when real data sets show trimodality. [ABSTRACT FROM AUTHOR]

The accelerated hazards model is one of the most commonly used models for regression analysis of failure time data and this is especially the case when, for example, the hazard functions may have monotonicity property. Correspondingly a large literature has been established for its estimation or inference when right-censored data are observed. Although several methods have also been developed for its inference based on interval-censored data, they apply only to limited situations or rely on some assumptions such as independent censoring. In this paper, we consider the situation where one observes case K interval-censored data, the type of failure time data that occur most in, for example, medical research such as clinical trials or periodical follow-up studies. For inference, we propose a sieve borrow-strength method and in particular, it allows for informative censoring. The asymptotic properties of the proposed estimators are established. Simulation studies demonstrate that the proposed inference procedure performs well. The method is applied to a set of real data set arising from an AIDS clinical trial. [ABSTRACT FROM AUTHOR]

The main concern of this paper is providing a flexible discrete model that captures every kind of dispersion (equi-, over- and under-dispersion). Based on the balanced discretization method, a new discrete version of Burr–Hatke distribution is introduced with the partial moment-preserving property. Some statistical properties of the new distribution are introduced, and the applicability of proposed model is evaluated by considering counting series. A new integer-valued autoregressive (INAR) process based on the mixing Pegram and binomial thinning operators with discrete Burr–Hatke innovations is introduced, which can model contagious data properly. The different estimation approaches of parameters of the new process are provided and compared through the Monte Carlo simulation scheme. The performance of the proposed process is evaluated by four data sets of the daily death counts of the COVID-19 in Austria, Switzerland, Nigeria and Slovenia in comparison with some competitor INAR(1) models, along with the Pearson residual analysis of the assessing model. The goodness of fit measures affirm the adequacy of the proposed process in modeling all COVID-19 data sets. The fundamental prediction procedures are considered for new process by classic, modified Sieve bootstrap and Bayesian forecasting methods for all COVID-19 data sets, which is concluded that the Bayesian forecasting approach provides more reliable results. [ABSTRACT FROM AUTHOR]

The Hosmer-Lemeshow (HL) test is a commonly used global goodness-of-fit (GOF) test that assesses the quality of the overall fit of a logistic regression model. In this paper, we give results from simulations showing that the type I error rate (and hence power) of the HL test decreases as model complexity grows, provided that the sample size remains fixed and binary replicates (multiple Bernoulli trials) are present in the data. We demonstrate that a generalized version of the HL test (GHL) presented in previous work can offer some protection against this power loss. These results are also supported by application of both the HL and GHL test to a real-life data set. We conclude with a brief discussion explaining the behavior of the HL test, along with some guidance on how to choose between the two tests. In particular, we suggest the GHL test to be used when there are binary replicates or clusters in the covariate space, provided that the sample size is sufficiently large. [ABSTRACT FROM AUTHOR]

This paper introduces a new family of quantile regression models whose response variable follows a reparameterized Marshall-Olkin distribution indexed by quantile, scale, and asymmetry parameters. The family has arisen by applying the Marshall-Olkin approach to distributions belonging to the location-scale family. Models of higher flexibility and whose structure is similar to generalized linear models were generated by quantile reparameterization. The maximum likelihood (ML) method is presented for the estimation of the model parameters, and simulation studies evaluated the performance of the ML estimators. The advantages of the family are illustrated through an application to a set of nutritional data, whose results indicate it is a good alternative for modeling slightly asymmetric response variables with support on the real line. [ABSTRACT FROM AUTHOR]

Methods for designing and analyzing multiple arms survival trials that incorporate patient's treatment choice are needed. In these trials, patients are randomized into two groups, random and choice. Participants in the choice group choose their treatment, which is not a current standard practice in randomized clinical trials. In this paper, we propose a new method based on the likelihood function to design and analyze these trials with time to event outcomes in the presence of non-informative right censoring. We use simulations to evaluate the methods for Weibull outcomes, complete and censored. Finally, we provide an illustration for designing a study in which we discuss some design considerations and demonstrate the methods. [ABSTRACT FROM AUTHOR]

Ensuring products maintain their normal function during the warranty period is essential for companies to control warranty costs and sustain their brand reputation. A variables acceptance sampling (VAS) plan based on the lifetime performance index (LPI) is a practical technology for verifying product lifetimes because it can provide the number of failures required for life testing and the acceptance criterion for lot disposition. Recently, an LPI-based variables quick-switch sampling (VQSS) system has been developed to improve the sampling efficiency. The VQSS system manipulates the normal and tightened VAS plans for a series of lot dispositions. However, the existing LPI-based VQSS system only alters acceptance criteria to build normal and tightened VAS plans, which may limit its applicability. In this paper, we propose an LPI-based VQSS system with a mechanism that can alter the required number of failures. Compared to the existing LPI-based VQSS system, the proposed method has superior discriminatory power and moderate cost-effectiveness. Finally, we illustrate a lifetime verification of an electric vehicle’s battery as a case study to demonstrate the practicality of the proposed method. [ABSTRACT FROM AUTHOR]

The advance of high-throughput biomedical technology makes it possible to access massive measurements of gene expression levels. An important statistical issue is identifying both under-expressed and over-expressed genes for a disease. Most existing multiple-testing procedures focus on selecting only the non-null or significant genes without further identifying their expression type. Only limited methods are designed for the directional problem, and yet they fail to separately control the numbers of falsely discovered over-expressed and under-expressed genes with only a unified index combining all the false discoveries. In this paper, based on a three-classification multiple testing framework, we propose a practical data-driven procedure to control separately the two directions of false discoveries. The proposed procedure is theoretically valid and optimal in the sense that it maximizes the expected number of true discoveries while controlling the false discovery rates for under-expressed and over-expressed genes simultaneously. The procedure allows different nominal levels for the two directions, exhibiting high flexibility in practice. Extensive numerical results and analysis of two large-scale genomic datasets show the effectiveness of our procedure. [ABSTRACT FROM AUTHOR]

It has long been an open problem to provide a unified test for the intercept of autoregressive (AR) models. In this paper, we use the empirical likelihood method to solve this issue. It turns out that the resulting test statistic always converges in distribution to a standard chi-squared distribution under the null hypothesis, whether the AR process is stationary or nonstationary, and with or without an intercept. The asymptotic distribution under the local alternative hypothesis is also derived under some mild conditions. Several simulations as well as a real data example are used to show how well the suggested test performs in terms of size and power on a finite sample. [ABSTRACT FROM AUTHOR]

Understanding wildfire spread in Canada is critical to promoting forest health and protecting human life and infrastructure. Quantifying fire spread from noisy images, where change-point boundaries separate regions of fire, is critical to accurately estimating fire spread rates. The challenge lies in denoising the fire images and accurately identifying highly non-linear fire lines without smoothing over boundaries. In this paper, we develop an iterative smoothing algorithm for change-point data that utilizes oversmoothed estimates of the underlying data generating process to inform re-smoothing. We demonstrate its effectiveness on simulated one- and two-dimensional change-point data, and robustness to response outliers. Then, we apply the methodology to fire spread images from laboratory micro-fire experiments and show that the regions fuel, burning and burnt-out are smoothed while boundaries are preserved. [ABSTRACT FROM AUTHOR]

The Mann–Whitney–Wilcoxon rank sum test (MWWRST) is a widely used method for comparing two treatment groups in randomized control trials, particularly when dealing with highly skewed data. However, when applied to observational study data, the MWWRST often yields invalid results for causal inference. To address this limitation, Wu et al. (Causal inference for Mann–Whitney–Wilcoxon rank sum and other nonparametric statistics, Stat. Med. 33 (2014), pp. 1261–1271) introduced an approach that incorporates inverse probability weighting (IPW) into this rank-based statistic to mitigate confounding effects. Subsequently, Mao (On causal estimation using U-statistics, Biometrika 105 (2018), pp. 215–220), Zhang et al. (Estimating Mann Whitney-type causal effects, J. Causal Inference 7 (2019), ARTICLE ID 20180010), and Ai et al. (A Mann–Whitney test of distributional effects in a multivalued treatment, J. Stat. Plan. Inference 209 (2020), pp. 85–100) extended this IPW estimator to develop doubly robust estimators. Nevertheless, each of these approaches has notable limitations. Mao's method imposes stringent assumptions that may not align with real-world study data. Zhang et al.'s (Estimating Mann Whitney-type causal effects, J. Causal Inference 7 (2019), ARTICLE ID 20180010) estimators rely on bootstrap inference, which suffers from computational inefficiency and lacks known asymptotic properties. Meanwhile, Ai et al. (A Mann–Whitney test of distributional effects in a multivalued treatment, J. Stat. Plan. Inference 209 (2020), pp. 85–100) primarily focus on testing the null hypothesis of equal distributions between two groups, which is a more stringent assumption that may not be well-suited to the primary practical application of MWWRST. In this paper, we aim to address these limitations by leveraging functional response models (FRM) to develop doubly robust estimators. We demonstrate the performance of our proposed approach using both simulated and real study data. [ABSTRACT FROM AUTHOR]

Survival data analysis often uses the Cox proportional hazards (PH) model. This model is widely applied due to its straightforward interpretation of the hazard ratio under the assumption that the hazard rates for two subjects remain constant over time. However, in several randomized clinical trials with long-term survival data comparing two new treatments, it is frequently observed that Kaplan-Meier plots exhibit crossing survival curves. This violation of the PH assumption of the Cox PH model can not be applied to evaluate the treatment's effect on survival. This paper introduces a novel long-term survival model with non-PH that incorporates a frailty term into the hazard function. This model allows us to examine the effect of prognostic factors on survival and quantify the degree of unobservable heterogeneity. The model parameters are estimated using the maximum likelihood estimation procedure, and we evaluate the performance of the proposed models through simulation studies. Additionally, we demonstrate the applicability of our approach by fitting the models to a real skin cancer dataset. [ABSTRACT FROM AUTHOR]

Software testing is an important step in software development where inputs are administered repeatedly to detect bugs present in the software. In this paper, we have considered the estimation of total number of bugs and software reliability as a size-biased sampling problem by introducing the concept of eventual bug size as a latent variable. We have developed a Bayesian generalised linear mixed model (GLMM) using software testing detection data to estimate software reliability and stopping phase. The model uses size-biased approach where the probability of detecting a bug is an increasing function of eventual size of the bug which is as an index for the potential number of inputs that may eventually pass through the bug. We have tested the sensitivity of the reliability estimates by varying the number of inputs and detection probability via a simulation study and have found that the key parameters could be accurately estimated. Further, we have applied our model to two empirical data sets – one from a commercial software and the other from ISRO launch mission software testing data set. The hierarchical modelling approach provides a unified modelling framework that may find applications in other fields (e.g. hydrocarbon explorations) apart from software management. [ABSTRACT FROM AUTHOR]

This paper investigates the estimation of a linear regression model with a possible structural break at a known point. We analytically derive the exact formulae of the MSE for the restricted SR and PSR estimators, which shrinks the OLS estimator toward the restriction of no structural break. We compare the MSE performance of restricted/unrestricted SR, PSR, and least squared estimators. We analytically show that the unrestricted SR estimator can have a smaller MSE than the restricted SR estimator even when the restriction is correct. Further, our numerical results show that the unrestricted PSR estimator has the best MSE performance over a wide region of parameter space. These results indicate that the use of the unrestricted PSR estimator is recommended even when a structural break may not exist, although the unrestricted PSR estimator does not take the possibility of no structural break into consideration. [ABSTRACT FROM AUTHOR]

We examine the use of time series data, derived from Electric Cell-substrate Impedance Sensing (ECIS), to differentiate between standard mammalian cell cultures and those infected with a mycoplasma organism. With the goal of easy visualization and interpretation, we perform low-dimensional feature-based classification, extracting application-relevant features from the ECIS time courses. We can achieve very high classification accuracy using only two features, which depend on the cell line under examination. Initial results also show the existence of experimental variation between plates and suggest types of features that may prove more robust to such variation. Our paper is the first to perform a broad examination of ECIS time course features in the context of detecting contamination; to combine different types of features to achieve classification accuracy while preserving interpretability; and to describe and suggest possibilities for ameliorating plate-to-plate variation. [ABSTRACT FROM AUTHOR]

Distribution-free or nonparametric control charts are used for monitoring the process parameters when there is a lack of knowledge about the underlying distribution. In this paper, we investigate a single distribution-free triple exponentially weighted moving average control chart based on the Lepage statistic (referred as TL chart) for simultaneously monitoring shifts in the unknown location and scale parameters of a univariate continuous distribution. The design and implementation of the proposed chart are discussed using time-varying and steady-state control limits for the zero-state case. The run-length distribution of the TL chart is evaluated by performing Monte Carlo simulations. The performance of the proposed chart is compared to those of the existing EWMA-Lepage (EL) and DEWMA-Lepage (DL) charts. It is observed that the TL chart with a time-varying control limit is superior to its competitors, especially for small to moderate shifts in the process parameters. We also provide a real example from a manufacturing process to illustrate the application of the proposed chart. [ABSTRACT FROM AUTHOR]

The growing popularity of personalized medicine motivates people to explore individualized treatment regimes according to heterogeneous characteristics of the patients. For the large-scale data analysis, however, the data are collected at different times and different locations, i.e. subjects are usually from a heterogeneous population, which causes that the optimal treatment regimes also vary for patients across different subgroups. In this paper, we mainly focus on the estimation of optimal treatment regimes for subjects come from a heterogeneous population with high-dimensional data. We first remove the main effects of the covariates for each subgroup to eliminate non-ignorable residual confounding. Based on the centralized outcome, we propose a penalized robust learning that estimates the coefficient matrix of the interactions between covariates and treatment by penalizing pairwise differences of the coefficients of any two subgroups for the same covariate, which can automatically identify the latent complex structure of the coefficient matrix with heterogeneous and homogeneous columns. At the same time, the penalized robust learning can also select the important variables that truly contribute to the individualized treatment decisions with commonly used sparsity structure penalty. Extensive simulation studies show that our proposed method outperforms current popular methods, and it is further illustrated in the real analysis of the Tamoxifen breast cancer data. [ABSTRACT FROM AUTHOR]

In this paper, we consider the structural change in a class of discrete valued time series, where the true conditional distribution of the observations is assumed to be unknown. The conditional mean of the process depends on a parameter $ \theta ^* $ θ ∗ which may change over time. We provide sufficient conditions for the consistency and the asymptotic normality of the Poisson quasi-maximum likelihood estimator (QMLE) of the model. We consider an epidemic change-point detection and propose a test statistic based on the QMLE of the parameter. Under the null hypothesis of a constant parameter (no change), the test statistic converges to a distribution obtained from increments of a Browninan bridge. The test statistic diverges to infinity under the epidemic alternative, which establishes that the proposed procedure is consistent in power. The effectiveness of the proposed procedure is illustrated by simulated and real data examples. [ABSTRACT FROM AUTHOR]

Traffic pattern identification and accident evaluation are essential for improving traffic planning, road safety, and traffic management. In this paper, we establish classification and regression models to characterize the relationship between traffic flows and different time points and identify different patterns of traffic flows by a negative binomial model with smoothing splines. It provides mean response curves and Bayesian credible bands for traffic flows, a single index, and the log-likelihood difference, for traffic flow pattern recognition. We further propose an impact measure for evaluating the influence of accidents on traffic flows based on the fitted negative binomial model. The proposed method has been successfully applied to real-world traffic flows, and it can be used for improving traffic management. [ABSTRACT FROM AUTHOR]

The formulation of variable selection has been widely developed in the Bayesian literature by linking a random binary indicator to each variable. This Bayesian inference has the advantage of stochastically exploring the set of possible sub-models, whatever their dimension. Bayesian selection approaches, appropriate for categorical predictors, are generally beyond the scope of the standard Bayesian selection of regressors in the linear model since all levels of a categorical variable should be jointly handled in the selection procedure. For categorical covariates, new strategies have been developed to detect the effect of grouped covariates rather than the single effect of a quantitative regressor. In this paper, we review three Bayesian selection methods for categorical predictors: Bayesian Group Lasso with Spike and Slab priors, Bayesian Sparse Group Selection and Bayesian Effect Fusion using model-based clustering. The motivation behind this paper is to provide detailed information about the implementation of the three Bayesian selection methods mentioned above, appropriate for categorical predictors, using the JAGS software. Selection performance and sensitivity analysis of the hyperparameters tuning for prior specifications are assessed under various simulated scenarios. JAGS helps user implement these three Bayesian selection methods for more complex model structures such as hierarchical ones with latent layers. [ABSTRACT FROM AUTHOR]

To study the effect of exposure mixture on the continuous health outcomes, one can use the linear model with a weighted sum of multiple standardized exposure variables as an index predictor and its coefficient for the overall effect. The unknown weights typically range between zero and one, indicating contributions of individual exposures to the overall effect. Because the weight parameters present only when the parameter for overall effect is non-zero, testing hypotheses on the overall effect can be challenging, especially when the number of exposure variables is above two. This paper presents a working model based approach to estimate the parameter for overall effect and to test specific hypotheses, including two tests for detecting the overall effect and one test for detecting unequal weights when the overall effect is evident. The statistics are computationally easy and one can apply existing statistical software to perform the analysis. A simulation study shows that the proposed estimators for the parameters of interest may have better finite sample performance than some other estimators. [ABSTRACT FROM AUTHOR]

Clustered current status data are frequently encountered in biomedical research and other areas that require survival analysis. This paper proposes graphical and formal model assessment procedures to evaluate the goodness of fit of the additive hazards model to clustered current status data. The test statistics proposed are based on sums of martingale-based residuals. Relevant asymptotic properties are established, and empirical distributions of the test statistics can be simulated utilizing Gaussian multipliers. Extensive simulation studies confirmed that the proposed test procedures work well for practical scenarios. This proposed method applies when failure times within the same cluster are correlated, and in particular, when cluster sizes can be informative about intra-cluster correlations. The method is applied to analyze clustered current status data from a lung tumorigenicity study. [ABSTRACT FROM AUTHOR]

Distributed interval estimation in linear regression may be computationally infeasible in the presence of big data that are normally stored in different computer servers or in cloud. The existing challenge represents the results from the distributed estimation may still contain redundant information about the population characteristics of the data. To tackle this computing challenge, we develop an optimization procedure to select the best subset from the collection of data subsets, based on which we perform interval estimation in the context of linear regression. The procedure is derived based on minimizing the length of the final interval estimator and maximizing the information remained in the selected data subset, thus is named as the LIC criterion. Theoretical performance of the LIC criterion is studied in this paper together with a simulation study and real data analysis. [ABSTRACT FROM AUTHOR]

In a large-scale environmental health population study that is composed of subprojects, often different fractions of participants out of the total enrolled have measures of specific outcomes. It's conceptually reasonable to assume the association study would benefit from utilizing additional exposure information from those with a specific outcome not measured. Partial least squares regression is a practical approach to determine the exposure-outcome associations for mixture data. Like a typical regression approach, however, the partial least squares regression requires that each data observation must have both complete covariate and outcome for model fitting. In this paper, we propose novel adjustments to the general partial least squares regression to estimate and examine the association effects of individual environmental exposure to an outcome within a more complete context of the study population's environmental mixture exposures. The proposed framework takes advantage of the bilinear model structure. It allows information from all participants, with or without the outcome values, to contribute to the model fitting and the assessment of association effects. Using this proposed framework, incorporation of additional information will lead to smaller root mean square errors in the estimation of association effects, and improve the ability to assess the significance of the effects. [ABSTRACT FROM AUTHOR]