Showing 1–25 of 25 results for author: Woodroofe, M

Quenched central limit theorems for a stationary linear process

Authors: Dalibor Volny, Michael Woodroofe

Abstract: We find a sufficient condition under which a central limit theorem for a stationary linear process is quenched. We find a stationary linear process szatisfying the Maxwell-Woodroofe condition for which the variances of partial sums are o(n), there is a CLT with a convergence towards N(0,1) when dividing by standard deviation of the partial sums, and the CLT is not quenched. The weak invariance pri… ▽ More We find a sufficient condition under which a central limit theorem for a stationary linear process is quenched. We find a stationary linear process szatisfying the Maxwell-Woodroofe condition for which the variances of partial sums are o(n), there is a CLT with a convergence towards N(0,1) when dividing by standard deviation of the partial sums, and the CLT is not quenched. The weak invariance principle does not hold. △ Less

Submitted 21 May, 2015; v1 submitted 9 April, 2015; originally announced April 2015.

MSC Class: 60F05; 60G10; 60G42; 28D05

Baseline zone estimation in two dimensions

Authors: Atul Mallik, Moulinath Banerjee, Michael Woodroofe

Abstract: We consider the problem of estimating the region on which a non-parametric regression function is at its baseline level in two dimensions. The baseline level typically corresponds to the minimum/maximum of the function and estimating such regions or their complements is pertinent to several problems arising in edge estimation, environmental statistics, fMRI and related fields. We assume the baseli… ▽ More We consider the problem of estimating the region on which a non-parametric regression function is at its baseline level in two dimensions. The baseline level typically corresponds to the minimum/maximum of the function and estimating such regions or their complements is pertinent to several problems arising in edge estimation, environmental statistics, fMRI and related fields. We assume the baseline region to be convex and estimate it via fitting a `stump' function to approximate $p$-values obtained from tests for deviation of the regression function from its baseline level. The estimates, obtained using an algorithm originally developed for constructing convex contours of a density, are studied in two different sampling settings, one where several responses can be obtained at a number of different covariate-levels (dose-response) and the other involving limited number of response values per covariate (standard regression). The shape of the baseline region and the smoothness of the regression function at its boundary play a critical role in determining the rate of convergence of our estimate: for a regression function which is `p-regular' at the boundary of the convex baseline region, our estimate converges at a rate $N^{2/(4p+3)}$ in the dose-response setting, $N$ being the total budget, and its analogue in the standard regression setting converges at a rate of $N^{1/(2p+2)}$. Extensions to non-convex baseline regions are explored as well. △ Less

Submitted 22 December, 2013; originally announced December 2013.

Bootstrap confidence intervals for isotonic estimators in a stereological problem

Authors: Bodhisattva Sen, Michael Woodroofe

Abstract: Let $\mathbf{X}=(X_1,X_2,X_3)$ be a spherically symmetric random vector of which only $(X_1,X_2)$ can be observed. We focus attention on estimating F, the distribution function of the squared radius $Z:=X_1^2+X_2^2+X_3^2$, from a random sample of $(X_1,X_2)$. Such a problem arises in astronomy where $(X_1,X_2,X_3)$ denotes the three dimensional position of a star in a galaxy but we can only observ… ▽ More Let $\mathbf{X}=(X_1,X_2,X_3)$ be a spherically symmetric random vector of which only $(X_1,X_2)$ can be observed. We focus attention on estimating F, the distribution function of the squared radius $Z:=X_1^2+X_2^2+X_3^2$, from a random sample of $(X_1,X_2)$. Such a problem arises in astronomy where $(X_1,X_2,X_3)$ denotes the three dimensional position of a star in a galaxy but we can only observe the projected stellar positions $(X_1,X_2)$. We consider isotonic estimators of F and derive their limit distributions. The results are nonstandard with a rate of convergence $\sqrt{n/{\log n}}$. The isotonized estimators of F have exactly half the limiting variance when compared to naive estimators, which do not incorporate the shape constraint. We consider the problem of constructing point-wise confidence intervals for F, state sufficient conditions for the consistency of a bootstrap procedure, and show that the conditions are met by the conventional bootstrap method (generating samples from the empirical distribution function). △ Less

Submitted 23 November, 2012; originally announced November 2012.

Comments: Published in at http://dx.doi.org/10.3150/12-BEJ378 the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm)

Report number: IMS-BEJ-BEJ378

Journal ref: Bernoulli 2012, Vol. 18, No. 4, 1249-1266

On the asymptotic normality of kernel density estimators for linear random fields

Authors: Yizao Wang, Michael Woodroofe

