Dataset Information

Optimizing experimental design for comparing models of brain function.

ABSTRACT: This article presents the first attempt to formalize the optimization of experimental design with the aim of comparing models of brain function based on neuroimaging data. We demonstrate our approach in the context of Dynamic Causal Modelling (DCM), which relates experimental manipulations to observed network dynamics (via hidden neuronal states) and provides an inference framework for selecting among candidate models. Here, we show how to optimize the sensitivity of model selection by choosing among experimental designs according to their respective model selection accuracy. Using Bayesian decision theory, we (i) derive the Laplace-Chernoff risk for model selection, (ii) disclose its relationship with classical design optimality criteria and (iii) assess its sensitivity to basic modelling assumptions. We then evaluate the approach when identifying brain networks using DCM. Monte-Carlo simulations and empirical analyses of fMRI data from a simple bimanual motor task in humans serve to demonstrate the relationship between network identification and the optimal experimental design. For example, we show that deciding whether there is a feedback connection requires shorter epoch durations, relative to asking whether there is experimentally induced change in a connection that is known to be present. Finally, we discuss limitations and potential extensions of this work.

SUBMITTER: Daunizeau J

PROVIDER: S-EPMC3219623 | biostudies-literature | 2011 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Optimizing experimental design for comparing models of brain function.

Daunizeau Jean J Preuschoff Kerstin K Friston Karl K Stephan Klaas K

PLoS computational biology 20111117 11

This article presents the first attempt to formalize the optimization of experimental design with the aim of comparing models of brain function based on neuroimaging data. We demonstrate our approach in the context of Dynamic Causal Modelling (DCM), which relates experimental manipulations to observed network dynamics (via hidden neuronal states) and provides an inference framework for selecting among candidate models. Here, we show how to optimize the sensitivity of model selection by choosing ...[more]

PMID: 22125485

Similar Datasets

Project description:MotivationNumerous sequencing studies, including transcriptomics of host-pathogen systems, sequencing of hybrid genomes, xenografts, mixed species systems, metagenomics and meta-transcriptomics, involve samples containing genetic material from divergent organisms. A crucial step in these studies is identifying from which organism each sequencing read originated, and the experimental design should be directed to minimize biases caused by cross-mapping of reads to incorrect source genomes. Additionally, pooling of sufficiently different genetic material into a single sequencing library could significantly reduce experimental costs but requires careful planning and assessment of the impact of cross-mapping. Having these applications in mind we designed Crossmapper, the first to our knowledge tool able to assess cross-mapping prior to sequencing, therefore allowing optimization of experimental design.ResultsUsing any combination of reference genomes, Crossmapper performs read simulation and back-mapping of those reads to the pool of references, quantifies and reports the cross-mapping rates for each organism. Crossmapper performs these analyses with numerous user-specified parameters, including, among others, read length, read layout, coverage, mapping parameters, genomic or transcriptomic data. Additionally, it outputs the results in highly interactive and publication-ready reports. This allows the user to perform multiple comparisons at once and choose the experimental setup minimizing cross-mapping rates. Moreover, Crossmapper can be used for resource optimization in sequencing facilities by pooling different samples into one sequencing library.Availability and implementationCrossmapper is a command line tool implemented in Python 3.6 and available as a conda package, allowing effortless installation. The source code, detailed information and a step-by-step tutorial is available at our GitHub page https://github.com/Gabaldonlab/crossmapper.Supplementary informationSupplementary data are available at Bioinformatics online.

Project description:ObjectiveDeep brain stimulation is a treatment for medically refractory essential tremor. To improve the therapy, closed-loop approaches are designed to deliver stimulation according to the system's state, which is constantly monitored by recording a pathological signal associated with symptoms (e.g. brain signal or limb tremor). Since the space of possible closed-loop stimulation strategies is vast and cannot be fully explored experimentally, how to stimulate according to the state should be informed by modeling. A typical modeling goal is to design a stimulation strategy that aims to maximally reduce the Hilbert amplitude of the pathological signal in order to minimize symptoms. Isostables provide a notion of amplitude related to convergence time to the attractor, which can be beneficial in model-based control problems. However, how isostable and Hilbert amplitudes compare when optimizing the amplitude response to stimulation in models constrained by data is unknown.ApproachWe formulate a simple closed-loop stimulation strategy based on models previously fitted to phase-locked deep brain stimulation data from essential tremor patients. We compare the performance of this strategy in suppressing oscillatory power when based on Hilbert amplitude and when based on isostable amplitude. We also compare performance to phase-locked stimulation and open-loop high-frequency stimulation.Main resultsFor our closed-loop phase space stimulation strategy, stimulation based on isostable amplitude is significantly more effective than stimulation based on Hilbert amplitude when amplitude field computation time is limited to minutes. Performance is similar when there are no constraints, however constraints on computation time are expected in clinical applications. Even when computation time is limited to minutes, closed-loop phase space stimulation based on isostable amplitude is advantageous compared to phase-locked stimulation, and is more efficient than high-frequency stimulation.SignificanceOur results suggest a potential benefit to using isostable amplitude more broadly for model-based optimization of stimulation in neurological disorders.

Project description:Progress continues in the field of cancer biology, yet much remains to be unveiled regarding the mechanisms of cancer invasion. In particular, complex biophysical mechanisms enable a tumor to remodel the surrounding extracellular matrix (ECM), allowing cells to invade alone or collectively. Tumor spheroids cultured in collagen represent a simplified, reproducible 3D model system, which is sufficiently complex to recapitulate the evolving organization of cells and interaction with the ECM that occur during invasion. Recent experimental approaches enable high resolution imaging and quantification of the internal structure of invading tumor spheroids. Concurrently, computational modeling enables simulations of complex multicellular aggregates based on first principles. The comparison between real and simulated spheroids represents a way to fully exploit both data sources, but remains a challenge. We hypothesize that comparing any two spheroids requires first the extraction of basic features from the raw data, and second the definition of key metrics to match such features. Here, we present a novel method to compare spatial features of spheroids in 3D. To do so, we define and extract features from spheroid point cloud data, which we simulated using Cells in Silico (CiS), a high-performance framework for large-scale tissue modeling previously developed by us. We then define metrics to compare features between individual spheroids, and combine all metrics into an overall deviation score. Finally, we use our features to compare experimental data on invading spheroids in increasing collagen densities. We propose that our approach represents the basis for defining improved metrics to compare large 3D data sets. Moving forward, this approach will enable the detailed analysis of spheroids of any origin, one application of which is informing in silico spheroids based on their in vitro counterparts. This will enable both basic and applied researchers to close the loop between modeling and experiments in cancer research.

Dataset Information

Optimizing experimental design for comparing models of brain function.

Publications

Optimizing experimental design for comparing models of brain function.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets