Adaptive nonparametric regression
Successful statistical modeling and data analysis often requires the use of smoothing procedures that can adapt to the unknown underlying structure. This includes an intrinsic detection of edges or discontinuities in imaging or signal processing as well as nonstationarities in time series. Established parametric and nonparametric methods often provide results that do not sufficiently meet demands from applications. The general principle of structural adaptation forms the basis of most of the new smoothing algorithms developed at WIAS. This approach enables an efficient description and analysis of complex discontinuous models. Within the recent years a series of innovative approaches and algorithms has been developed. Examples are adaptive weights smoothing (AWS) (Polzehl & Spokoiny 2000), , propagation separation methods (PS), (Polzehl & Spokoiny 2006, Becker & Mathé 2013), Pointwise Adaptive Methods (PA), (Spokoiny 1998, Polzehl & Spokoiny 2003, Spokoiny & Vial 2009), local model selection (LMS), stagewise aggregation (SA) and (Belomestny & Spokoiny 2007) local change point detection (LCP) (Mercurio & Spokoiny 2004, Spokoiny 2009). All these procedures are based on a new local parametric approach and have been justified within a series of publications on theoretical results. In this context important problems and questions in the field of modern parametric and nonparametric statistics have been solved. This includes the problem of automatic selection of smoothing parameters in case of qualitative structural assumptions. Structural adaptive methods have been developed and investigated for several classes of statistical models. Apart from regression this especially includes nonstationary time series and volatility models but also adaptive classification methods and models for extremes. Strong methodological connections exist to adaptive methods for dimension reduction and to multiscale methods. Several algorithms have been implemented as contributed packages ( adaptive weights smoothing aws, adaptive image processing adimpro, functional magnetic resonance imaging fmri and diffusion weighted magnetic resonance imaging dti) (Image processing/Neuroimaging Packages) of the R Environment for Statistical Computing . For applications in neuroimaging see Computational neuroscience and medical imaging.Software Adaptive weights clustering
Software description can be found at a separate page.Dimension reduction
The problem of dimension reduction arises in various forms with modeling or analysis of complex, usually highdimensional data structures. Examples are relations between variables in regression models where single or multiindex models (Hristache et al. 2001, Dalalyan et al. 2008) allow for a more efficient description and methods that infer on linear or nonlinear subspaces that contain most of the relevant information on underlying structures. Such models and corresponding procedures are often a prerequisite for an efficient statistical modeling and data analysis. Classical procedures like principal component analysis (PCA) or independent component analysis (ICA) are based on strict assumptions and are of limited use in this context. Based on the principle of structural adaptation novel methods and procedures for dimension reduction have been developed at WIAS. The proposed dimension reduction methods for regression can be interpreted as nonlinear PCA. These methods allow to efficiently estimate a dimension reduction space after a finite number of iterations. NonGaussian componet analysis (NGCA) (Blanchard et al. 2006 , Diederichs et al. 2013) provides a new promising approach to the analysis of highdimensional data structures. It is based on the assumption that structure lives in lowdimensional subspaces and projections onto these subspaces are characterized by nonGaussianity. This assumption is then used to identify these informative subspaces.Simultaneous Statistical Inference
The more questions you ask, the more wrong answers you are expected to receive  even if every single source of your information is quite trustworthy. We consider cases in which the sources of information are data, and the questions are formalized by statistical hypothesisalternative pairs. From the mathematical point of view, this leads to multiple test problems. Criteria and methods (multiple tests) are developed which ensure that with high probability not too many wrong decisions are made, even if many hypotheses are of interest under the scope of one and the same statistical model, i. e., regarding one and the same dataset. Highthroughput technologies in different fields of modern life sciences have led to massive multiplicity and given rise to multiple test problems with more hypotheses than observations. Driven by these developments, also new statistical paradigms have arisen. In particular, it is fair to say that a new era of multiple testing began when Yoav Benjamini and Yosef Hochberg formally introduced the false discovery rate (FDR) and the linear stepup test for FDR control in 1995. At WIAS, FDR methodology and multivariate multiple tests are developed for temporally, spatially, or spatiotemporally dependent data. Furthermore, multiple test theory can also be used in related fields of simultaneous statistical inference and statistical learning, for example in the context of binary classification and model selection. Our contributions to simultaneous statistical inference are made with particular emphasis on life science applications in (behavioral) genetics, epigenetics, braincomputer interfacing, functional magnetic resonance imaging, and neuroeconomics. Main publications: Blanchard & Roquain 2009, Blanchard et al. 2010, Dickhaus 2013, Dickhaus et al 2013, Bodnar & Dickhaus 2017, Schildknecht et al. 2015, Schildknecht et al. 2016.Publications
Monographs

J. Polzehl, K. Tabelow, Magnetic Resonance Brain Imaging: Modeling and Data Analysis using R, Series: Use R!, Springer International Publishing, Cham, 2019, 231 pages, (Monograph Published), DOI 10.1007/9783030291846 .
Abstract
This book discusses the modeling and analysis of magnetic resonance imaging (MRI) data acquired from the human brain. The data processing pipelines described rely on R. The book is intended for readers from two communities: Statisticians who are interested in neuroimaging and looking for an introduction to the acquired data and typical scientific problems in the field; and neuroimaging students wanting to learn about the statistical modeling and analysis of MRI data. Offering a practical introduction to the field, the book focuses on those problems in data analysis for which implementations within R are available. It also includes fully worked examples and as such serves as a tutorial on MRI analysis with R, from which the readers can derive their own data processing scripts. The book starts with a short introduction to MRI and then examines the process of reading and writing common neuroimaging data formats to and from the R session. The main chapters cover three common MR imaging modalities and their data modeling and analysis problems: functional MRI, diffusion MRI, and MultiParameter Mapping. The book concludes with extended appendices providing details of the nonparametric statistics used and the resources for R and MRI data.The book also addresses the issues of reproducibility and topics like data organization and description, as well as open data and open science. It relies solely on a dynamic report generation with knitr and uses neuroimaging data publicly available in data repositories. The PDF was created executing the R code in the chunks and then running LaTeX, which means that almost all figures, numbers, and results were generated while producing the PDF from the sources. 
H.G. Bartel, H.J. Mucha, Chapter 2: Incomparability/Inequality Measures and Clustering, in: Partial Order Concepts in Applied Sciences, M. Fattore, R. Brüggemann, eds., Springer International Publishing, New York, 2017, pp. 2134, (Chapter Published).

V. Spokoiny, Th. Dickhaus, Basics of Modern Mathematical Statistics, 18 of Springer Texts in Statistics, Springer, Berlin et al., 2015, 296 pages, (Monograph Published).

H.G. Bartel, H.J. Mucha, Chapter 3: Measures of Incomparability and of Inequality and Their Applications, in: Multiindicator Systems and Modelling in Partial Order, R. Brüggemann, L. Carlsen, J. Wittmann, eds., Springer, New York et al., 2014, pp. 4767, (Chapter Published).

TH. Dickhaus, Simultaneous Statistical Inference, Springer, Berlin et al., 2014, 180 pages, (Monograph Published).

J. Polzehl, K. Tabelow, Chapter 4: Structural Adaptive Smoothing: Principles and Applications in Imaging, in: Mathematical Methods for Signal and Image Analysis and Representation, L. Florack, R. Duits, G. Jongbloed, M.C. VAN Lieshout, L. Davies, eds., 41 of Computational Imaging and Vision, Springer, London et al., 2012, pp. 6581, (Chapter Published).
Articles in Refereed Journals

N. Puchkin, V. Spokoiny, Structureadaptive manifold estimation, Journal of Machine Learning Research (JMLR). MIT Press, Cambridge, MA. English, English abstracts., 23 (2022), pp. 162.
Abstract
We consider a problem of manifold estimation from noisy observations. Many manifold learning procedures locally approximate a manifold by a weighted average over a small neighborhood. However, in the presence of large noise, the assigned weights become so corrupted that the averaged estimate shows very poor performance. We suggest a novel computationally efficient structureadaptive procedure, which simultaneously reconstructs a smooth manifold and estimates projections of the point cloud onto this manifold. The proposed approach iteratively refines the weights on each step, using the structural information obtained at previous steps. After several iterations, we obtain nearly öracle" weights, so that the final estimates are nearly efficient even in the presence of relatively large noise. In our theoretical study we establish tight lower and upper bounds proving asymptotic optimality of the method for manifold estimation under the Hausdorff loss. Our finite sample study confirms a very reasonable performance of the procedure in comparison with the other methods of manifold estimation. 
A. Kroshnin, V. Spokoiny, A. Suvorikova, Statistical inference for BuresWasserstein barycenters, The Annals of Applied Probability, 31 (2021), pp. 12641298, DOI 10.1214/20AAP1618 .

L. Antoine, P. Pigato, Maximum likelihood drift estimation for a threshold diffusion, , published online on 23.10.2019, urlhttps://doi.org/10.1111/sjos.12417, DOI 10.1111/sjos.12417 .
Abstract
We study the maximum likelihood estimator of the drift parameters of a stochastic differential equation, with both drift and diffusion coefficients constant on the positive and negative axis, yet discontinuous at zero. This threshold diffusion is called the drifted Oscillating Brownian motion. The asymptotic behaviors of the positive and negative occupation times rule the ones of the estimators. Differently from most known results in the literature, we do not restrict ourselves to the ergodic framework: indeed, depending on the signs of the drift, the process may be ergodic, transient or null recurrent. For each regime, we establish whether or not the estimators are consistent; if they are, we prove the convergence in long time of the properly rescaled difference of the estimators towards a normal or mixed normal distribution. These theoretical results are backed by numerical simulations. 
N. Alia, V. John, S. Ollila, Revisiting the singlephase flow model for liquid steel ladle stirred by gas, Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems. Elsevier Science Inc., New York, NY. English, English abstracts., 67 (2019), pp. 549556 (published online on 21.11.2018), DOI 10.1016/j.apm.2018.11.005 .

A. Lejay, P. Pigato, Statistical estimation of the oscillating Brownian motion, Bernoulli. Official Journal of the Bernoulli Society for Mathematical Statistics and Probability, 24 (2018), pp. 35683602, DOI 10.3150/17BEJ969 .
Abstract
We study the asymptotic behavior of estimators of a twovalued, discontinuous diffusion coefficient in a Stochastic Differential Equation, called an Oscillating Brownian Motion. Using the relation of the latter process with the Skew Brownian Motion, we propose two natural consistent estimators, which are variants of the integrated volatility estimator and take the occupation times into account. We show the stable convergence of the renormalized errors? estimations toward some Gaussian mixture, possibly corrected by a term that depends on the local time. These limits stem from the lack of ergodicity as well as the behavior of the local time at zero of the process. We test both estimators on simulated processes, finding a complete agreement with the theoretical predictions. 
D. Belomestny, H. Mai, J.G.M. Schoenmakers, Generalized PostWidder inversion formula with application to statistics, Journal of Mathematical Analysis and Applications, 455 (2017), pp. 89104.
Abstract
In this work we derive an inversion formula for the Laplace transform of a density observed on a curve in the complex domain, which generalizes the well known PostWidder formula. We establish convergence of our inversion method and derive the corresponding convergence rates for the case of a Laplace transform of a smooth density. As an application we consider the problem of statistical inference for variancemean mixture models. We construct a nonparametric estimator for the mixing density based on the generalized PostWidder formula, derive bounds for its root mean square error and give a brief numerical example. 
Y. Nesterov, V. Spokoiny, Random gradientfree minimization of convex functions, Foundations of Computational Mathematics. The Journal of the Society for the Foundations of Computational Mathematics, 17 (2017), pp. 527566.
Abstract
Summary: In this paper, we prove new complexity bounds for methods of convex optimization based only on computation of the function value. The search directions of our schemes are normally distributed random Gaussian vectors. It appears that such methods usually need at most nn times more iterations than the standard gradient methods, where nn is the dimension of the space of variables. This conclusion is true for both nonsmooth and smooth problems. For the latter class, we present also an accelerated scheme with the expected rate of convergence O(n2k2)O(n2k2), where kk is the iteration counter. For stochastic optimization, we propose a zeroorder scheme and justify its expected rate of convergence O(nk1/2)O(nk1/2). We give also some bounds for the rate of convergence of the random gradientfree methods to stationary points of nonconvex functions, for both smooth and nonsmooth cases. Our theoretical results are supported by preliminary computational experiments. 
H.J. Mucha, Assessment of stability in partitional clustering using resampling techniques, Archives of Data Science Series A, 1 (2017), pp. 2135, DOI 10.5445/KSP/1000058747/02 .

V. Spokoiny, Penalized maximum likelihood estimation and effective dimension, Annales de l'Institut Henri Poincare. Probabilites et Statistiques, 53 (2017), pp. 389429, DOI 10.1214/15AIHP720 .

A. Kalinina, A. Suvorikova, V. Spokoiny, M. Gelfand, Detection of homologous recombination in closely related strains, Journal of Bioinformatics and Computational Biology, 14 (2016), pp. 1641001/11641001/12.

A. Andresen, V. Spokoiny, Convergence for an alternation maximization procedure, Journal of Machine Learning Research (JMLR). MIT Press, Cambridge, MA. English, English abstracts., 17 (2016), pp. 153.

D. Belomestny, J.G.M. Schoenmakers, Statistical inference for timechanged Lévy processes via Mellin transform approach, Stochastic Processes and their Applications, 126 (2016), pp. 20922122.

A. Gasnikov, P. Dvurechensky, V. Spokoiny, P. Stetsyuk, A. Suvorikova, Superposition of the balancing algorithm and the universal gradient method for search of the regularized Wasserstein barycenter and equilibria in multistage transport models, Proceedings of Moscow Institute of Physics and Technology, 8 (2016), pp. 524.

K. Schildknecht, K. Tabelow, Th. Dickhaus, More specific signal detection in functional magnetic resonance imaging by false discovery rate control for hierarchically structured systems of hypotheses, PLOS ONE, 11 (2016), pp. e0149016/1e0149016/21, DOI 10.1371/journal.pone.0149016 .

CH. Bayer, A. Moraes, R. Tempone, P. Villanova, An efficient forwardreverse expectationmaximization algorithm for statistical inference in stochastic reaction networks, Stochastic Analysis and Applications, 34 (2016), pp. 193231.

J. Polzehl, K. Tabelow, Low SNR in diffusion MRI models, Journal of the American Statistical Association, 111 (2016), pp. 14801490, DOI 10.1080/01621459.2016.1222284 .
Abstract
Noise is a common issue for all magnetic resonance imaging (MRI) techniques such as diffusion MRI and obviously leads to variability of the estimates in any model describing the data. Increasing spatial resolution in MR experiments further diminish the signaltonoise ratio (SNR). However, with low SNR the expected signal deviates from the true value. Common modeling approaches therefore lead to a bias in estimated model parameters. Adjustments require an analysis of the data generating process and a characterization of the resulting distribution of the imaging data. We provide an adequate quasilikelihood approach that employs these characteristics. We elaborate on the effects of typical data preprocessing and analyze the bias effects related to low SNR for the example of the diffusion tensor model in diffusion MRI. We then demonstrate the relevance of the problem using data from the Human Connectome Project. 
R. Hildebrand, Centroaffine hypersurface immersions with parallel cubic form, Beitrage zur Algebra und Geometrie. Contributions to Algebra and Geometry, 56 (2015), pp. 593640.
Abstract
We consider nondegenerate centroaffine hypersurface immersions in Rn whose cubic form is parallel with respect to the LeviCivita connection of the affine metric. There exists a bijective correspondence between homothetic families of proper affine hyperspheres with center in the origin and with parallel cubic form, and Köchers conic ? domains, which are the maximal connected sets consisting of invertible elements in a real semisimple Jordan algebra. Every level surface of the ? function in an ? domain is an affine complete, Euclidean complete proper affine hypersphere with parallel cubic form and with center in the origin. On the other hand, every proper affine hypersphere with parallel cubic form and with center in the origin can be represented as such a level surface. We provide a complete classification of proper affine hyperspheres with parallel cubic form based on the classification of semisimple real Jordan algebras. Centroaffine hypersurface immersions with parallel cubic form are related to the wider class of real unital Jordan algebras. Every such immersion can be extended to an affine complete one, whose conic hull is the connected component of the unit element in the set of invertible elements in a real unital Jordan algebra. Our approach can be used to study also other classes of hypersurfaces with parallel cubic form. 
K. Schildknecht, S. Olek, Th. Dickhaus, Simultaneous statistical inference for epigenetic data, PLOS ONE, 10 (2015), pp. e0125587/1e0125587/15.
Abstract
Epigenetic research leads to complex data structures. Since parametric model assumptions for the distribution of epigenetic data are hard to verify we introduce in the present work a nonparametric statistical framework for twogroup comparisons. Furthermore, epigenetic analyses are often performed at various genetic loci simultaneously. Hence, in order to be able to draw valid conclusions for specific loci, an appropriate multiple testing correction is necessary. Finally, with technologies available for the simultaneous assessment of many interrelated biological parameters (such as gene arrays), statistical approaches also need to deal with a possibly unknown dependency structure in the data. Our statistical approach to the nonparametric comparison of two samples with independent multivariate observables is based on recently developed multivariate multiple permutation tests. We adapt their theory in order to cope with families of hypotheses regarding relative effects. Our results indicate that the multivariate multiple permutation test keeps the preassigned type I error level for the global null hypothesis. In combination with the closure principle, the familywise error rate for the simultaneous test of the corresponding locus/parameterspecific null hypotheses can be controlled. In applications we demonstrate that group differences in epigenetic data can be detected reliably with our methodology. 
K. Tabelow, S. Mohammadi, N. Weiskopf, J. Polzehl, POAS4SPM  A toolbox for SPM to denoise diffusion MRI data, Neuroinformatics, 13 (2015), pp. 1929.
Abstract
We present an implementation of a recently developed noise reduction algorithm for dMRI data, called multishell position orientation adaptive smoothing (msPOAS), as a toolbox for SPM. The method intrinsically adapts to the structures of different size and shape in dMRI and hence avoids blurring typically observed in nonadaptive smoothing. We give examples for the usage of the toolbox and explain the determination of experimentdependent parameters for an optimal performance of msPOAS. 
K. Tabelow, H.U. Voss, J. Polzehl, Local estimation of the noise level in MRI using structural adaptation, Medical Image Analysis, 20 (2015), pp. 7686.
Abstract
We present a method for local estimation of the signaldependent noise level in magnetic resonance images. The procedure uses a multiscale approach to adaptively infer on local neighborhoods with similar data distribution. It exploits a maximumlikelihood estimator for the local noise level. The validity of the method was evaluated on repeated diffusion data of a phantom and simulated data using T1data corrupted with artificial noise. Simulation results are compared with a recently proposed estimate. The method was applied to a highresolution diffusion dataset to obtain improved diffusion model estimation results and to demonstrate its usefulness in methods for enhancing diffusion data. 
Y. Chen, V. Spokoiny, Modeling nonstationary and leptokurtic financial time series, Econometric Theory, (2015), pp. 703728.