Abstract: We establish sufficient conditions for the asymptotic normality of kernel density estimators, applied to causal linear random fields. Our conditions on the coefficients of linear random fields are weaker than known results, although our assumption on the bandwidth is not minimal. The proof is based on the $m$-approximation method. As a key step, we prove a central limit theorem for triangular arra… ▽ More We establish sufficient conditions for the asymptotic normality of kernel density estimators, applied to causal linear random fields. Our conditions on the coefficients of linear random fields are weaker than known results, although our assumption on the bandwidth is not minimal. The proof is based on the $m$-approximation method. As a key step, we prove a central limit theorem for triangular arrays of stationary $m$-dependent random fields with unbounded $m$. We also apply a moment inequality recently established for stationary random fields. △ Less

Submitted 31 December, 2011; originally announced January 2012.

A New Condition for the Invariance Principle for Stationary Random Fields

Authors: Yizao Wang, Michael Woodroofe

Abstract: We establish a central limit theorem and an invariance principle for stationary random fields, with projective-type conditions. Our result is obtained via an m-dependent approximation method. As applications, we establish invariance principles for orthomartingales and functionals of linear random fields. We establish a central limit theorem and an invariance principle for stationary random fields, with projective-type conditions. Our result is obtained via an m-dependent approximation method. As applications, we establish invariance principles for orthomartingales and functionals of linear random fields. △ Less

Submitted 10 April, 2012; v1 submitted 26 January, 2011; originally announced January 2011.

Inconsistency of bootstrap: The Grenander estimator

Authors: Bodhisattva Sen, Moulinath Banerjee, Michael Woodroofe

Abstract: In this paper, we investigate the (in)-consistency of different bootstrap methods for constructing confidence intervals in the class of estimators that converge at rate $n^{1/3}$. The Grenander estimator, the nonparametric maximum likelihood estimator of an unknown nonincreasing density function $f$ on $[0,\infty)$, is a prototypical example. We focus on this example and explore different approach… ▽ More In this paper, we investigate the (in)-consistency of different bootstrap methods for constructing confidence intervals in the class of estimators that converge at rate $n^{1/3}$. The Grenander estimator, the nonparametric maximum likelihood estimator of an unknown nonincreasing density function $f$ on $[0,\infty)$, is a prototypical example. We focus on this example and explore different approaches to constructing bootstrap confidence intervals for $f(t_0)$, where $t_0\in(0,\infty)$ is an interior point. We find that the bootstrap estimate, when generating bootstrap samples from the empirical distribution function $\mathbb{F}_n$ or its least concave majorant $\tilde{F}_n$, does not have any weak limit in probability. We provide a set of sufficient conditions for the consistency of any bootstrap method in this example and show that bootstrapping from a smoothed version of $\tilde{F}_n$ leads to strongly consistent estimators. The $m$ out of $n$ bootstrap method is also shown to be consistent while generating samples from $\mathbb{F}_n$ and $\tilde{F}_n$. △ Less

Submitted 19 October, 2010; originally announced October 2010.

Comments: Published in at http://dx.doi.org/10.1214/09-AOS777 the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOS-AOS777

Journal ref: Annals of Statistics 2010, Vol. 38, No. 4, 1953-1977

A Central Limit Theorem For Linear Random Fields

Authors: Atul Mallik, Michael Woodroofe

Abstract: A Central Limit Theorem is proved for linear random fields when sums are taken over finite disjoint union of rectangles. The approach does not rely upon the use of Beveridge Nelson decomposition and the conditions needed are similar to those given by Ibragimov for linear processes. When specializing this result to the case when sums are being taken over rectangles, a complete analogue of Ibragimov… ▽ More A Central Limit Theorem is proved for linear random fields when sums are taken over finite disjoint union of rectangles. The approach does not rely upon the use of Beveridge Nelson decomposition and the conditions needed are similar to those given by Ibragimov for linear processes. When specializing this result to the case when sums are being taken over rectangles, a complete analogue of Ibragimov result is obtained with a lot of uniformity. △ Less

Submitted 13 July, 2010; v1 submitted 8 July, 2010; originally announced July 2010.

Quenched Central Limit Theorems for Sums of Stationary Processes

Authors: Dalibor Volný, Michael Woodroofe

Abstract: It is shown that the existence of an L^1 co boundary does not imply the quenched version of the central limit theorem. In another result it is shown that Hannan's condition does imply quenched convergence for an appropriately centered version of the sum. It is shown that the existence of an L^1 co boundary does not imply the quenched version of the central limit theorem. In another result it is shown that Hannan's condition does imply quenched convergence for an appropriately centered version of the sum. △ Less

Submitted 9 June, 2010; originally announced June 2010.

