Search

Closed testing has recently been shown to be optimal for simultaneous true discovery proportion control. It is, however, challenging to construct true discovery guarantee procedures in such a way that it focuses power on some feature sets chosen by users based on their specific interest or expertise. We propose a procedure that allows users to target power on prespecified feature sets, that is, "focus sets." Still, the method also allows inference for feature sets chosen post hoc, that is, "nonfocus sets," for which we deduce a true discovery lower confidence bound by interpolation. Our procedure is built from partial true discovery guarantee procedures combined with Holm's procedure and is a conservative shortcut to the closed testing procedure. A simulation study confirms that the statistical power of our method is relatively high for focus sets, at the cost of power for nonfocus sets, as desired. In addition, we investigate its power property for sets with specific structures, for example, trees and directed acyclic graphs. We also compare our method with AdaFilter in the context of replicability analysis. The application of our method is illustrated with a gene ontology analysis in gene expression data. [ABSTRACT FROM AUTHOR]

Objective: Impaired amine metabolism has been associated with the etiology of migraine, that is, why patients continue to get migraine attacks. However, evidence from cerebrospinal fluid (CSF) is lacking. Here, we evaluated individual amine levels, global amine profiles, and amine pathways in CSF and plasma of interictal migraine patients and healthy controls. Methods: CSF and plasma were sampled between 8:30 am and 1:00 pm, randomly and interchangeably over the time span to avoid any diurnal and seasonal influences, from healthy volunteers and interictal migraine patients, matched for age, sex, and sampling time. The study was approved by the local medical ethics committee. Individual amines (n = 31), global amine profiles, and specific amine pathways were analyzed using a validated ultraperformance liquid chromatography mass spectrometry platform. Results: We analyzed n = 99 participants with migraine with aura, n = 98 with migraine without aura, and n = 96 healthy volunteers. Univariate analysis with Bonferroni correction indicated that CSF L‐arginine was reduced in migraine with aura (10.4%, p < 0.001) and without aura (5.0%, p = 0.03). False discovery rate‐corrected CSF L‐phenylalanine was also lower in migraine with aura (6.9%, p = 0.011) and without aura (8.1%, p = 0.001), p = 0.088 after Bonferroni correction. Multivariate analysis revealed that CSF global amine profiles were similar for both types of migraine (p = 0.64), but distinct from controls (p = 0.009). Global profile analyses were similar in plasma. The strongest associated pathways with migraine were related to L‐arginine metabolism. Interpretation: L‐Arginine was decreased in the CSF (but not in plasma) of interictal patients with migraine with or without aura, and associated pathways were altered. This suggests that dysfunction of nitric oxide signaling is involved in susceptibility to getting migraine attacks. ANN NEUROL 2023;93:715–728 [ABSTRACT FROM AUTHOR]

Purpose: To assess errors associated with EPI‐accelerated intracardiac 4D flow MRI (4DEPI) with EPI factor 5, compared with non‐EPI gradient echo (4DGRE). Methods: Three 3T MRI experiments were performed comparing 4DEPI to 4DGRE: steady flow through straight tubes, pulsatile flow in a left‐ventricle phantom, and intracardiac flow in 10 healthy volunteers. For each experiment, 4DEPI was repeated with readout and blip phase‐encoding gradient in different orientations, parallel or perpendicular to the flow direction. In vitro flow rates were compared with timed volumetric collection. In the left‐ventricle phantom and in vivo, voxel‐based speed and spatio‐temporal median speed were compared between sequences, as well as mitral and aortic transvalvular net forward volume. Results: In steady‐flow phantoms, the flow rate error was largest (12%) for high velocity (>2 m/s) with 4DEPI readout gradient parallel to the flow. Voxel‐based speed and median speed in the left‐ventricle phantom were ≤5.5% different between sequences. In vivo, mean net forward volume inconsistency was largest (6.4 ± 8.5%) for 4DEPI with nonblip phase‐encoding gradient parallel to the main flow. The difference in median speed for 4DEPI versus 4DGRE was largest (9%) when the 4DEPI readout gradient was parallel to the flow. Conclusions: Velocity and flow rate are inaccurate for 4DEPI with EPI factor 5 when flow is parallel to the readout or blip phase‐encoding gradient. However, mean differences in flow rate, voxel‐based speed, and spatio‐temporal median speed were acceptable (≤10%) when comparing 4DEPI to 4DGRE for intracardiac flow in healthy volunteers. [ABSTRACT FROM AUTHOR]