T. Bodnar, Th. Dickhaus, On the Simes inequality in elliptical models, Annals of the Institute of Statistical Mathematics, (published online on Sept. 4, 2015), pp. 116, DOI 10.1007/s1046301505394 .
Abstract
We provide necessary and sufficient conditions for the validity of the inequality of Simes (1986) in models with elliptical dependencies. Necessary conditions are presented in terms of sufficient conditions for the reverse Simes inequality. One application of our main results concerns the problem of model misspecification, in particular the case that the assumption of Gaussianity of test statistics is violated. Since our sufficient conditions require nonnegativity of correlation coefficients between test statistics, we also develop exact tests for vectors of correlation coefficients. 
A. Gasnikov, Y. Nesterov, V. Spokoiny, On the efficiency of a randomized mirror descent algorithm in online optimization problems, Computational Mathematics and Mathematical Physics, 55 (2015), pp. 580596.

V.G. Gitis, A.B. Derendyaev, S.A. Pirogov, V. Spokoiny, E.F. Yurkov, Adaptive estimation of seismic parameter fields from earthquakes catalogs, Journal of Communications Technology and Electronics, 60 (2015), pp. 14591465.

S. Mohammadi, K. Tabelow, L. Ruthotto, Th. Feiweier, J. Polzehl, N. Weiskopf, Highresolution diffusion kurtosis imaging at 3T enabled by advanced postprocessing, Frontiers in Neuroscience, 8 (2015), pp. 427/1427/14.

TH. Dickhaus, Th. Royen, On multivariate chisquare distributions and their applications in testing multiple hypotheses, Statistics. A Journal of Theoretical and Applied Statistics, (2015), pp. 427454.
Abstract
We are considered with three different types of multivariate chisquare distributions. Their members play important roles as limiting distributions of vectors of test statistics in several applications of multiple hypotheses testing. We explain these applications and provide formulas for computing multiplicityadjusted $p$values under the respective global hypothesis. 
V. Spokoiny, M. Zhilova, Bootstrap confidence sets under a model misspecification, The Annals of Statistics, 43 (2015), pp. 26532675.
Abstract
A multiplier bootstrap procedure for construction of likelihoodbased confidence sets is considered for finite samples and possible model misspecification. Theoretical results justify the bootstrap consistency for small or moderate sample size and allow to control the impact of the parameter dimension: the bootstrap approximation works if the ratio of cube of the parameter dimension to the sample size is small. The main result about bootstrap consistency continues to apply even if the underlying parametric model is misspecified under the so called Small Modeling Bias condition. In the case when the true model deviates significantly from the considered parametric family, the bootstrap procedure is still applicable but it becomes a bit conservative: the size of the constructed confidence sets is increased by the modeling bias. We illustrate the results with numerical examples of misspecified constant and logistic regressions. 
A. Andresen, V. Spokoiny, Critical dimension in profile semiparametric estimation, Electronic Journal of Statistics, 8 (2014), pp. 30773125.

S. Becker, K. Tabelow, S. Mohammadi, N. Weiskopf, J. Polzehl, Adaptive smoothing of multishell diffusionweighted magnetic resonance data by msPOAS, NeuroImage, 95 (2014), pp. 90105.
Abstract
In this article we present a noise reduction method (msPOAS) for multishell diffusionweighted magnetic resonance data. To our knowledge, this is the first smoothing method which allows simultaneous smoothing of all qshells. It is applied directly to the diffusion weighted data and consequently allows subsequent analysis by any model. Due to its adaptivity, the procedure avoids blurring of the inherent structures and preserves discontinuities. MsPOAS extends the recently developed positionorientation adaptive smoothing (POAS) procedure to multishell experiments. At the same time it considerably simplifies and accelerates the calculations. The behavior of the algorithm msPOAS is evaluated on diffusionweighted data measured on a single shell and on multiple shells. 
D. Belomestny, V. Spokoiny, Concentration inequalities for smooth random fields, Theory of Probability and its Applications, 58 (2014), pp. 314323.
Abstract
In this note we derive a sharp concentration inequality for the supremum of a smooth random field over a finite dimensional set. It is shown that this supremum can be bounded with high probability by the value of the field at some deterministic point plus an intrinsic dimension of the optimisation problem. As an application we prove the exponential inequality for a function of the maximal eigenvalue of a random matrix is proved. 
G. Blanchard, Th. Dickhaus, E. Roquain, F. Villers, On least favorable configurations for stepupdown tests, Statistica Sinica, 24 (2014), pp. 123.

T. Bodnar, Th. Dickhaus, False discovery rate control under Archimedean copula, Electronic Journal of Statistics, 8 (2014), pp. 22072241.

S. Becker, P. Mathé, A different perspective on the PropagationSeparation approach, Electronic Journal of Statistics, 7 (2013), pp. 27022736.
Abstract
The PropagationSeparation approach is an iterative procedure for pointwise estimation of local constant and local polynomial functions. The estimator is defined as a weighted mean of the observations with datadriven weights. Within homogeneous regions it ensures a similar behavior as nonadaptive smoothing (propagation), while avoiding smoothing among distinct regions (separation). In order to enable a proof of stability of estimates, the authors of the original study introduced an additional memory step aggregating the estimators of the successive iteration steps. Here, we study theoretical properties of the simplified algorithm, where the memory step is omitted. In particular, we introduce a new strategy for the choice of the adaptation parameter yielding propagation and stability for local constant functions with sharp discontinuities. 
M. Zhilova, V. Spokoiny, Uniform properties of the local maximum likelihood estimate, Automation and Remote Control, 74 (2013), pp. 16561669.

M. Welvaert, K. Tabelow, R. Seurinck, Y. Rosseel, Adaptive smoothing as inference strategy: More specificity for unequally sized or neighboring regions, Neuroinformatics, 11 (2013), pp. 435445.
Abstract
Although spatial smoothing of fMRI data can serve multiple purposes, increasing the sensitivity of activation detection is probably its greatest benefit. However, this increased detection power comes with a loss of specificity when nonadaptive smoothing (i.e. the standard in most software packages) is used. Simulation studies and analysis of experimental data was performed using the R packages neuRosim and fmri. In these studies, we systematically investigated the effect of spatial smoothing on the power and number of false positives in two particular cases that are often encountered in fMRI research: (1) Single condition activation detection for regions that differ in size, and (2) multiple condition activation detection for neighbouring regions. Our results demonstrate that adaptive smoothing is superior in both cases because less false positives are introduced by the spatial smoothing process compared to standard Gaussian smoothing or FDR inference of unsmoothed data. 
D. Belomestny, V. Spokoiny, Concentration inequalities for smooth random fields, Theory of Probability and its Applications, (2013), pp. 401410.

R. Bruggemann, H.J. Mucha, H.G. Bartel, Ranking of polluted regions in south west Germany based on a multiindicator system, MATCH. Communications in Mathematical and in Computer Chemistry, 69 (2013), pp. 433462.

E. Diederichs, A. Juditsky, A. Nemirovski, V. Spokoiny, Sparse non Gaussian component analysis by semidefinite programming, Journal of Machine Learning Research (JMLR). MIT Press, Cambridge, MA. English, English abstracts., (2013), pp. 128.

F. Gach, R. Nickl, V. Spokoiny, Spatially adaptive density estimation by localised Haar projections, Annales de l'Institut Henri Poincare. Probabilites et Statistiques, 49 (2013), pp. 900914.

TH. Dickhaus, B. Blankertz, F.C. Meinecke, Binary classification with pFDRpFNR losses, Biometrical Journal, 55 (2013), pp. 463477.

TH. Dickhaus, J. Stange, Multiple point hypothesis test problems and effective numbers of tests for control of the familywise error rate, Calcutta Statistical Association Bulletin, 65 (2013), pp. 123144.

TH. Dickhaus, Randomized pvalues for multiple testing of composite null hypotheses, Journal of Statistical Planning and Inference, 143 (2013), pp. 19681979.

V. Spokoiny, M. Zhilova, Sharp deviation bounds for quadratic forms, Mathematical Methods of Statistics, 22 (2013), pp. 100113.

V. Spokoiny, W. Wang, W. Härdle, Local quantile regression (with rejoinder), Journal of Statistical Planning and Inference, 143 (2013), pp. 11091129.

S. Becker, K. Tabelow, H.U. Voss, A. Anwander, R.M. Heidemann, J. Polzehl, Positionorientation adaptive smoothing of diffusion weighted magnetic resonance data (POAS), Medical Image Analysis, 16 (2012), pp. 11421155.
Abstract
We introduce an algorithm for diffusion weighted magnetic resonance imaging data enhancement based on structural adaptive smoothing in both space and diffusion direction. The method, called POAS, does not refer to a specific model for the data, like the diffusion tensor or higher order models. It works by embedding the measurement space into a space with defined metric and group operations, in this case the Lie group of threedimensional Euclidean motion SE(3). Subsequently, pairwise comparisons of the values of the diffusion weighted signal are used for adaptation. The positionorientation adaptive smoothing preserves the edges of the observed fine and anisotropic structures. The POASalgorithm is designed to reduce noise directly in the diffusion weighted images and consequently also to reduce bias and variability of quantities derived from the data for specific models. We evaluate the algorithm on simulated and experimental data and demonstrate that it can be used to reduce the number of applied diffusion gradients and hence acquisition time while achieving similar quality of data, or to improve the quality of data acquired in a clinically feasible scan time setting. 
K. Tabelow, H.U. Voss, J. Polzehl, Modeling the orientation distribution function by mixtures of angular central Gaussian distributions, Journal of Neuroscience Methods, 203 (2012), pp. 200211.
Abstract
In this paper we develop a tensor mixture model for diffusion weighted imaging data using an automatic model selection criterion for the order of tensor components in a voxel. We show that the weighted orientation distribution function for this model can be expanded into a mixture of angular central Gaussian distributions. We show properties of this model in extensive simulations and in a high angular resolution experimental data set. The results suggest that the model may improve imaging of cerebral fiber tracts. We demonstrate how inference on canonical model parameters may give rise to new clinical applications. 
H.G. Bartel, H.J. Mucha, Finding incomparable pairs of subsets by using formal concept analysis, Statistica & Applicazioni, Special Issue (2011), pp. 6179.

S. Arlot, G. Blanchard, E. Roquain, Some nonasymptotic results on resampling in high dimension, I: Confidence regions, The Annals of Statistics, 38 (2010), pp. 5182.

S. Arlot, G. Blanchard, E. Roquain, Some nonasymptotic results on resampling in high dimension, II: Multiple tests, The Annals of Statistics, 38 (2010), pp. 8399.

Y. Chen, W. Härdle, V. Spokoiny, GHICA  Risk analysis with GH distributions and independent components, Journal of Empirical Finance, 17 (2010), pp. 255269.

D. Belomestny, Spectral estimation of the fractional order of a Lévy process, The Annals of Statistics, 38 (2010), pp. 317351.

G. Blanchard, G. Lee, C. Scott, Semisupervised novelty detection, Journal of Machine Learning Research (JMLR). MIT Press, Cambridge, MA. English, English abstracts., 11 (2010), pp. 29733009.
Abstract
A common setting for novelty detection assumes that labeled examples from the nominal class are available, but that labeled examples of novelties are unavailable. The standard (inductive) approach is to declare novelties where the nominal density is low, which reduces the problem to density level set estimation. In this paper, we consider the setting where an unlabeled and possibly contaminated sample is also available at learning time. We argue that novelty detection in this semisupervised setting is naturally solved by a general reduction to a binary classification problem. In particular, a detector with a desired false positive rate can be achieved through a reduction to NeymanPearson classification. Unlike the inductive approach, semisupervised novelty detection (SSND) yields detectors that are optimal (e.g., statistically consistent) regardless of the distribution on novelties. Therefore, in novelty detection, unlabeled data have a substantial impact on the theoretical properties of the decision rule. We validate the practical utility of SSND with an extensive experimental study. We also show that SSND provides distributionfree, learningtheoretic solutions to two well known problems in hypothesis testing. First, our results provide a general solution to the general twosample problem, that is, the problem of determining whether two random samples arise from the same distribution. Second, a specialization of SSND coincides with the standard $p$value approach to multiple testing under the socalled random effects model. Unlike standard rejection regions based on thresholded $p$values, the general SSND framework allows for adaptation to arbitrary alternative distributions. 
G. Blanchard, P. Mathé, Conjugate gradient regularization under general smoothness and noise assumptions, Journal of Inverse and IllPosed Problems, 18 (2010), pp. 701726.

G. Blanchard, Adaptive FDR control under independence and dependence, Journal of Machine Learning Research (JMLR). MIT Press, Cambridge, MA. English, English abstracts., 10 (2010), pp. 28372871.

E. Diederichs, A. Juditsky, V. Spokoiny, Ch. Schütte, Sparse nonGaussian component analysis, IEEE Transactions on Information Theory, 56 (2010), pp. 30333047.

J. Polzehl, H.U. Voss, K. Tabelow, Structural adaptive segmentation for statistical parametric mapping, NeuroImage, 52 (2010), pp. 515523.
Abstract
Functional Magnetic Resonance Imaging inherently involves noisy measurements and a severe multiple test problem. Smoothing is usually used to reduce the effective number of multiple comparisons and to locally integrate the signal and hence increase the signaltonoise ratio. Here, we provide a new structural adaptive segmentation algorithm (AS) that naturally combines the signal detection with noise reduction in one procedure. Moreover, the new method is closely related to a recently proposed structural adaptive smoothing algorithm and preserves shape and spatial extent of activation areas without blurring the borders. 
N. Serdyukova, Dependence on the dimension for complexity of approximation of random fields, Theory of Probability and its Applications, 54 (2010), pp. 272284.

K. Tabelow, V. Piëch, J. Polzehl, H.U. Voss, Highresolution fMRI: Overcoming the signaltonoise problem, Journal of Neuroscience Methods, 178 (2009), pp. 357365.
Abstract
Increasing the spatial resolution in functional Magnetic Resonance Imaging (fMRI) inherently lowers the signaltonoise ratio (SNR). In order to still detect functionally significant activations in highresolution images, spatial smoothing of the data is required. However, conventional nonadaptive smoothing comes with a reduced effective resolution, foiling the benefit of the higher acquisition resolution. We show how our recently proposed structural adaptive smoothing procedure for functional MRI data can improve signal detection of highresolution fMRI experiments regardless of the lower SNR. The procedure is evaluated on human visual and sensorymotor mapping experiments. In these applications, the higher resolution could be fully utilized and highresolution experiments were outperforming normal resolution experiments by means of both statistical significance and information content. 
H.G. Bartel, J. Dolata, H.J. Mucha, Klassifikation von 613 Proben als Referenzen für die Herstellungsprovenienzen römischer Baukeramik im nördlichen Obergermanien, Mainzer Archäologische Zeitschrift, 8 (2009), pp. 5176.

P. Čížek, W. Härdle, V. Spokoiny, Adaptive pointwise estimation in timeinhomogeneous conditional heteroscedasticity models, The Econometrics Journal, 12 (2009), pp. 248271.

E. Giacomini, W. Härdle, V. Krätschmer, Dynamic semiparametric factor models in risk neutral density estimation, AStA. Advances in Statistical Analysis. A Journal of the German Statistical Society, 93 (2009), pp. 387402.

E. Giacomini, W. Härdle, V. Spokoiny, Inhomogeneous dependency modelling with time varying copulae, Journal of Business & Economic Statistics, 27 (2009), pp. 224234.
Abstract
Measuring dependence in a multivariate time series is tantamount to modelling its dynamic structure in space and time. In the context of a multivariate normally distributed time series, the evolution of the covariance (or correlation) matrix over time describes this dynamic. A wide variety of applications, though, requires a modelling framework different from the multivariate normal. In risk management the nonnormal behaviour of most financial time series calls for nonlinear dependency. The correct modelling of nongaussian dependencies is therefore a key issue in the analysis of multivariate time series. In this paper we use copulae functions with adaptively estimated time varying parameters for modelling the distribution of returns, free from the usual normality assumptions. Further, we apply copulae to estimation of ValueatRisk (VaR) of a portfolio and show its better performance over the RiskMetrics approach, a widely used methodology for VaR estimation 
Y. Golubev, V. Spokoiny, Exponential bounds for minimum contrast estimators, Electronic Journal of Statistics, 3 (2009), pp. 712746.

A. Schwaighöfer, T. Schröter, S. Mika, G. Blanchard, How wrong can we get? A review of machine learning approaches and error bars, Combinatorial Chemistry & High Throughput Screening, 12 (2009), pp. 453468.