MSC Class: 60F05

A central limit theorem for reversible processes with non-linear growth of variance

Authors: Ou Zhao, Michael Woodroofe, Dalibor Volny

Abstract: Kipnis and Varadhan showed that for an additive functional, $S_n$ say, of a reversible Markov chain the condition $E(S_n^{2})/n \to κ\in (0,\infty)$ implies the convergence of the conditional distribution of $S_n/\sqrt{E(S_n^{2}})$, given the starting point, to the standard normal distribution. We revisit this question under the weaker condition, $E(S_n^{2}) = n\ell(n)$, where $\ell$ is a slowly v… ▽ More Kipnis and Varadhan showed that for an additive functional, $S_n$ say, of a reversible Markov chain the condition $E(S_n^{2})/n \to κ\in (0,\infty)$ implies the convergence of the conditional distribution of $S_n/\sqrt{E(S_n^{2}})$, given the starting point, to the standard normal distribution. We revisit this question under the weaker condition, $E(S_n^{2}) = n\ell(n)$, where $\ell$ is a slowly varying function. It is shown by example that the conditional distribution of $S_n/\sqrt{E(S_n^{2}})$ need not converge to the standard normal distribution in this case; and sufficient conditions for convergence to a (possibly non-standard) normal distribution are developed. △ Less

Submitted 22 May, 2010; originally announced May 2010.

Comments: 11 pages

MSC Class: 60F05; 60J05

Streaming motion in Leo I

Authors: Bodhisattva Sen, Moulinath Banerjee, Michael Woodroofe, Mario Mateo, Matthew Walker

Abstract: Whether a dwarf spheroidal galaxy is in equilibrium or being tidally disrupted by the Milky Way is an important question for the study of its dark matter content and distribution. This question is investigated using 328 recent observations from the dwarf spheroidal Leo I. For Leo I, tidal disruption is detected, at least for stars sufficiently far from the center, but the effect appears to be qu… ▽ More Whether a dwarf spheroidal galaxy is in equilibrium or being tidally disrupted by the Milky Way is an important question for the study of its dark matter content and distribution. This question is investigated using 328 recent observations from the dwarf spheroidal Leo I. For Leo I, tidal disruption is detected, at least for stars sufficiently far from the center, but the effect appears to be quite modest. Statistical tools include isotonic and split point estimators, asymptotic theory, and resampling methods. △ Less

Submitted 15 May, 2009; originally announced May 2009.

Comments: Published in at http://dx.doi.org/10.1214/08-AOAS211 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOAS-AOAS211

Journal ref: Annals of Applied Statistics 2009, Vol. 3, No. 1, 96-116

Estimating a monotone trend

Authors: Ou Zhao, Michael Woodroofe

Abstract: Motivated by global warming issues, we consider a time se- ries that consists of a nondecreasing trend observed with station- ary fluctuations, nonparametric estimation of the trend under monotonicity assumption is considered. The rescaled isotonic es- timators at an interior point are shown to converge to Chernoff's distribution under minimal conditions on the stationary errors. Since the isoto… ▽ More Motivated by global warming issues, we consider a time se- ries that consists of a nondecreasing trend observed with station- ary fluctuations, nonparametric estimation of the trend under monotonicity assumption is considered. The rescaled isotonic es- timators at an interior point are shown to converge to Chernoff's distribution under minimal conditions on the stationary errors. Since the isotonic estimators suffer from the spiking problem at the end point, two modifications are proposed. The estima- tion errors for both estimators of the boundary point are shown to have interesting limiting distributions. Approximation accu- racies are assessed through simulations. One highlight of our treatment is the proof of the weak convergence results which involve several recent techniques developed in the study of con- ditional central limit questions. These weak convergences can be shown to hold conditionally given the starting values. △ Less

Submitted 16 December, 2008; originally announced December 2008.

Comments: 20 pages

MSC Class: 62G05; 62E20; 62G08

Clean Kinematic Samples in Dwarf Spheroidals: An Algorithm for Evaluating Membership and Estimating Distribution Parameters When Contamination is Present

Authors: Matthew G. Walker, Mario Mateo, Edward W. Olszewski, Bodhisattva Sen, Michael Woodroofe

