Nonparametric Clustering Research Articles

Investigating a hybrid mixed population leads to recognizing a new species of Arctostaphylos (Ericaceae).

While investigating the potential for Arctostaphylos species to hybridize in the mixed populations of Point Sal and Burton Mesa in Santa Barbara County, California, we discovered that Arctostaphylos from the Nipomo Mesa (San Luis Obispo County), formerly considered a northern population of A.rudis, are genetically and morphologically distinct. We name this new taxon A.nipumu after the ytt (Northern Chumash language) word for the Nipomo Mesa region. For morphological and molecular analyses, we sampled 54 plants, focusing on A.purissima, A.rudis, and A.crustacea from multiple species and comparative single species populations. Parametric and nonparametric clustering analyses (STRUCTURE and PCA) of ddRADseq data show that Arctostaphylos from the Nipomo Mesa segregate from all other samples in the dataset. In mixed populations A.purissima and A.crustacea samples cluster with samples from other unmixed populations of the same species but A.rudis samples form two distinct clusters. One is composed of the mixed populations in Santa Barbara County, and the other consists of the Nipomo Mesa population. Additionally, the Santa Barbara County A.rudis samples are admixed in STRUCTURE analysis unlike the samples from the Nipomo Mesa. A principal component analysis of eight morphological characters shows that A.rudis individuals from Santa Barbara County tend to be phenotypically variable, occurring in a wide morphological cluster that overlaps with the tight clusters formed by A.purissima, A.crustacea, and Arctostaphylos from the Nipomo Mesa. Based on this evidence we describe the Nipomo Mesapopulation as a new species of Arctostaphylos. Given its limited and fragmented distribution we believe that A.nipumu is of critical conservation concern.

Bayesian subtyping for multi-state brain functional connectome with application on preadolescent brain cognition.

Converging evidence indicates that the heterogeneity of cognitive profiles may arise through detectable alternations in brain functional connectivity. Despite an unprecedented opportunity to uncover neurobiological subtypes through clustering or subtyping analyses on multi-state functional connectivity, few existing approaches are applicable to accommodate the network topology and unique biological architecture. To address this issue, we propose an innovative Bayesian nonparametric network-variate clustering analysis to uncover subgroups of individuals with homogeneous brain functional network patterns under multiple cognitive states. In light of the existing neuroscience literature, we assume there are unknown state-specific modular structures within functional connectivity. Concurrently, we identify informative network features essential for defining subtypes. To further facilitate practical use, we develop a computationally efficient variational inference algorithm to approximate posterior inference with satisfactory estimation accuracy. Extensive simulations show the superiority of our method. We apply the method to the Adolescent Brain Cognitive Development (ABCD) study, and identify neurodevelopmental subtypes and brain sub-network phenotypes under each state to signal neurobiological heterogeneity, suggesting promising directions for further exploration and investigation in neuroscience.

Longitudinal Data Clustering Methods: A Systematic ReviewLongitudinal Data Clustering Methods: A Systematic Review

Introduction: In the last few decades, in many research fields, different methods were introduced to discover groups with the same trends in longitudinal data. The clustering process is an unsupervised learning method, which classifies longitudinal data based on different criteria by performing algorithms. The current study was performed with the aim of reviewing various methods of longitudinal data clustering, including two general categories of non-parametric methods and model-based methods. Methods: In this research, to obtain related scientific articles, PubMed, Science Direct Scopus, ISI, and Google Scholar were searched between 2000 and 2021. Results: According to our systematic review, the non-parametric k-means Clustering Method utilizing Euclidean distance emerges as a leading approach for clustering longitudinal data. Conclusion: This research, with an overview of the studies done in the field of clustering, can help researchers as a toolbox to choose various methods of longitudinal data clustering in idea generation and choosing the appropriate method in the classification and analysis of longitudinal data.

Bayesian clustering with uncertain data.