D. Belomestny, G.N. Milstein, V. Spokoiny, Regression methods in pricing American and Bermudan options using consumption processes, Quantitative Finance, 9 (2009), pp. 315327.
Abstract
Here we develop methods for efficient pricing multidimensional discretetime American and Bermudan options by using regression based algorithms together with a new approach towards constructing upper bounds for the price of the option. Applying sample space with payoffs at the optimal stopping times, we propose sequential estimates for continuation values, values of the consumption process, and stopping times on the sample paths. The approach admits constructing both low and upper bounds for the price by Monte Carlo simulations. The methods are illustrated by pricing Bermudan swaptions and snowballs in the Libor market model. 
J. Polzehl, S. Sperlich, A note on structural adaptive dimension reduction, Journal of Statistical Computation and Simulation, 79 (2009), pp. 805818.

J. Polzehl, K. Tabelow, Structural adaptive smoothing in diffusion tensor imaging: The R package dti, Journal of Statistical Software, 31 (2009), pp. 124.
Abstract
Diffusion Weighted Imaging has become and will certainly continue to be an important tool in medical research and diagnostics. Data obtained with Diffusion Weighted Imaging are characterized by a high noise level. Thus, estimation of quantities like anisotropy indices or the main diffusion direction may be significantly compromised by noise in clinical or neuroscience applications. Here, we present a new package dti for R, which provides functions for the analysis of diffusion weighted data within the diffusion tensor model. This includes smoothing by a recently proposed structural adaptive smoothing procedure based on the PropagationSeparation approach in the context of the widely used Diffusion Tensor Model. We extend the procedure and show, how a correction for Rician bias can be incorporated. We use a heteroscedastic nonlinear regression model to estimate the diffusion tensor. The smoothing procedure naturally adapts to different structures of different size and thus avoids oversmoothing edges and fine structures. We illustrate the usage and capabilities of the package through some examples. 
V. Spokoiny, C. Vial, Parameter tuning in pointwise adaptation using a propagation approach, The Annals of Statistics, 37 (2009), pp. 27832807.
Abstract
This paper discusses the problem of adaptive estimating a univariate object like the value of a regression function at a given point or a linear functional in a linear inverse problem. We consider an adaptive procedure originated from Lepski (1990) which selects in a datadriven way one estimate out of a given class of estimates ordered by their variability. A serious problem with using this and similar procedures is the choice of some tuning parameters like thresholds. Numerical results show that the theoretically recommended proposals appear to be too conservative and lead to a strong oversmoothing effects. A careful choice of the parameters of the procedure is extremely important for getting the reasonable quality of estimation. The main contribution of this paper is the new approach for choosing the parameters of the procedure by providing the prescribed behavior of the resulting estimate in the simple parametric situation. We establish a nonasymptotical “oracle” bound which shows that the estimation risk is, up to a logarithmic multiplier, equal to the risk of the “oracle” estimate which is optimally selected from the given family. A numerical study demonstrates the nice performance of the resulting procedure in a number of simulated examples. 
V. Spokoiny, Multiscale local change point detection with applications to ValueatRisk, The Annals of Statistics, 37 (2009), pp. 14051436.
Abstract
This paper offers a new procedure for nonparametric estimation and forecasting of time series with applications to volatility modeling for financial data. The approach is based on the assumption of local homogeneity: for every time point there exists a historical emphinterval of homogeneity, in which the volatility parameter can be well approximated by a constant. The procedure recovers this interval from the data using the local change point (LCP) analysis. Afterwards the estimate of the volatility can be simply obtained by local averaging. The approach carefully addresses the question of choosing the tuning parameters of the procedure using the so called “propagation” condition. The main result claims a new “oracle” inequality in terms of the modeling bias which measures the quality of the local constant approximation. This result yields the optimal rate of estimation for smooth and piecewise constant volatility functions. Then the new procedure is applied to some data sets and a comparison with a standard GARCH model is also provided. Finally we discuss applications of the new method to the Value at Risk problem. The numerical results demonstrate a very reasonable performance of the new method. 
K. Tabelow, J. Polzehl, A.M. Uluğ, J.P. Dyke, R. Watts, L.A. Heier, H.U. Voss, Accurate localization of brain activity in presurgical fMRI by structure adaptive smoothing, IEEE Transactions on Medical Imaging, 27 (2008), pp. 531537.
Abstract
An important problem of the analysis of fMRI experiments is to achieve some noise reduction of the data without blurring the shape of the activation areas. As a novel solution to this problem, the PropagationSeparation approach (PS), a structure adaptive smoothing method, has been proposed recently. PS adapts to different shapes of activation areas by generating a spatial structure corresponding to similarities and differences between time series in adjacent locations. In this paper we demonstrate how this method results in more accurate localization of brain activity. First, it is shown in numerical simulations that PS is superior over Gaussian smoothing with respect to the accurate description of the shape of activation clusters and and results in less false detections. Second, in a study of 37 presurgical planning cases we found that PS and Gaussian smoothing often yield different results, and we present examples showing aspects of the superiority of PS as applied to presurgical planning. 
K. Tabelow, J. Polzehl, V. Spokoiny, H.U. Voss, Diffusion tensor imaging: Structural adaptive smoothing, NeuroImage, 39 (2008), pp. 17631773.
Abstract
Diffusion Tensor Imaging (DTI) data is characterized by a high noise level. Thus, estimation errors of quantities like anisotropy indices or the main diffusion direction used for fiber tracking are relatively large and may significantly confound the accuracy of DTI in clinical or neuroscience applications. Besides pulse sequence optimization, noise reduction by smoothing the data can be pursued as a complementary approach to increase the accuracy of DTI. Here, we suggest an anisotropic structural adaptive smoothing procedure, which is based on the PropagationSeparation method and preserves the structures seen in DTI and their different sizes and shapes. It is applied to artificial phantom data and a brain scan. We show that this method significantly improves the quality of the estimate of the diffusion tensor and hence enables one either to reduce the number of scans or to enhance the input for subsequent analysis such as fiber tracking. 
H.G. Bartel, H.J. Mucha, J. Dolata, Über Identifikationsmethoden, dargestellt am Beispiel römischer Baukeramik, Berliner Beiträge zur Archäometrie, 21 (2008), pp. 115132.

A. Dalalyan, A. Juditsky, V. Spokoiny, A new algorithm for estimating the effective dimensionreduction subspace, Journal of Machine Learning Research (JMLR). MIT Press, Cambridge, MA. English, English abstracts., 9 (2008), pp. 16471678.

D. Divine, J. Polzehl, F. Godtliebsen, A propagationseparation approach to estimate the autocorrelation in a timeseries, Nonlinear Processes in Geophysics, 15 (2008), pp. 591599.

I.G. Grama, V. Spokoiny, Statistics of extremes by oracle estimation, The Annals of Statistics, 36 (2008), pp. 16191648.
Abstract
We use the fitted Pareto law to construct an accompanying approximation of the excess distribution function. A selection rule of the location of the excess distribution function is proposed based on a stagewise lackoffit testing procedure. Our main result is an oracle type inequality for the KullbackLeibler loss of the obtained adaptive estimator. 
V. Katkovnik, V. Spokoiny, Spatially adaptive estimation via fitted local likelihood techniques, IEEE Transactions on Signal Processing, 56 (2008), pp. 873886.
Abstract
This paper offers a new technique for spatially adaptive estimation. The local likelihood is exploited for nonparametric modelling of observations and estimated signals. The approach is based on the assumption of a local homogeneity of the signal: for every point there exists a neighborhood in which the signal can be well approximated by a constant. The fitted local likelihood statistics is used for selection of an adaptive size of this neighborhood. The algorithm is developed for quite a general class of observations subject to the exponential distribution. The estimated signal can be uni and multivariable. We demonstrate a good performance of the new algorithm for Poissonian image denoising and compare of the new method versus the intersection of confidence interval (ICI) technique that also exploits a selection of an adaptive neighborhood for estimation. 
Y. Chen, W. Härdle, V. Spokoiny, Portfolio value at risk based on independent components analysis, Journal of Computational and Applied Mathematics, 205 (2007), pp. 594607.
Abstract
Risk management technology applied to highdimensional portfolios needs simple and fast methods for calculation of value at risk (VaR). The multivariate normal framework provides a simple offtheshelf methodology but lacks the heavytailed distributional properties that are observed in data. A principle componentbased method (tied closely to the elliptical structure of the distribution) is therefore expected to be unsatisfactory. Here, we propose and analyze a technology that is based on independent component analysis (ICA). We study the proposed ICVaR methodology in an extensive simulation study and apply it to a highdimensional portfolio situation. Our analysis yields very accurate VaRs. 
H.U. Voss, K. Tabelow, J. Polzehl, O. Tchernichovski, K. Maul, D. SalgadoCommissariat, D. Ballon, S.A. Helekar, Functional MRI of the zebra finch brain during song stimulation suggests a lateralized response topography, Proceedings of the National Academy of Sciences of the United States of America, 104 (2007), pp. 1066710672.
Abstract
Electrophysiological and activitydependent gene expression studies of birdsong have contributed to the understanding of the neural representation of natural sounds. However, we have limited knowledge about the overall spatial topography of song representation in the avian brain. Here, we adapt the noninvasive functional MRI method in mildly sedated zebra finches (Taeniopygia guttata) to localize and characterize song driven brain activation. Based on the blood oxygenation leveldependent signal, we observed a differential topographic responsiveness to playback of bird's own song, tutor song, conspecific song, and a pure tone as a nonsong stimulus. The bird's own song caused a stronger response than the tutor song or tone in higher auditory areas. This effect was more pronounced in the medial parts of the forebrain. We found leftright hemispheric asymmetry in sensory responses to songs, with significant discrimination between stimuli observed only in the right hemisphere. This finding suggests that perceptual responses might be lateralized in zebra finches. In addition to establishing the feasibility of functional MRI in sedated songbirds, our results demonstrate spatial coding of song in the zebra finch forebrain, based on developmental familiarity and experience. 
D. Belomestny, V. Spokoiny, Spatial aggregation of local likelihood estimates with applications to classification, The Annals of Statistics, 35 (2007), pp. 22872311.
Abstract
This paper presents a new method for spatially adaptive local (constant) likelihood estimation which applies to a broad class of nonparametric models, including the Gaussian, Poisson and binary response models. The main idea of the method is given a sequence of local likelihood estimates (”weak” estimates), to construct a new aggregated estimate whose pointwise risk is of order of the smallest risk among all “weak” estimates. We also propose a new approach towards selecting the parameters of the procedure by providing the prescribed behavior of the resulting estimate in the simple parametric situation. We establish a number of important theoretical results concerning the optimality of the aggregated estimate. In particular, our “oracle” results claims that its risk is up to some logarithmic multiplier equal to the smallest risk for the given family of estimates. The performance of the procedure is illustrated by application to the classification problem. A numerical study demonstrates its nice performance in simulated and real life examples. 
J. Polzehl, K. Tabelow, Adaptive smoothing of digital images: The R package adimpro, Journal of Statistical Software, 19 (2007), pp. 117.
Abstract
Digital imaging has become omnipresent in the past years with a bulk of applications ranging from medical imaging to photography. When pushing the limits of resolution and sensitivity noise has ever been a major issue. However, commonly used nonadaptive filters can do noise reduction at the cost of a reduced effective spatial resolution only. Here we present a new package adimpro for R, which implements the PropagationSeparation approach by Polzehl and Spokoiny (2006) for smoothing digital images. This method naturally adapts to different structures of different size in the image and thus avoids oversmoothing edges and fine structures. We extend the method for imaging data with spatial correlation. Furthermore we show how the estimation of the dependence between variance and mean value can be included. We illustrate the use of the package through some examples. 
J. Polzehl, K. Tabelow, fmri: A package for analyzing fmri data, Newsletter of the R Project for Statistical Computing, 7 (2007), pp. 1317.

K. Tabelow, J. Polzehl, H.U. Voss, V. Spokoiny, Analyzing fMRI experiments with structural adaptive smoothing procedures, NeuroImage, 33 (2006), pp. 5562.
Abstract
Data from functional magnetic resonance imaging (fMRI) consists of time series of brain images which are characterized by a low signaltonoise ratio. In order to reduce noise and to improve signal detection the fMRI data is spatially smoothed. However, the common application of a Gaussian filter does this at the cost of loss of information on spatial extent and shape of the activation area. We suggest to use the propagationseparation procedures introduced by Polzehl and Spokoiny (2006) instead. We show that this significantly improves the information on the spatial extent and shape of the activation region with similar results for the noise reduction. To complete the statistical analysis, signal detection is based on thresholds defined by random field theory. Effects of ad aptive and nonadaptive smoothing are illustrated by artificial examples and an analysis of experimental data. 
G. Blanchard, M. Kawanabe, M. Sugiyama, V. Spokoiny, K.R. Müller, In search of nonGaussian components of a highdimensional distribution, Journal of Machine Learning Research (JMLR). MIT Press, Cambridge, MA. English, English abstracts., 7 (2006), pp. 247282.
Abstract
Finding nonGaussian components of highdimensional data is an important preprocessing step for efficient information processing. This article proposes a new em linear method to identify the “nonGaussian subspace” within a very general semiparametric framework. Our proposed method, called NGCA (NonGaussian Component Analysis), is essentially based on the fact that we can construct a linear operator which, to any arbitrary nonlinear (smooth) function, associates a vector which belongs to the low dimensional nonGaussian target subspace up to an estimation error. By applying this operator to a family of different nonlinear functions, one obtains a family of different vectors lying in a vicinity of the target space. As a final step, the target space itself is estimated by applying PCA to this family of vectors. We show that this procedure is consistent in the sense that the estimaton error tends to zero at a parametric rate, uniformly over the family. Numerical examples demonstrate the usefulness of our method. 
A. Goldenshluger, V. Spokoiny, Recovering convex edges of image from noisy tomographic data, IEEE Transactions on Information Theory, 52 (2006), pp. 13221334.

J. Polzehl, V. Spokoiny, Propagationseparation approach for local likelihood estimation, Probability Theory and Related Fields, 135 (2006), pp. 335362.
Abstract
The paper presents a unified approach to local likelihood estimation for a broad class of nonparametric models, including, e.g., regression, density, Poisson and binary response models. The method extends the adaptive weights smoothing (AWS) procedure introduced by the authors [Adaptive weights smoothing with applications to image sequentation. J. R. Stat. Soc., Ser. B 62, 335354 (2000)] in the context of image denoising. The main idea of the method is to describe a greatest possible local neighborhood of every design point in which the local parametric assumption is justified by the data. The method is especially powerful for model functions having large homogeneous regions and sharp discontinuities. The performance of the proposed procedure is illustrated by numerical examples for density estimation and classification. We also establish some remarkable theoretical nonasymptotic results on properties of the new algorithm. This includes the “propagation” property which particularly yields the root$n$ consistency of the resulting estimate in the homogeneous case. We also state an “oracle” result which implies rate optimality of the estimate under usual smoothness conditions and a “separation” result which explains the sensitivity of the method to structural changes. 
A. Samarov, V. Spokoiny, C. Vial, Component identification and estimation in nonlinear highdimensional regression models by structural adaptation, Journal of the American Statistical Association, 100 (2005), pp. 429445.

M. Giurcanu, V. Spokoiny, Confidence estimation of the covariance function of stationary and locally stationary processes, Statistics & Decisions. International Journal for Statistical Theory and Related Fields, 22 (2004), pp. 283300.

A. Goldenshluger, V. Spokoiny, On the shapefrommoments problem and recovering edges from noisy Radon data, Probability Theory and Related Fields, 128 (2004), pp. 123140.

G.N. Milstein, J.G.M. Schoenmakers, V. Spokoiny, Transition density estimation for stochastic differential equations via forwardreverse representations, Bernoulli. Official Journal of the Bernoulli Society for Mathematical Statistics and Probability, 10 (2004), pp. 281312.
Abstract
The general reverse diffusion equations are derived and applied to the problem of transition density estimation of diffusion processes between two fixed states. For this problem we propose density estimation based on forward?reverse representations and show that this method allows essentially better results to be achieved than the usual kernel or projection estimation based on forward representations only. 
J. Polzehl, S. Zwanzig, On a symmetrized simulation extrapolation estimator in linear errorsinvariables models, Computational Statistics & Data Analysis, 47 (2004), pp. 675688.
Abstract
Application of naive regression estimates in errorsinvariables models suffers from a severe bias. The simulation extrapolation estimator (SIMEX) was introduced by Cook and Stefanski as a correction method modeling the dependence of error variance in the regressors and bias of the regression method. Our symmetrized simulation extrapolation estimator (SYMEX), a generalization of SIMEX, allows to employ the symmetric structure of errorsinvariables models. Relations of both SIMEX and SYMEX to total least squares are investigated. 
V. Spokoiny, D. Mercurio, Statistical inference for timeinhomogeneous volatility models, The Annals of Statistics, 32 (2004), pp. 577602.

W. Härdle, H. Herwatz, V. Spokoiny, Time inhomogeneous multiple volatility modelling, Journal of Financial Econometrics, 1 (2003), pp. 5595.