Abstract: (abridged) We develop an algorithm for estimating parameters of a distribution sampled with contamination, employing a statistical technique known as ``expectation maximization'' (EM). Given models for both member and contaminant populations, the EM algorithm iteratively evaluates the membership probability of each discrete data point, then uses those probabilities to update parameter estimates… ▽ More (abridged) We develop an algorithm for estimating parameters of a distribution sampled with contamination, employing a statistical technique known as ``expectation maximization'' (EM). Given models for both member and contaminant populations, the EM algorithm iteratively evaluates the membership probability of each discrete data point, then uses those probabilities to update parameter estimates for member and contaminant distributions. The EM approach has wide applicability to the analysis of astronomical data. Here we tailor an EM algorithm to operate on spectroscopic samples obtained with the Michigan-MIKE Fiber System (MMFS) as part of our Magellan survey of stellar radial velocities in nearby dwarf spheroidal (dSph) galaxies. These samples are presented in a companion paper and contain discrete measurements of line-of-sight velocity, projected position, and Mg index for ~1000 - 2500 stars per dSph, including some fraction of contamination by foreground Milky Way stars. The EM algorithm quantifies both dSph and contaminant distributions, returning maximum-likelihood estimates of the means and variances, as well as the probability that each star is a dSph member. Applied to our MMFS data, the EM algorithm identifies more than 5000 probable dSph members. We test the performance of the EM algorithm on simulated data sets that represent a range of sample size, level of contamination, and amount of overlap between dSph and contaminant velocity distributions. The simulations establish that for samples ranging from large (N ~3000) to small (N~30), the EM algorithm distinguishes members from contaminants and returns accurate parameter estimates much more reliably than conventional methods of contaminant removal (e.g., sigma clipping). △ Less

Submitted 12 November, 2008; originally announced November 2008.

Comments: Accepted for publication in The Astronomical Journal. Download pdf with full-resolution figures from http://www.ast.cam.ac.uk/~walker/dsph_em.pdf

On martingale approximations

Authors: Ou Zhao, Michael Woodroofe

Abstract: Consider additive functionals of a Markov chain $W_k$, with stationary (marginal) distribution and transition function denoted by $π$ and $Q$, say $S_n=g(W_1)+...+g(W_n)$, where $g$ is square integrable and has mean 0 with respect to $π$. If $S_n$ has the form $S_n=M_n+R_n$, where $M_n$ is a square integrable martingale with stationary increments and $E(R_n^2)=o(n)$, then $g$ is said to admit a… ▽ More Consider additive functionals of a Markov chain $W_k$, with stationary (marginal) distribution and transition function denoted by $π$ and $Q$, say $S_n=g(W_1)+...+g(W_n)$, where $g$ is square integrable and has mean 0 with respect to $π$. If $S_n$ has the form $S_n=M_n+R_n$, where $M_n$ is a square integrable martingale with stationary increments and $E(R_n^2)=o(n)$, then $g$ is said to admit a martingale approximation. Necessary and sufficient conditions for such an approximation are developed. Two obvious necessary conditions are $E[E(S_n|W_1)^2]=o(n)$ and $\lim_{n\to \infty}E(S_n^2)/n<\infty$. Assuming the first of these, let $\Vert g\Vert^2_+=\limsup_{n\to \infty}E(S_n^2)/n$; then $\Vert\cdot\Vert_+$ defines a pseudo norm on the subspace of $L^2(π)$ where it is finite. In one main result, a simple necessary and sufficient condition for a martingale approximation is developed in terms of $\Vert\cdot\Vert_+$. Let $Q^*$ denote the adjoint operator to $Q$, regarded as a linear operator from $L^2(π)$ into itself, and consider co-isometries ($QQ^*=I$), an important special case that includes shift processes. In another main result a convenient orthonormal basis for $L_0^2(π)$ is identified along with a simple necessary and sufficient condition for the existence of a martingale approximation in terms of the coefficients of the expansion of $g$ with respect to this basis. △ Less

Submitted 14 November, 2008; v1 submitted 30 August, 2007; originally announced August 2007.

Comments: Published in at http://dx.doi.org/10.1214/07-AAP505 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AAP-AAP505 MSC Class: 60F05 (Primary) 60J10 (Secondary)

Journal ref: Annals of Applied Probability 2008, Vol. 18, No. 5, 1831-1847

Estimating a Polya frequency function_2

Authors: Jayanta Kumar Pal, Michael Woodroofe, Mary Meyer

Abstract: We consider the non-parametric maximum likelihood estimation in the class of Polya frequency functions of order two, viz. the densities with a concave logarithm. This is a subclass of unimodal densities and fairly rich in general. The NPMLE is shown to be the solution to a convex programming problem in the Euclidean space and an algorithm is devised similar to the iterative convex minorant algor… ▽ More We consider the non-parametric maximum likelihood estimation in the class of Polya frequency functions of order two, viz. the densities with a concave logarithm. This is a subclass of unimodal densities and fairly rich in general. The NPMLE is shown to be the solution to a convex programming problem in the Euclidean space and an algorithm is devised similar to the iterative convex minorant algorithm by Jongbleod (1999). The estimator achieves Hellinger consistency when the true density is a PFF_2 itself. △ Less

