Program on Mathematical and Statistical Methods for Climate and the Earth System (CLIM)
|Working Group||Group Leaders|
||Working Group Leaders: Amy Braverman (JPL); Jessica Matthews (NCEI/CICS)|
||Working Group Leader: Peter Challenor, (Exeter)|
||Working Group Leaders: Chris Jones (UNC-Chapel Hill); Erik Van Vleck (Kansas)|
||Working Group Leaders: Doug Nychka (NCAR); Vipin Kumar (Minnesota)|
||Working Group Leader: Leonard Smith (London School of Economics)|
||Working Group Leaders: Dan Cooley (Colorado State); Richard Smith (UNC-Chapel Hill)|
||Working Group Leader: Adam Monahan (University of Victoria, Canada)|
||Working Group Leader: Brian Reich (NCSU)|
||Working Group Leaders: Hans Kaper (Georgetown); Mary Lou Zeeman (Bowdoin College)|
||Working Group Leaders: Chris Jones (UNC-Chapel Hill); Alberto Carrassi (NERSC) and Murali Haran (Penn State)|
|Detection and Attribution
||Working Group Leaders: Dorit Hammerling (NCAR) and Matthias Katzfuss (Texas A&M)|
|Risk & Coastal Hazards
||Working Group Leaders: Brian Blanton (RENCI); Slava Lyubchich (Maryland) and Richard Smith (UNC-Chapel Hill)|
||Working Group Leaders: Michael Stein (University of Chicago) and Mikael Kuusela (Carnegie Mellon University)|
|Course: Fall 2017: Statistics for Climate Research||Description: This course is being offered in conjunction with the SAMSI year-long research program on Mathematical and Statistical Methods for Climate and the Earth System. The course will cover statistical and computational methods for the analysis of data arising in climate research. Specific topics will include:
Schedule: August 29 – December 5, 2017 / Tuesdays at 4:30PM – 7:00PM
|Course: Spring 2018: Data Assimilation in Dynamical Systems||Description: This course is being offered in conjunction with the SAMSI year-long research program on Mathematical and Statistical Methods for Climate and the Earth System. The course will cover statistical and computational methods for the analysis of data arising in climate research. Specific topics will include:
Schedule: January 16 – April 24, 2018 / Tuesdays at 4:30PM – 7:00PM
Program on Quasi-Monte Carlo and High-Dimensional Sampling Methods for Applied Mathematics (QMC)
|QMC Working Group I: Parallel Monte Carlo Neutronics Simulation||Working Group Leaders: C.T. Kelley (NCSU); T.M. Evans (Oak Ridge National Laboratory) and S.P. Hamilton (Oak Ridge National Laboratory)
Description: We simulate neutron flux intensity far from the source or behind shielding. Applications include nuclear power plants where we model radiation in machine rooms near the core, but behind shielding. Monte Carlo simulation randomly generates particles using information on sources and scattering properties of the problem domain and then follows those particles until they either leave the domain or are absorbed. This is not a PDE-based approach and has no discretization error.
Fluxes are very small and standard MC would resolve the radiation field inefficiently. We address this problem with CADIS (Consistent Adjoint-Driven Importance Sampling) which eliminates particles that move away from a region of interest. The computation must be done in parallel and the research focus of this project is to optimize the domain decomposition that is central to parallelization.
|QMC Working Group II: Probabilistic Numbers||Working Group Leaders: Chris Oates (University of Newcastle Upon Tyne and Alan Turing Institute, UK) and Tim Sullivan (Free University of Berlin, Germany)
Description: The accuracy and robustness of numerical predictions that are based on mathematical models depend critically upon the construction of accurate discrete approximations to key quantities of interest. The exact error due to approximation will be unknown to the analyst, but worst-case upper bounds can often be obtained. This working group aims, instead, to develop Probabilistic Numerical Methods, which provide the analyst with a richer, probabilistic quantification of the numerical error in their output, thus providing better tools for reliable statistical inference.
|QMC Working Group III: Sampling by Interacting Particle Systems||Working Group Leaders: Jianfeng Lu (Duke University) and Jonathan Mattingly (Duke University)
Description: In recent years, several algorithms, based on interacting samplers, have been proposed and further advanced; such algorithms have been used in various settings, including high dimensional Bayesian sampling, particle filtering, computational physics, etc. This working group plans to explore the connections between various ideas proposed in different communities and try to advance the theoretical understanding of efficiency and design principle for such sampling algorithms.
|QMC Working Group IV: Representative Points for Small-data and Big-data Problems||Working Group Leaders: Roshan Vengazhiyil (Georgia Institute of Technology) and Simon Mak (Georgia Institute of Technology)
Description: Representative points, which compact a probability distribution into a finite point set, are useful for a wide array of “small-data” and “big-data” problems. Small-data problems arise naturally in many engineering applications, where a key challenge is to allocate limited experimental runs for performing functional approximation, uncertainty propagation or design optimization. Similarly, given the massive volume, variety and velocity of big-data (particularly in Bayesian problems), the reduction of such datasets using representative points allows for meaningful and timely analysis. This working group aims to investigate the theory and application of representative points to the aforementioned small-data and big-data problems, with an emphasis on engineering applications and Bayesian computation.
|QMC Working Group V: Sampling and Analysis in High Dimensions When Samples Are Expensive||Working Group Leaders: Fred J. Hickernell (Illinois Institute of Technology); Mac Hyman (Tulane University) and Paul Constantine (University of Colorado)
Description: For some problems, generating observations (multivariate function values) is expensive, such as running a time-consuming computer code or conducting a real-world experiment. This can result in having only a few samples for high-dimensional input variables. These samples can be noisy and may not be available for some important ranges of the parameters. Examples include medical trials, large-scale computer simulations, such as climate models, and financial market data. This working group will study how to best use the available small samples and how strategically guide picking future samples. The group will consider algorithms based on assumptions about the underlying input-output relationship to improve the effectiveness of the statistical analysis under these difficult constraints. The assumption will be that the cost of the observations will far exceed the computational cost of the post-processing algorithms.
|Working Group VI: Adaptive Choice of Sobolev Space Weights||Working Group Leader: Art Owen (Stanford University)
Description: High dimensional numerical integration problems can often be solve efficiently using algorithms devised for weighted Sobolev spaces. A given black box function may belong to all of the spaces people use. That raises the question of which weights to be used. This WG looks at methods to select weights using ideas from global sensitivity analysis and active subspaces.
|Working Group VII: Multivariate Decomposition Method (MDM) and Applications||Working Group Leaders: Dirk Nuyens (KU Leuven – Belgium) and Alec Gilbert (University of New South Wales Sydney, Australia)
Description: The multivariate decomposition method allows us to approximate infinite-dimensional problems by a sum of finite-dimensional subproblems in the “active set” whose size depends on the requested error demand. The weight structure imposed on these finite-dimensional subspaces often implies that the maximum number of variables in each subspace grows very slowly in terms of the requested error demand. Thus we only need to solve many low-dimensional problems (say, 1, 2, 3 dimensional) and a few medium-dimensional problems (say, up to 10 dimensions) to reach a reasonable error demand. These problems can be solved in parallel. The working group will develop efficient implementation of the MDM and consider its application to PDE problems with random coefficients as well as other potential applications.
|Working Group VIII: Application of QMC to PDEs with random coefficients||Working Group Leaders: Frances Kuo (University of New South Wales Sydney, Australia) and Alec Gilbert (University of New South Wales Sydney, Australia)
Description: Quasi-Monte Carlo methods have recently been successfully used in calculating expected values of quantities of interests which depend on the solutions of PDEs with random coefficients, for so-called “uniform” and “lognormal” random fields. The construction of “randomly shifted lattice rules” and “interlaced polynomial lattice rules” based on decay properties of the random field has been implemented in the QMC4PDE software package.
The working group will extend the theory and implementation of these methods in various directions: for example, to generate lognormal random fields by the circulant embedding technique, to stochastic wave propagation, to neutron diffusion, to Bayesian inverse problems in uncertainty quantification, and more.
|Working Group IX: Multivariate Integration and Approximation in the Context of IBC||Working Group Leader: Peter Kritzer (Johann Radon Institute, Austria)
Description: The working group studies multivariate algorithms for numerical integration and function approximation in the context of Information-Based Complexity (IBC). Apart from deriving error bounds for these algorithms, we are particularly interested in studying how much information about a given problem is needed to solve it to within a given error bound.
|Working Group X: Simpson Plays Billiards?||Working Group Leader: Dirk Nuyens (KU Leuven – Belgium)
Description: Lattice rules are normally applied in the context of integrating periodic functions. There are however several results known in which lattice rules can be used in the context of non-periodic function spaces. We are interested in obtaining higher-order convergence for such non-periodic function spaces.
As a first wild idea, inspired by an interpretation of tent-transormed lattice rules as trapezoidal rules over a billiard ball trajectory, we might look at a fractal-transform which generates a Simpson’s rule. We are then interested in deriving worst-case error bounds which show higher-order convergence.
Program on Optimization (OPT)
|Statistical Inverse Problems
Working Group Leaders: Alen Alexandrian, Arvind Saibaba
|EM & MM Algorithms
||Working Group Leaders: Hua Zhou, Eric Chi, Ekkehard Sachs, Dirk Lorenz
SAMSI Webmaster: Liuyu Hu
|Sums of Squares and Semidefinite Programming||Contact: email@example.com
Working Group Leaders: David Papp
SAMSI Webmaster: Sercan Yildiz
|Bayesian Optimization & Decision Analysis||Contact: firstname.lastname@example.org
Working Group Leaders: Mike West
SAMSI Webmaster: Mike Lindon
|Mixed Interger – PDE Constrained Optimization||Contact: email@example.com
Working Group Leaders: Sven Leyffer, Ekkehard Sachs and Bart van Bloemen Wannders
SAMSI Webmaster: Joey Hart
|Modeling & Computation of Equilibrium Problems||Contact: firstname.lastname@example.org
Working Group Leaders: Shu Lu
SAMSI Webmaster: Hongsheng Liu
|Applications to Energy and the Environment||Contact: email@example.com
Working Group Leaders: Mihai Anitescu and Jianfeng Lu
SAMSI Webmaster: Peter Diao
|Probabilistic Numerics||Working Group Leaders: David Bortz and Vanja Dukic|
|Small Area Estimation||Contact: firstname.lastname@example.org|
|Radiotherapy Optimization||Contact: email@example.com
Working Group Leaders: David Papp
SAMSI Webmaster: Melissa Gaddy
|Faster Statistics||Contact: firstname.lastname@example.org
Working Group Leaders: Xiaoming Huo
SAMSI Webmaster: Cheng Huang
|Optimization for Electronic Structure Models||Working Group Leaders: Jianfeng Lu|
|Minimum action principle for multi-physics problems||Working Group Leaders: Jianfeng Lu
SAMSI Webmaster: Jeff LaComb
|Course: Fall 2016: Numerical Optimization and Applications – Part I||Description: In this course we consider the numerical solution of constrained and unconstrained optimization problems. We review the classical schemes with a special emphasis on new and modern developments in this area like e.g. nonlinear trust region methods, inexact Newton methods or conditional gradient methods.
We also analyze how these methods can be especially adapted to optimization problem with constraints that come from ordinary, partial or even stochastic differential equations. The computation of Lagrange multipliers play an important role in this context from a traditional point of view but also from recent research using second order adjoints.
Another important topic in this course will be applications in optimization which come from problems in statistics and which often exhibit special structure that can be exploited to reduced computation time and storage.
Schedule: Tuesdays 4:30 pm at SAMSI, Research Triangle Park, NC, beginning Tuesday, September 6, 2016
No class during week of Thanksgiving: Tuesday November 22, 2016
Last class: Tuesday, November 29, 2016
|Course: Numerical Optimization and Applications – Part II||Description: The main topic of this class is the numerical solution of constrained optimization problems. We review the classical schemes with a special emphasis on new and modern developments in this area. We also analyze how the algorithms can be especially adapted to optimization problem with constraints that come from ordinary, partial or even stochastic differential equations. Another important topic in this course will be applications in optimization which come from problems in statistics and which often exhibit special structure that can be exploited to reduced computation time and storage.
Schedule: Tuesdays 4:30 pm at SAMSI, Research Triangle Park, NC, beginning Tuesday, January 17, 2017
(No class March 14, 2017 for Spring Break)
Last class: Tuesday, April 25, 2017
Program on Statistical, Mathematical and Computational Methods for Astronomy (ASTRO)
|Uncertainty Quantification and Astrophysical Emulation (UQAE)||Contact: email@example.com
Working Group Leaders: Derek Bingham (Simon Fraser) and Earl Lawrence (LANL)
SAMSI Webmaster: David Stenning
|Synoptic Time Domain Surveys (TDA)||Contact: firstname.lastname@example.org
Working Group Leaders: Ashish Mahabal (Astro, Caltech) and G. Jogesh Babu (Stat, PSU)
SAMSI Webmaster: David Jones
|Multivariate and Irregularly Sampled Time Series (MISTS)||Contact: email@example.com
Working Group Leaders: Ben Farr (Astro, U. Chicago) and Soumen Lahiri (Stat, NCSU)
SAMSI Webmaster: Hyungsuk Tak/David Jones
|Astrophysical Populations (AP)||Contact: firstname.lastname@example.org
Working Group Leaders: Jessi Cisewski (Stat, Yale) and Eric Ford (Astro, Penn State)
SAMSI Webmaster: David Stenning
|Statistics, Computation, and Modeling in Cosmology (COSMO)||Contact: email@example.com
Working Group Leaders: Jeff Jewell (Astro, JPL) and Joe Guinness (Stat, NCSU)
SAMSI Webmaster: Hyungsuk Tak
|Course: Fall 2016: Analytical Methods and Applications to Astrophysics and Astronomy||Description: With the advance of digital imaging techniques, astronomy has become a data science in which knowledge creation depends on applying and developing sophisticated statistical methodology to large and/or complex data sets. This course will cover common types of data in astronomy such as light curves, spectra, and images as well as statistical methods used for analyzing these data sets, such as functional data analysis, measurement error models, hierarchical models, survival analysis, and machine learning techniques. An emphasis will be placed on the complexity of the inference tasks faced by astronomers and the propagation of uncertainty across several levels of inference. Guest lecturers will discuss topical issues in the analysis of astronomy data. Students will complete a final project in groups based on data or statistical methodology presented during the course.
The course will be aimed at a wide audience in an effort to appeal to students with either an astronomy or statistics background. While there are no formal prerequisites students will benefit from having some past experience with medium to large data sets and a familiarity with statistical methods such as maximum likelihood and regression. Theory will be kept to a minimum with an emphasis instead on astronomy data and statistical methodology.
Schedule: Wednesdays 4:30 pm at SAMSI, Research Triangle Park, NC, beginning Wednesday, September 7, 2016
No class during week of Thanksgiving: Wednesday November 23, 2016
Last class: Wednesday, November 30, 2016
|Course: Spring 2017: Time Series Methods for Astronomy||Description: The course starts with an overview of variable cosmic phenomena and characteristics of astronomical time series. Classical time series analysis in the time and frequency domain for evenly spaced data will be reviewed. This includes Gaussian and Poisson processes, smoothing and interpolation, autocorrelation and autoregressive modeling, Fourier analysis, and wavelet analysis. The class then proceeds to treatments of unevenly spaced time series commonly found in astronomical datasets, again in both the time and frequency domain. Guest lectures by expert SAMSI scholars developing advanced techniques for unevenly spaced data will be featured. Throughout the course, methods will be exercised using the public domain R statistical software environment using contemporary astronomical datasets. Students will complete R-based homeworks and a personal project in time series analysis involving a dataset of their choice.
Schedule: January 11 – April 26, 2017 / Wednesdays at 4:30PM – 7:00PM
|Multiple Sources of Bias||Contact: firstname.lastname@example.org
Working Group Leaders: Sandy Zabell, Cliff Spiegelman
|Pattern Evidence||Contact: email@example.com
Working Group Leaders: Karen Kafadar, Anil Jain
|Forensic Experiments||Contact: firstname.lastname@example.org
Working Group Leader: Dennis Lin
|Statistical Evidence||Contact: email@example.com
Working Group Leaders: Colin Aitken, Anjali Mazumder
|Possible Matches||Contact: firstname.lastname@example.org
Working Group Leader: Len Stefanski
|Ballistic Images||Contact: email@example.com
Working Group Leaders: Nell Sedransk, Cliff Spiegelman, Sarena Wiesner
|Forensic Evidence||Contact: firstname.lastname@example.org
Working Group Leader: Cedric Neumann
|Clinical Brain Imaging||Contact: email@example.com
Working Group Leader: Ciprian Crainiceanu, Johns Hopkins University
|Computational Approaches to Large-scale
Inverse Problems with Applications to Neuroscience
Working Group Leader: Arvind Saibaba, North Carolina State University
|Understanding Neuromechanical Processes in
Locomotion with Physical Modeling and Network Analysis
Working Group Leaders: Laura Miller, UNC and Katie Newhall, UNC
|Mathematical and Statistical Approaches to Modeling Brain Networks||Contact: firstname.lastname@example.org
Working Group Leaders: Rob Kass, Carnegie Mellon University; Uri Eden, Boston U; Mark Kramer, Boston U
|Theory of neural networks: structure and dynamics||Contact: email@example.com
Working Group Leaders: Carina Curto, PSU; Brent Doiron, U. of Pittsburgh; Chris Hillar, MSRI
|Acquisition, Reconstruction, and Processing of MRI Data||Contact: firstname.lastname@example.org
Working Group Leader: Daniel Rowe, Marquette University
|Imaging Genetics||Contact: email@example.com
Working Group Leader: Hongtu Zhu, UNC
|Structural Connectivity||Contact: firstname.lastname@example.org
Working Group Leaders: David Dunson, Duke University; Hongtu Zhu, UNC
|Functional Imaging Methods and Functional Connectivity||Contact: email@example.com
Working Group Leaders: Hernando Ombao, UCI; John Aston, University of Cambridge
|Big Data Integration in Neuroimaging||Contact: firstname.lastname@example.org
Working Group Leaders: Martin Lindquist and Timothy Johnson
Login: Login: https://sakai.duke.edu/sakai-login-tool/container
|Analysis of Optical Imaging Data||Contact: email@example.com
Working group leader: Mark Reimers
Login: Login: https://sakai.duke.edu/sakai-login-tool/container
|Fall Course 2014-2015: Stochastic Process Modeling for Ecological Processes||Schedule: Wednesdays 4:30 pm at SAMSI, Research Triangle Park, NC, beginning September 3, 2014 No class during week of Thanksgiving: Wednesday Nov 26, 2014 Last class: December 3, 2014 Instructors: L. Miller, A. Lloyd, D. Adalsteinsson, J. Clark, A. Gelfand The central theme of this proposed course is the use of stochastic process models to introduce desired behaviors into models for complex ecological processes. Objective include both emulation and data analysis. Specifically, there will be a focus on (stochastic) pde’s for dynamical systems with application to invasive and emerging infections (Lloyd) and to dispersal and spatial effects in competition models (Miller). There can be demanding coding aspects to this work with the possibility of some flipped classroom “lectures” (Adalsteinsson, Miller). On the more statistical side, the emphasis would be on hierarchical modeling with examples of hierarchical models for inference on seed dispersal and pde’s for population growth (Clark) supplemented with space-time diffusions for invasive species (Gelfand). Model development here would be complemented with some computer labs in R. Registration for this course is being processed through your respective university: UNC-CH: STOR 930 Section 001 (cross-listed with MATH 892 Section 001) Duke: STA 790.04 NCSU: MA 810.002 For additional information about this course, send e-mail to firstname.lastname@example.org|
|Kepler Working Group||This group is specifically geared toward starting discussions to help people prepare for the actual meeting in June 2013. There are already discussion threads where you are encouraged to: introduce yourself, discuss ideas for focused working groups, and suggest background reading that may be helpful for other participants. Please contribute to via comments section of each discussion. Also, feel free to start a discussion topic of your own. (Membership is limited to invited participants.)|
|Online Streaming and Sketching||Purpose: Methodology and fast algorithms for computing leverage scores, with application to astronomy and genomics. Preliminary list of topics: * Leverage scores: Computation, fast approximation, sensitivity, numerical stability of algorithms, behavior under sketching * Randomized low-rank approximations: Subset selection, CUR, Nyström, PCA, robust PCA, subsampled regression, regression on manifolds, construction of robust linear models, windowed and online streaming approaches * Randomized sketching and importance sampling strategies: Methodology and numerical computation * Local vs global: Eigenvector and invariant subspace localization, eigen-analysis of data connectivity matrices, numerical stability of streaming and updating methods, relation to generalized eigenvalue problems * Data fusion/integration: Robust and fast/streaming methods with application to galaxy formation and evolution|
|2013-14 Course: Geometric and Topological Summaries of Data and Inference||The course will focus on geometric and topological summaries computed from data that are routinely generated across science and engineering. The focus is on modeling objects that have geometric or topological structure. Examples include curves, or surfaces such as bones or teeth, or objects of higher dimension such as positive definite matrices, or subspaces that describe variation in phenotypic traits due to genetic variation, or the geometry of multivariate trajectories generated from cellular processes. Specific topics will include the following. (1) Geometry in statistical inference — Material covered will include recent work in machine learning and statistics on the topics of manifold learning, subspace inference, factor models, and inferring covariance/positive definite matrices. Applications will be used to highlight methodologies. The focus will be on methods used to reduce high-dimensional data to low-dimensional summaries using geometric ideas. (2) Topology in statistical inference — Material covered will focus on probabilistic perspectives on topological summaries such as persistence homology and on inference of topological summaries based on the Hodge operator and the Laplacian on forms. Again, applications will be used to highlight methodologies. (3) Random geometry and topology — Material will cover the geometry and topology induced by random processes. Topics include the topology of random clique complexes, random geometric complexes, limit theorems of Betti numbers of random simplicial complexes. (4) Applications of the Laplacian operator in data analysis — Material will cover the various uses of the Laplacian in data analysis, including manifold learning, spectral clustering, and Cheeger inequalities. More advanced topics will include the Hodge operator or combinatorial Laplacian and applications to data analysis including decomposing ranked data into consistent and inconsistent components, inference of structure in social networks, and decomposing games into parts that have Nash equilibria and parts that cycle. Prerequisites: Background in calculus and linear algebra and some reasonable foundation in statistics and probability. Course Format: The main instructor will be Sayan Mukherjee but there will be several guest lecturers, with material and instructors paralleling certain of the major themes in the 2013-2014 year-long SAMSI program on Low-Dimensional Structure in High-Dimensional Systems (LDHD). All course updates including example projects, reading material, and lecture slides will be posted at http://www.stat.duke.edu/~sayan/SAMSI/.|
|Fall 2013 Course: Computational Methods for Social Sciences||Coordinating Instructor: Richard L. Smith (Department of STOR, UNC) Lead Instructors: David Banks, Department of Statistical Science, Duke; Thomas Carsey, Department of Political Science and Odum Institute, UNC; Peter Mucha, Departments of Mathematics and Applied Physical Sciences, UNC; Jerry Reiter, Department of Statistical Science, Duke Time and Place: Wednesdays from 4:30pm-7:00pm beginning August 28, 2013; Statistical and Applied Mathematical Sciences Institute, 19 T.W. Alexander Drive, RTP, N.C. First Class is Wednesday, August 28, 2013 Course description: The Statistical and Applied Mathematical Sciences Institute (SAMSI) is hosting a year-long research program on Computational Methods for Social Sciences. As part of this program, there will be an advanced graduate course on the topics of the program. The syllabus will cover the three main research themes in the program: (a) Social Networks; (b) Statistical Methods for Censuses and Surveys; (c) Agent-based models. The course will meet once a week at SAMSI and will consist primarily of lectures by senior researchers. The course will be suitable for advanced graduate students in Mathematics, Statistics, Biostatistics or quantitative social sciences (e.g. Sociology, Political Science, Psychology). There are no specific prerequisites. Assessment will be by class presentations or a written project, the exact format to be determined partly based on the number of students participating. For additional information about the course, send e-mail to email@example.com.|
|CMSS: Social Networks|
|CMSS: Causal Inference|
|CMSS: Censuses and Surveys||This workng group is investigating several methods related to computation in surveys and censuses. The group is particularly interested in combining information from multiple sources. One specific application is record linkage. How does one account for uncertainty in imperfect linkage? How effective is creating a joint distribution from available information compared to imperfect, multi-way record linkage? What should one do with survey weights in linked files? Another specific application is merging big data with probability samples. Can we use information from surveys to help generalize analyses from large-scale administrative/private or organic data?|
|CMSS: Weighting in Surveys|
|CMSS: Agent-based Models|
|High-dimensional Graphical Models||Graphical models are now a standard tool in statistics, and high-dimensional graphical models have been studied extensively. One of the ideas behind graphical models is to break down a high-dimensional problem into several low-dimensional ones. The difficult problem is model selection: one way to do that is to identify the low-dimensional components and see how they fit together. Buhlmann and Meinhausen (2006) were among the first to introduce the local method of neighborhood regression for model selection among graphical Gaussian models. More recent approaches include local l_1 regularized logistic regression for model selection among discrete Ising models, the concept of sparse local separators, and the usage of neighborhood structure. One of the activities of the working group will be review of the current iterature on local methods for structure estimation in high dimension. Another will be the exploration of geometric and topological local methods for the identification of structure in graphical models. Suggestions for other activities are, of course, welcome.|
|Data Analysis on Hilbert Manifolds and their Applications||The theoretical focus of this Working Group is nonparametric statistics on Hilbert manifolds and dimensionality reduction from infinity to low dimension as small as 1. This will include an extension of CLT for iid variables to infinite dimensional Hilbert manifolds and beyond — to infinite dimensional stratified spaces — as well as the neighborhood hypothesis testing methodology, extending recent results to two or multiple samples on Hilbert manifolds or even to infinite dimensional stratified spaces. Applications to Hilbert manifolds data analysis that will be potentially discussed, depending on the WG participants, could be to any of the following: i. MRI imaging including MRI, DTI, or f-MRI; ii. CT imaging; iii. eye medical imaging, including stereo imaging; iv. 2D and 3D scene recognition from similarity shapes of curves and surfaces; v. projective shape of 3D scenes reconstructed from digital camera imaging; vi. color and texture imaging data; vii. spatial and temporal data on the geoid; viii. plate tectonics data and continental drift; ix. paths of eyes of storms on planet Earth; x. volcanic activity, earthquakes, and other Earth Sciences data; xi. astronomy data; xii. solar system data; xiii. DNA based data.|
|Nonlinear Low-dimensional Structures in High-dimensions for Biological Data||Many problems in biology, particularly at the molecular level, involve very high dimensions. Traditional methods, such as principal component analysis provide, an important first step toward understanding the underlying lower dimensional structure. Nevertheless, finding even lower-dimensional and possibly nonlinear structures requires more qualitative procedures. Some possible candidates, among others, include geodesic principal component analysis, SiZER analysis, and persistent homology. The aim of this Working Group will be to describe and understand the nature of certain biological data coming from areas such as brain artery tree networks, gene networks, biomechanical motion data, and metagenomics, to name a few. We will explore and formulate strategies to best apply these qualitative procedures. In addition, a recent insight into a connection between SiZER analysis and persistent homology will be discussed.|
|Inference: Dimension Reduction||Dimension reduction methods search for low dimensional structure in high dimensional data. Although unsupervised methods such as PCA have been used for over 80 years, new issues in robustness, very high dimensionality p >> n, nonlinear structure, and so on continue to provide rich sources of statistical and computational issues. Supervised methods, such as Sliced Inverse Regression, which look for low dimensional association between predictors and response have been developed within the past 35 years, and many useful extensions continue to be developed. Performing valid inference to correctly account for the impact of dimension reduction including variable selection also poses problems that have not been thoroughly explored. This working group will work on problems of extending unsupervised and supervised methods to some of the interesting new areas such as machine learning, functional data, discrete data, and complex data structures such as networks. As well, we will work on inferential issues.|
|SAMSI Fall Course – Operations Research Methods in Healthcare||Principal Instructor: V. Kulkarni Course Day and Time: Course will be held at SAMSI (driving directions) in RTP on Wednesdays, 4:30-7:00 p.m. in Room 150. Schedule: First class Wednesday, September 5, 2012 ; last class day, Thursday, December 6, 2012 Course Description: This is a seminar-style course treating application of operations research methods such as stochastic modeling, queuing theory (including fluid models), optimization and simulation to problems in healthcare. Potential problems to be studied are data-based design of healthcare operations, patient flow, scheduling of facilities and personnel, management of transplant lists, mass casualty events and comparative effectiveness research. Students will be expected to read and make presentations of material from the relevant literature. Registration for this course is being processed through your respective university: Duke: STA 790-02 NCSU: MA 810.002 UNC: STOR 892.1 Questions about the course or the Healthcare program should be emailed to firstname.lastname@example.org|
|MD Imaging||The focus of the Imaging Working Group will be on methodological and computational questions of statistics, mathematics and computer science posed by imaging science and technology, with applications to either astronomy, high energy physics, the environment, health sciences or other areas. This will be a venue to discuss the challenges, motivate and develop innovative approaches towards computing environments, analyses, methods, algorithms and tools, in relation to imaging science and innovative applications.|
|MD Online Streaming & Sketching||Working Group leaders: Petros Drineas, Ilse Ipsen, and Michael Mahoney Webmasters: John Holodnak and Kevin Penner Purpose: Development and analysis of fast randomized algorithms for computing leverage scores, and their application Preliminary, Partial List of Topics: Approximating leverage scores for L2 and other regression problems: Online, streaming, incremental streaming algorithms Numerical analysis: Sensitivity of leverage scores, numerical stability of algorithms Applications in astronomy: Characterization of streaming & time dependent aspects of low rank approximations Incremental computation of leverage scores Applications in feature selection: How to distinguish among almost identical columns with high leverage scores (RRQR factorization, clustering) Derivation of formal bounds|
|Statistics of Extremes – Climate and Methodology UQ||This group will examine the characterization of extreme events from a statistical point-of-view. The overall group will be composed of several project groups, each with a specific focus. The project groups may be application-driven, may work develop new statistical methodologies, or may work to further the theory on which extreme value analyses rely. The overall group will provide a structure for the project groups as well as provide an environment for investigating aspects which have general interest.|
|Parallel Computing Issues – Climate UQ||– adaptive design of experiments, and resolution issues – embedding of emulated sub-models to resolve sub-processes that are now computationally prohibitive – python open source software open platform|
|Simulation of Rare Events – Methodology UQ||This group will focus on methods for simulating rare events in high-dimensional physical systems, especially PDE models. We will explore the use of importance sampling and large deviation theory in order to identify important mechanisms or configurations of the parameters that lead to rare events. We also will consider the role of asymptotic analysis in constructing effective sampling weights for such computations.|
|Data Assimilation – Methodology UQ||Data assimilation is the process of fusing information from imperfect models, noisy measurements, and priors, to produce an optimal representation of the state of a physical system. Data assimilation can be interpreted and carried out in a Bayesian framework. Practical methods for large-scale systems include suboptimal and the ensemble Kalman filter approaches, optimal interpolation, and three and four dimensional variational methods. This working group will focus on emerging problems that include, but are not limited to: new computational algorithms, modeling of model errors, impact of observations, and quantification of posterior uncertainties. There will be strong ties between theory and applications investigated throughout the program.|
|Stochastic to Deterministic Models and Back Again – Methodology UQ||Models of complex multiscale and/or multiphysics phenomena often require combining stochastic and deterministic models. Direct coupling of stochastic and deterministic models, e.g. molecular dynamics with a continuum model Stochastic parameterization, with parameters determined by a stochastic model simulation or other statistical models Such models are often used to predict “engineering scale” questions from limited microscale information. Of course, this is a classic analysis/modeling problem. The working group will focus on computational issues, including: Rigorous formulation and analysis of coupling mechanisms and their discretizations Numerical treatment of averaging and computed expectations and the effect of approximations A posteriori error analysis, resolutions required in different components, adaptive computation Rigorous treatment of feedback between stochastic and deterministic models, e.g. nonlinear iterative methods, convergence|
|Model Validation – Methodology UQ||Model validation refers to the process of assessing the accuracy with which mathematical models can predict physical events, or, more specifically, quantities of interest observed in physical phenomena. Validation should be a prerequisite for predictive modeling, which often forms the basis for decision-making. This working group will study the principles, merits, and limitations of various probabilistic approaches to model validation. Special emphasis will be laid on methods for splitting datasets for calibration and validation purposes, on the analysis of model discrepancies, on the development of rejection metrics, and on any other issues of interest raised during the working group meetings. The working group is organized and managed by Serge Prudhomme (ICES, UT Austin), Sujit Ghosh (Statistics, NCSU), and Jan Hannig (STOR, UNC).|
|Multiphysics – Methodology UQ||Multiphysics models comprising compositions of models of several physical processes, often at different scales, dominate many areas of science and engineering. The working group will study UQ topics for MP models (including both forward and inverse topics) Research issues: Complex feedback between physical processes, highly nonlinear responses Complex and unresolved coupling mechanisms Different kinds and representations of uncertainty for different components, and complex interactions between sources of uncertainty and error Complex, high dimension parameter space Bifurcations and discontinuous model changes High performance computational issues|
|Approximating Computationally Intensive Functions and Sampling Design in High Dimensions – Methodology UQ||This is a core problem in constructing surrogates, with application to uncertainty propagation, inference, prediction, and design. Relevant research questions: Cross-examination of different methods: Projection, regression, interpolation, L1 minimization, Gaussian process/kriging, etc Appropriate measures of performance/accuracy and their dependence on the intended use of the surrogate model. Error analysis and convergence properties Sparse representations: l1 minimization, pursuit algorithms, low-rank approximation “Optimal” choices of nodes/design points for different methods Adaptive approaches: a posteriori error estimates; derivative properties; dimension reduction; additional optimization of nodes; ANOVA; relation to sequential design of experiments Interpreting and combining uncertainty information from stochastic surrogates (e.g., Gaussian process variance) and deterministic error bounds Deriving optimality criteria and search algorithms that are good for high dimensions Borrow existing theoretical results in high dimensional statistics (Donoho, etc.) to shed light on the structure of “optimal” designs in high dimension. Additional issues— Lack of regularity: discovering and approximating discontinuities in high dimensions Incorporating gradient information Enforcing constraints on output and input domains Fault tolerance and missing samples|
|Inverse Function-based Inference – Methodology UQ||Jan Hannig (lead), D. Estep, Troy Butler (U. Texas), Simon Tavener (CSU) The working group will study the use of set-valued inversion of models for inference Research issues: Approximation of set-valued inverses in complex spaces Computation of inverse measures in parameter space Convergence and accuracy of computed inverse measures Theoretical issues regarding inversion of multiple observations Relation to fiducial inference and Dempster and Shafer calculus Intrusive and non-intrusive algorithms, dimension-benign computational algorithms|
|Surrogate Models – Methodology UQ||This working group focuses on the exploration of properties, utility, and performance of two classes of model surrogates, namely Polynomial Chaos and Gaussian Process surrogates. The study will be done in the context of specific model problems with a range of difficulty involving nonlinearity and dimensionality. Test problems will include both algebraic functions as well as simple ODE/PDE problems.|
|Engineered Systems – Engineering UQ|
|Sustainability – Engineering UQ|
|Materials – Engineering UQ|
|Renewable Energy – Engineering UQ||The group will consider uncertainty quantification issues arising in specific applications linked to renewable energy. In particular, we will study biofuels and wind farms. Other aspects involves the inclusion, in a UQ framework, of factors such as technological advances and/or regulations.|
|Nuclear Energy – Engineering UQ|
|Geosciences – Geosciences UQ|
|Data Assimilation in IPCC Level Models – Climate UQ||In the IPCC AR4, the results were based on runs from ca. 24 models. These were built and run at climate research centers around the world and are each integrated Earth system models that comprise many components, including atmosphere, ocean, ice and land. The so-called dynamical core of such models is a computational model covering both the atmosphere and ocean and based on the primitive equations of GFD. While this is essentially computational, data come into the process of forming the final model. This incorporation occurs at a number of stages of the model development, including parametrization of sub-grid scale effects and model tuning. The process is not, however, done systematically and current practice is not thought of as “data assimilation.” There seems to be a growing realization that DA will have a significant role to play in future climate model development. This is, in part, driven by the need to quantify uncertainty in the model predictions. Nevertheless, there is not a consensus as to how DA should be used in these large-scale climate models. This working group will consider the issues involved in formulating a plan for DA in such models. The first step will be to understand how such a model is put together and uncover all the steps where data is currently used in the model formation. For this purpose, we will look at the latest CESM from NCAR.|
|Numerical Methods for Uncertainty Quantification – Spring Part 2||Principal Instructors: Various Course Day and Time: Course will be held at SAMSI (driving directions) in RTP on Wednesdays, 4:30-7:00 p.m. in Room 150. Schedule: First class Wednesday, January 18, 2012 ; last class Wednesday, April 25, 2012 This course focuses on numerical methods for stochastic computation and uncertainty quantification (UQ). It is a two-semester course, where the first semester focuses on fundamental materials for UQ computing and the second semester on more advanced research materials. The main topics covered in the second semester include advanced numerical techniques for SPDE: adaptive methods and compressive sensing propagation of probability distributions Bayesian inference data assimilation model calibration Prerequisites: Numerical linear algebra, numerical methods for ordinary and partial differential equations; programming skill in one language, e.g., C/C++, FORTRAN, or Matlab. The course is open to all levels of graduate students in Mathematics, Statistics, as well as to those in other departments of sciences and engineering. Senior level undergraduate students with outstanding background are also considered. Registration for this course is being processed through your university. Duke: STA 294-01 NCSU: MA 810.002 UNC: STOR 891-001 Questions about the course or the UQ program should be emailed to email@example.com.|
|Dynamics OF Networks||Complex Networks Program working group.
The dynamics of networks working group is exploring a variety of mathematical and statistical approaches for describing and understanding the changing connection topology of networks over time, the interplay of these network dynamics with other dynamic processes on the network, and the connections between these different mathematical and statistical methodologies.
|Sampling / Modeling / Inference||CN Program working group.
The Working Group on Sampling/Modeling/Inference in networks aims to work towards moving the current state of knowledge on these inter-related tasks — in the specific context of networks — to rest on a more principled and integrated mathematical and statistical foundation. We are pursuing this goal by focusing on a handful of specific prototype problems in the context of certain application areas, ranging from information networks to animal communities to neuroscience.
|Dynamics ON Networks||CN Program working group.
Random graphs are useful models of social and technological networks. To date most of the research in this area has concerned geometric properties of the graphs. This working group will focus on processes taking place ON the network. In particular we are interested in how their behavior on networks differs from that in homogeneously mixing populations or on regular lattices of the type commonly used in ecology and physics.
|Geometrical / Spectral Analysis||CN Program working group
This working group is concerned with the following topics: detection of communities in networks, multiscale spectral methods for the analysis of the geometry of networks, algorithms that simplify graphs into simpler graphs in order to speed up certain optimization problems, metrics for comparing graphs, and multiscale homogenization of random walks. These topics have applications biology and to spread of “epidemics” in financial networks.