D. Belomestny, V. Jentsch, M. Schreckenberg, Completion and continuation of nonlinear traffic time series: A probabilistic approach, Journal of Physics. A. Mathematical and General, 36 (2003), pp. 1136911383.
Abstract
When dealing with nonlinear time series of car traffic on highways, one of the outstanding problems to be solved is completion and continuation of data in space and time. To this end, the underlying process is decomposed into stochastic and deterministic components. The former is approximated by Gaussian white noise, while the latter refers, apart from always existing trends, to the space and timedependent jam propagation process. Jams are modelled in terms of dynamical Bayesian networks with radial basis functions involved. The models developed are used to tackle travel time estimation and prediction. Results are obtained for one of the most crowded traffic areas of Europe, namely the ringlike highway around Cologne. 
J. Polzehl, V. Spokoiny, Image denoising: Pointwise adaptive approach, The Annals of Statistics, 31 (2003), pp. 3057.
Abstract
A new method of pointwise adaptation has been proposed and studied in Spokoiny (1998) in context of estimation of piecewise smooth univariate functions. The present paper extends that method to estimation of bivariate greyscale images composed of large homogeneous regions with smooth edges and observed with noise on a gridded design. The proposed estimator $, hatf(x) ,$ at a point $, x ,$ is simply the average of observations over a window $, hatU(x) ,$ selected in a datadriven way. The theoretical properties of the procedure are studied for the case of piecewise constant images. We present a nonasymptotic bound for the accuracy of estimation at a specific grid point $, x ,$ as a function of the number of pixel $n$, of the distance from the point of estimation to the closest boundary and of smoothness properties and orientation of this boundary. It is also shown that the proposed method provides a near optimal rate of estimation near edges and inside homogeneous regions. We briefly discuss algorithmic aspects and the complexity of the procedure. The numerical examples demonstrate a reasonable performance of the method and they are in agreement with the theoretical issues. An example from satellite (SAR) imaging illustrates the applicability of the method. 
J.L. Horowitz, V. Spokoiny, An adaptive, rateoptimal test of linearity for median regression models, Journal of the American Statistical Association, 97 (2002), pp. 822835.

R. Liptser, A.Y. Veretennikov, V. Spokoiny, FreidlinWentzell type moderate deviations for smooth processes, Markov Processes and Related Fields, 8 (2002), pp. 611636.

V. Spokoiny, Variance estimation for highdimensional regression models, Journal of Multivariate Analysis, 82 (2002), pp. 111133.

L. Dümbgen, V. Spokoiny, Multiscale testing of qualitative hypotheses, The Annals of Statistics, 29 (2001), pp. 124152.

W. Härdle, S. Sperlich, V. Spokoiny, Structural tests for additive regression, Journal of the American Statistical Association, 96 (2001), pp. 13331347.

J.L. Horowitz, V. Spokoiny, An adaptive, rateoptimal test of a parametric meanregression model against a nonparametric alternative, Econometrica. Journal of the Econometric Society, 69 (2001), pp. 599631.

M. Hristache, A. Juditsky, J. Polzehl, V. Spokoiny, Structure adaptive approach for dimension reduction, The Annals of Statistics, 29 (2001), pp. 15371566.
Abstract
We propose a new method of effective dimension reduction for a multiindex model which is based on iterative improvement of the family of average derivative estimates. The procedure is computationally straightforward and does not require any prior information about the structure of the underlying model. We show that in the case when the effective dimension $m$ of the index space does not exceed $3$, this space can be estimated with the rate $n^1/2$ under rather mild assumptions on the model. 
M. Hristache, A. Juditsky, V. Spokoiny, Direct estimation of the index coefficient in a singleindex model, The Annals of Statistics, 29 (2001), pp. 595623.

S. Jaschke, U. Küchler, Coherent risk measures and gooddeal bounds, Statistics & Probability Letters, 5 (2001), pp. 181200.

J. Polzehl, B. Grund, Semiparametric lackoffit tests in an additive hazard regression model, Statistics and Computing, 11 (2001), pp. 323335.
Abstract
In the semiparametric additive hazard regression model of McKeague and Sasieni (Biometrika 81: 501514), the hazard contributions of some covariates are allowed to change over time, without parametric restrictions (Aalen model), while the contributions of other covariates are assumed to be constant. In this paper, we develop tests that help to decide which of the covariate contributions indeed change over time. The remaining covariates may be modelled with constant hazard coefficients, thus reducing the number of curves that have to be estimated nonparametrically. Several bootstrap tests are proposed. The behavior of the tests is investigated in a simulation study. In a practical example, the tests consistently identify covariates with constant and with changing hazard contributions. 
J. Polzehl, V. Spokoiny, Functional and dynamic Magnetic Resonance Imaging using vector adaptive weights smoothing, Journal of the Royal Statistical Society. Series C. Applied Statistics, 50 (2001), pp. 485501.
Abstract
We consider the problem of statistical inference for functional and dynamic Magnetic Resonance Imaging (MRI). A new approach is proposed which extends the adaptive weights smoothing (AWS) procedure from Polzehl and Spokoiny (2000) originally designed for image denoising. We demonstrate how the AWS method can be applied for time series of images, which typically occur in functional and dynamic MRI. It is shown how signal detection in functional MRI and analysis of dynamic MRI can benefit from spatially adaptive smoothing. The performance of the procedure is illustrated using real and simulated data. 
V. Spokoiny, Data driven testing the fit of linear models, Mathematical Methods of Statistics, 10 (2001), pp. 465497.

R. Liptser, V. Spokoiny, Deviation probability bound for martingales with applications to statistical estimation, Statistics & Probability Letters, 46 (2000), pp. 347357.

R. Liptser, V. Spokoiny, On estimating a dynamic function of a stochastic system with averaging, Statistical Inference for Stochastic Processes. An International Journal Devoted to Time Series Analysis and the Statistics of Continuous Time Processes and Dynamical Systems, 3 (2000), pp. 225249.

J. Polzehl, V. Spokoiny, Adaptive Weights Smoothing with applications to image restoration, Journal of the Royal Statistical Society. Series B. Statistical Methodology, 62 (2000), pp. 335354.
Abstract
We propose a new method of nonparametric estimation which is based on locally constant smoothing with an adaptive choice of weights for every pair of datapoints. Some theoretical properties of the procedure are investigated. Then we demonstrate the performance of the method on some simulated univariate and bivariate examples and compare it with other nonparametric methods. Finally we discuss applications of this procedure to magnetic resonance and satellite imaging. 
V. Spokoiny, Adaptive drift estimation for nonparametric diffusion model, The Annals of Statistics, 28 (2000), pp. 815836.
Contributions to Collected Editions

M. Kantner, L. Mertenskötter, DataDriven modeling of NonMarkovian noise in semiconductor lasers, in: 22nd International Conference on Numerical Simulation of Optoelectronic Devices (NUSOD), J. Piprek, P. Bardella, eds., IEEE xplorer, 2022, pp. 5758, DOI 10.1109/NUSOD54938.2022.9894788 .
Abstract
NonMarkovian noise degrades the coherence properties of semiconductor lasers and contributes significantly to broadening of the linewidth. Since modeling of such colored noise systems from first principles is not accessible, we aim for a datadriven modeling approach in which a system of stochastic rate equations shall be reconstructed from time series data. 
H.J. Mucha, H.G. Bartel, Distance and data transformation, in: The Encyclopedia of Archaeological Sciences, S.L. López Varela, ed., John Wiley & Sons, Inc., pp. published online on 26.11.2018, urlhttps://doi.org/10.1002/9781119188230.saseas0194, DOI 10.1002/9781119188230.saseas0194 .

H.J. Mucha, T.M. Gluhak, Finding groups in compositional data  Some experiments, in: Big data clustering: Data preprocessing, variable selection, and dimension reduction, H.J. Mucha, ed., WIAS Reports, 2017, pp. 97105.
Abstract
The talk is concerned with finding groups (clusters) in compositional data, that is nonnegative data with row sums (or column sums, respectively) equal to a constant, usually 1 in case of proportions or 100 in case of percentages. Without loss of generality, the cluster analysis of observations (row points) of compositional data is considered here, where the row profiles contains parts of some whole. Special distance functions between the profiles are proposed. Finally, applications to archaeometry are presented. 
N. Buzun, A. Suvorikova, V. Spokoiny, Multiscale parametric approach for change point detection, in: Proceedings of Information Technology and Systems 2016  The 40th Interdisciplinary Conference & School, Institute for Information Transmission Problems (Kharkevich Institute), Moscow, pp. 979996.

K. Tabelow, J. Polzehl, SHOWCASE 21  Towards invivo histology, in: MATHEON  Mathematics for Key Technologies, M. Grötschel, D. Hömberg, J. Sprekels, V. Mehrmann ET AL., eds., 1 of EMS Series in Industrial and Applied Mathematics, European Mathematical Society Publishing House, Zurich, 2014, pp. 378379.

H. Lamecker, H.Ch. Hege, K. Tabelow, J. Polzehl, F2  Image processing, in: MATHEON  Mathematics for Key Technologies, M. Grötschel, D. Hömberg, J. Sprekels, V. Mehrmann ET AL., eds., 1 of EMS Series in Industrial and Applied Mathematics, European Mathematical Society Publishing House, Zurich, 2014, pp. 359376.

H.J. Mucha, H.G. Bartel, J. Dolata, Dual scaling classification and its application to archaeometry, in: Data Analysis, Machine Learning and Knowledge Discovery, M. Spiliopoulou, L. SchmidtThieme, R. Janning, eds., 48 of Studies in Classification, Data Analysis, and Knowledge Organization, Springer, Heidelberg et al., 2014, pp. 105113.

H.J. Mucha, H.G. Bartel, Soft bootstrapping in cluster analysis and its comparison with other resampling methods, in: Data Analysis, Machine Learning and Knowledge Discovery, M. Spiliopoulou, L. SchmidtThieme, R. Janning, eds., 48 of Studies in Classification, Data Analysis, and Knowledge Organization, Springer, Heidelberg et al., 2014, pp. 97104.
Abstract
The bootstrap approach is resampling taken with replacement from the original data. Here we consider sampling from the empirical distribution of a given data set in order to investigate the stability of results of cluster analysis. Concretely, the original bootstrap technique can be formulated by choosing the following weights of observations: m i = n, if the corresponding object i is drawn n times, and m i = 0, otherwise. We call the weights of observations masses. In this paper, we present another bootstrap method, called soft bootstrapping, that consists of random change of the ?bootstrap masses? to some degree. Soft bootstrapping can be applied to any cluster analysis method that makes (directly or indirectly) use of weights of observations. This resampling scheme is especially appropriate for small sample sizes because no object is totally excluded from the soft bootstrap sample. At the end we compare different resampling techniques with respect to cluster analysis. 
B. Maldonado, C. MoralesMerino, H.J. Mucha, Ancient technology, modern science: Archaeometallurgy in northern Chile, in: Archäometrie und Denkmalpflege 2013, A. Hauptmann, O. Mecking, M. Prange, eds., Metalla, Sonderheft 6, Deutsches BergbauMuseum Bochum, 2013, pp. 7074.

H.J. Mucha, H.G. Bartel, C. MoralesMerino, Visualisation of cluster analysis results, in: Classification and Data Mining, A. Giusti, G. Ritter, M. Vichi, eds., Studies in Classification, Data Analysis, and Knowledge Organization, Springer, Berlin, 2013, pp. 261270.

TH. Dickhaus, J. Gierl, Simultaneous test procedures in terms of pvalue copulae, in: Proceedings of the 2nd Annual International Conference on Computational Mathematics, Computational Geometry & Statistics (CMCGS 2013), K. Chen, C. Raju, eds., Global Science and Technology Forum, Singapore, 2013, pp. 7580.

H.J. Mucha, J. Dolata, H.G. Bartel, Classification of Roman tiles with stamp PARDALIVS, in: Algorithms from & for Nature and Life, B. Lausen, D. VAN DEN Poel, A. Ultsch, eds., Studies in Classification, Data Analysis, and Knowledge Organization, Springer, Heidelberg, 2013, pp. 309317.

H.J. Mucha, R. Naumann, G. Schneider, H.G. Bartel, M. Daszkiewicz, Diskriminanzanalytische Überprüfung einer Klassifikation von Keramikfunden aus Uruk, in: Archäometrie und Denkmalpflege 2013, A. Hauptmann, O. Mecking, M. Prange, eds., Metalla, Sonderheft 6, Deutsches BergbauMuseum Bochum, 2013, pp. 159163.

S. Behrendt, H.J. Mucha, H.G. Bartel, D.P. Mielke, Zur Problematik der multivariaten statistischen Analyse umfangreicher PRFADatenmengen phönizischer Keramik, in: Archäometrie und Denkmalpflege 2012, F. Schlütter, S. Greif, M. Prange, eds., Metalla, Sonderheft 5, Deutsches BergbauMuseum Bochum, 2012, pp. 157159.

J. Dolata, H.J. Mucha, H.G. Bartel, Eine neue Referenzgruppe für die römische Ziegelproduktion der I. Thrakerkohorte am Mittelrhein, in: Archäometrie und Denkmalpflege 2012, F. Schlütter, S. Greif, M. Prange, eds., Metalla, Sonderheft 5, Deutsches BergbauMuseum Bochum, 2012, pp. 171173.

C. MoralesMerino, H.J. Mucha, H.G. Bartel, Multivariate statistical analysis of clay and ceramicdata for provenance of Bronze Age pottery from Troia, in: Archäometrie und Denkmalpflege 2012, F. Schlütter, S. Greif, M. Prange, eds., Metalla, Sonderheft 5, Deutsches BergbauMuseum Bochum, 2012, pp. 174176.

H.J. Mucha, H.G. Bartel, An accelerated Kmeans algorithm based on adaptive distances, in: Challenges at the Interface of Data Analysis, Computer Science, and Optimization  Proceedings of the 34th Annual Conference of the Gesellschaft für Klassifikation e. V., Karlsruhe, July 2123, 2010, W. Gaul, A. GeyerSchulz, L. SchmidtThieme, J. Kunze, eds., Studies in Classification, Data Analysis, and Knowledge Organization, Springer, Heidelberg, 2012, pp. 3747.

H.J. Mucha, Method selection in cluster analysis followed by builtin validation, in: Data Analysis Methods and its Applications, J. Pociecha, R. Decker, eds., Wydawnictwo C. H. Beck, Warszawa, 2012, pp. 125142.

H.G. Bartel, H.J. Mucha, J. Dolata, Geochemical and statistical investigation of Roman stamped tiles of the Legio XXI Rapax, in: Classification as a Tool for Research. 11th IFCS Biennial Conference and 33rd Annual Conference of the Gesellschaft für Klassifikation e. V., Dresden, March 1318, 2009, H. LocarekJunge, C. Weihs, eds., Studies in Classification, Data Analysis, and Knowledge Organization, Springer, Berlin/Heidelberg, 2010, pp. 427434.

J. Dolata, D. Lück, H.G. Bartel, H.J. Mucha, C. Wenzel, Analysen römischer Ziegel aus dem Kastellbad von GroßGerau: ihre statistische Auswertung und archäologische Bewertung, in: Archäometrie und Denkmalpflege 2010, O. Hahn, A. Hauptmann, D. ModarressiTehrani, M. Prange, eds., Metalla, Sonderheft 3, Deutsches BergbauMuseum Bochum, 2010, pp. 129131.

J. Dolata, H.J. Mucha, H.G. Bartel, Mathematischarchäometrische Auswertung der GISgestützten Fundkartierung römischer Ziegel aus Mainz, in: Archäometrie und Denkmalpflege 2010, O. Hahn, A. Hauptmann, D. ModarressiTehrani, M. Prange, eds., Metalla, Sonderheft 3, Deutsches BergbauMuseum Bochum, 2010, pp. 119121.

J. Dolata, H.J. Mucha, H.G. Bartel, Mapping findspots of Roman military brickstamps in em Mogontiacum (Mainz) and archaeometrical analysis, in: Advances in Data Analysis, Data Handling and Business Intelligence, A. Fink, B. Lausen, W. Seidel, A. Ultsch, eds., Studies in Classification, Data Analysis, and Knowledge Organization, Springer, Heidelberg, 2010, pp. 595603.

E. Haimerl, H.J. Mucha, Comparing the stability of clustering results of dialect data based on several distance matrices, in: Classification as a Tool for Research. Proceedings of the 11th IFCS Biennial Conference and 33rd Annual Conference of the Gesellschaft für Klassifikation e. V., Dresden, March 1318, 2009, H. LocarekJunge, C. Weihs, eds., Studies in Classification, Data Analysis, and Knowledge Organization, Springer, Berlin/Heidelberg, 2010, pp. 665672.

C. MoralesMerino, H.J. Mucha, H.G. Bartel, E. Pernicka, Clay sediments analysis in the Troad and its segmentation, in: Archäometrie und Denkmalpflege 2010, O. Hahn, A. Hauptmann, D. ModarressiTehrani, M. Prange, eds., Metalla, Sonderheft 3, Deutsches BergbauMuseum Bochum, 2010, pp. 122124.

G. Blanchard, N. Krämer, Kernel partial least squares is universally consistent, in: Proceedings of the 13th International Conference on Artificial Intelligence and Statistics (AISTATS 2010), Y.W. Teh, M. Titterington, eds., 9 of JMLR Workshop and Conference Proceedings, Journal of Machine Learning Research, Cambridge, MA, USA, 2010, pp. 5764.

P. Čížek, V. Spokoiny, Varying coefficient GARCH models, in: Handbook of Financial Time Series, T.G. Andersen, R.A. Davis, J.P. Kreiss, T. Mikosch, eds., Springer, Berlin, 2009, pp. 169185.

J. Dolata, H.G. Bartel, H.J. Mucha, Geochemische und statistische Erkundung der Herstellungsorte von Ziegeln der em Legio XXI Rapax, in: Oedenburg Fouilles Françaises, Allemandes et Suisses à Biesheim et Kunheim, HautRhin, France. Volume 1: Les Camps Militaires JulioClaudiens, M. Reddé, ed., 79 of Monographien des RömischGermanischen Zentralmuseums, Verlag des RömischGermanischen Zentralmuseums, Mainz, 2009, pp. 355364.

C. Scott, G. Blanchard, Novelty detection: Unlabeled data definitely help, in: Proceedings of the Twelfth International Conference on Artificial Intelligence and Statistics, April 1618, 2009, Clearwater Beach, Florida, USA, D. VAN Dyke, M. Welling, eds., 5 of JMLR Workshop and Conference Proceedings, JMLR, 2009, pp. 464471.