Submitted 8 August, 2007; originally announced August 2007.

Comments: Published at http://dx.doi.org/10.1214/074921707000000184 in the IMS Lecture Notes Monograph Series (http://www.imstat.org/publications/lecnotes.htm) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-LNMS54-LNMS5418 MSC Class: 62G07; 62G08 (Primary) 90C25 (Secondary)

Journal ref: IMS Lecture Notes Monograph Series 2007, Vol. 54, 239-249

Velocity Dispersion Profiles of Seven Dwarf Spheroidal Galaxies

Authors: Matthew G. Walker, Mario Mateo, Edward W. Olszewski, Oleg Y. Gnedin, Xiao Wang, Bodhisattva Sen, Michael Woodroofe

Abstract: We present stellar velocity dispersion profiles for seven Milky Way dwarf spheroidal (dSph) satellite galaxies. We have measured 8394 line-of-sight velocities (+/- 2.5 km/s) for 6804 stars from high-resolution spectra obtained at the Magellan and MMT telescopes. We combine these new data with previously published velocities to obtain the largest available kinematic samples, which include more th… ▽ More We present stellar velocity dispersion profiles for seven Milky Way dwarf spheroidal (dSph) satellite galaxies. We have measured 8394 line-of-sight velocities (+/- 2.5 km/s) for 6804 stars from high-resolution spectra obtained at the Magellan and MMT telescopes. We combine these new data with previously published velocities to obtain the largest available kinematic samples, which include more than 5500 dSph members. All the measured dSphs have stellar velocity dispersion of order 10 km/s that remains approximately constant with distance from the dSph center, out to and in some cases beyond the radius at which the mean surface brightness falls to the background level. Assuming dSphs reside within dark matter halos characterized by the NFW density profile, we obtain reasonable fits to the empirical velocity dispersion profiles. These fits imply that, among the seven dSphs, M_vir ~ 10^[8-9] M_sun. The mass enclosed at a radius of 600 pc, the region common to all data sets, ranges from (2-7) x 10^7 M_sun . △ Less

Submitted 31 July, 2007; originally announced August 2007.

Comments: Accepted for publication in the Astrophysical Journal Letters

The Michigan/MIKE Fiber System Survey of Stellar Radial Velocities in Dwarf Spheroidal Galaxies: Acquisition and Reduction of Data

Authors: Matthew G. Walker, Mario Mateo, Edward W. Olszewski, Rebecca Bernstein, Bodhisattva Sen, Michael Woodroofe

Abstract: We introduce a stellar velocity survey of dwarf spheroidal galaxies, undertaken using the Michigan/MIKE Fiber System (MMFS) at the Magellan/Clay 6.5 m telescope at Las Campanas Observatory. As of 2006 November we have used MMFS to collect 6415 high-resolution (R= 20000-25000) spectra from 5180 stars in four dwarf spheroidal galaxies: Carina, Fornax, Sculptor and Sextans. Spectra sample the range… ▽ More We introduce a stellar velocity survey of dwarf spheroidal galaxies, undertaken using the Michigan/MIKE Fiber System (MMFS) at the Magellan/Clay 6.5 m telescope at Las Campanas Observatory. As of 2006 November we have used MMFS to collect 6415 high-resolution (R= 20000-25000) spectra from 5180 stars in four dwarf spheroidal galaxies: Carina, Fornax, Sculptor and Sextans. Spectra sample the range 5140-5180 Angstroms, which includes the prominent magnesium triplet absorption feature. We measure radial velocity (RV) to a median precision of 2.0 km/s for stars as faint as V ~ 20.5. From the spectra we also are able to measure the strength of iron and magnesium absorption features using spectral indices that correlate with effective temperature, surface gravity and chemical abundance. Measurement of line strength allows us to identify interloping foreground stars independently of velocity, and to examine the metallicity distribution among dSph members. Here we present detailed descriptions of MMFS, our target selection and spectroscopic observations, the data reduction procedure, and error analysis. We compare our RV results to previously published measurements for individual stars. In some cases we find evidence for a mild, velocity-dependent offset between the RVs we measure using the magnesium triplet and previously published RV measurements derived from the infrared calcium triplet. In companion papers we will present the complete data sets and kinematic analyses of these new observations. △ Less

Submitted 14 March, 2007; v1 submitted 12 March, 2007; originally announced March 2007.

Comments: Accepted for publication in the Astrophysical Journal Supplemement Series. 30 pages, 24 figures (6 color). Dowload pdf including high-resolution figures and complete Table 1 from http://www.astro.lsa.umich.edu/~mgwalker/mmfsfigs/ms.pdf

