21,611 research outputs found
Bayesian Semiparametric Multivariate Density Deconvolution
We consider the problem of multivariate density deconvolution when the
interest lies in estimating the distribution of a vector-valued random variable
but precise measurements of the variable of interest are not available,
observations being contaminated with additive measurement errors. The existing
sparse literature on the problem assumes the density of the measurement errors
to be completely known. We propose robust Bayesian semiparametric multivariate
deconvolution approaches when the measurement error density is not known but
replicated proxies are available for each unobserved value of the random
vector. Additionally, we allow the variability of the measurement errors to
depend on the associated unobserved value of the vector of interest through
unknown relationships which also automatically includes the case of
multivariate multiplicative measurement errors. Basic properties of finite
mixture models, multivariate normal kernels and exchangeable priors are
exploited in many novel ways to meet the modeling and computational challenges.
Theoretical results that show the flexibility of the proposed methods are
provided. We illustrate the efficiency of the proposed methods in recovering
the true density of interest through simulation experiments. The methodology is
applied to estimate the joint consumption pattern of different dietary
components from contaminated 24 hour recalls
The Magnetic Topology of the Weak-Lined T Tauri Star V410 - A Simultaneous Temperature and Magnetic Field Inversion
We present a detailed temperature and magnetic investigation of the T Tauri
star V410 Tau by means of a simultaneous Doppler- and Zeeman-Doppler Imaging.
Moreover we introduce a new line profile reconstruction method based on a
singular value decomposition (SVD) to extract the weak polarized line profiles.
One of the key features of the line profile reconstruction is that the SVD line
profiles are amenable to radiative transfer modeling within our Zeeman-Doppler
Imaging code iMap. The code also utilizes a new iterative regularization scheme
which is independent of any additional surface constraints. To provide more
stability a vital part of our inversion strategy is the inversion of both
Stokes I and Stokes V profiles to simultaneously reconstruct the temperature
and magnetic field surface distribution of V410 Tau. A new image-shear analysis
is also implemented to allow the search for image and line profile distortions
induced by a differential rotation of the star. The magnetic field structure we
obtain for V410 Tau shows a good spatial correlation with the surface
temperature and is dominated by a strong field within the cool polar spot. The
Zeeman-Doppler maps exhibit a large-scale organization of both polarities
around the polar cap in the form of a twisted bipolar structure. The magnetic
field reaches a value of almost 2 kG within the polar region but smaller fields
are also present down to lower latitudes. The pronounced non-axisymmetric field
structure and the non-detection of a differential rotation for V410 Tau
supports the idea of an underlying -type dynamo, which is predicted
for weak-lined T Tauri stars.Comment: Accepted for A&A, 18 pages, 10 figure
Neutron Star Properties with Hyperons
In the light of the recent discovery of a neutron star with a mass accurately
determined to be almost two solar masses, it has been suggested that hyperons
cannot play a role in the equation of state of dense matter in
-equilibrium. We re-examine this issue in the most recent development of
the quark-meson coupling model. Within a relativistic Hartree-Fock approach and
including the full tensor structure at the vector-meson-baryon vertices, we
find that not only must hyperons appear in matter at the densities relevant to
such a massive star but that the maximum mass predicted is completely
consistent with the observation.Comment: Minor correction
Gravitational Wave Burst Source Direction Estimation using Time and Amplitude Information
In this article we study two problems that arise when using timing and
amplitude estimates from a network of interferometers (IFOs) to evaluate the
direction of an incident gravitational wave burst (GWB). First, we discuss an
angular bias in the least squares timing-based approach that becomes
increasingly relevant for moderate to low signal-to-noise ratios. We show how
estimates of the arrival time uncertainties in each detector can be used to
correct this bias. We also introduce a stand alone parameter estimation
algorithm that can improve the arrival time estimation and provide
root-sum-squared strain amplitude (hrss) values for each site. In the second
part of the paper we discuss how to resolve the directional ambiguity that
arises from observations in three non co-located interferometers between the
true source location and its mirror image across the plane containing the
detectors. We introduce a new, exact relationship among the hrss values at the
three sites that, for sufficiently large signal amplitudes, determines the true
source direction regardless of whether or not the signal is linearly polarized.
Both the algorithm estimating arrival times, arrival time uncertainties, and
hrss values and the directional follow-up can be applied to any set of
gravitational wave candidates observed in a network of three non co-located
interferometers. As a case study we test the methods on simulated waveforms
embedded in simulations of the noise of the LIGO and Virgo detectors at design
sensitivity.Comment: 10 pages, 14 figures, submitted to PR
A study of the application of singular perturbation theory
A hierarchical real time algorithm for optimal three dimensional control of aircraft is described. Systematic methods are developed for real time computation of nonlinear feedback controls by means of singular perturbation theory. The results are applied to a six state, three control variable, point mass model of an F-4 aircraft. Nonlinear feedback laws are presented for computing the optimal control of throttle, bank angle, and angle of attack. Real Time capability is assessed on a TI 9900 microcomputer. The breakdown of the singular perturbation approximation near the terminal point is examined Continuation methods are examined to obtain exact optimal trajectories starting from the singular perturbation solutions
- …