Search

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]

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]

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, 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 the empirical performance of the proposed methods. We further apply the proposed methodology to analyze different levels of COVID-19 data in the United States. [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]

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]

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]

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]

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]

Applied researchers often confront two issues when using the fixed effect-two-stage least squares (FE-2SLS) estimator for panel data models. One is that it may lose its consistency due to too many instruments. The other is that the gain of using FE-2SLS may not exceed its loss when the endogeneity is weak. In this paper, an $ L_{2} $ L 2 Boosting regularization procedure for panel data models is proposed to tackle the many instruments issue. We then construct a Stein-like model-averaging estimator to take advantage of FE and FE-2SLS-Boosting estimators. Finite sample properties are examined in Monte Carlo and an empirical application is presented. [ABSTRACT FROM AUTHOR]

The Sharpe ratio function is a commonly used risk/return measure in financial econometrics. To estimate this function, most existing methods take a two-step procedure that first estimates the mean and volatility functions separately and then applies the plug-in method. In this paper, we propose a direct method via local maximum likelihood to simultaneously estimate the Sharpe ratio function and the negative log-volatility function as well as their derivatives. We establish the joint limiting distribution of the proposed estimators, and moreover extend the proposed method to estimate the multivariate Sharpe ratio function. We also evaluate the numerical performance of the proposed estimators through simulation studies, and compare them with existing methods. Finally, we apply the proposed method to the three-month US Treasury bill data and that captures a well-known covariate-dependent effect on the Sharpe ratio. [ABSTRACT FROM AUTHOR]

The present communication develops the tools for estimation and prediction of the Burr-III distribution under unified progressive hybrid censoring scheme. The maximum likelihood estimates of model parameters are obtained. It is shown that the maximum likelihood estimates exist uniquely. Expectation maximization and stochastic expectation maximization methods are employed to compute the point estimates of unknown parameters. Based on the asymptotic distribution of the maximum likelihood estimators, approximate confidence intervals are proposed. In addition, the bootstrap confidence intervals are constructed. Furthermore, the Bayes estimates are derived with respect to squared error and LINEX loss functions. To compute the approximate Bayes estimates, Metropolis–Hastings algorithm is adopted. The highest posterior density credible intervals are obtained. Further, maximum a posteriori estimates of the model parameters are computed. The Bayesian predictive point, as well as interval estimates, are proposed. A Monte Carlo simulation study is employed in order to evaluate the performance of the proposed statistical procedures. Finally, two real data sets are considered and analysed to illustrate the methodologies established in this paper. [ABSTRACT FROM AUTHOR]

This paper presents a class of regression models with cumulative logistic functions that are chiefly designed to analyse spatially dependent ordinal data. In contrast to previous works, the proposed model requires neither the sites to be regularly spaced nor the assumption of an underlying continuous variable. It belongs to a more general class of Markov random field models, and can be considered an extension of the ordinal regression model with the proportional odds link function. Our proposed model allows practitioners to interpret the model parameters using odds ratios. Apart from the theoretical developments, this work also highlights the practical aspects of model fitting, including parameterisation, selection of neighbourhood, and calculation of standard errors. Simulation studies with regularly and irregularly spaced sites were conducted. Modelling strategies including pseudo-likelihood methods were found to be useful in both settings. The proposed model and the non-spatial counterpart were applied to the daily air quality index measured in the United Kingdom. The results indicate the presence of spatial effects and the incorporation of spatial effects led to better model performance in terms of various goodness-of-fit measures. [ABSTRACT FROM AUTHOR]

Applications of circular regression models are ubiquitous in many disciplines, particularly in meteorology, biology and geology. In circular regression models, variable selection problem continues to be a remarkable open question. In this paper, we address variable selection in linear-circular regression models where uni-variate linear dependent and a mixed set of circular and linear independent variables constitute the data set. We consider Bayesian lasso which is a popular choice for variable selection in classical linear regression models. We show that Bayesian lasso in linear-circular regression models is not able to produce robust inference as the coefficient estimates are sensitive to the choice of hyper-prior setting for the tuning parameter. To eradicate the problem, we propose a robustified Bayesian lasso that is based on an empirical Bayes (EB) type methodology to construct a hyper-prior for the tuning parameter while using Gibbs Sampling. This hyper-prior construction is computationally more feasible than the hyper-priors that are based on correlation measures. We show in a comprehensive simulation study that Bayesian lasso with EB-GS hyper-prior leads to a more robust inference. Overall, the method offers an efficient Bayesian lasso for variable selection in linear-circular regression while reducing model complexity. [ABSTRACT FROM AUTHOR]