H.J. Mucha, H.G. Bartel, J. Dolata, Finding Roman brickyards in Germania Superior by modelbased cluster analysis of archaeometric data, in: Proceedings of the 35th International Conference on Computer Applications and Quantitative Methods in Archaeology (CAA 2007), A. Posluschny, K. Lambers, I. Herzog, eds., 10 of Kolloquien zur Vor und Frühgeschichte, Habelt Verlag, Bonn, 2009, pp. 360365.

H.J. Mucha, H.G. Bartel, J. Dolata, Zur Klassifikation römischer Ziegel von Fundorten im südlichen Obergermanien. I. Multivariate statistische Analyse der archäometrischen Daten, in: Archäometrie und Denkmalpflege 2009, Jahrestagung in der Pinakothek der Moderne, München, 25.28. März 2009, A. Hauptmann, H. Stege, eds., Metalla, Sonderheft 2, Deutsches BergbauMuseum Bochum, 2009, pp. 144146.

M. Elagin, V. Spokoiny, Locally time homogeneous time series modelling, in: Applied Quantitative Finance (2nd edition), W. Härdle, N. Hautsch, L. Overbeck, eds., Springer, Berlin, 2008, pp. 345362.
Preprints, Reports, Technical Reports

N. Nüsken, D.R.M. Renger, Stein variational gradient descent: Manyparticle and longtime asymptotics, Preprint no. 2819, WIAS, Berlin, 2021, DOI 10.20347/WIAS.PREPRINT.2819 .
Abstract, PDF (430 kByte)
Stein variational gradient descent (SVGD) refers to a class of methods for Bayesian inference based on interacting particle systems. In this paper, we consider the originally proposed deterministic dynamics as well as a stochastic variant, each of which represent one of the two main paradigms in Bayesian computational statistics: emphvariational inference and emphMarkov chain Monte Carlo. As it turns out, these are tightly linked through a correspondence between gradient flow structures and largedeviation principles rooted in statistical physics. To expose this relationship, we develop the cotangent space construction for the Stein geometry, prove its basic properties, and determine the largedeviation functional governing the manyparticle limit for the empirical measure. Moreover, we identify the emphSteinFisher information (or emphkernelised Stein discrepancy) as its leading order contribution in the longtime and manyparticle regime in the sense of $Gamma$convergence, shedding some light on the finiteparticle properties of SVGD. Finally, we establish a comparison principle between the SteinFisher information and RKHSnorms that might be of independent interest. 
V. Avanesov, Datadriven confidence bands for distributed nonparametric regression, Preprint no. 2729, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2729 .
Abstract, PDF (316 kByte)
Gaussian Process Regression and Kernel Ridge Regression are popular nonparametric regression approaches. Unfortunately, they suffer from high computational complexity rendering them inapplicable to the modern massive datasets. To that end a number of approximations have been suggested, some of them allowing for a distributed implementation. One of them is the divide and conquer approach, splitting the data into a number of partitions, obtaining the local estimates and finally averaging them. In this paper we suggest a novel computationally efficient fully datadriven algorithm, quantifying uncertainty of this method, yielding frequentist $L_2$confidence bands. We rigorously demonstrate validity of the algorithm. Another contribution of the paper is a minimaxoptimal highprobability bound for the averaged estimator, complementing and generalizing the known risk bounds. 
V. Avanesov, Nonparametric change point detection in regression, Preprint no. 2687, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2687 .
Abstract, PDF (329 kByte)
This paper considers the prominent problem of changepoint detection in regression. The study suggests a novel testing procedure featuring a fully datadriven calibration scheme. The method is essentially a black box, requiring no tuning from the practitioner. The approach is investigated from both theoretical and practical points of view. The theoretical study demonstrates proper control of firsttype error rate under H0 and power approaching 1 under H1. The experiments conducted on synthetic data fully support the theoretical claims. In conclusion, the method is applied to financial data, where it detects sensible changepoints. Techniques for changepoint localization are also suggested and investigated 
V. Avanesov, How to gamble with nonstationary Xarmed bandits and have no regrets, Preprint no. 2686, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2686 .
Abstract, PDF (287 kByte)
In Xarmed bandit problem an agent sequentially interacts with environment which yields a reward based on the vector input the agent provides. The agent's goal is to maximise the sum of these rewards across some number of time steps. The problem and its variations have been a subject of numerous studies, suggesting sublinear and sometimes optimal strategies. The given paper introduces a new variation of the problem. We consider an environment, which can abruptly change its behaviour an unknown number of times. To that end we propose a novel strategy and prove it attains sublinear cumulative regret. Moreover, the obtained regret bound matches the best known bound for GPUCB for a stationary case, and approaches the minimax lower bound in case of highly smooth relation between an action and the corresponding reward. The theoretical result is supported by experimental study.
Talks, Poster

L. Mertenskötter, M. Kantner, U. Bandelow, H. Wenzel, NonMarkovian noise in semiconductor lasers, MATH+ Day 2022, Berlin, November 18, 2022.

H. Kremer, J.J. Zhu, K. Muandet, B. Schölkopf, Functional generalized empirical likelihood estimation for conditional moment restrictions (spotlight, online talk), ICML 2022: 39th International Conference on Machine Learning (Online Event), July 18  23, 2022, Baltimore, USA, July 19, 2022.

A. Maltsi, Modelbased geometry reconstruction of TEM images, MATH+ Day 2021 (Online Event), Technische Universität Berlin, November 5, 2021.

A. Suvorikova, Optimal transport in machine learning (onlne talk), Seminar of the International Laboratory of Stochastic Algorithms and Highdimensional Inference (HDI Lab), National Research University  Higher School of Economics (HSE), Faculty of Computer Science (Online Event), Moscow, Russian Federation, April 29, 2021.

A. Suvorikova, Statistics for nonstatisticians (online talk), ITaS Interdisciplinary Conference 2021 (Online Event), November 15  17, 2021, Russian Academy of Sciences, Institute for Information Transmission Problems, Moscow, Russian Federation, November 16, 2021.

A. Suvorikova, Survey of methods of kmeans clustering with optimal transport (online talk), Third HSEYandex Autumn School on Generative Models (Hybrid Event), November 23  26, 2021, National Research University  Higher School of Economics (HSE), Yandex SDA Campus, Moscow, Russian Federation, November 26, 2021.

V. Avanesov, Datadriven confidence bands for distributed nonparametric regression (online talk), The 33rd Annual Conference on Learning Theory (COLT 2020) (Online Event), July 9  12, 2020, Graz, Austria, July 10, 2020.

A. Suvorikova, Change point detection in highdimensional data (online talk), Joint AramcoHSE Reserach Seminar, Higher School of Economics, Faculty of Computer Science, Moscow, Russian Federation, April 15, 2020.

P. Dvurechensky, Distributed optimization for Wasserstein barycenters (online talk), 15th INFORMS Telecommunication and Network Analytics Conference 2020 (Online Event), October 20  21, 2020, INFORMS Telecommunications & Networks Analytics, Catonsville, USA, October 21, 2020.

TH. Koprucki, K. Tabelow, T. Streckenbach, T. Niermann, A. Maltsi, Modelbased geometry reconstruction of TEM images, MATH+ Day 2020 (Online Event), Berlin, November 6, 2020.

F. Besold, Manifold clustering with adaptive weights, Structural Inference in HighDimensional Models 2, National Research University Higher School of Economics, HDILab, St. Petersburg, Russian Federation, August 26  30, 2019.

F. Besold, Manifold clustering with adaptive weights, Joint Workshop of BBDC, BZML and RIKEN AIP, Fraunhofer Institute HHI, September 9  10, 2019.

D. Dvinskikh, Decentralized and parallelized primal and dual accelerated methods, Structural Inference in HighDimensional Models 2, National Research University Higher School of Economics, HDILab, St. Petersburg, Russian Federation, August 26  30, 2019.

TH. Koprucki, On a database of simulated TEM images for In(Ga)As/GaAs quantum dots with various shapes, 19th International Conference on Numerical Simulation of Optoelectronic Devices (NUSOD 2019) , Session ``Nanostructures", July 8  12, 2019, University of Ottawa, Canada, July 8, 2019.

V. Spokoiny, Advanced statistical methods, April 9  11, 2019, Higher School of Economics, National Research University, Moscow, Russian Federation.

V. Spokoiny, Optimal stopping and control via reinforced regression, Optimization and Statistical Learning, March 25  28, 2019, Les Houches School of Physics, France, March 26, 2019.

A. Suvorikova, Statistical inference with optimal transport, Spring School ``Structural Inference 2018'' and Closing Workshop, FOR 1735 ``Structural Inference in Statistics'', Lübbenau, March 4  9, 2018.

P. Pigato, Estimation of piecewiseconstant coefficients in a stochastic differential equation, The 40th Conference on Stochastic Processes and their Applications (SPA 2018), University of Gothenburg, Göteborg, Sweden, June 13, 2018.

V. Spokoiny, Adaptive nonparametric clustering, Multiscale Problems in Materials Science and Biology: Analysis and Computation, January 8  12, 2018, Tsinghua University, Yau Mathematical Sciences Center, Sanya, Hainan, China, January 10, 2018.

V. Spokoiny, Bootstrap confidence sets for spectral projectors of sample covariance, 12th International Vilnius Conference on Probability Theory and Mathematical Statistics and 2018 IMS Annual Meeting on Probability and Statistics, July 2  6, 2018, Vilnius University, Lithuanian Mathematical Society and the Institute of Mathematical Statistics, Lithuania, July 5, 2018.

V. Spokoiny, Gaussian process forecast with multidimensional distributional input, 4th Conference of the International Society for Nonparametric Statistics, June 11  15, 2018, University of Salerno, Italy, June 15, 2018.

V. Spokoiny, Inference for spectral projectors, Workshop ``Statistical Inference for Structured Highdimensional Models'', March 11  16, 2018, Mathematisches Forschungsinstitut Oberwolfach, March 14, 2018.

V. Spokoiny, Large ball probability with applications in statistics, Mathematical Workshop of the School of Applied Mathematics and Computer Science, Moscow Institute of Physics and Technology, Dolgoprudny, Russian Federation, November 30, 2018.

V. Spokoiny, Large ball probability with applications to statistical inference, 6th Princeton Day of Statistics, Princeton University, Department of Operations Research and Financial Engineering, USA, November 9, 2018.

V. Spokoiny, Prior impact in Bayesian inference, International Seminar & Workshop: ``Stochastic Dynamics on Large Networks: Prediction and Inference'', October 15  16, 2018, MaxPlanckInstitut für Physik komplexer Systeme, Dresden, October 15, 2018.

V. Spokoiny, Structured nonparametric and highdimensional statistics, 2018 Joint Meeting of the Korean Mathematical Society and the German Mathematical Society, October 3  6, 2018, Korean Mathematical Society, Seoul, Korea (Republic of), October 5, 2018.

A. Suvorikova, Construction of confidence sets in 2Wasserstein space, Haindorf Seminar 2017, January 24  28, 2017, HumboldtUniversität zu Berlin, Wirtschaftswissenschaftliche Fakultät, Hejnice, Czech Republic, January 26, 2017.

A. Suvorikova, Construction of confidence sets in 2Wasserstein space, Machine Learning Seminar, Université PaulSabatier, Institut de Mathématiques de Toulouse, France, December 1, 2017.

A. Suvorikova, Construction of nonasymptotic confidence sets in 2Wasserstein space, Spring School ``Structural Inference'' 2017, Bad Malente, March 5  10, 2017.

P. Pigato, Estimation of the parameters of a diffusion with discontinuous coefficients, 7th Annual ERC BerlinOxford Young Researchers Meeting on Applied Stochastic Analysis, May 18  20, 2017, WIASBerlin, May 20, 2017.

A. Koziuk, Bootstrap for the regression problem with instrumental variables, Haindorf Seminar 2017, January 24  28, 2017, HumboldtUniversität zu Berlin, Wirtschaftswissenschaftliche Fakultät, Hejnice, Czech Republic, January 26, 2017.

H.J. Mucha, About the influence of multiple points of bootstrapping on the validation of clustering results, European Classification and Data Analysis Conference 2017 (ECDA2017), September 27  29, 2017, Wrocław University of Economics, Poland, September 27, 2017.

H.J. Mucha, Big data clustering: Comparison of the performance of a new fast preclustering and subsampling, German Polish Seminar on Data Analysis and Applications 2017, September 25  26, 2017, Wrocław University of Economics, Poland, September 26, 2017.

H.J. Mucha, Hierarchical clustering of big data using a special bootstrapping version, AG DANK Herbsttagung, November 17  18, 2017, Gesellschaft für KlassifikationGesellschaft für Klassifikation, Data Science Society, November 18, 2017.

A. Suvorikova, Bootstrap confidence sets for Wasserstein barycenters, Meeting in Mathematical Statistics 2016 ``Advances in Nonparametric and Highdimensional Statistics'', December 12  16, 2016, Fréjus, France, December 16, 2016.

A. Suvorikova, Bootstrap procedure in the space of Gaussian measures, Information Technology and Systems 2016, September 25  30, 2016, Russian Academy of Sciences, Institute for Information Transmission Problems, St. Petersburg.

A. Suvorikova, Multiscale change point detection, GeorgAugustUniversität Göttingen, Institut für Mathematische Stochastik, November 9, 2016.

CH. Bayer, The forwardreverse method for conditional Markov processes, Bayes in Paris, École Nationale de la Statistique et de l'Administration Économique, Paris, France, January 28, 2016.

V. Spokoiny, Adaptive weights clustering, Mathematisches Kolloquium, Universität Ulm, Institut für Analysis, November 18, 2016.

V. Spokoiny, Clustering using adaptive weights, Yandex, Moscow, Russian Federation, October 28, 2016.

V. Spokoiny, Inference for Structural Nonparametrics, Open Access: urlhttp://www.mathnet.ru/php/conference.phtml?option_lang=eng&eventID=25&confid=872, February 15  March 1, 2016, Independent University of Moscow (IUM), Russian Federation.

V. Spokoiny, Inference for structured regression, Meeting in Mathematical Statistics 2016 ``Advances in Nonparametric and Highdimensional Statistics'', December 12  16, 2016, Fréjus, France, December 15, 2016.

N. Buzun, Multiplier bootstrap with misspecified models, Statistics Mathematics and Applications, August 31  September 4, 2015, La Villa Clythia, Fréjus, France, September 1, 2015.

N. Buzun, Multiscale parametric approach for change point detection, Information Technologies and Systems 2015, September 6  11, 2015, Russian Academy of Sciences, Institute for Information Transmission Problems, Sochi, Russian Federation, September 9, 2015.

K. Schildknecht, Simultaneous statistical inference for epigenetic data, 61. Biometrischen Kolloquium in Dortmund, March 16  18, 2015, Technische Universität Dortmund, Fakultät Statistik, March 18, 2015.

J. Stange, On the correlation structure of test statistics in genetic casecontrol association studies, 61. Biometrischen Kolloquium in Dortmund, March 16  18, 2015, Technische Universität Dortmund, Fakultät Statistik, March 17, 2015.

CH. Bayer, The forwardreverse method for conditional Markov processes, MInisymposium ``Multilevel Monte Carlo Methods and Applications'' of the 8th International Congress on Industrial and Applied Mathematics (ICIAM2015), August 10  14, 2015, International Council for Industrial and Applied Mathematics, Beijing, China, August 10, 2015.

CH. Bayer, The forwardreverse method for conditional diffusion processes, Advances in Uncertainty Quantification Methods, Algorithms and Applications (UQAW 2015), January 6  9, 2015, King Abdullah University of Science and Technology (KAUST), Jeddah, Saudi Arabia, January 6, 2015.

CH. Bayer, The forwardreverse method for conditional diffusion processes, Oberseminar Finanz und Versicherungsmathematik, Technische Universität München, Fakultät für Mathematik, January 19, 2015.

CH. Bayer , The forwardreverse method for conditional Markov processes, International Conference on Scientific Computation and Differential Equations, Universität Potsdam, Institut für Mathematik, September 16, 2015.

V. Spokoiny, Bootstrap tuning in model selection, Lomonosov State University, Faculty of Computer Sciences, Moscow, Russian Federation, September 18, 2015.

V. Spokoiny, Predictive modeling: Methods and applications, Skolkov Institute of Science and Technology, Moscow, Russian Federation, August 27, 2015.

M. Zhilova, Nonasymptotic confidence bounds via multiplier bootstrap, CRC 649 "Economic Risk" Conference, July 10  12, 2014, HumboldtUniversität zu Berlin, Sonderforschungsbereich 649: Ökonomisches Risiko, July 11, 2014.

M. Zhilova, Parametric confidence sets using multiplier bootstrap, Premoday, February 21, 2014, Russian Academy of Sciences, Institute of Information Transmission Problems (Kharkevich Institute), Moscow, Russian Federation, February 21, 2014.

M. Zhilova, Uniform confidence bands for generalized regression via multiplier bootstrap, International Symposium on Financial Engineering and Risk Management 2014, June 27  28, 2014, Central University of Finance and Economics, Beijing, China, June 27, 2014.

M. Zhilova, Uniform confidence bands in generalized regression via multiplier bootstrap, The 3rd Institute of Mathematical Statistics Asia Pacific Rim Meeting, June 29  July 3, 2014, Taipei, Taiwan, Province Of China, July 1, 2014.

M. Zhilova, Uniform confidence sets via multiplier bootstrap, New Horizons in Statistical Decision Theory, September 7  13, 2014, Mathematisches Forschungsinstitut Oberwolfach, September 9, 2014.

A. Andresen, Finite sample analysis of Mestimators, International Workshop ``Advances in Optimization and Statistics'', May 15  16, 2014, Russian Academy of Sciences, Institute of Information Transmission Problems (Kharkevich Institute), Moscow, May 16, 2014.