Clustering is widely used in bioinformatics and many other fields, with applications from exploratory analysis to prediction. Many types of data have associated uncertainty or measurement error, but this is rarely used to inform the clustering. We present Dirichlet Process Mixtures with Uncertainty (DPMUnc), an extension of a Bayesian nonparametric clustering algorithm which makes use of the uncertainty associated with data points. We show that DPMUnc out-performs existing methods on simulated data. We cluster immune-mediated diseases (IMD) using GWAS summary statistics, which have uncertainty linked with the sample size of the study. DPMUnc separates autoimmune from autoinflammatory diseases and isolates other subgroups such as adult-onset arthritis. We additionally consider how DPMUnc can be used to cluster gene expression datasets that have been summarised using gene signatures. We first introduce a novel procedure for generating a summary of a gene signature on a dataset different to the one where it was discovered, which incorporates a measure of the variability in expression across signature genes within each individual. We summarise three public gene expression datasets containing patients with a range of IMD, using three relevant gene signatures. We find association between disease and the clusters returned by DPMUnc, with clustering structure replicated across the datasets. The significance of this work is two-fold. Firstly, we demonstrate that when data has associated uncertainty, this uncertainty should be used to inform clustering and we present a method which does this, DPMUnc. Secondly, we present a procedure for using gene signatures in datasets other than where they were originally defined. We show the value of this procedure by summarising gene expression data from patients with immune-mediated diseases using relevant gene signatures, and clustering these patients using DPMUnc.

Characterizing the clinical heterogeneity of early symptomatic Alzheimer's disease: a data-driven machine learning approach.

Alzheimer's disease (AD) is highly heterogeneous, with substantial individual variabilities in clinical progression and neurobiology. Amyloid deposition has been thought to drive cognitive decline and thus a major contributor to the variations in cognitive deterioration in AD. However, the clinical heterogeneity of patients with early symptomatic AD (mild cognitive impairment or mild dementia due to AD) already with evidence of amyloid abnormality in the brain is still unknown. Participants with a baseline diagnosis of mild cognitive impairment or mild dementia, a positive amyloid-PET scan, and more than one follow-up Alzheimer's Disease Assessment Scale-Cognitive Subscale-13 (ADAS-Cog-13) administration within a period of 5-year follow-up were selected from the Alzheimer's Disease Neuroimaging Initiative database (n = 421; age = 73±7; years of education = 16 ± 3; percentage of female gender = 43%; distribution of APOE4 carriers = 68%). A non-parametric k-means longitudinal clustering analysis in the context of the ADAS-Cog-13 data was performed to identify cognitive subtypes. We found a highly variable profile of cognitive decline among patients with early AD and identified 4 clusters characterized by distinct rates of cognitive progression. Among the groups there were significant differences in the magnitude of rates of changes in other cognitive and functional outcomes, clinical progression from mild cognitive impairment to dementia, and changes in markers presumed to reflect neurodegeneration and neuronal injury. A nomogram based on a simplified logistic regression model predicted steep cognitive trajectory with an AUC of 0.912 (95% CI: 0.88 - 0.94). Simulation of clinical trials suggested that the incorporation of the nomogram into enrichment strategies would reduce the required sample sizes from 926.8 (95% CI: 822.6 - 1057.5) to 400.9 (95% CI: 306.9 - 516.8). Our findings show usefulness in the stratification of patients in early AD and may thus increase the chances of finding a treatment for future AD clinical trials.

Streaming variational inference-empowered Bayesian nonparametric clustering for online structural damage detection with transmissibility function

Transmissibility function (TF) is widely applied in damage detection due to its sensitivity to damage and robustness to external excitations, but its application in online damage detection is rarely reported due to challenges in handling data streams. This study proposes a new TF-based online damage detection method that integrates a truncation-free variational inference-based full Dirichlet process Gaussian mixture model (VI-FDPGMM) within a streaming variational inference (SVI) paradigm. As an improved Bayesian nonparametric approach, the truncation-free VI-FDPGMM addresses the issue of truncating mixing components in traditional VI-DPGMM for online learning with increasing data by strategically setting the variational distributions of parameters for the components without assigned data (i.e., inactivated components) to their prior distributions based on the Bayesian viewpoint, which enables computing the probabilities to assign data points to these components and determining the creation of new components. As a result, the truncation-free VI-FDPGMM allows dynamically adding components to the mixture model, providing the flexibility to automatically adapt the number of components for arbitrary amounts of data. This characteristic enables its intuitive integration into the SVI paradigm featured as the variational posterior conditioned on the previous data as the prior when new data are observed, facilitating continuous refinement of the mixture model without repeatedly making inference of previous data. Therefore, the proposed method is highly efficient and well-suited for online damage detection. The proposed method is validated using two case studies, demonstrating its capability to dynamically generate new clusters as new data are available online to indicate the emergence of new damage patterns during the monitoring process, which enables it to perform structural anomaly detection tasks in a semi-supervised manner. Furthermore, the method outperforms some state-of-the-art methods due to its capability for continuous model refinement and robustness in interpreting and capturing uncertainties.