Law of the iterated logarithm for stationary processes

Authors: Ou Zhao, Michael Woodroofe

Abstract: There has been recent interest in the conditional central limit question for (strictly) stationary, ergodic processes $...,X_{-1},X_0,X_1,...$ whose partial sums $S_n=X_1+...+X_n$ are of the form $S_n=M_n+R_n$, where $M_n$ is a square integrable martingale with stationary increments and $R_n$ is a remainder term for which $E(R_n^2)=o(n)$. Here we explore the law of the iterated logarithm (LIL) f… ▽ More There has been recent interest in the conditional central limit question for (strictly) stationary, ergodic processes $...,X_{-1},X_0,X_1,...$ whose partial sums $S_n=X_1+...+X_n$ are of the form $S_n=M_n+R_n$, where $M_n$ is a square integrable martingale with stationary increments and $R_n$ is a remainder term for which $E(R_n^2)=o(n)$. Here we explore the law of the iterated logarithm (LIL) for the same class of processes. Letting $\Vert\cdot\Vert$ denote the norm in $L^2(P)$, a sufficient condition for the partial sums of a stationary process to have the form $S_n=M_n+R_n$ is that $n^{-3/2}\Vert E(S_n|X_0,X_{-1},...)\Vert$ be summable. A sufficient condition for the LIL is only slightly stronger, requiring $n^{-3/2}\log^{3/2}(n)\Vert E(S_n|X_0,X_{-1},...)\Vert$ to be summable. As a by-product of our main result, we obtain an improved statement of the conditional central limit theorem. Invariance principles are obtained as well. △ Less

Submitted 11 December, 2007; v1 submitted 23 December, 2006; originally announced December 2006.

Comments: Published in at http://dx.doi.org/10.1214/009117907000000079 the Annals of Probability (http://www.imstat.org/aop/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOP-AOP321 MSC Class: 60F15 (Primary); 60F05 (Secondary)

Journal ref: Annals of Probability 2008, Vol. 36, No. 1, 127-142

A non-linear Renewal Theorem with stationary and slowly changing perturbations

Authors: Dong-Yun Kim, Michael Woodroofe

Abstract: Non-linear renewal theory is extended to include random walks perturbed by both a slowly changing sequence and a stationary one. Main results include a version of the Key Renewal Theorem, a derivation of the limiting distribution of the excess over a boundary, and an expansion for the expected first passage time. The formulation is motivated by problems in sequential analysis with staggered entr… ▽ More Non-linear renewal theory is extended to include random walks perturbed by both a slowly changing sequence and a stationary one. Main results include a version of the Key Renewal Theorem, a derivation of the limiting distribution of the excess over a boundary, and an expansion for the expected first passage time. The formulation is motivated by problems in sequential analysis with staggered entry, where subjects enter a study at random times. △ Less

Submitted 22 November, 2006; originally announced November 2006.

Comments: Published at http://dx.doi.org/10.1214/074921706000000680 in the IMS Lecture Notes--Monograph Series (http://www.imstat.org/publications/lecnotes.htm) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-LNMS50-LNMS5013 MSC Class: 60K05 (Primary)

Journal ref: IMS Lecture Notes--Monograph Series 2006, Vol. 50, 76-189

A Kiefer--Wolfowitz comparison theorem for Wicksell's problem

Authors: Xiao Wang, Michael Woodroofe

Abstract: We extend the isotonic analysis for Wicksell's problem to estimate a regression function, which is motivated by the problem of estimating dark matter distribution in astronomy. The main result is a version of the Kiefer--Wolfowitz theorem comparing the empirical distribution to its least concave majorant, but with a convergence rate $n^{-1}\log n$ faster than $n^{-2/3}\log n$. The main result is… ▽ More We extend the isotonic analysis for Wicksell's problem to estimate a regression function, which is motivated by the problem of estimating dark matter distribution in astronomy. The main result is a version of the Kiefer--Wolfowitz theorem comparing the empirical distribution to its least concave majorant, but with a convergence rate $n^{-1}\log n$ faster than $n^{-2/3}\log n$. The main result is useful in obtaining asymptotic distributions for estimators, such as isotonic and smooth estimators. △ Less

Submitted 18 October, 2007; v1 submitted 3 November, 2006; originally announced November 2006.

Comments: Published in at http://dx.doi.org/10.1214/009053606000001604 the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOS-AOS0226 MSC Class: 62G05; 62G08; 62G20 (Primary)

Journal ref: Annals of Statistics 2007, Vol. 35, No. 4, 1559-1575