Brain imaging research is a very challenging topic due to complex structure and lack of explicitly identifiable features in the image. With the advancement of magnetic resonance imaging (MRI) technologies, such as diffusion tensor imaging (DTI), developing classification methods to improve clinical diagnosis is crucial. This paper proposes a classification method for DTI data based on a novel neural network strategy that combines a convolutional neural network (CNN) with a multilayer neural network using central-peripheral deviation (CPD), which reflects diffusion dynamics in the white matter by spatially evaluating the deviation of diffusion coefficients between the inner and outer parts of the brain. In our method, a multilayer perceptron (MLP) using CPD is combined with the final layers for classification after reducing the dimensions of images in the convolutional layers of the neural network architecture. In terms of training loss and the classification error, the proposed classification method improves the existing image classification with CNN. For real data analysis, we demonstrate how to process raw DTI image data sets obtained from a traumatic brain injury study (MagNeTS) and a brain atlas construction study (ICBM), and apply the proposed approach to the data, successfully improving classification performance with two age groups. [ABSTRACT FROM AUTHOR]

This paper characterizes the legislators voting behavior in the Colombian Senate 2010-2014, by implementing a one-dimensional standard Bayesian ideal point estimator via Markov chain Monte Carlo algorithms. Our main goal is to retrieve the political preferences of legislators from their roll-call voting records, which individualizes the electoral behavior of the legislative chamber. Furthermore, we conclude about the nature of the latent trait underlying the deputies voting decisions and the legislators locations in political space. Finally, we also offer several methodological and theoretical tools to guide the analysis of nominal voting data in the context of unbalanced parliaments (multi-party systems), taking as reference the particular case of the Colombian Senate. [ABSTRACT FROM AUTHOR]

Calculating the entropy for complex systems is a significant problem in science and engineering problems. However, this calculation is usually computationally expensive when the entropy is computed directly. This paper introduces three classes of information sources that for all members of each class, the entropy value is the same. These classes are characterized according to special dynamics created by three kinds of self-mappings on Ω, and A, where Ω is a probability space and A is a finite set. An approximation of rank variables of the product of information sources is made, and it is proved that the topological entropy of the product of two information sources is equal to the summation of their topological entropies. [ABSTRACT FROM AUTHOR]

Although there have been a number of available tests of bivariate exchangeability, i.e. bivariate symmetry for bivariate distributions, the literature is void of tests whether a multivariate distribution with more than two dimensions is exchangeable or not. In this paper, multivariate permutation tests of exchangeability of multivariate distributions are proposed, which are based on the non-parametric combination methodology, i.e. on combining non-parametric bivariate exchangeability tests. Numerical experiments on real as well as simulated multivariate data with more than two dimensions are presented here. The multivariate permutation test turns out to be typically more powerful than a bivariate exchangeability test performed only over a single pair of variables, and also more suitable compared to tests exploiting the approaches of Benjamini–Yekutieli or Bonferroni. [ABSTRACT FROM AUTHOR]

This paper builds on the recently proposed prediction test for muliple endpoints. The prediction test combines information across multiple endpoints while accounting for the correlation between them. The test performs well with small samples relative to the number of endpoints of interest and is flexible in the hypotheses across the individual endpoints that can be combined. The prediction test addresses a global hypothesis that is of particular interest in early-stage studies and can be used as justification for continuing on to a larger trial. However, the prediction test has several limitations which we seek to address. First, the prediction test is overly conservative when both the effect sizes across all endpoints and the number of endpoints are small. By using a parametric bootstrap to estimate the null distribution, we show that the test achieves the nominal error rate in this situation and increases the power of the test. Second, we provide a framework to allow for predictions of a difference on one or more endpoints. Finally, we extend the test with a composite null hypothesis that allows for different null hypothesized predictive abilities across the endpoints which can be especially useful if the study contains both familiar and novel endpoints. We use an example from a physical activity trial to illustrate these extensions. [ABSTRACT FROM AUTHOR]