A. Andresen, Finite sample analysis of maximum likelihood estimators, Kolloquium über Mathematische Statistik und Stochastische Prozesse, Universität Hamburg, April 22, 2014.

R. Hildebrand, Spectrahedral cones generated by rank 1 matrices, Workshop ``Real Algebraic Geometry With A View Toward Systems Control and Free Positivity'', April 6  12, 2014, Mathematisches Forschungsinstitut Oberwolfach, April 11, 2014.

K. Schildknecht, A twostage hierarchical multiple test procedure based on the asymptotically optimal rejection curve, 60. Biometrisches Kolloquium, Bremen, March 10  13, 2014.

K. Schildknecht, A twostage hierarchical multiple test procedure based on the asymptotically optimal rejection curve, Workshop Adaptive Designs and Multiple Testing Procedures, Internationalen Biometrischen Gesellschaft, Basel, Switzerland, June 5, 2014.

J. Stange, Computation of an effective number of simultaneous $chi^2$ tests, International Workshop ``Advances in Optimization and Statistics'', May 15  16, 2014, Russian Academy of Sciences, Institute of Information Transmission Problems (Kharkevich Institute), Moscow, May 15, 2014.

J. Stange, Multiple testing adjustments with copula models, 27th International Biometric Conference (IBC 2014), Florence, Italy, July 6  11, 2014.

J. Stange , Multiple testing adjustments with copula models, Haindorf Seminar 2014, February 6  8, 2014, HumboldtUniversitát zu Berlin, Wirtschaftswissenschaftliche Fakultät, Hejnice, Czech Republic, February 6, 2014.

N. Willrich, Model selection by Lepski's method for penalized likelihood, International Workshop ``Advances in Optimization and Statistics'', May 15  16, 2014, Russian Academy of Sciences, Institute of Information Transmission Problems (Kharkevich Institute), Moscow, May 16, 2014.

S. Mohammadi, L. Ruthotto, K. Tabelow, T. Feiweier, J. Polzehl, N. Weiskopf, ACID  A postprocessing toolbox for advanced diffusion MRI, 20th Annual Meeting of the Organization for Human Brain Mapping, Hamburg, June 8  12, 2014.

TH. Dickhaus, T. Bodnar, False discovery rate control under Archimedean copula, 27th International Biometric Conference (IBC 2014), July 6  11, 2014, International Biometric Society, Florence, Italy, July 7, 2014.

TH. Dickhaus, Randomized pvalues for multiple testing of composite null hypotheses, 60. Biometrisches Kolloquium, March 11  13, 2014, Universität Bremen, March 12, 2014.

TH. Dickhaus, Simultaneous Bayesian analysis of contingency tables in genetic association studies, Bayesian Biostatistics 2014, July 2  5, 2014, University of Zurich, Switzerland, July 2, 2014.

TH. Dickhaus, Simultaneous Bayesian analysis of contingency tables in genetic association studies, International Workshop ``Advances in Optimization and Statistics'', May 15  16, 2014, Russian Academy of Sciences, Institute of Information Transmission Problems (Kharkevich Institute), Moscow, May 15, 2014.

P. Friz, Implied, local and stochastic volatility: Analytical results in Heston and Stein Stein, Thematic Cycle on Robust Management in Finance, Workshop ``Model Approximation and Numerical Methods'' (February  March 2014), Institut Louis Bachelier, Paris, France, March 14, 2014.

H.J. Mucha, Bottomup variable selection in cluster analysis using bootstrapping: A proposal, European Conference on Data Analysis (ECDA 2014), July 2  4, 2014, Jacobs University Bremen gGmbH, July 4, 2014.

H.J. Mucha, From univariate to multivariate clustering with application to PEDXRF data, Workshop ``Application of Portable Energydispersive Xray Fluorescence to the Analysis of Archaeological Ceramics'', June 20  21, 2014, Excellence Cluster TOPOI, Berlin, June 21, 2014.

H.J. Mucha, Variable selection in cluster analysis by means of bootstrapping, International Workshop ``Advances in Optimization and Statistics'', May 15  16, 2014, Russian Academy of Sciences, Institute of Information Transmission Problems (Kharkevich Institute), Moscow, May 16, 2014.

J. Polzehl, Estimation of sparse precision matrices, MMSWorkshop ``large p small n'', WIASBerlin, April 15, 2014.

J. Polzehl, Quantification of noise in MR experiments, Statistical Challenges in Neuroscience, September 3  5, 2014, University of Warwick, Centre for Research in Statistical Methodology, UK, September 4, 2014.

J. Polzehl, Statistical problems in diffusion weighted MR, CoSy Seminar, University of Uppsala, Department of Mathematics, Sweden, November 11, 2014.

V. Spokoiny, Bernsteinvon Mises theorem for a quasiposterior, International Symposium on Financial Engineering and Risk Management 2014, June 27  28, 2014, Central University of Finance and Economics, Beijing, China, June 27, 2014.

V. Spokoiny, Bootstrap confidence sets under a model misspecification, Workshop ``New Horizons in Statistical Decision Theory'', September 7  13, 2014, Mathematisches Forschungsinstitut Oberwolfach, September 8, 2014.

V. Spokoiny, Construction of the sharp confidence bands using multiplier bootstrap, Advances in Stochastic Analysis, September 3  5, 2014, National Research University  Higher School of Economics, Laboratory of Stochastic Analysis and its Applications, Moscow, Russian Federation, September 5, 2014.

V. Spokoiny, Modern Parametric Statistics, February 10  25, 2014, Independent University of Moscow, Russian Federation.

K. Tabelow, S. Mohammadi, N. Weiskopf, J. Polzehl, Adaptive noise reduction in multishell dMRI data with SPM by POAS4SPM, 20th Annual Meeting of the Organization for Human Brain Mapping, Hamburg, June 8  12, 2014.

K. Tabelow, H.U. Voss, J. Polzehl, Local estimation of noise standard deviation in MRI images using propagation separation, 20th Annual Meeting of the Organization for Human Brain Mapping, Hamburg, June 8  12, 2014.

K. Tabelow, H.U. Voss, J. Polzehl, Local estimation of the noise level in MRI images using structural adaptation, 5th UltraHighfield MRI Scientific Symposium, Max Delbrück Center, Berlin, June 20, 2014.

K. Tabelow, Highresolution diffusion MRI by msPOAS, Statistical Challenges in Neuroscience, September 3  5, 2014, University of Warwick, Centre for Research in Statistical Methodology, UK, September 4, 2014.

S. Becker, The propagationseparation approach: Choice of the adaptation bandwidth and its (in)dependence of the unknown parameter, PreMoLab Workshop on: Advances in predictive modeling and optimization, May 16  17, 2013, WIASBerlin, May 17, 2013.

M. Zhilova, Resampling based uniform confidence bands in quantile regression, Tagung des SFB 649 "Ökonomisches Risiko" in Motzen, June 20  22, 2013, HumboldtUniversität zu Berlin, School of Business and Economics, June 22, 2013.

M. Zhilova, Uniform confidence bands in local estimation, 29th European Meeting of Statisticians (EMS), July 20  25, 2013, Eötvös Loránd University, Budapest, Hungary, July 20, 2013.

M. Zhilova, Uniform confidence bands in local estimation, Premoday February 2013, Moscow Institute of Physics and Technology State University, Russian Federation, February 22, 2013.

M. Zhilova, Uniform confidence bands in local parametric estimation, HaindorfSeminar 2013, February 7  9, 2013, HumboldtUniversität zu Berlin, Czech Republic, February 9, 2013.

R. Hildebrand, Convex projective programming, International Conference on Continuous Optimization (ICCOPT2013), July 29  August 1, 2013, Universidade Nove de Lisboa, Portugal, July 31, 2013.

R. Hildebrand, Scaling relationships between the copositive cone and the cone of sum of squares polynomials, International Conference on Continuous Optimization (ICCOPT2013), July 29  August 1, 2013, Universidade Nove de Lisboa, Portugal, July 30, 2013.

K. Tabelow, S. Becker, S. Mohammadi, N. Weiskopf, J. Polzehl, Multishell positionorientation adaptive smoothing (msPOAS), 19th Annual Meeting of the Organization for Human Brain Mapping, Seattle, USA, June 16  20, 2013.

M. Welvaert, K. Tabelow, R. Seurinck, Y. Rosseel, Defining ROIs based on localizer studies: More specific localization using adaptive smoothing, 19th Annual Meeting of the Organization for Human Brain Mapping, Seattle, USA, June 16  20, 2013.

TH. Dickhaus, Bootstrap methods for simultaneous statistical inference in dynamic factor models, PrincetonHumboldt Conference 2013, November 1  2, 2013, Princeton University, Department of Operations Research and Financial Engineering, USA, November 1, 2013.

H.J. Mucha, Assessment of stability in partitional clustering using resampling techniques, Third Bilateral GermanPolish Symposium on Data Analysis and Its Applications (GPSDAA 2013), September 26  28, 2013, Technische Universität Dresden, September 27, 2013.

H.J. Mucha, Resampling techniques in cluster analysis: Is subsampling better than bootstrapping?, European Conference on Data Analysis, July 10  12, 2013, University of Luxembourg, Luxembourg, July 11, 2013.

H.J. Mucha, Selection of variables in cluster analysis using resampling techniques: A proposal, Conference of the International Federation of Classification Societies (IFCS), July 14  17, 2013, Tilburg University, Netherlands, July 16, 2013.

H.J. Mucha, Variables election in cluster analysis using resampling techniques, British Classification Society Meeting 2013, November 8  9, 2013, University College London, Department of Statistical Science, UK, November 8, 2013.

V. Spokoiny, Identification and critical dimension in semiparametric estimation, Workshop ``Mathematical Statistics of Partially Identified Objects'', April 21  27, 2013, Mathematisches Forschungsinstitut Oberwolfach, April 21, 2013.

V. Spokoiny, Robust clustering using adaptive weights, International Workshop on Statistical Learning, June 26  28, 2013, Moscow Institute of Physics and Technology, PreMoLab, Russian Federation, June 27, 2013.

S. Becker, Diffusion weighted imaging: Modeling and analysis beyond the diffusion tensor, Methodological Workshop: Structural Brain Connectivity: Diffusion ImagingState of the Art and Beyond, October 30  November 2, 2012, HumboldtUniversität zu Berlin, November 2, 2012.

S. Becker, Image processing via orientation scores, Workshop ``Computational Inverse Problems'', October 23  26, 2012, Mathematisches Forschungsinstitut Oberwolfach, October 25, 2012.

S. Becker, Revisiting: Propagationseparation approach for local likelihood estimation, PreMoLab: MoscowBerlinStochastic and Predictive Modeling, May 29  June 1, 2012, Russian Academy of Sciences, Institute for Information Transmission Problems (Kharkevich Institute), Moscow, May 31, 2012.

M. Zhilova, Uniform properties of local maximum likelihood estimate, 55th Scientific conference of Moscow Institute of Physics and Technology, November 19  24, 2012, Moscow Institute of Physics and Technology, Dolgoprudny, Russian Federation, November 24, 2012.

M. Becker, CPdetection, outlook, ECONS Projekttreffen, PIK Potsdam, December 12, 2012.

M. Becker, Change point detection via universal compressors, PreMoLab: MoscowBerlin Stochastic and Predictive Modeling, May 30  June 2, 2012, Russian Academy of Sciences, Institute for Information Transmission Problems (Kharkevich Institute), Moscow, May 31, 2012.

M. Becker, Changepoint detection via universal compressors, ECONS  Projekttreffen, WIAS  Berlin, June 7, 2012.

K. Tabelow, Positionorientation adaptive smoothing (POAS) diffusion weighted imaging data, Workshop on Neurogeometry, November 15  17, 2012, Masaryk University, Department of Mathematics and Statistics, Brno, Czech Republic, November 16, 2012.

J. Dolata, H.J. Mucha, H.G. Bartel, Eine neue Referenzgruppe für die römische Ziegelproduktion der I. Thrakerkohorte am Mittelrhein, Archäometrie und Denkmalpflege 2012, Tübingen, March 27  31, 2012.

C. MoralesMerino, H.J. Mucha, H.G. Bartel, Multivariate statistical analysis of clay and ceramicdata for provenance of Bronze age pottery from Troia, Archäometrie und Denkmalpflege 2012, Tübingen, March 27  31, 2012.

M. Zhilova, Quasi MLE in Instrumental variable model, PreMoDay II, March 16, 2012, Russian Academy of Sciences, Institute for Information Transmission Problems (Kharkevich Institute), Russian Federation, March 16, 2012.

H.J. Mucha, Clustering and visualisation based on dual scaling, 4th JapaneseGerman Symposium on Classification (JGSC2012), March 8  10, 2012, Doshisha University, Kyoto, Japan, March 9, 2012.

H.J. Mucha, Dual scaling classification and its application to archaeometry, 36th Annual Conference of the German Classification Society, August 1  3, 2012, Universität Hildesheim, Information Systems and Machine Learning Lab (ISMLL), August 2, 2012.

H.J. Mucha, Measures of incomparability and of inequality and their application, 10th Workshop on Partial Order, Theory and Application: MultiIndicator Systems and Modelling in Partial Order, September 27  28, 2012, Hochschule für Technik und Wirtschaft Berlin, September 28, 2012.

H.J. Mucha, Soft bootstrapping and its comparison with other resampling methods, 36th Annual Conference of the German Classification Society, August 1  3, 2012, Universität Hildesheim, Information Systems and Machine Learning Lab (ISMLL), August 2, 2012.

H.J. Mucha, Visualisation and cluster analysis, Herbsttagung der AG Datenanalyse und numerische Klassifikation, October 5  6, 2012, Landschaftsverband Rheinland, LVRAmt für Bodendenkmalpflege im Rheinland, Bonn, October 6, 2012.

J. Polzehl, Adaptive methods for noise reduction in diffusion weighted MR, BRIC Seminar Series, University of North Carolina, School of Medicine, Chapel Hill, NC, USA, July 10, 2012.

V. Spokoiny, Basics of modern parametric statistics, February 13  28, 2012, Independent University of Moscow, Center for Continuous Mathematical Education, Russian Federation.

V. Spokoiny, Bernsteinvon Mises theorem for quasi posteriors, International Workshop on Recent Advances in Time Series Analysis (RATS 2012), June 8  12, 2012, University of Cyprus, Department of Mathematics and Statistics, Protaras, June 9, 2012.

V. Spokoiny, Bernsteinvon Mises theorem for quasi posteriors, Workshop ``Frontiers in Nonparametric Statistics'', March 11  17, 2012, Mathematisches Forschungsinstitut Oberwolfach, March 12, 2012.

V. Spokoiny, Bernsteinvon Mises theorem for quasi posteriors, Workshop II on Financial Time Series Analysis: Highdimensionality, Nonstationarity and the Financial Crisis, June 19  22, 2012, National University of Singapore, Institute for Mathematical Sciences, June 21, 2012.

V. Spokoiny, Bernsteinvon Mises theorem for quasi posteriors, Workshop on Recent Developments in Statistical Multiscale Methods, July 16  18, 2012, GeorgAugustUniversität Göttingen, Institut für Mathematische Stochastik, July 17, 2012.

V. Spokoiny, Bernsteinvon Mises theorem for quasi posteriors, PreMoLab Seminar, Russian Academy of Sciences, Institute for Information Transmission Problems (Kharkevich Institute), Moscow, March 15, 2012.

V. Spokoiny, Parametric estimation: Modern view, PreMoDay I, February 24, 2012, Russian Academy of Sciences, Institute for Information Transmission Problems (Kharkevich Institute), Moscow, February 24, 2012.

V. Spokoiny, Some methods of modern statistics, Information Technology and Systems 2, August 19  25, 2012, Russian Academy of Sciences, Institute for Information Transmission Problems (Kharkevich Institute), Petrozavodsk, August 20, 2012.

M. Becker, Change point detection via universal compressors, ECONS (Evolving Complex Networks) Spring School and Workshop, March 28  31, 2011, Potsdam Institut für Klimafolgenforschung, Wandlitz, March 29, 2011.

S. Becker, Structural adaptive smoothing in R3xS2 with application to diffusion weighted imaging, Statistics Mathematics and Applications (StatMathAppli 2011), August 29  September 2, 2011, Fréjus, France, September 1, 2011.

S. Becker, Structural adaptive smoothing in R3xS2 with application to diffusion weighted imaging, Chemnitz Symposium on Inverse Problems 2011, September 22  23, 2011, Technische Universität, Fakultät für Mathematik, September 22, 2011.

E. Diederichs, Modellselektion durch Semidefinite Relaxation, Jahrestagung der Deutschen MathematikerVereinigung (DMV) 2011, September 19  22, 2011, Universität zu Köln, Mathematisches Institut, September 21, 2011.

E. Diederichs, Recent trends in large scale optimization, ECONS (Evolving Complex Networks) Spring School and Workshop, March 28  31, 2011, Potsdam Institut für Klimafolgenforschung, Wandlitz, March 31, 2011.

V. Panov, NonGaussian component analysis, ECONS Spring School and Workshop, March 28  31, 2011, Potsdam Institut für Klimaforschung, March 31, 2011.

S. Becker, Modelfree structural adaptive smoothing of diffusion weighted images, Workshop on Statistics and Neuroimaging 2011, November 23  25, 2011, WIAS, November 24, 2011.

M. Becker, Current state of research on SILT and CPdetection, ECONSprojektmeeting, November 21, 2011, LeibnizInstitut für Gewässerökologie und Binnenfischerei (IGB), November 21, 2011.

E. Diederichs, Sparse nonGaussian component analysis, Structural Methods of Data Analysis and Optimization, December 15  17, 2011, Russian Academy of Sciences, Institute for Information Transmission Problems (Kharkevich Institute), Moscow, December 16, 2011.