Summary: The statistical literature is known to be inconsistent in the use of the terms 'permutation test' and 'randomisation test'. Several authors successfully argue that these terms should be used to refer to two distinct classes of tests and that there are major conceptual differences between these classes. The present paper explains an important difference in mathematical reasoning between these classes: a permutation test fundamentally requires that the set of permutations has a group structure, in the algebraic sense; the reasoning behind a randomisation test is not based on such a group structure, and it is possible to use an experimental design that does not correspond to a group. In particular, we can use a randomisation scheme where the number of possible treatment patterns is larger than in standard experimental designs. This leads to exact p values of improved resolution, providing increased power for very small significance levels, at the cost of decreased power for larger significance levels. We discuss applications in randomised trials and elsewhere. Further, we explain that Fisher's famous Lady Tasting Tea experiment, which is commonly referred to as the first permutation test, is in fact a randomisation test. This distinction is important to avoid confusion and invalid tests. [ABSTRACT FROM AUTHOR]

We propose a top-down approach for pathway analysis of longitudinal metabolite data. We apply a score test based on a shared latent process mixed model which can identify pathways with differentially progressing metabolites. The strength of our approach is that it can handle unbalanced designs, deals with potential missing values in the longitudinal markers, and gives valid results even with small sample sizes. Contrary to bottom-up approaches, correlations between metabolites are explicitly modeled leveraging power gains. For large pathway sizes, a computationally efficient solution is proposed based on pseudo-likelihood methodology. We demonstrate the advantages of the proposed method in identification of differentially expressed pathways through simulation studies. Finally, longitudinal metabolite data from a mice experiment is analyzed to demonstrate our methodology. [ABSTRACT FROM AUTHOR]

It is well known that the maximum likelihood estimator (MLE) is inadmissible when estimating the multidimensional Gaussian location parameter. We show that the verdict is much more subtle for the binary location parameter. We consider this problem in a regression framework by considering a ridge logistic regression (RR) with three alternative ways of shrinking the estimates of the event probabilities. While it is shown that all three variants reduce the mean squared error (MSE) of the MLE, there is at the same time, for every amount of shrinkage, a true value of the location parameter for which we are overshrinking, thus implying the minimaxity of the MLE in this family of estimators. Little shrinkage also always reduces the MSE of individual predictions for all three RR estimators; however, only the naive estimator that shrinks toward 1/2 retains this property for any generalized MSE (GMSE). In contrast, for the two RR estimators that shrink toward the common mean probability, there is always a GMSE for which even a minute amount of shrinkage increases the error. These theoretical results are illustrated on a numerical example. The estimators are also applied to a real data set, and practical implications of our results are discussed. [ABSTRACT FROM AUTHOR]

Hommel's and Hochberg's procedures for familywise error control are both derived as shortcuts in a closed testing procedure with the Simes local test. Hommel's shortcut is exact but takes quadratic time in the number of hypotheses. Hochberg's shortcut takes only linear time after the P‐values are sorted, but is conservative. In this paper, we present an exact shortcut in linear time on sorted P‐values, combining the strengths of both procedures. The novel shortcut also applies to a robust variant of Hommel's procedure that does not require the assumption of the Simes inequality. [ABSTRACT FROM AUTHOR]

Abstract: Muscular dystrophies are characterized by a progressive loss of muscle tissue and/or muscle function. While metabolic alterations have been described in patients’‐derived muscle biopsies, non‐invasive readouts able to describe these alterations are needed in order to objectively monitor muscle condition and response to treatment targeting metabolic abnormalities. We used a metabolomic approach to study metabolites concentration in serum of patients affected by multiple forms of muscular dystrophy such as Duchenne and Becker muscular dystrophies, limb‐girdle muscular dystrophies type 2A and 2B, myotonic dystrophy type 1 and facioscapulohumeral muscular dystrophy. We show that 15 metabolites involved in energy production, amino acid metabolism, testosterone metabolism and response to treatment with glucocorticoids were differentially expressed between healthy controls and Duchenne patients. Five metabolites were also able to discriminate other forms of muscular dystrophy. In particular, creatinine and the creatine/creatinine ratio were significantly associated with Duchenne patients performance as assessed by the 6‐minute walk test and north star ambulatory assessment. The obtained results provide evidence that metabolomics analysis of serum samples can provide useful information regarding muscle condition and response to treatment, such as to glucocorticoids treatment. [ABSTRACT FROM AUTHOR]