The joint models for longitudinal and survival data have recently received significant attention in medical and epidemiological studies. Joint models typically combine linear mixed effects models for repeated measurement data and Cox models for survival time. When we are jointly modeling the longitudinal and survival data, variable selection and efficient estimation of parameters are especially important for performing reliable statistical analyzes, both of which are currently lacking in the literature. In this paper we discuss the pretest and shrinkage estimation methods for jointly modeling longitudinal data and survival time data when some of the covariates in both longitudinal and survival components may not be relevant for predicting survival times. In this situation, we fit two models: the full model that contains all the covariates and the subset model that contains a reduced number of covariates. We combine the full model estimators and the estimators that are restricted by a linear hypothesis to define pretest and shrinkage estimators. We provide their numerical mean squared errors (MSE) and relative MSE. We show that if the shrinkage dimension exceeds two, the risk of the shrinkage estimators is strictly less than that of the full model estimators. Our proposed methods are illustrated by extensive simulation studies and a real-data example. [ABSTRACT FROM AUTHOR]

The C statistic, also known as the Cash statistic, is often used in astronomy for the analysis of low-count Poisson data. The main advantage of this statistic, compared to the more commonly used χ 2 statistic, is its applicability without the need to combine data points. This feature has made the C statistic a very useful method to analyze Poisson data that have small (or even null) counts in each resolution element. One of the challenges of the C statistic is that its probability distribution, under the null hypothesis that the data follow a parent model, is not known exactly. This paper presents an effort towards improving our understanding of the C statistic by studying (a) the distribution of C statistic for a fully specified model, (b) the distribution of Cmin resulting from a maximum-likelihood fit to a simple one-parameter constant model, i.e. a model that represents the sample mean of N Poisson measurements, and (c) the distribution of the associated Δ C statistic that is used for parameter estimation. The results confirm the expectation that, in the high-count limit, both C statistic and Cmin have the same mean and variance as a χ 2 statistic with same number of degrees of freedom. It is also found that, in the low-count regime, the expectation of the C statistic and Cmin can be substantially lower than for a χ 2 distribution. The paper makes use of recent X-ray observations of the astronomical source PG 1116+215 to illustrate the application of the C statistic to Poisson data. [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]

Multistage sequential decision-making occurs in many real-world applications such as healthcare diagnosis and treatment. One concrete example is when the doctors need to decide to collect which kind of information from subjects so as to make the good medical decision cost-effectively. In this paper, an active learning-based method is developed to model the doctors' decision-making process that actively collects necessary information from each subject in a sequential manner. The effectiveness of the proposed model, especially its two-stage version, is validated on both simulation studies and a case study of common bile duct stone evaluation for pediatric patients. [ABSTRACT FROM AUTHOR]

High-throughput plant phenotyping (HTPP) has become an emerging technique to study plant traits due to its fast, labor-saving, accurate and non-destructive nature. It has wide applications in plant breeding and crop management. However, the resulting massive image data has raised a challenge associated with efficient plant traits prediction and anomaly detection. In this paper, we propose a two-step image-based online detection framework for monitoring and quick change detection of the individual plant leaf area via real-time imaging data. Our proposed method is able to achieve a smaller detection delay compared with some baseline methods under some predefined false alarm rate constraint. Moreover, it does not need to store all past image information and can be implemented in real time. The efficiency of the proposed framework is validated by a real data analysis. [ABSTRACT FROM AUTHOR]

In this paper, we present an efficient statistical method (denoted as 'Adaptive Resources Allocation CUSUM') to robustly and efficiently detect the hotspot with limited sampling resources. Our main idea is to combine the multi-arm bandit (MAB) and change-point detection methods to balance the exploration and exploitation of resource allocation for hotspot detection. Further, a Bayesian weighted update is used to update the posterior distribution of the infection rate. Then, the upper confidence bound (UCB) is used for resource allocation and planning. Finally, CUSUM monitoring statistics to detect the change point as well as the change location. For performance evaluation, we compare the performance of the proposed method with several benchmark methods in the literature and showed the proposed algorithm is able to achieve a lower detection delay and higher detection precision. Finally, this method is applied to hotspot detection in a real case study of county-level daily positive COVID-19 cases in Washington State WA) and demonstrates the effectiveness with very limited distributed samples. [ABSTRACT FROM AUTHOR]