Understanding non‐stationarity patterns in basin‐scale hydroclimatic extremes

AbstractStationarity, a cornerstone in hydraulic design, is now under scrutiny due to anthropogenic activities and climate change. Numerous studies have sought to identify non‐stationarity (NS); however, a comprehensive assessment of time invariance in all statistical properties of a time series is less explored. This study presents a non‐overlapping block‐stratified random sampling (NBRS) framework leveraging the strengths of several nonparametric tests to assess NS. The NBRS approach exclusively detects NS and distinguishes between various forms of stationarity, including weak and strict. A variant of NBRS is proposed in this study to identify the underlying stochastic process(es) influencing NS in hydroclimatic extremes. Furthermore, a nonparametric clustering approach is used to unveil spatial clusters showcasing NS due to shifts in mean, variance, distribution of time series or a combination of these factors. A comparative assessment of the modified NBRS approach with traditional trend and change point methods is also performed. The proposed methodology is applied to assess the presence of NS in 28 hydroclimatic indices derived for the west‐central river basins of India, exhibiting diverse physio‐climatic settings, for the study period 1973–2021. The modified NBRS approach rigorously explores NS within extreme hydroclimatic indices, conclusively pinpointing its root causes and profound implications for hydrologic design. The applicability of the modified NBRS approach to gridded and point datasets is also demonstrated. The findings highlight the limitations of conventional trend and change point tests in capturing time‐invariant characteristics in heteroscedastic variables (such as streamflow and rainfall extremes) compared to the NBRS approach. The research reveals that NS in rainfall and streamflow extremes primarily results from distributional shifts, whilst temperature extremes are influenced by changes in mean and distribution properties. This research deepens our understanding of the evolving patterns in hydroclimatic extremes in a changing climate.

Efficient Feature Clustering for High-Dimensional Datasets: A Non-Parametric Approach

Clustering methods have become widely popular in the field of data analysis, as they enable the grouping of similar data points. However, the challenge of clustering high-dimensional data remains a significant obstacle due to the ”curse of dimensionality.” Traditional methods like Principal Component Analysis (PCA) and hierarchy-based feature agglomeration have been used to overcome this challenge by reducing the dimensionality of the dataset. However, there is currently no efficient method for clustering features of a high-dimensional dataset. To address this issue, we present a non-parametric clustering method in this paper, specifically designed for identifying features in a high-dimensional dataset. Our proposed feature clustering algorithm utilizes a partitioning algorithm that creates clusters of features. Through experimental results using publicly available datasets, we demonstrate the effectiveness of our method. Our approach offers a valuable tool for researchers and practitioners to better understand and analyze high-dimensional data.

WhARIO: whole-slide-image-based survival analysis for patients treated with immunotherapy.

Immune checkpoint inhibitors (ICIs) are now one of the standards of care for patients with lung cancer and have greatly improved both progression-free and overall survival, although of the patients respond to the treatment, and some face acute adverse events. Although a few predictive biomarkers have integrated the clinical workflow, they require additional modalities on top of whole-slide images and lack efficiency or robustness. In this work, we propose a biomarker of immunotherapy outcome derived solely from the analysis of histology slides. We develop a three-step framework, combining contrastive learning and nonparametric clustering to distinguish tissue patterns within the slides, before exploiting the adjacencies of previously defined regions to derive features and train a proportional hazards model for survival analysis. We test our approach on an in-house dataset of 193 patients from 5 medical centers and compare it with the gold standard tumor proportion score (TPS) biomarker. On a fivefold cross-validation (CV) of the entire dataset, the whole-slide image-based survival analysis for patients treated with immunotherapy (WhARIO) features are able to separate a low- and a high-risk group of patients with a hazard ratio (HR) of 2.29 ( to 3.56), whereas the TPS 1% reference threshold only reaches a HR of 1.81 ( to 2.69). Combining the two yields a higher HR of 2.60 ( to 3.94). Additional experiments on the same dataset, where one out of five centers is excluded from the CV and used as a test set, confirm these trends. Our uniquely designed WhARIO features are an efficient predictor of survival for lung cancer patients who received ICI treatment. We achieve similar performance to the current gold standard biomarker, without the need to access other imaging modalities, and show that both can be used together to reach even better results.