Aims POLE exonuclease domain mutations identify a subset of endometrial cancer ( EC) patients with an excellent prognosis. The use of this biomarker has been suggested to refine adjuvant treatment decisions, but the necessary sequencing is not widely performed and is relatively expensive. Therefore, we aimed to identify histopathological and immunohistochemical characteristics to aid in the detection of POLE-mutant ECs. Methods and results Fifty-one POLE-mutant endometrioid, 67 POLE-wild-type endometrioid and 15 POLE-wild-type serous ECs were included (total N = 133). An expert gynaecopathologist, blinded to molecular features, evaluated each case (two or more slides) for 16 morphological characteristics. Immunohistochemistry was performed for p53, p16, MLH1, MSH2, MSH6, and PMS2. POLE-mutant ECs were characterised by a prominent immune infiltrate: 80% showed peritumoral lymphocytes and 59% showed tumour-infiltrating lymphocytes, as compared with 43% and 28% of POLE-wild-type endometrioid ECs, and 27% and 13% of their serous counterparts ( P < 0.01, all comparisons). Of POLE-mutant ECs, 33% contained tumour giant cells; this proportion was significantly higher than that in POLE-wild-type endometrioid ECs (10%; P = 0.003), but not significantly different from that in serous ECs (53%). Serous-like features were as often (focally) present in POLE-mutant as in POLE-wild-type endometrioid ECs (6-24%, depending on the feature). The majority of POLE-mutant ECs showed wild-type p53 (86%), negative/focal p16 (82%) and normal mismatch repair protein expression (90%). Conclusions A simple combination of morphological and immunohistochemical characteristics (tumour type, grade, peritumoral lymphocytes, MLH1, and p53 expression) can assist in prescreening for POLE exonuclease domain mutations in EC, increasing the probability of a mutation being detected from 7% to 33%. This facilitates the use of this important prognostic biomarker in routine pathology. [ABSTRACT FROM AUTHOR]

We define an adaptive procedure for control of the false discovery rate that is uniformly more powerful than the procedure of Benjamini and Hochberg. The power gain is tiny, however, and only appreciable for small numbers of hypotheses. We illustrate the new method with the case of two hypotheses, for which so far no procedure was known that controls false discovery rate but not also familywise error rate under positive dependence. [ABSTRACT FROM AUTHOR]

In high-dimensional omics studies where multiple molecular profiles are obtained for each set of patients, there is often interest in identifying complex multivariate associations, for example, copy number regulated expression levels in a certain pathway or in a genomic region. To detect such associations, we present a novel approach to test for association between two sets of variables. Our approach generalizes the global test, which tests for association between a group of covariates and a single univariate response, to allow high-dimensional multivariate response. We apply the method to several simulated datasets as well as two publicly available datasets, where we compare the performance of multivariate global test (G2) with univariate global test. The method is implemented in R and will be available as a part of the globaltest package in R. [ABSTRACT FROM AUTHOR]

We propose a Cochran-Armitage-type and a score-free global test that can be used to assess the presence of an association between a set of ordinally scaled covariates and an outcome variable within the range of generalized linear models. Both tests are developed within the framework of the well-established 'global test' methodology and as such are feasible in high-dimensional data situations under any correlation and enable adjustment for covariates. The Cochran-Armitage-type test, for which an intimate connection with the traditional score-based Cochran-Armitage test is shown, rests upon explicit assumptions on the distances between the covariates' ordered categories. The score-free test, in contrast, parametrizes these distances and thus keeps them flexible, rendering it ideally suited for covariates measured on an ordinal scale. As confirmed by means of simulations, the Cochran-Armitage-type test focuses its power on set-outcome relationships where the distances between the covariates' categories are equal or close to those assumed, whereas the score-free test spreads its power over a wide range of possible set-outcome relationships, putting more emphasis on monotonic than on non-monotonic ones. Based on the tests' power properties, it is discussed when to favour one or the other, and the practical merits of both of them are illustrated by an application in the field of rehabilitation medicine. Our proposed tests are implemented in the R package globaltest. Copyright © 2016 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