Count data occur widely in many bio-surveillance and healthcare applications, e.g. the numbers of new patients of different types of infectious diseases from different cities/counties/states repeatedly over time, say, daily/weekly/monthly. For this type of count data, one important task is the quick detection and localization of hot-spots in terms of unusual infectious rates so that we can respond appropriately. In this paper, we develop a method called Poisson assisted Smooth Sparse Tensor Decomposition (PoSSTenD), which not only detect when hot-spots occur but also localize where hot-spots occur. The main idea of our proposed PoSSTenD method is articulated as follows. First, we represent the observed count data as a three-dimensional tensor including (1) a spatial dimension for location patterns, e.g. different cities/countries/states; (2) a temporal domain for time patterns, e.g. daily/weekly/monthly; (3) a categorical dimension for different types of data sources, e.g. different types of diseases. Second, we fit this tensor into a Poisson regression model, and then we further decompose the infectious rate into two components: smooth global trend and local hot-spots. Third, we detect when hot-spots occur by building a cumulative sum (CUSUM) control chart and localize where hot-spots occur by their LASSO-type sparse estimation. The usefulness of our proposed methodology is validated through numerical simulation studies and a real-world dataset, which records the annual number of 10 different infectious diseases from 1993 to 2018 for 49 mainland states in the United States. [ABSTRACT FROM AUTHOR]

The main aim of this paper is to propose a novel method (RMD-MRCD-PCA) of identification of High Leverage Points (HLPs) in high-dimensional sparse data. It is to address the weakness of the Robust Mahalanobis Distance (RMD) method which is based on the Minimum Regularized Covariance Determinant (RMD-MRCD), which indicates a decrease in its performance as the number of independent variables (p) increases. The RMD-MRCD-PCA is developed by incorporating the Principal Component Analysis (PCA) in the MRCD algorithm whereby this robust approach shrinks the covariance matrix to make it invertible and thus, can be employed to compute the RMD for high dimensional data. A simulation study and two real data sets are used to illustrate the merit of our proposed method compared to the RMD-MRCD and Robust PCA (ROBPCA) methods. Findings show that the performance of the RMD-MRCD is similar to the performance of the RMD-MRCD-PCA for p close to 200. However, its performance tends to decrease when the number of p is more than 200 and worsens at p equals 700 and larger. On the other hand, the ROBPCA is not effective for less than 20% contamination as it suffers from serious swamping problems. [ABSTRACT FROM AUTHOR]

While there is wide agreement that physical activity is an important component of a healthy lifestyle, it is unclear how many people adhere to public health recommendations on physical activity. The Physical Activity Guidelines (PAG), published by the CDC, provides guidelines to American adults, but it is difficult to assess compliance with these guidelines. The PAG further complicates adherence assessment by recommending activity to occur in at least 10 min bouts. To better understand the measurement capabilities of various instruments to quantify activity, and to propose an approach to evaluate activity relative to the PAG, researchers at Iowa State University administered the Physical Activity Measurement Survey (PAMS) to over 1000 participants in four different Iowa counties. In this paper, we develop a two-part Bayesian measurement error model and apply it to the PAMS data in order to assess compliance with the PAG in the Iowa adult population. The model accurately accounts for the 10 min bout requirement put forth in the PAG. The measurement error model corrects biased estimates and accounts for day-to-day variation in activity. The model is also applied to the nationally representative National Health and Nutrition Examination Survey. [ABSTRACT FROM AUTHOR]

Joint modelling of longitudinal and time-to-event data is usually described by a joint model which uses shared or correlated latent effects to capture associations between the two processes. Under this framework, the joint distribution of the two processes can be derived straightforwardly by assuming conditional independence given the random effects. Alternative approaches to induce interdependency into sub-models have also been considered in the literature and one such approach is using copulas to introduce non-linear correlation between the marginal distributions of the longitudinal and time-to-event processes. The multivariate Gaussian copula joint model has been proposed in the literature to fit joint data by applying a Monte Carlo expectation-maximisation algorithm. In this paper, we propose an exact likelihood estimation approach to replace the more computationally expensive Monte Carlo expectation-maximisation algorithm and we consider results based on using both the multivariate Gaussian and t copula functions. We also provide a straightforward way to compute dynamic predictions of survival probabilities, showing that our proposed model is comparable in prediction performance to the shared random effects joint model. [ABSTRACT FROM AUTHOR]