A Prototype-Based Neural Network for Image Anomaly Detection and Localization

Image anomaly detection and localization perform not only image-level anomaly classification but also locate pixel-level anomaly regions. Recently, it has received much research attention due to its wide application in various fields. This paper proposes ProtoAD, a prototype-based neural network for image anomaly detection and localization. First, the patch features of normal images are extracted by a deep network pre-trained on nature images. Then, the prototypes of the normal patch features are learned by non-parametric clustering. Finally, we construct an image anomaly localization network (ProtoAD) by appending the feature extraction network with L2 feature normalization, a 1×1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1\ imes 1$$\\end{document} convolutional layer, a channel max-pooling, and a subtraction operation. We use the prototypes as the kernels of the 1×1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1\ imes 1$$\\end{document} convolutional layer; therefore, our neural network does not need a training phase and can conduct anomaly detection and localization in an end-to-end manner. Extensive experiments on two challenging industrial anomaly detection datasets, MVTec AD and BTAD, demonstrate that ProtoAD achieves competitive performance compared to the state-of-the-art methods with a higher inference speed. The code and pre-trained models are publicly available at https://github.com/98chao/ProtoAD.

Spectral decomposition of ground motions in New Zealand using the generalized inversion technique

SUMMARY To gain new insights into ground-motion phenomena in New Zealand (NZ), we apply the non-parametric generalized inversion technique (GIT) in the Fourier domain to isolate the systematic source, path, and site effects from 20 813 seismograms, recorded by 693 sensors at 439 unique locations, from 1200 shallow crustal events (Mw > 3) during the period 2000–2021. From the inverted source spectra, we derive Brune's stress parameter, ∆σ, which is found to follow a lognormal distribution with a log10 standard deviation of 0.36 or equivalently 0.83 in natural log unit. ∆σ slightly increases with focal depth and is practically independent of earthquake size (i.e. self-similar), but displays a statistically significant spatial clustering. Based on the inverted attenuation, a trilinear geometric-spreading function, and a distance-dependent quality-factor Q(f) model are found to well describe the attenuation in NZ; though a single $Q( f )$ model is also obtained for the whole distance range: $Q( f ) = 149.1{f}^{0.62}$. Using the site response decomposed from GIT, we find that the soil classification scheme specified in NZ seismic code, NZS1170.5, has a limited capability in discerning the site-specific frequency-dependent amplification functions in comparison to a non-parametric clustering with the same number of discrete classes. The potential use of the spatial variation in source parameters from this GIT analysis in region-specific physics-based simulations is discussed.

Bergen metrics: composite error metrics for assessing performance of climate models using EURO-CORDEX simulations

Abstract. Error metrics are useful for evaluating model performance and have been used extensively in climate change studies. Despite the abundance of error metrics in the literature, most studies use only one or two metrics. Since each metric evaluates a specific aspect of the relationship between the reference data and model data, restricting the comparison to just one or two metrics limits the range of insights derived from the analysis. This study proposes a new framework and composite error metrics called Bergen metrics to summarize the overall performance of climate models and to ease interpretation of results from multiple error metrics. The framework of Bergen metrics are based on the p norm, and the first norm is selected to evaluate the climate models. The framework includes the application of a non-parametric clustering technique to multiple error metrics to reduce the number of error metrics with minimum information loss. An example of Bergen metrics is provided through its application to the large ensemble of regional climate simulations available from the EURO-CORDEX initiative. This study calculates 38 different error metrics to assess the performance of 89 regional climate simulations of precipitation and temperature over Europe. The non-parametric clustering technique is applied to these 38 metrics to reduce the number of metrics to be used in Bergen metrics for eight different sub-regions in Europe. These provide useful information about the performance of the error metrics in different regions. Results show it is possible to observe contradictory behaviour among error metrics when examining a single model. Therefore, the study also underscores the significance of employing multiple error metrics depending on the specific use case to achieve a thorough understanding of the model behaviour.

