Department of Statistics, UCLA

A brief overview of the papers and commentaries in this special edition.

Ambulance response times in Santa Barbara County for 2006 are analyzed using point process techniques, including kernel intensity estimates and K-functions. Clusters of calls result in significantly higher response times, and this effect is quantified. In particular, calls preceded by other calls within 20 km and within the previous hour are significantly more likely to result in violations. This effect appears to be especially pronounced within semi-rural neighborhoods.

[WestJEM. 2009;10:42-47.]

Introductory remarks

Two correspondences raised concerns or comments about our analyses regarding exaggerated false positives found by differential expression (DE) methods. Here, we discuss the points they raise and explain why we agree or disagree with these points. We add new analysis to confirm that the Wilcoxon rank-sum test remains the most robust method compared to the other five DE methods (DESeq2, edgeR, limma-voom, dearseq, and NOISeq) in two-condition DE analyses after considering normalization and winsorization, the data preprocessing steps discussed in the two correspondences.

Summer temperature extremes can have large impacts on humans and the biosphere, and an increase in heat extremes is one of the most visible symptoms of climate change. Multiple mechanisms have been proposed that would predict faster warming of heat extremes than typical summer days, but it is unclear whether this is occurring. Here, we show that, in both observations and historical climate model simulations, the hottest summer days have warmed at the same pace as the median globally, in each hemisphere, and in the tropics from 1959 to 2023. In contrast, the coldest summer days have warmed more slowly than the median in the global average, a signal that is not simulated in any of 262 simulations across 28 CMIP6 models. The observed stretching of the cold tail indicates that observed summertime temperatures have become more variable despite the lack of hot day amplification. The interannual variability and trend in the warming of both hot and cold extremes compared to the median can be explained from a surface energy balance perspective based on changes in net surface radiation and evaporative fraction. Tropical hot day amplification is projected to emerge in the future (2024-2099, SSP3-7.0 scenario), while Northern Hemisphere heat extremes are expected to continue to follow the median.

Abstract: Antarctic sea ice, a key component in the complex Antarctic climate system, is an important driver and indicator of the global climate. In the relatively short satellite‐observed period from 1979 to 2022 the sea ice extent has continuously increased (contrasting a major decrease in Arctic sea ice) up to a dramatic decrease between 2014 and 2017. Recent years have seen record sea ice lows in February 2022–February 2023. We use a statistical ensemble reconstruction of Antarctic sea ice to put the observed changes into the historical context of the entire 20th century. We propose a seasonal Vector Auto‐Regressive Moving Average (VARMA) model fit in a Bayesian framework using regularized horseshoe priors on the regression coefficients to create a stochastic ensemble reconstruction of monthly Antarctic Sea ice extent from 1900 to 1979. This novel model produces a set of 2,500 plausible sea ice extent reconstructions for the sea ice by sector that incorporate the autocorrelation structure of sea ice over time as well as the dependence of sea ice between the sectors. These fully observed reconstructions exhibit plausible month‐to‐month changes in reconstructed sea ice as well as plausible interactions between the sectors and the total. We reconstruct an overall higher sea ice extent earlier in the 20th century with a relatively sharp decline in the 1970s. These trends agree well with previous reconstructions of Antarctic sea ice based on ice core data, whaling locations, and climatological data, as well as early satellite observations in the reconstruction period.

Dialysis patients experience frequent hospitalizations and a higher mortality rate compared to other Medicare populations, in whom hospitalizations are a major contributor to morbidity, mortality, and healthcare costs. Patients also typically remain on dialysis for the duration of their lives or until kidney transplantation. Hence, there is growing interest in studying the spatiotemporal trends in the correlated outcomes of hospitalization and mortality among dialysis patients as a function of time starting from transition to dialysis across the United States Utilizing national data from the United States Renal Data System (USRDS), we propose a novel multivariate spatiotemporal functional principal component analysis model to study the joint spatiotemporal patterns of hospitalization and mortality rates among dialysis patients. The proposal is based on a multivariate Karhunen-Loéve expansion that describes leading directions of variation across time and induces spatial correlations among region-specific scores. An efficient estimation procedure is proposed using only univariate principal components decompositions and a Markov Chain Monte Carlo framework for targeting the spatial correlations. The finite sample performance of the proposed method is studied through simulations. Novel applications to the USRDS data highlight hot spots across the United States with higher hospitalization and/or mortality rates and time periods of elevated risk.

Abstract: The cool season (November–March) of 2022–2023 was exceptional in the western United States (US), with the highest precipitation totals in ≥128 years in some areas. Recent precipitation extremes and expectations based on thermodynamics motivate us to evaluate the evidence for an anthropogenic intensification of western US cool‐season precipitation to date. Over cool seasons 1951–2023, trends in precipitation totals on the wettest cool‐season days were neutral or negative across the western US, and significantly negative in northern California and parts of the Pacific Northwest, counter to the expected net intensification effect from anthropogenic forcing. Multiple reanalysis data sets indicate a corresponding lack of increase in moisture transports into the western US, suggesting that atmospheric circulation trends over the North Pacific have counteracted the increases in atmospheric moisture expected from warming alone. The lack of precipitation intensification to date is generally consistent with climate model simulations. A large ensemble of 648 simulations from 35 climate models suggests it is too soon to detect anthropogenic intensification of precipitation across much of the western US. In California, the 35‐model median time of emergence for intensification of the wettest days is 2080 under a mid‐level emissions scenario. On the other hand, observed reductions of precipitation extremes in California and the Pacific Northwest are near the lower edge of the large ensemble of simulated trends, calling into question model representation of western US precipitation variability.