In this paper we introduce a discrete-time and continuous state-space Markov stationary process { X n ; n = 1 , 2 , ... } , where X n has a two-parameter Weibull distribution, X n 's are dependent and there is a positive probability that X n = X n + 1 . The motivation came from the gold price data where there are several instances for which X n = X n + 1 . Hence, the existing methods cannot be used to analyze this data. We derive different properties of the proposed Weibull process. It is observed that the joint cumulative distribution function of X n and X n + 1 has a very convenient copula structure. Hence, different dependence properties and dependence measures can be obtained. The maximum likelihood estimators cannot be obtained in explicit forms, we have proposed a simple profile likelihood method to compute these estimators. We have used this model to analyze two synthetic data sets and one gold price data set of the Indian market, and it is observed that the proposed model fits quite well with the data set. [ABSTRACT FROM AUTHOR]

Identifying and locating areas – hot spots – that present high concentration of observations in a high-dimensional data set is crucial in many data processing and analysis methods and techniques, since observations that belong to the same hot spot share information and behave in a similar way. A useful tool towards that aim is the reduction of the data dimensionality and the graphical representation of them. In the present paper, a new method to identify and locate hot spots is proposed, based on the Andrews curves. Simulations results demonstrate the performance of the proposed method, which is also applied to a high-dimensional data set, regarding caregiver distress related to symptoms of people with neurocognitive disorder and to the mental effects of the recent outbreak of the COVID-19 pandemic. [ABSTRACT FROM AUTHOR]

In this paper, we develop a mixture of autoregressive (MoAR) process model with time varying and freely indexed covariates under the flexible class of two–piece distributions using the scale mixtures of normal (TP-SMN) family. This novel family of time series (TP-SMN-MoAR) models was used to examine flexible and robust clustering of reported cases of Covid-19 across 313 counties in the U.S. The TP-SMN distributions allow for symmetrical/ asymmetrical distributions as well as heavy-tailed distributions providing for flexibility to handle outliers and complex data. Developing a suitable hierarchical representation of the TP-SMN family enabled the construction of a pseudo-likelihood function to derive the maximum pseudo-likelihood estimates via an EM-type algorithm. [ABSTRACT FROM AUTHOR]

To make informative public policy decisions in battling the ongoing COVID-19 pandemic, it is important to know the disease prevalence in a population. There are two intertwined difficulties in estimating this prevalence based on testing results from a group of subjects. First, the test is prone to measurement error with unknown sensitivity and specificity. Second, the prevalence tends to be low at the initial stage of the pandemic and we may not be able to determine if a positive test result is a false positive due to the imperfect test specificity. The statistical inference based on a large sample approximation or conventional bootstrap may not be valid in such cases. In this paper, we have proposed a set of confidence intervals, whose validity doesn't depend on the sample size in the unweighted setting. For the weighted setting, the proposed inference is equivalent to hybrid bootstrap methods, whose performance is also more robust than those based on asymptotic approximations. The methods are used to reanalyze data from a study investigating the antibody prevalence in Santa Clara County, California in addition to several other seroprevalence studies. Simulation studies have been conducted to examine the finite-sample performance of the proposed method. [ABSTRACT FROM AUTHOR]