H.J. Mucha, Geochemical and statistical investigation of clay deposits in the troad and its implication for provenance of bronze age fine pottery from Troia, Joint Conference of the German Classification Society (GfKl), German Association for Pattern Recognition (DAGM), Symposium of the International Federation of Classification Societies (IFCS), August 30  September 2, 2011, Goethe Universität Frankfurt am Main, September 1, 2011.

H.J. Mucha, Method selection in cluster analysis followed by builtin validation, Second Bilateral GermanPolish Symposium on Data Analysis and Its Applications (GPSDAA 2011), April 14  16, 2011, Cracow University of Economics, Poland, April 16, 2011.

J. Polzehl, Statistical issues in modeling diffusion weighted magnetic resonance data, 3rd International Conference on Statistics and Probability 2011 (IMSChina), July 8  11, 2011, Institute of Mathematical Statistics, Xian, China, July 10, 2011.

V. Spokoiny, Alternating and semiparametric efficiency, Conference in Honour of Joel Horowitz, June 23  24, 2011, Centre for Microdata Methods and Practice, London, UK, June 24, 2011.

V. Spokoiny, Alternating and semiparametric efficiency, École Nationale de la Statistique et de l'Analyse de l'Information (ENSAI), Rennes, France, September 16, 2011.

V. Spokoiny, Methods of dimension reduction, ECONS (Evolving Complex Networks) Spring School and Workshop, March 28  30, 2011, Potsdam Institut für Klimafolgenforschung, March 29, 2011.

V. Spokoiny, Modern parametric theory, September 13  16, 2011, École Nationale de la Statistique et de l'Analyse de l'Information (ENSAI), Rennes, France.