On Kinematic Substructure in the Sextans Dwarf Spheroidal Galaxy

Authors: Matthew G. Walker, Mario Mateo, Edward W. Olszewski, Jayanta Kumar Pal, Bodhisattva Sen, Michael Woodroofe

Abstract: We present multifiber echelle radial velocity results for 551 stars in the Sextans dwarf spheroidal galaxy and identify 294 stars as probable Sextans members. The projected velocity dispersion profile of the binned data remains flat to a maximum angular radius of $30^{\prime}$. We introduce a nonparametric technique for estimating the projected velocity dispersion surface, and use this to search… ▽ More We present multifiber echelle radial velocity results for 551 stars in the Sextans dwarf spheroidal galaxy and identify 294 stars as probable Sextans members. The projected velocity dispersion profile of the binned data remains flat to a maximum angular radius of $30^{\prime}$. We introduce a nonparametric technique for estimating the projected velocity dispersion surface, and use this to search for kinematic substructure. Our data do not confirm previous reports of a kinematically distinct stellar population at the Sextans center. Instead we detect a region near the Sextans core radius that is kinematically colder than the overall Sextans sample with 95% confidence. △ Less

Submitted 24 March, 2006; originally announced March 2006.

Comments: accepted for publication in ApJ Letters; 4 figures (2 color)

Journal ref: Astrophys.J.642:L41-L44,2006

Internal Kinematics of the Fornax Dwarf Spheroidal Galaxy

Authors: Matthew G. Walker, Mario Mateo, Edward W. Olszewski, Rebecca A. Bernstein, Xiao Wang, Michael Woodroofe

Abstract: We present new radial velocity results for 176 stars in the Fornax dwarf spheroidal galaxy, of which at least 156 are probable Fornax members. We combine with previously published data to obtain a radial velocity sample with 206 stars, of which at least 176 are probable Fornax members. We detect the hint of rotation about an axis near Fornax's morphological minor axis, though the significance of… ▽ More We present new radial velocity results for 176 stars in the Fornax dwarf spheroidal galaxy, of which at least 156 are probable Fornax members. We combine with previously published data to obtain a radial velocity sample with 206 stars, of which at least 176 are probable Fornax members. We detect the hint of rotation about an axis near Fornax's morphological minor axis, though the significance of the rotation signal in the galactic rest frame is sensitive to the adopted value of Fornax's proper motion. Regardless, the observed stellar kinematics are dominated by random motions, and we do not find kinematic evidence of tidal disruption. The projected velocity dispersion profile of the binned dataset remains flat over the sampled region, which reaches a maximum angular radius of 65 arcmin. Single-component King models in which mass follows light fail to reproduce the observed flatness of the velocity dispersion profile. Two-component (luminous plus dark matter) models can reproduce the data, provided the dark component extends sufficiently beyond the luminous component, and the central dark matter density is of the same order as the central luminous density. These requirements suggest a more massive, darker Fornax than standard core-fitting analyses have previously concluded, with M/L_V over the sampled region reaching 10 to 40 times the M/L_V of the luminous component. We also apply a non-parametric mass estimation technique, introduced in a companion paper. Though it is designed to operate on datasets containing velocities for $>$1000 stars, the estimation yields preliminary results suggesting M/L_V ~ 15 inside r < 1.5 kpc. △ Less

Submitted 26 February, 2006; v1 submitted 16 November, 2005; originally announced November 2005.

Comments: 30 pages, 12 figures, 6 tables, AJ in press. Update of Table 1 to identify 3 repeat measurements; modification of equation 6; minor corrections of text (no qualitative change in results)

Journal ref: Astron.J.131:2114-2139,2006; Erratum-ibid.132:968-968,2006

Estimating Dark Matter Distributions

Authors: Xiao Wang, Michael Woodroofe, Matthew G. Walker, Mario Mateo, Edward Olszewski