In this paper, a statistical model for COVID-19 infection dynamics is described, using only the observed daily statistics of infected individuals. For this purpose, two special classes of branching processes without or with an immigration component are considered. These models are intended to estimate the main parameter of the infection and to give a prediction of the mean value of the non-observed population of the infected individuals. This is a serious advantage in comparison with other more complicated models where the officially reported data are not sufficient for estimation of the model parameters. The model is applied for different regions in the world and the corresponding parameters of the infection dynamics are estimated. [ABSTRACT FROM AUTHOR]

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 dynamic infectious disease model, we tested the impact of temperature, precipitation, state emergency order and stay home order on the spread of COVID-19 based on the United States county-wise daily counts data. Our numerical results show that warm and humid weather can significantly slow the spread of COVID-19, and the state emergency and stay home orders also help to slow it. This finding provides guidance and support to future policies or acts for mitigating the community transmission and lowering the mortality rate of COVID-19. [ABSTRACT FROM AUTHOR]

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 many. Yet, due to the exponential nature of epidemic growth, the error due to this simplification can be rapidly amplified, and its extent remains unknown. How can this intrinsic variability be percolated into epidemic models, and its impact, better quantified? We study this question here through a Bayesian perspective that captures at scale the heterogeneity of a population and environmental conditions, creating a bridge between the traditional agent-based and compartmental approaches. We use our model to simulate the spread as well as the impact of different social distancing strategies on real COVID-19 data, and highlight the significant impact of the heterogeneity. We emphasize that the contribution of this paper focuses on discussing the importance of the impact of R's heterogeneity on uncertainty quantification from a statistical viewpoint, rather than developing new predictive models. [ABSTRACT FROM AUTHOR]

Skip-lot sampling plan is often applied in industries for reducing the cost and effort of the inspection of the product having excellent quality history. Consequence of skip-lot sampling plans is to reduce the cost of inspection so which are more attractive in economical aspect. In this paper, we develop a sampling plan by incorporating the idea of resampling in two-level skip lot sampling plan and the new plan is designated as SkSP-2L.1-R. This paper presents the Markov chain formulation of the proposed plan along with the derivation of performance measures of the plan. We also provide the designing methodology to determine the optimal parameters of the SkSP-2L.1-R plan so as to minimize the average sample number by using two points on the operating characteristic curve approach. By contemplating various combinations of producer and consumer quality levels along with respective risks, a table is constructed to determine the optimal parameters. An industrial application of the proposed SkSP-2L.1-R plan is discussed. The SkSP-2L.1-R with single sampling plan as a reference plan is compared with the conventional single sampling plan, SkSP-2 plan and SkSP-2-R plan and proved that the proposed SkSP-2L.1-R plan outperforms these plans. [ABSTRACT FROM AUTHOR]

The Cash statistic, also known as the C statistic, is commonly used for the analysis of low-count Poisson data, including data with null counts for certain values of the independent variable. The use of this statistic is especially attractive for low-count data that cannot be combined, or re-binned, without loss of resolution. This paper presents a new maximum-likelihood solution for the best-fit parameters of a linear model using the Poisson-based Cash statistic. The solution presented in this paper provides a new and simple method to measure the best-fit parameters of a linear model for any Poisson-based data, including data with null counts. In particular, the method enforces the requirement that the best-fit linear model be non-negative throughout the support of the independent variable. The method is summarized in a simple algorithm to fit Poisson counting data of any size and counting rate with a linear model, by-passing entirely the use of the traditional χ 2 statistic. [ABSTRACT FROM AUTHOR]

Ecological momentary assessment (EMA) studies investigate intensive repeated observations of the current behavior and experiences of subjects in real time. In particular, such studies aim to minimize recall bias and maximize ecological validity, thereby strengthening the investigation and inference of microprocesses that influence behavior in real-world contexts by gathering intensive information on the temporal patterning of behavior of study subjects. Throughout this paper, we focus on the data analysis of an EMA study that examined behavior of intermittent smokers (ITS). Specifically, we sought to explore the pattern of clustered smoking behavior of ITS, or smoking 'bouts', as well as the covariates that predict such smoking behavior. To do this, in this paper we introduce a framework for characterizing the temporal behavior of ITS via the functions of event gap time to distinguish the smoking bouts. We used the time-varying coefficient models for the cumulative log gap time and to characterize the temporal patterns of smoking behavior, while simultaneously adjusting for behavioral covariates, and incorporated the inverse probability weighting into the models to accommodate missing data. Simulation studies showed that irrespective of whether missing by design or missing at random, the model was able to reliably determine prespecified time-varying functional forms of a given covariate coefficient, provided the the within-subject level was small. [ABSTRACT FROM AUTHOR]