Spatiotemporal Variation and Causes of Typical Extreme Precipitation Events in Shandong Province over the Last 50 Years

In this study, based on hourly ERA5 reanalysis data from July to September, from 1971 to 2020, for Shandong Province, we used mathematical statistical analysis, the Mann–Kendall nonparametric statistical test, cluster analysis, and other methods to extract and analyze the spatiotemporal variation characteristics and causes of typical extreme precipitation events. The results indicated the following: (1) The total number and duration of precipitation events show a nonsignificant upward trend, while the average and extreme rainfall intensities show a nonsignificant downward trend. (2) Extreme precipitation events are primarily concentrated in Qingdao, Jinan, Heze, and Binzhou, with fewer events occurring in central Shandong Province. (3) Extreme precipitation events are classified into four types (namely, patterns I, II, III, and IV). Pattern I exhibits two rain peaks, with the primary rain peak occurring after the secondary rain peak. Similarly, pattern II also displays two rain peaks, with equivalent rainfall amounts for both peaks. In contrast, pattern III has multiple, evenly distributed rain peaks. Finally, pattern IV shows a rain peak during the first half of the precipitation event. Pattern I has the highest occurrence probability (46%), while pattern IV has the lowest (7%). (4) The spatial distributions of the different rain patterns are similar, with most being found in the eastern coastal and western regions. (5) Extreme precipitation events result from interactions between large-scale circulation configurations and mesoscale convective systems. The strong blocking situation and significant circulation transport at middle and low latitudes in East Asia, along with strong convergent uplift, abnormally high specific humidity, and high-water-vapor convergence centers, play crucial roles in supporting large-scale circulation systems and triggering mesoscale convective systems.

Persistence-based clustering with outlier-removing filtration

This article describes a non-parametric clustering algorithm with an outlier removal step. Our method is based on tools from topological data analysis: we define a new filtration on metric spaces which is a variant of the Vietoris–Rips filtration that adds information about the points' nearest neighbor to the persistence diagram. We prove a stability theorem for this filtration, and evaluate our method on point cloud and graph datasets, showing that it can compete with state-of-the-art methods while being non-parametric.

A Non-parametric Graph Clustering Framework for Multi-View Data

Multi-view graph clustering (MVGC) derives encouraging grouping results by seamlessly integrating abundant information inside heterogeneous data, and has captured surging focus recently. Nevertheless, the majority of current MVGC works involve at least one hyper-parameter, which not only requires additional efforts for tuning, but also leads to a complicated solving procedure, largely harming the flexibility and scalability of corresponding algorithms. To this end, in the article we are devoted to getting rid of hyper-parameters, and devise a non-parametric graph clustering (NpGC) framework to more practically partition multi-view data. To be specific, we hold that hyper-parameters play a role in balancing error item and regularization item so as to form high-quality clustering representations. Therefore, under without the assistance of hyper-parameters, how to acquire high-quality representations becomes the key. Inspired by this, we adopt two types of anchors, view-related and view-unrelated, to concurrently mine exclusive characteristics and common characteristics among views. Then, all anchors' information is gathered together via a consensus bipartite graph. By such ways, NpGC extracts both complementary and consistent multi-view features, thereby obtaining superior clustering results. Also, linear complexities enable it to handle datasets with over 120000 samples. Numerous experiments reveal NpGC's strong points compared to lots of classical approaches.

Nonparametric clustering of discrete probability distributions with generalized Shannon’s entropy and heatmap

This article proposes a non-parametric clustering heatmap for non-ordinal samples based on their underlying discrete probability distributions. The proposed heatmap is a preliminary tool to extract immediate information from non-ordinal data without making any parametric assumption.

Functional Data Clustering Method Based on Shape Information and Functional Mahalanobis Distance