V. Spokoiny, NonGaussian component analysis using semidefinite relaxation, MASCOT NUM (Méthodes d'Analyse Stochastique pour les COdes et Traitements NUMériques) 2011 Workshop, March 23  25, 2011, Villard de Lans, France, March 25, 2011.

V. Spokoiny, Parametric inference. Revisited, 5èmes Journées Statistiques du Sud, June 14  16, 2011, Université Nice Sophia Antipolis, Faculté des Sciences, June 15, 2011.

V. Spokoiny, Structure adaptive estimation by alternating, Workshop ``Stuctural Inference Day'', April 18, 2011, Universität Hamburg, April 18, 2011.

H.J. Mucha, Classification based on dual scaling, Herbsttagung der AG DANK der GfKl, November 11  12, 2011, Heinrich Heine Universität Düsseldorf, November 11, 2011.

H.J. Mucha, Classification of roman tiles with stamp Pardalis, Joint Conference of the German Classification Society (GfKl), German Association for Pattern Recognition (DAGM), Symposium of the International Federation of Classification Societies (IFCS), August 30  September 2, 2011, Goethe Universität Frankfurt am Main, September 1, 2011.

V. Spokoiny, Alternating and semiparametric efficiency, Workshop ``Very High Dimensional Semiparametric Models'', October 3  7, 2011, Mathematisches Forschungsinstitut Oberwolfach, October 4, 2011.

V. Spokoiny, Semiparametric estimation alternating and efficiency, École Nationale de la Statistique et de l'Administration Économique (ENSAE), Paris, France, December 12, 2011.

A. Beinrucker, An application of Gaussian processes to Bayesian statistics, Seminar ``Gaussian Processe'', Universität Potsdam, Fachbereich Mathematik, November 25, 2010.

A. Beinrucker, Classification using FISTA, MASH Meeting, November 8  9, 2010, EUProjekt MASH, Prag, Czech Republic, November 9, 2010.

A. Beinrucker, MASH  Massive sets of heuristics, Machine Learning Summer School, Canberra, Australia, September 27  October 6, 2010.

N. Krämer, Conjugate gradient regularization  A statistical framework for partial least squares regression, 4th Workshop on Partial Least Squares and Related Methods for Cuttingedge Research in Experimental Sciences, May 10  11, 2010, École Supérieure d'Électricité, Department of Information Systems and Decision Sciences, GifsurYvette, France, May 11, 2010.

N. Krämer, The degrees of freedom of partial least squares regression, DAGStat 2010, March 23  26, 2010, Technische Universität Dortmund, Fakultät Statistik, March 23, 2010.

V. Panov, NonGaussian component classification with applications to American Option Pricing, Haindorf Seminar 2010 (Klausurtagung des SFB 649), February 11  14, 2010, HumboldtUniversität zu Berlin, Wirtschaftswissenschaftliche Fakultät, Hejnice, Czech Republic, February 13, 2010.

K. Tabelow, J.D. Clayden, P. Lafaye DE Micheaux, J. Polzehl, V.J. Schmid, B. Whitcher, Image analysis and statistical inference in NeuroImaging with R., Human Brain Mapping 2010, Barcelona, Spain, June 6  10, 2010.

K. Tabelow, J. Polzehl, S. Mohammadi, M. Deppe, Impact of smoothing on the interpretation of FA maps, Human Brain Mapping 2010, Barcelona, Spain, June 6  10, 2010.

G. Blanchard, N. Krämer, Kernel partial least squares is universally consistent, AI & Statistics 2010, Sardinia, Italy, May 13  15, 2010.

G. Blanchard, N. Krämer, Optimal rates for conjugate gradient regularization, AI & Statistics 2010, Sardinia, Italy, May 13  15, 2010.

G. Blanchard, On optimal rates for kernel conjugate gradient regularization under random design and noise, Mini Special Semester on Inverse Problems, Part II: ``Inverse Problems in Data Driven Modelling'', July 19  22, 2010, Johann Radon Institute for Computational and Applied Mathematics (RICAM), Linz, Austria, July 20, 2010.

G. Blanchard, Optimal rate for kernel conjugate gradient regularization, Workshop "Modern Nonparametric Statistics: Going Beyond Asymptotic Minimax", March 29  April 2, 2010, Mathematisches Forschungsinstitut Oberwolfach, March 30, 2010.

N. Krämer, Modeling of highdimensional data with partial least squares, University of Leeds, Faculty of Medicine and Health, Leeds Institute of Genetics, Health and Therapeutics, UK, September 2, 2010.

N. Krämer, Regularized estimation of largescale gene associaton networks, 55. GMDSJahrestagung: Effiziente und wirtschaftliche Gesundheitsversorgung von heute und morgen  nur mit Medizinischer Dokumentation, Medizinischer Informatik, Medizinischer Biometrie und Epidemiologie, September 6  7, 2010, Hochschule Mannheim, Institut für Medizinische Informatik, September 6, 2010.

N. Krämer, The consistency of partial least spuares and conjugate gradient regularization, 73rd Annual Meeting of the Institute of Mathematical Statistics, August 9  13, 2010, Institute of Mathematical Statistics (IMS), Sweden, August 13, 2010.

P. Mathé, The Lepskii balancing principle for conjugate gradient regularization, Workshop on Inverse Problems: Developments in Theory and Applications (IPTA 2010), February 8  12, 2010, Stefan Banach International Mathematical Center, Warsaw, Poland, February 9, 2010.

H.J. Mucha, An accelerated Kmeans algorithm based on adaptive distances, 34th Annual Conference of the German Classification Society (GfKl) and International Symposium on the Data Analysis Interface, July 21  23, 2010, Karlsruhe Institute of Technology (KIT), July 22, 2010.

H.J. Mucha, Automatische Klassifikation (Clusteranalyse), IBM SPSS Kundentag 2010, June 8, 2010, SPSS GmbH Software, München, June 8, 2010.

H.J. Mucha, Pairwise data clustering accompanied by validation and visualisation, 3rd GermanJapanese Workshop ``Advances in Data Analysis and Related New Techniques and Applications'', July 20  21, 2010, Karlsruhe Institute of Technology (KIT), July 20, 2010.

J. Polzehl, K. Tabelow, Image and signal processing in the biomedical sciences: Diffusionweighted imaging modeling and beyond, 1st Annual Scientific Symposium ``Ultrahigh Field Magnetic Resonance'', Max Delbrück Center, Berlin, April 16, 2010.

J. Polzehl, Statistical problems in functional and diffusion weighted magnetic resonance, Uppsala University, Dept. of Mathematics, Graduate School in Mathematics and Computing, Sweden, May 27, 2010.

V. Spokoiny, Local parametric estimation, October 18  22, 2010, École Nationale de la Statistique et de l'Analyse de l'Information (ENSAI), Rennes, France.

V. Spokoiny, Semidefinite nonGaussian component analysis, Bivariate Penalty Choice in Model Selection, Deutsches Diabetes Zentrum Düsseldorf, June 17, 2010.

V. Spokoiny, Semidefinite optimization in instrumental variable estimation with shape constraints, Conference on ``Shape Restrictions in Non and Semiparametric Estimation of Econometric Models'', November 1  7, 2010, Northwestern University, Department of Statisics, Chicago, USA, November 6, 2010.

V. Spokoiny, Sparse nonGaussian component analysis, Journées STAR 2010, October 21  22, 2010, Université Rennes II, France, October 22, 2010.

V. Spokoiny, Sparse nonGaussian component analysis, School of Mathematics Colloquium, November 19, 2010, University of Edinburgh, School of Mathematics, UK, November 19, 2010.

V. Spokoiny, Sparse nonGaussian component analysis, University of Cambridge, Faculty of Mathematics, Statistical Laboratory, UK, May 7, 2010.

V. Spokoiny, Statistical estimation: Modern view, Rencontres de Statistiques Mathématiques 10, December 13  17, 2010, Centre International de Rencontres Mathématiques (CIRM), Luminy, France, December 16, 2010.

V. Spokoiny, Statistical estimation: Modern view, Rencontres de Statistiques Mathématiques 10, December 13  17, 2010, Centre International de Rencontres Mathématiques (CIRM), Luminy, France, December 15, 2010.

V. Panov, Pricing Bermudan options via dimension reduction, Klausurtagung des SFB 649, June 4  6, 2009, HumboldtUniversität zu Berlin, Motzen, June 5, 2009.

V. Panov, Solving classification task using detection of nonGaussian components, Haindorf Seminar 2009, February 12  15, 2009, HumboldtUniversität zu Berlin, CASE  Center for Applied Statistics and Economics, Hejnice, Czech Republic, February 14, 2009.

D. Peschka, Dewetting of thin liquid films on viscoelastic substrates, European Coating Symposium, September 7  9, 2009, Karlsruhe, September 7, 2009.

D. Belomestny, Estimation of the jump activity of a Lévy process from low frequency data, Haindorf Seminar 2009, February 12  15, 2009, HumboldtUniversität zu Berlin, CASE  Center for Applied Statistics and Economics, Hejnice, Czech Republic, February 12, 2009.

D. Belomestny, Regression methods for stochastic control problems and their convergence analysis, Fourth General Conference on Advanced Mathematical Methods in Finance, May 4  10, 2009, University of Oslo, Norway, May 9, 2009.

D. Belomestny, Spectral estimation of the fractional order of a Lévy process, Workshop ``Statistical Inference for Lévy Processes with Applications to Finance'', July 15  17, 2009, EURANDOM, Eindhoven, Netherlands, July 16, 2009.

G. Blanchard, Convergence du gradient conjugué fonctionnel pour l'apprentissage statistique, École Normale Supérieure, Paris, France, March 16, 2009.

G. Blanchard, Nonasymptotic adaptive control of the familywise error rate in multiple testing, Journées Statistiques du Sud 2009, June 17  19, 2009, Université de Provence AixMarseille I, Porquerolles, France, June 17, 2009.

H.J. Mucha, Clustering accompanied by visualization and validation, 11th Conference of the International Federation of Classification Societies (IFCS 2009), March 13  18, 2009, Technische Universität Dresden, March 13, 2009.

H.J. Mucha, Explorative statistische Methoden in der Archäometrie: von Daten zu Distanzen und weiter zu Gruppierungen, CurtEngelhornZentrum Archäometrie Mannheim, July 9, 2009.

H.J. Mucha, Geochemical and statistical investigation of Roman stamped tiles of the Legio XXI Rapax, 11th Conference of the International Federation of Classification Societies (IFCS 2009), March 13  18, 2009, Technische Universität Dresden, March 14, 2009.

H.J. Mucha, Remarks on Validation of Hierarchical Clustering, First Bilateral GermanPolish Symposium on Data Analysis and its Applications, October 8  10, 2009, RWTH Aachen, October 10, 2009.

H.J. Mucha, Zur Klassifikation römischer Ziegel von Fundorten im südlichen Obergermanien, Tagung ``Archäometrie und Denkmalpflege 2009'', March 25  28, 2009, Pinakothek der Moderne, München, March 26, 2009.

N. Serdyukova, Local parametric estimation under noise misspecification in regression problem, Workshop on structure adapting methods, November 6  8, 2009, WIAS, November 7, 2009.

V. Spokoiny, Adaptive local parametric estimation, Université Joseph Fourier Grenoble I, Équipe de Statistique et Modélisation Stochastique, Laboratoire Jean Kuntzmann, France, February 26, 2009.

V. Spokoiny, Adaptive local parametric methods in imaging, Technische Universität Kaiserslautern, Fachbereich Mathematik, January 23, 2009.

V. Spokoiny, Cross concentration in semiparametric estimation, Rencontres de Statistiques Mathématiques 9, December 14  18, 2009, International Center for Mathematical Meetings (CIRM), Marseille, France, December 15, 2009.

V. Spokoiny, Modern nonparametric statistics (block lecture), October 2  13, 2009, École Nationale de la Statistique et de l'Analyse de l'Information (ENSAI), Rennes, France.

V. Spokoiny, Modern nonparametric statistics (block lecture), October 18  29, 2009, Yale University, New Haven, USA.

V. Spokoiny, Modern nonparametric statistics (block lecture), January 13  16, 2009, École Nationale de la Statistique et de l'Analyse de l'Information (ENSAI), Rennes, France.

V. Spokoiny, NonGaussian component analysis, Workshop ``Semiparametric and Nonparametric Methods in Econometrics'', April 4  11, 2009, Banff International Research Station for Mathematical Innovation and Discovery, Canada, April 6, 2009.

V. Spokoiny, Parameter tuning in statistical inverse problem, European Meeting of Statisticians (EMS2009), July 20  22, 2009, Université Paul Sabatier, Toulouse, France, July 21, 2009.

V. Spokoiny, Saddle point model selection, 2nd HumboldtPrinceton Conference: Perceiving and Measuring Financial Risk: Credit, Energy and Illiquidity, October 30  31, 2009, Princeton University, USA, October 30, 2009.

V. Spokoiny, Saddle point model selection, Université Toulouse 1 Capitole, Toulouse School of Economics, France, November 24, 2009.

V. Spokoiny, Saddle point model selection, Workshop on structure adapting methods, November 6  8, 2009, WIAS, November 7, 2009.

V. Spokoiny, Saddle point model selection for inverse problems, Workshop ``Challenges in Statistical Theory: Complex Data Structures and Algorithmic Optimization'', August 23  29, 2009, Mathematisches Forschungsinstitut Oberwolfach, August 26, 2009.

V. Spokoiny, Sparse nonGaussian component analysis, Workshop ``Sparse Recovery Problems in High Dimensions: Statistical Inference and Learning Theory'', March 15  21, 2009, Mathematisches Forschungsinstitut Oberwolfach, March 16, 2009.

M. Elagin, Comparison of some methods for the adaptive volatility estimation, Haindorf Seminar 2008 (Klausurtagung des SFB 649), February 7  10, 2008, Hejnice, Czech Republic, February 9, 2008.

K. Tabelow, Structure adaptive smoothing medical images, 22. Treffpunkt Medizintechnik: Fortschritte in der medizinischen Bildgebung, Charité, Campus Virchow Klinikum Berlin, May 22, 2008.

K. Tabelow, Strukturadaptive Bild und Signalverarbeitung, Workshop of scshape Matheon with Siemens AG (Health Care Sector) in cooperation with Center of Knowledge Interchange (CKI) of Technische Universität (TU) Berlin and Siemens AG, TU Berlin, July 8, 2008.

E. Giacomini, Dynamic semiparametric factor modus in risk preferences estimation, 2008 North American Winter Meetings of the Econometric Society, January 4  6, 2008, North American Econometric Society, New Orleans, USA, January 6, 2008.

H.J. Mucha, Clustering, multivariate visualization and validation, 30th Fall Meeting of the German Classification Society, Berlin, Gesellschaft für Klassifikation e.V., Berlin, November 14, 2008.

H.J. Mucha, Mapping find spots of Roman military brickstamps in Mogontiacum (Mainz) and archeometrical analysis, Annual Meeting of Gesellschaft für Klassifikation, July 16  18, 2008, HelmutSchmidtUniversität, Hamburg, August 17, 2008.

J. Polzehl, New developments in structural adaptive smoothing: Images, fMRI and DWI, University of Tromsoe, Norway, May 27, 2008.

J. Polzehl, Smoothing fMRI and DWI data using the propagationseparation approach, University of Utah, Computing and Scientific Imaging Institute, Salt Lake City, USA, September 11, 2008.

J. Polzehl, Structural adaptive smoothing using the propagationseparation approach, University of Chicago, Department of Statistics, USA, September 3, 2008.

A. Rohde, Spatially adaptive comparison of multivariate samples via randomization, Université Paris VI ``Pierre et Marie Curie'', France, March 10, 2008.

A. Rohde, Spatially adaptive comparison of multivariate samples via randomization, Universität Bern, Switzerland, February 20, 2008.

N. Serdyukova, Dependence on the dimension for complexity of approximation of random fields, 4. Doktorandenkonferenz Stochastik, September 4  6, 2008, Technische Universität Berlin, September 6, 2008.

V. Spokoiny, Foundations and applications of modern nonparametric statistics (block lecture), May 15  16, 2008, Graduate School ``Mathematics and Practice'', Technische Universität Kaiserslautern/Fraunhofer Institut für Techno und Wirtschaftsmathematik.

V. Spokoiny, Local parametric modeling of nonstationary time series, International Workshop on Recent Advances in Time Series Analysis, June 8  11, 2008, University of Cyprus, Department of Mathematics and Statistics, Protaras, June 11, 2008.

K. Tabelow, J. Polzehl, H.U. Voss, Increasing SNR in high resolution fMRI by spatially adaptive smoothing, Human Brain Mapping Conference 2007, Chicago, USA, June 10  14, 2007.

K. Tabelow, A3: Image and signal processing in medicine and biosciences, ADay des sc Matheon, KonradZuseZentrum für Informationstechnik Berlin (ZIB), December 5, 2007.

K. Tabelow, Improving data quality in fMRI and DTI by structural adaptive smoothing, Cornell University, Weill Medical College, New York, USA, June 18, 2007.

K. Tabelow, Structural adaptive signal detection in fMRI and structure enhancement in DTI, International Workshop on Image Analysis in the Life Sciences, Theory and Applications, February 28  March 2, 2007, Johannes Kepler Universität Linz, Austria, March 2, 2007.

K. Tabelow, Structural adaptive smoothing in medical imagine, WIAS Day, WIAS Berlin, February 22, 2007.

K. Tabelow, Structural adaptive smoothing in medical imaging, Seminar ``Visualisierung und Datenanalyse'', KonradZuseZentrum für Informationstechnik Berlin (ZIB), January 30, 2007.

J. Polzehl, Propagationseparation procedures for image processing, International Workshop on Image Analysis in the Life Sciences, Theory and Applications, February 28  March 2, 2007, Johannes Kepler Universität Linz, Austria, March 2, 2007.

J. Polzehl, Structural adaptive smoothing in imaging problems, Spring Seminar Series, University of Minnesota, School of Statistics, College of Liberal Arts, USA, May 24, 2007.

J. Polzehl, Structural adaptive smoothing procedures by propagationseparation methods, Final meeting of the DFG Priority Program 1114, November 7  9, 2007, Freiburg, November 7, 2007.

K. Tabelow, J. Polzehl, H.U. Voss, V. Spokoiny, Analyzing fMRI experiments with structural adaptive smoothing methods, Human Brain Mapping Conference, Florence, Italy, June 12  15, 2006.

K. Tabelow, Analyzing fMRI experiments with structural adaptive smoothing methods, BCCN PhD Symposium 2006, June 7  8, 2006, Bernstein Center for Computational Neuroscience Berlin, Bad Liebenwalde, June 8, 2006.

K. Tabelow, Image and signal processing in medicine and biosciences, Evaluation Colloquium of the DFG Research Center sc Matheon, Berlin, January 24  25, 2006.

J. Polzehl, Structural adaptive smoothing by propagationseparation, 69th Annual Meeting of the IMS and 5th International Symposium on Probability and its Applications, July 30  August 4, 2006, Rio de Janeiro, Brazil, July 30, 2006.

K. Tabelow, Adaptive weights smoothing in the analysis of fMRI data, LudwigMaximiliansUniversität München, SFB 386, December 8, 2005.

K. Tabelow, Detecting shape and borders of activation areas infMRI data, Forschungsseminar ''Mathematische Statistik'', WIAS, Berlin, November 23, 2005.

K. Tabelow, Spatially adaptive smoothing infMRI analysis, Neuroimaging Center, Cahrité, Berlin, November 10, 2005.

J. Polzehl, Adaptive smoothing by propagationseparation, Australian National University, Center of Mathematics and its Applications, Canberra, March 31, 2005.

J. Polzehl, Image reconstruction and edge enhancement by structural adaptive smoothing, 55th Session of the International Statistical Institute (ISI), April 5  12, 2005, Sydney, Australia, April 8, 2005.

J. Polzehl, Propagationseparation at work: Main ideas and applications, National University of Singapore, Department of Probability Theory and Statistics, March 24, 2005.

J. Polzehl, Spatially adaptive smoothing: A propagationseparation approach for imaging problems, Joint Statistical Meetings, August 7  11, 2005, Minneapolis, USA, August 11, 2005.

J. Polzehl, Structural adaptive smoothing by propagationseparation methods, LudwigMaximiliansUniversität München, SFB 386, December 7, 2005.

J.G.M. Schoenmakers, Robust Libor modelling and pricing of derivative products, Delft University of Technology, Netherlands, June 9, 2005.

J. Polzehl, Local likelihood modeling by structural adaptive smoothing, University of Minnesota, School of Statistics, Minneapolis, USA, September 9, 2004.

J. Polzehl, On a nonstationary structural adaptive approach to volatility estimation, University of Gothenburg, Centre for Finance, Sweden, May 5, 2004.

J. Polzehl, Smoothing by adaptive weights: An overview, Chalmers University of Technology, Department of Mathematical Statistics, Gothenburg, Sweden, May 11, 2004.

J. Polzehl, Structural adaptive smoothing methods, GeorgAugustUniversität Göttingen, Institut für Mathematische Stochastik, January 14, 2004.

J. Polzehl, Structural adaptive smoothing methods, TandemWorkshop on Nonlinear Optimization at the Crossover of Discrete Geometry and Numerical Analysis, July 15  16, 2004, Technische Universität Berlin, Institut für Mathematik, July 15, 2004.

J. Polzehl, Structural adaptive smoothing methods and possible applications in imaging, Charité Berlin, NeuroImaging Center, Berlin, July 1, 2004.

J. Polzehl, Structural adaptive smoothing methods for imaging problems, Annual Conference of Deutsche MathematikerVereinigung (DMV), September 13  17, 2004, Heidelberg, September 14, 2004.

J. Polzehl, Structural adaptive smoothing methods for imaging problems, GermanIsraeli Binational Workshop, October 20  22, 2004, Ollendorff Minerva Center for Vision and Image Sciences, Technion, Haifa, Israel, October 21, 2004.

J. Polzehl, Adaptive smoothing procedures for image processing, Workshop on Nonlinear Analysis of Multidimensional Signals, February 25  28, 2003, Teistungenburg, February 25, 2003.

J. Polzehl, Image processing using Adaptive Weights Smoothing, Uppsala University, Department of Mathematics, Sweden, May 7, 2003.

J. Polzehl, Local likelihood modeling by Adaptive Weights Smoothing, Joint Statistical Meetings, August 3  7, 2003, San Francisco, USA, August 6, 2003.

J. Polzehl, Structural adaptation I: Pointwise adaptive smoothing and imaging, University of Tromso, Department of Mathematics, Norway, April 11, 2002.

J. Polzehl, Structural adaptation I: Varying coefficient regression modeling by adaptive weights smoothing, Workshop on Nonparametric Smoothing in Complex Statistical Models, April 27  May 4, 2002, Ascona, Switzerland, April 30, 2002.

J. Polzehl, Structural adaptation II: Time series and estimation of dimension reduction spaces, University of Tromso, Department of Mathematics, Norway, April 17, 2002.

J. Polzehl, Structural adaptation methods in imaging, Joint Statistical Meetings 2002, August 11  15, 2002, New York, USA, August 12, 2002.

J. Polzehl, Structural adaptive smoothing and its applications in imaging and time series, Uppsala University, Department of Mathematics, Sweden, May 2, 2002.

J. Polzehl, Structuraladaptive smoothing methods, FrenchGerman Seminar, Universität Potsdam, April 6, 2002.

J. Polzehl, Can structural assumptions be used to improve nonparametric estimates?, University of Minnesota, School of Statistics, Minneapolis, USA, May 31, 2001.

J. Polzehl, Structural adaptation in nonparametric regression, Workshop on HighDimensional Nonlinear Statistical Modelling, September 15  19, 2001, Wulkow, September 16, 2001.

J. Polzehl, Structural adaption  A method to estimate the effective dimension reduction space, Closed Meeting of Sfb 373, May 17  19, 2001, Wulkow, May 18, 2001.

J. Polzehl, Structural adaption in nonparametric smoothing, Departamento de Estadistica y Econometria, Universidad Carlos III de Madrid, Spain, March 2, 2001.

J. Polzehl, Structural adaptive estimation, Bayer AG, Leverkusen, November 29, 2001.

J. Polzehl, Adaptive weights smoothing with applications in imaging, Universität Essen, Fachbereich Mathematik, Sfb 475, November 6, 2000.

J. Polzehl, Adaptive weights smoothing with applications to image denoising and signal detection, Université Catholique de LouvainlaNeuve, Institut de Statistique, Belgium, September 29, 2000.

J. Polzehl, Functional and dynamic Magnet Resonance Imaging using adaptive weights smoothing, Workshop "`Mathematical Methods in Brain Mapping"', Université de Montréal, Centre de Recherches Mathématiques, Canada, December 11, 2000.

J. Polzehl, Spatially adaptive procedures for signal detection in fMRI, Tagung "`Controlling Complexity for Strong Stochastic Dependencies"', September 10  16, 2000, Mathematisches Forschungsinstitut Oberwolfach, September 11, 2000.

J. Polzehl, Spatially adaptive smoothing techniques for signal detection in functional and dynamic Magnet Resonance Imaging, Human Brain Mapping 2000, San Antonio, Texas, USA, June 12  16, 2000.
External Preprints

A. Daneshmand, G. Scutari, P. Dvurechensky, A. Gasnikov, Newton method over networks is fast up to the statistical precision, Preprint no. arXiv:2102.06780, Cornell University Library, arXiv.org, 2021.

A. Kroshnin, V. Spokoiny, A. Suvorikova, Multiplier bootstrap for BuresWasserstein barycenters, Preprint no. arXiv:2111.12612, Cornell University Library, arXiv.org, 2021.
Abstract
BuresWasserstein barycenter is a popular and promising tool in analysis of complex data like graphs, images etc. In many applications the input data are random with an unknown distribution, and uncertainty quantification becomes a crucial issue. This paper offers an approach based on multiplier bootstrap to quantify the error of approximating the true BuresWasserstein barycenter Q? by its empirical counterpart Qn. The main results state the bootstrap validity under general assumptions on the data generating distribution P and specifies the approximation rates for the case of subexponential P. The performance of the method is illustrated on synthetic data generated from the weighted stochastic block model. 
A. Rogozin, A. Beznosikov, D. Dvinskikh, D. Kovalev, P. Dvurechensky, A. Gasnikov, Decentralized distributed optimization for saddle point problems, Preprint no. arXiv:2102.07758, Cornell University Library, arXiv.org, 2021.

P. Dvurechensky, D. Kamzolov, A. Lukashevich, S. Lee, E. Ordentlich, C.A. Uribe, A. Gasnikov, Hyperfast secondorder local solvers for efficient statistically preconditioned distributed optimization, Preprint no. arXiv:2102.08246, Cornell University Library, arXiv.org, 2021.

Y.W. Sun, K. Papagiannouli, V. Spokoiny, Online graphbased changepoint detection for high dimensional data, Preprint no. arXiv:1906.03001, Cornell University Library, arXiv.org, 2019.
Abstract
Online changepoint detection (OCPD) is important for application in various areas such as finance, biology, and the Internet of Things (IoT). However, OCPD faces major challenges due to highdimensionality, and it is still rarely studied in literature. In this paper, we propose a novel, online, graphbased, changepoint detection algorithm to detect change of distribution in low to highdimensional data. We introduce a similarity measure, which is derived from the graphspanning ratio, to test statistically if a change occurs. Through numerical study using artificial online datasets, our datadriven approach demonstrates high detection power for highdimensional data, while the false alarm rate (type I error) is controlled at a nominal significant level. In particular, our graphspanning approach has desirable power with small and multiple scanning window, which allows timely detection of changepoint in the online setting. 
A. Kroshnin, V. Spokoiny, A. Suvorikova, Statistical inference for BuresWasserstein barycenters, Preprint no. arXiv:1901.00226, Cornell University Library, arXiv.org, 2019.

N. Puchkin, V. Spokoiny, Structureadaptive manifold estimation, Preprint no. arXiv:1906.05014, Cornell University Library, arXiv.org, 2019.
Abstract
We consider a problem of manifold estimation from noisy observations. Many manifold learning procedures locally approximate a manifold by a weighted average over a small neighborhood. However, in the presence of large noise, the assigned weights become so corrupted that the averaged estimate shows very poor performance. We suggest a novel computationally efficient structureadaptive procedure, which simultaneously reconstructs a smooth manifold and estimates projections of the point cloud onto this manifold. The proposed approach iteratively refines the weights on each step, using the structural information obtained at previous steps. After several iterations, we obtain nearly öracle" weights, so that the final estimates are nearly efficient even in the presence of relatively large noise. In our theoretical study we establish tight lower and upper bounds proving asymptotic optimality of the method for manifold estimation under the Hausdorff loss. Our finite sample study confirms a very reasonable performance of the procedure in comparison with the other methods of manifold estimation. 
N. Alia, V. John, S. Ollila, Revisiting the singlephase flow model for liquid steel ladle stirred by gas, Preprint no. arXiv.1811.11535, Cornell University Library, arXiv.org, 2018, DOI 10.1016/j.apm.2018.11.005 .
Abstract
Ladle stirring is an important step of the steelmaking process to homogenize the temperature and the chemical composition of the liquid steel and to remove inclusions before casting. Gas is injected from the bottom of the bath to induce a turbulent flow of the liquid steel. Multiphase modeling of ladle stirring can become computationally expensive, especially when used within optimal flow control problems. This paper focuses therefore on singlephase flow models. It aims at improving the existing models from the literature. Simulations in a 2d axialsymmetrical configuration, as well as, in a real 3d laboratoryscale ladle, are performed. The results obtained with the present model are in a relative good agreement with experimental data and suggest that it can be used as an efficient model in optimal flow control problems. 
A. Koziuk, V. Spokoiny, Instrumental variables regression, Preprint no. arXiv:1806.06111v1, Cornell University Library, arXiv.org, 2018.

N. Buzun, V. Avanesov, Bootstrap for change point detection, Preprint no. arXiv:1710.07285, Cornell University Library, arXiv.org, 2017.

V. Avanesov, N. Buzun , Changepoint detection in highdimensional covariance structure, Preprint no. arXiv:1610.03783, Cornell University Library, arXiv.org, 2016.
Abstract
In this paper we introduce a novel approach for an important problem of change point detection. Specifically, we are interested in detection of an abrupt change in the covariance structure of a highdimensional random process  a problem, which has applications many areas e.g., neuroimaging and finance. The developed approach is essentially a testing procedure requiring a proper choice of a critical level. To that end a nonstandard bootstrap scheme is proposed and theoretically justified under mild assumptions. Multiscale nature of the approach allows for a tradeoff between sensitivity of changepoint detection and localization of it. The approach can be naturally used in an online setting. A simulation study demonstrates that the approach matches the nominal level of false alarm probability and exhibits high power, outperforming competing approaches. 
A. Andresen, V. Spokoiny, Critical dimension in profile semiparametric estimation, Preprint no. arXiv:1303.4640, Cornell University Library, arXiv.org, 2014.

A. Andresen, A note on critical dimensions in profile semiparametric estimation, Preprint no. arXiv:1410.4709, Cornell University Library, arXiv.org, 2014.

A. Andresen, A result on the bias of sieve profile estimators, Preprint no. arXiv:1406.4045, Cornell University Library, arXiv.org, 2014.

A. Andresen, Finite sample analysis of profile Mestimation in the single index model, Preprint no. arXiv:1406.4052, Cornell University Library, arXiv.org, 2014.

T. Bodnar, Th. Dickhaus, False discovery rate control under Archimedean copula, Preprint no. arXiv:1305.3897, Cornell University Library, arXiv.org, 2013.

V. Spokoiny, W. Wang, W.K. Härdle, Local quantile regression, Preprint no. arXiv:1208.5384, Cornell University Library, arXiv.org, 2012.

V. Spokoiny, Parametric estimation. Finite sample theory, Preprint no. arXiv:1111.3029, Cornell University Library, arXiv.org, 2012.

V. Spokoiny, Roughness penalty, Wilks phenomenon, and Bernsteinvon Mises theorem, Preprint no. arXiv:1205.0498, Cornell University Library, arXiv.org, 2012.

W. Härdle, V. Spokoiny, Local quantile regression, Discussion paper no. 2011005, HumboldtUniversität zu Berlin, SFB 649, 2011.

J. Kampen, On the multivariate Burgers equation and the incompressible NavierStokes equation (electronic only), Preprint no. arXiv:0910.5672, Cornell University Library, arXiv.org, 2009.

S. Arlot, G. Blanchard, E. Roquain, Some nonasymptotic results on resampling in high dimension, I: Confidence regions, II: Multiple tests (electronic only), Preprint no. arXiv:0712.0775, Cornell University Library, arXiv.org, 2009.
Abstract
We study generalized bootstrap confidence regions for the mean of a random vector whose coordinates have an unknown dependency structure. The random vector is supposed to be either Gaussian or to have a symmetric and bounded distribution. The dimensionality of the vector can possibly be much larger than the number of observations and we focus on a nonasymptotic control of the confidence level, following ideas inspired by recent results in learning theory. We consider two approaches, the first based on a concentration principle (valid for a large class of resampling weights) and the second on a direct resampled quantile, specifically using Rademacher weights. Several intermediate results established in the approach based on concentration principles are of selfinterest. We also discuss the question of accuracy when using MonteCarlo approximations of the resampled quantities. We present an application of these results to the onesided and twosided multiple testing problem, in which we derive several resamplingbased stepdown procedures providing a nonasymptotic FWER control. We compare our different procedures in a simulation study, and we show that they can outperform Bonferroni's or Holm's procedures as soon as the observed vector has sufficiently correlated coordinates.