Comparative lifetime experiments are of particular importance in production processes when one wishes to determine the relative merits of several competing products with regard to their reliability. This paper confines itself to the data obtained by running a joint progressive Type-II censoring plan on samples in a combined manner. The problem of Bayesian predicting failure times of surviving units is discussed in details when parent populations are exponential. Two real data sets are analyzed in order to illustrate all the inferential procedures developed here. When destructive experiments under a censoring scheme finished, the researchers are usually interested to estimate remaining lifetimes of surviving units for sequel experiments. Findings of this paper are useful for these purposes specially when samples are non-homogeneous such as those taken from industrial storages. [ABSTRACT FROM AUTHOR]

In this paper, we propose two multimodal circular distributions which are suitable for modeling circular data sets with two or more modes. Both distributions belong to the regular exponential family of distributions and are considered as extensions of the von Mises distribution. Hence, they possess the highly desirable properties, such as the existence of non-trivial sufficient statistics and optimal inferences for their parameters. Fine particulates (PM2.5) are generally emitted from activities such as industrial and residential combustion and from vehicle exhaust. We illustrate the utility of our proposed models using a real data set consisting of fine particulates (PM2.5) pollutant levels in Houston region during Fall season in 2019. Our results provide a strong evidence that its diurnal pattern exhibits four modes; two peaks during morning and evening rush hours and two peaks in between. [ABSTRACT FROM AUTHOR]

In this paper, two new general families of distributions supported on the unit interval are introduced. The proposed families include several known models as special cases and define at least twenty (each one) new special models. Since the list of well-being indicators may include several double bounded random variables, the applicability for modeling those is the major practical motivation for introducing the distributions on those families. We propose a parametrization of the new families in terms of the median and develop a shiny application to provide interactive density shape illustrations for some special cases. Various properties of the introduced families are studied. Some special models in the new families are discussed. In particular, the complementary unit Weibull distribution is studied in some detail. The method of maximum likelihood for estimating the model parameters is discussed. An extensive Monte Carlo experiment is conducted to evaluate the performances of these estimators in finite samples. Applications to the literacy rate in Brazilian and Colombian municipalities illustrate the usefulness of the two new families for modeling well-being indicators. [ABSTRACT FROM AUTHOR]

In this paper, we propose a hierarchical Bayesian approach for modeling the evolution of the 7-day moving average for the number of deaths due to COVID-19 in a country, state or city. The proposed approach is based on a Gaussian process regression model. The main advantage of this model is that it assumes that a nonlinear function f used for modeling the observed data is an unknown random parameter in opposite to usual approaches that set up f as being a known mathematical function. This assumption allows the development of a Bayesian approach with a Gaussian process prior over f. In order to estimate the parameters of interest, we develop an MCMC algorithm based on the Metropolis-within-Gibbs sampling algorithm. We also present a procedure for making predictions. The proposed method is illustrated in a case study, in which, we model the 7-day moving average for the number of deaths recorded in the state of São Paulo, Brazil. Results obtained show that the proposed method is very effective in modeling and predicting the values of the 7-day moving average. [ABSTRACT FROM AUTHOR]

This paper develops an extension of spatiotemporal models that handle count data using nonhomogeneous Poisson processes. In this new proposal, we incorporate a seasonal cycle component in the definition of the intensity function to control possible effects produced by the occurrence of the event of interest in regular periods. The seasonal cycle can cause problems in estimating the shape parameter of the Weibull and generalized Goel intensity functions. This shape parameter serves to confront the research hypothesis that seeks to identify a trend in the occurrence rate of an event of interest. In the case of the Weibull intensity function, a value significantly equal to one of the shape parameters indicates a constant rate of occurrence, less than one indicates a decreasing rate, and greater than one indicates an increasing rate. In the case of the Goel intensity function, parameter values less than or equal to one indicate a decreasing occurrence rate, and values greater than one indicate the presence of a change point. We also built a spatial model using the Musa-Okumoto intensity function as an alternative to approximate counting processes for which there is a decreasing trend in the occurrence rate of the event of interest. We estimated the parameters of the proposed method from a Bayesian perspective. Finally, we fitted the proposed model and compared it with other approximations to analyze the frequency of extreme rainfall in the northern region of the states of Maranhão and Piauí in northeastern Brazil over ten years. Among the main results, we found that (1) the proposed method has proven superior in terms of fit and prediction performance than the other models, and (2) unlike other approximations, the proposed model does not detect changes in the rate of extreme rainfall occurrences. [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]

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]