We present a novel multiple testing method for testing null hypotheses that are structured in a directed acyclic graph (DAG). The method is a top-down method that strongly controls the familywise error rate and can be seen as a generalization of Meinshausen's procedure for tree-structured hypotheses. Just as Meinshausen's procedure, our proposed method can be used to test for variable importance, only the corresponding variable clusters can be chosen more freely, because the method allows for multiple parent nodes and partially overlapping hypotheses. An important application of our method is in gene set analysis, in which one often wants to test multiple gene sets as well as individual genes for their association with a clinical outcome. By considering the genes and gene sets as nodes in a DAG, our method enables us to test both for significant gene sets as well as for significant individual genes within the same multiple testing procedure. The method will be illustrated by testing Gene Ontology terms for evidence of differential expression in a survival setting and is implemented in the [ABSTRACT FROM AUTHOR]

This paper presents an overview of the current state of the art in multiple testing in genomics data from a user's perspective. We describe methods for familywise error control, false discovery rate control and false discovery proportion estimation and confidence, both conceptually and practically, and explain when to use which type of error rate. We elaborate on the assumptions underlying the methods and discuss pitfalls in the interpretation of results. In our discussion, we take into account the exploratory nature of genomics experiments, looking at selection of genes before or after testing, and at the role of validation experiments. Copyright © 2014 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

This paper presents an efficient algorithm based on the combination of Newton Raphson and Gradient Ascent, for using the fused lasso regression method to construct a genome-based classifier. The characteristic structure of copy number data suggests that feature selection should take genomic location into account for producing more interpretable results for genome-based classifiers. The fused lasso penalty, an extension of the lasso penalty, encourages sparsity of the coefficients and their differences by penalizing the L1-norm for both of them at the same time, thus using genomic location. The major advantage of the algorithm over other existing fused lasso optimization techniques is its ability to predict binomial as well as survival response efficiently. We apply our algorithm to two publicly available datasets in order to predict survival and binary outcomes. [ABSTRACT FROM AUTHOR]

The bodily decline that occurs with advancing age strongly impacts on the prospects for future health and life expectancy. Despite the profound role of age in disease etiology, knowledge about the molecular mechanisms driving the process of aging in humans is limited. Here, we used an integrative network-based approach for combining multiple large-scale expression studies in blood (2539 individuals) with protein-protein Interaction ( PPI) data for the detection of consistently coexpressed PPI modules that may reflect key processes that change throughout the course of normative aging. Module detection followed by a meta-analysis on chronological age identified fifteen consistently coexpressed PPI modules associated with chronological age, including a highly significant module ( P = 3.5 × 10−38) enriched for ' T-cell activation' marking age-associated shifts in lymphocyte blood cell counts ( R2 = 0.603; P = 1.9 × 10−10). Adjusting the analysis in the compendium for the ' T-cell activation' module showed five consistently coexpressed PPI modules that robustly associated with chronological age and included modules enriched for ' Translational elongation', ' Cytolysis' and ' DNA metabolic process'. In an independent study of 3535 individuals, four of five modules consistently associated with chronological age, underpinning the robustness of the approach. We found three of five modules to be significantly enriched with aging-related genes, as defined by the GenAge database, and association with prospective survival at high ages for one of the modules including ASF1A. The hereby-detected age-associated and consistently coexpressed PPI modules therefore may provide a molecular basis for future research into mechanisms underlying human aging. [ABSTRACT FROM AUTHOR]