Function-type data can provide insights into the internal structure of the data and facilitate the extraction of data features from the perspective of interactive derivative functions. The paper proposes a nonparametric clustering method for function-based data that incorporates first-order and second-order derivative function information into the Marginal distance for clustering function-based data. The method is based on the traditional K-means algorithm and is designed to cluster function-based data that contains rich shape information. The paper evaluated the performance of the algorithm of the text by comparing its purity and adjusted Rand coefficients against six other clustering algorithms on three different datasets. The results show that the algorithm of the text outperforms the other algorithms, demonstrating its outstanding performance, wide applicability, and practical significance in solving real-world problems.

Nonparametric Bayesian functional clustering with applications to racial disparities in breast cancer

AbstractAs we have easier access to massive data sets, functional analyses have gained more interest. However, such data sets often contain large heterogeneities, noises, and dimensionalities. When generalizing the analyses from vectors to functions, classical methods might not work directly. This paper considers noisy information reduction in functional analyses from two perspectives: functional clustering to group similar observations and thus reduce the sample size and functional variable selection to reduce the dimensionality. The complicated data structures and relations can be easily modeled by a Bayesian hierarchical model due to its flexibility. Hence, this paper proposes a nonparametric Bayesian functional clustering and peak point selection method via weighted Dirichlet process mixture (WDPM) modeling that automatically clusters and provides accurate estimations, together with conditional Laplace prior, which is a conjugate variable selection prior. The proposed method is named WDPM‐VS for short, and is able to simultaneously perform the following tasks: (1) Automatic cluster without specifying the number of clusters or cluster centers beforehand; (2) Cluster for heterogeneously behaved functions; (3) Select vibrational peak points; and (4) Reduce noisy information from the two perspectives: sample size and dimensionality. The method will greatly outperform its comparison methods in root mean squared errors. Based on this proposed method, we are able to identify biological factors that can explain the breast cancer racial disparities.

Facial feedback manipulation influences the automatic detection of unexpected emotional body expressions

Unexpected or changing facial expressions are known to be able to engage more automatic processing than frequently occurring facial expressions, thereby inducing a neural differential wave response known as expression mismatch negativity (EMMN). Recent studies have shown that EMMN can be modulated by the observer's facial feedback (i.e., feedback from their own facial movements). A similar EMMN activity has been discovered for body expressions, but thus far only a few emotion types have been investigated. It is unknown whether the EMMNs evoked by body expressions can be influenced by facial feedback. To explore this question, we recorded EEG activity of 29 participants in the reverse oddball paradigm. Here two unexamined categories of body expressions were presented, happy and sad, placed in two paired stimulus sequences: in one the happy body was presented with a probability of 80% (standards) while the sad body was presented with a probability of 20% (deviants), and in the other the probabilities were reversed. The facial feedback was manipulated by different pen holding conditions (i.e., participants holding the pen with the teeth, lips, or nondominant hand). The nonparametric cluster permutation test revealed significant happy and sad body-related EMMN (bEMMN) activities. The happy-bEMMN were more negative than sad-bEMMN within the range of 100–150 ms. Additionally, the bEMMN amplitude of both emotions is modulated by the facial feedback conditions. These results expand the range of emotional types applicable to bEMMN and provide evidence for the validity of the facial feedback hypothesis across emotional carriers.

Proteome trait regulation of marine Synechococcus elemental stoichiometry under global change.

Recent studies have demonstrated regional differences in marine ecosystem C:N:P with implications for carbon and nutrient cycles. Due to strong co-variance, temperature and nutrient stress explain variability in C:N:P equally well. A reductionistic approach can link changes in individual environmental drivers with changes in biochemical traits and cell C:N:P. Thus, we quantified effects of temperature and nutrient stress on Synechococcus chemistry using laboratory chemostats, chemical analyses, and data-independent acquisition mass spectrometry proteomics. Nutrient supply accounted for most C:N:Pcell variability and induced tradeoffs between nutrient acquisition and ribosomal proteins. High temperature prompted heat-shock, whereas thermal effects via the "translation-compensation hypothesis" were only seen under P-stress. A Nonparametric Bayesian Local Clustering algorithm suggested that changes in lipopolysaccharides, peptidoglycans, and C-rich compatible solutes may also contribute to C:N:P regulation. Physiological responses match field-based trends in ecosystem stoichiometry and suggest a hierarchical environmental regulation of current and future ocean C:N:P.