Statistical Methods and Applications最新文献

The first cluster of coronavirus cases in Europe was officially detected on 21st February 2020 in Northern Italy, even if recent evidence showed sporadic first cases in Europe since the end of 2019. In this study, we have tested the presence of coronavirus in Italy and, even more importantly, we have assessed whether the virus had already spread sooner than 21st February. We use a counterfactual approach and certified daily data on the number of deaths (deaths from any cause, not only related to coronavirus) at the municipality level. Our estimates confirm that coronavirus began spreading in Northern Italy in mid-January.

Graphical models are powerful tools that are regularly used to investigate complex dependence structures in high-throughput biomedical datasets. They allow for holistic, systems-level view of the various biological processes, for intuitive and rigorous understanding and interpretations. In the context of large networks, Bayesian approaches are particularly suitable because it encourages sparsity of the graphs, incorporate prior information, and most importantly account for uncertainty in the graph structure. These features are particularly important in applications with limited sample size, including genomics and imaging studies. In this paper, we review several recently developed techniques for the analysis of large networks under non-standard settings, including but not limited to, multiple graphs for data observed from multiple related subgroups, graphical regression approaches used for the analysis of networks that change with covariates, and other complex sampling and structural settings. We also illustrate the practical utility of some of these methods using examples in cancer genomics and neuroimaging.

Detecting changes in COVID-19 disease transmission over time is a key indicator of epidemic growth. Near real-time monitoring of the pandemic growth is crucial for policy makers and public health officials who need to make informed decisions about whether to enforce lockdowns or allow certain activities. The effective reproduction number $R_t$ is the standard index used in many countries for this goal. However, it is known that due to the delays between infection and case registration, its use for decision making is somewhat limited. In this paper a near real-time COVINDEX is proposed for monitoring the evolution of the pandemic. The index is computed from predictions obtained from a GAM beta regression for modelling the test positive rate as a function of time. The proposal is illustrated using data on COVID-19 pandemic in Italy and compared with $R_t$ . A simple chart is also proposed for monitoring local and national outbreaks by policy makers and public health officials.

Providing support outside the household can be considered an actual sign of an active social life for the elderly. Adopting an ego-network perspective, we study support Italian elders provide to kin or non-kin. More specifically, using Italian survey data, we build the ego-centered networks of social contacts elders entertain and the ego-networks of support elders provide to other non-cohabitant kin or non-kin. Since ego-network data are inherently multilevel, we use Bayesian multilevel models to analyze variation in support ties, controlling for the characteristics of elders and their contacts. This modeling strategy enables dealing with sparseness and alter-alter overlap in the ego support network data and to disentangle the effects related to the ego (the elder), the dyad ego-alter, the kind of support provided, as well as social contacts and contextual variables. The results suggest that the elderly in Italy who provide support outside their household - compared to all elders in the sample - are younger, healthier, more educated, and embedded in a more diversified ego-network of social contacts. The latter also conveys both the type and the recipient of the support, with the elderly who entertain few relationships with kin being more prone to provide aid to non-kin. Further, a "peer homophily" effect in directing elder support to a non-kin is also found.

The purpose of this paper is to show in regression clustering how to choose the most relevant solutions, analyze their stability, and provide information about best combinations of optimal number of groups, restriction factor among the error variance across groups and level of trimming. The procedure is based on two steps. First we generalize the information criteria of constrained robust multivariate clustering to the case of clustering weighted models. Differently from the traditional approaches which are based on the choice of the best solution found minimizing an information criterion (i.e. BIC), we concentrate our attention on the so called optimal stable solutions. In the second step, using the monitoring approach, we select the best value of the trimming factor. Finally, we validate the solution using a confirmatory forward search approach. A motivating example based on a novel dataset concerning the European Union trade of face masks shows the limitations of the current existing procedures. The suggested approach is initially applied to a set of well known datasets in the literature of robust regression clustering. Then, we focus our attention on a set of international trade datasets and we provide a novel informative way of updating the subset in the random start approach. The Supplementary material, in the spirit of the Special Issue, deepens the analysis of trade data and compares the suggested approach with the existing ones available in the literature.

Motivated by the ongoing COVID-19 pandemic, this article introduces Bayesian dynamic network actor models for the analysis of infected individuals' movements in South Korea during the first three months of 2020. The relational event data modelling framework makes use of network statistics capturing the structure of movement events from and to several country's municipalities. The fully probabilistic Bayesian approach allows to quantify the uncertainty associated to the relational tendencies explaining where and when movement events are established and where they are directed. The observed patient movements' patterns at an early stage of the pandemic can provide interesting insights about the spread of the disease in the Asian country.

An important problem in network analysis is the online detection of anomalous behaviour. In this paper, we introduce a network surveillance method bringing together network modelling and statistical process control. Our approach is to apply multivariate control charts based on exponential smoothing and cumulative sums in order to monitor networks generated by temporal exponential random graph models (TERGM). The latter allows us to account for temporal dependence while simultaneously reducing the number of parameters to be monitored. The performance of the considered charts is evaluated by calculating the average run length and the conditional expected delay for both simulated and real data. To justify the decision of using the TERGM to describe network data, some measures of goodness of fit are inspected. We demonstrate the effectiveness of the proposed approach by an empirical application, monitoring daily flights in the United States to detect anomalous patterns.

The special issue on Statistical Analysis of Networks aspires to convey the breadth and depth of statistical learning with networks, ranging from networks that are observed to networks that are unobserved and learned from data. It includes ten select papers with methodological and theoretical advances, and demonstrates the usefulness of the network paradigm by applications to current problems.

We propose a weighted stochastic block model (WSBM) which extends the stochastic block model to the important case in which edges are weighted. We address the parameter estimation of the WSBM by use of maximum likelihood and variational approaches, and establish the consistency of these estimators. The problem of choosing the number of classes in a WSBM is addressed. The proposed model is applied to simulated data and an illustrative data set.