The involvement of epigenetic alterations in the pathogenesis of melanoma is increasingly recognized. Here, we performed genome-wide DNA methylation analysis of primary cutaneous melanoma and benign melanocytic nevus interrogating 14 495 genes using BeadChip technology. This genome-wide view of promoter methylation in primary cutaneous melanoma revealed an array of recurrent DNA methylation alterations with potential diagnostic applications. Among 106 frequently hypermethylated genes, there were many novel methylation targets and tumor suppressor genes. Highly recurrent methylation of the HOXA9, MAPK13, CDH11, PLEKHG6, PPP1R3C, and CLDN11 genes was established. Promoter methylation of MAPK13, encoding p38δ, was present in 67% of primary and 85% of metastatic melanomas. Restoration of MAPK13 expression in melanoma cells exhibiting epigenetic silencing of this gene reduced proliferation, indicative of tumor suppressive functions. This study demonstrates that DNA methylation alterations are widespread in melanoma and suggests that epigenetic silencing of MAPK13 contributes to melanoma progression. [ABSTRACT FROM AUTHOR]

In model building and model evaluation, cross-validation is a frequently used resampling method. Unfortunately, this method can be quite time consuming. In this article, we discuss an approximation method that is much faster and can be used in generalized linear models and Cox' proportional hazards model with a ridge penalty term. Our approximation method is based on a Taylor expansion around the estimate of the full model. In this way, all cross-validated estimates are approximated without refitting the model. The tuning parameter can now be chosen based on these approximations and can be optimized in less time. The method is most accurate when approximating leave-one-out cross-validation results for large data sets which is originally the most computationally demanding situation. In order to demonstrate the method's performance, it will be applied to several microarray data sets. An R package penalized, which implements the method, is available on CRAN. [ABSTRACT FROM AUTHOR]

When fitting generalized linear models or the Cox proportional hazards model, it is important to have tools to test for lack of fit. Because lack of fit comes in all shapes and sizes, distinguishing among different types of lack of fit is of practical importance. We argue that an adequate diagnosis of lack of fit requires a specified alternative model. Such specification identifies the type of lack of fit the test is directed against so that if we reject the null hypothesis, we know the direction of the departure from the model. The goodness-of-fit approach of this paper allows to treat different types of lack of fit within a unified general framework and to consider many existing tests as special cases. Connections with penalized likelihood and random effects are discussed, and the application of the proposed approach is illustrated with medical examples. Tailored functions for goodness-of-fit testing have been implemented in the R package globaltest. Copyright © 2012 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

We propose three-sided testing, a testing framework for simultaneous testing of inferiority, equivalence and superiority in clinical trials, controlling for multiple testing using the partitioning principle. Like the usual two-sided testing approach, this approach is completely symmetric in the two treatments compared. Still, because the hypotheses of inferiority and superiority are tested with one-sided tests, the proposed approach has more power than the two-sided approach to infer non-inferiority or non-superiority. Applied to the classical point null hypothesis of equivalence, the three-sided testing approach shows that it is sometimes possible to make an inference on the sign of the parameter of interest, even when the null hypothesis itself could not be rejected. Relationships with confidence intervals are explored, and the effectiveness of the three-sided testing approach is demonstrated in a number of recent clinical trials. Copyright © 2010 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

The bodily decline that occurs with advancing age strongly impacts on the prospects for future health and life expectancy. Despite the profound role of age in disease etiology, knowledge about the molecular mechanisms driving the process of aging in humans is limited. Here, we used an integrative network-based approach for combining multiple large-scale expression studies in blood (2539 individuals) with protein-protein Interaction (PPI) data for the detection of consistently coexpressed PPI modules that may reflect key processes that change throughout the course of normative aging. Module detection followed by a meta-analysis on chronological age identified fifteen consistently coexpressed PPI modules associated with chronological age, including a highly significant module (P = 3.5 × 10(-38)) enriched for 'T-cell activation' marking age-associated shifts in lymphocyte blood cell counts (R(2) = 0.603; P = 1.9 × 10(-10)). Adjusting the analysis in the compendium for the 'T-cell activation' module showed five consistently coexpressed PPI modules that robustly associated with chronological age and included modules enriched for 'Translational elongation', 'Cytolysis' and 'DNA metabolic process'. In an independent study of 3535 individuals, four of five modules consistently associated with chronological age, underpinning the robustness of the approach. We found three of five modules to be significantly enriched with aging-related genes, as defined by the GenAge database, and association with prospective survival at high ages for one of the modules including ASF1A. The hereby-detected age-associated and consistently coexpressed PPI modules therefore may provide a molecular basis for future research into mechanisms underlying human aging., (© 2013 The Authors. Aging Cell published by the Anatomical Society and John Wiley & Sons Ltd.)