Abstract: Thanks to instrumental advances, new, very large kinematic datasets for nearby dwarf spheroidal (dSph) galaxies are on the horizon. A key aim of these datasets is to help determine the distribution of dark matter in these galaxies. Past analyses have generally relied on specific dynamical models or highly restrictive dynamical assumptions. We describe a new, non-parametric analysis of the kinema… ▽ More Thanks to instrumental advances, new, very large kinematic datasets for nearby dwarf spheroidal (dSph) galaxies are on the horizon. A key aim of these datasets is to help determine the distribution of dark matter in these galaxies. Past analyses have generally relied on specific dynamical models or highly restrictive dynamical assumptions. We describe a new, non-parametric analysis of the kinematics of nearby dSph galaxies designed to take full advantage of the future large datasets. The method takes as input the projected positions and radial velocities of stars known to be members of the galaxies, but does not use any parametric dynamical model, nor the assumption that the mass distribution follows that of the visible matter. The problem of estimating the radial mass distribution, M(r) (the mass interior to true radius r), is converted into a problem of estimating a regression function non-parametrically. From the Jeans Equation we show that the unknown regression function is subject to fundamental shape restrictions which we exploit in our analysis using statistical techniques borrowed from isotonic estimation and spline smoothing. Simulations indicate that M(r) can be estimated to within a factor of two or better with samples as small as 1000 stars over almost the entire radial range sampled by the kinematic data. The technique is applied to a sample of 181 stars in the Fornax dSph galaxy. We show that the galaxy contains a significant, extended dark halo some ten times more massive than its baryonic component. Though applied here to dSph kinematics, this approach can be used in the analysis of any kinematically hot stellar system in which the radial velocity field is discretely sampled. △ Less

Submitted 11 April, 2005; originally announced April 2005.

Comments: Accepted for publication in The Astrophysical Journal

Journal ref: Astrophys.J. 626 (2005) 145-158

Martingale approximations for sums of stationary processes

Authors: Wei Biao Wu, Michael Woodroofe

Abstract: Approximations to sums of stationary and ergodic sequences by martingales are investigated. Necessary and sufficient conditions for such sums to be asymptotically normal conditionally given the past up to time 0 are obtained. It is first shown that a martingale approximation is necessary for such normality and then that the sums are asymptotically normal if and only if the approximating martinga… ▽ More Approximations to sums of stationary and ergodic sequences by martingales are investigated. Necessary and sufficient conditions for such sums to be asymptotically normal conditionally given the past up to time 0 are obtained. It is first shown that a martingale approximation is necessary for such normality and then that the sums are asymptotically normal if and only if the approximating martingales satisfy a Lindeberg-Feller condition. Using the explicit construction of the approximating martingales, a central limit theorem is derived for the sample means of linear processes. The conditions are not sufficient for the functional version of the central limit theorem. This is shown by an example, and a slightly stronger sufficient condition is given. △ Less

Submitted 6 October, 2004; originally announced October 2004.

Comments: Published by the Institute of Mathematical Statistics (http://www.imstat.org) in the Annals of Probability (http://www.imstat.org/aop/) at http://dx.doi.org/10.1214/009117904000000351

Report number: IMS-AOP-AOP259 MSC Class: 60F17; 60G42; 60F05 (Primary) 60J10. (Secondary)

Journal ref: Annals of Probability 2004, Vol. 32, No. 2, 1674-1690

Setting Confidence Belts

Authors: Byron P. Roe, Michael B. Woodroofe

Abstract: We propose using a Bayes procedure with uniform improper prior to determine credible belts for the mean of a Poisson distribution in the presence of background and for the continuous problem of measuring a non-negative quantity $θ$ with a normally distributed measurement error. Within the Bayesian framework, these belts are optimal. The credible limits are then examined from a frequentist point… ▽ More We propose using a Bayes procedure with uniform improper prior to determine credible belts for the mean of a Poisson distribution in the presence of background and for the continuous problem of measuring a non-negative quantity $θ$ with a normally distributed measurement error. Within the Bayesian framework, these belts are optimal. The credible limits are then examined from a frequentist point of view and found to have good frequentist and conditional frequentist properties. △ Less

Submitted 13 October, 2000; v1 submitted 24 July, 2000; originally announced July 2000.

Comments: 26 pages, 8 figures, submitted to Phys. Rev. D

Journal ref: Phys.Rev. D63 (2001) 013009

Improved Probability Method for Estimating Signal in the Presence of Background

Authors: Byron P. Roe, Michael B. Woodroofe

Abstract: A suggestion is made for improving the Feldman Cousins method of estimating signal counts in the presence of background. The method concentrates on finding essential information about the signal and ignoring extraneous information about background. An appropriate method is found which uses the condition that the number of background events obtained does not exceed the total number of events obta… ▽ More A suggestion is made for improving the Feldman Cousins method of estimating signal counts in the presence of background. The method concentrates on finding essential information about the signal and ignoring extraneous information about background. An appropriate method is found which uses the condition that the number of background events obtained does not exceed the total number of events obtained. Several alternative approaches are explored. △ Less

Submitted 11 August, 1999; v1 submitted 18 December, 1998; originally announced December 1998.

Comments: Modified 12/21 for singlespace to save trees, 9 pages, 1 figure. Modified 8/11/99 to add small modifications made for the Phys. Rev. article

Journal ref: Phys.Rev. D60 (1999) 053009