ترغب بنشر مسار تعليمي؟ اضغط هنا

Estimation of Polarized Power Spectra by Gibbs sampling

162   0   0.0 ( 0 )
 نشر من قبل Hans Kristian Eriksen
 تاريخ النشر 2006
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Earlier papers introduced a method of accurately estimating the angular cosmic microwave background (CMB) temperature power spectrum based on Gibbs sampling. Here we extend this framework to polarized data. All advantages of the Gibbs sampler still apply, and exact analysis of mega-pixel polarized data sets is thus feasible. These advantages may be even more important for polarization measurements than for temperature measurements. While approximate methods can alias power from the larger E-mode spectrum into the weaker B-mode spectrum, the Gibbs sampler (or equivalently, exact likelihood evaluations) allows for a statistically optimal separation of these modes in terms of power spectra. To demonstrate the method, we analyze two simulated data sets: 1) a hypothetical future CMBPol mission, with the focus on B-mode estimation; and 2) a Planck-like mission, to highlight the computational feasibility of the method.



قيم البحث

اقرأ أيضاً

The application of the lasso is espoused in high-dimensional settings where only a small number of the regression coefficients are believed to be nonzero. Moreover, statistical properties of high-dimensional lasso estimators are often proved under th e assumption that the correlation between the predictors is bounded. In this vein, coordinatewise methods, the most common means of computing the lasso solution, work well in the presence of low to moderate multicollinearity. The computational speed of coordinatewise algorithms degrades however as sparsity decreases and multicollinearity increases. Motivated by these limitations, we propose the novel Deterministic Bayesian Lasso algorithm for computing the lasso solution. This algorithm is developed by considering a limiting version of the Bayesian lasso. The performance of the Deterministic Bayesian Lasso improves as sparsity decreases and multicollinearity increases, and can offer substantial increases in computational speed. A rigorous theoretical analysis demonstrates that (1) the Deterministic Bayesian Lasso algorithm converges to the lasso solution, and (2) it leads to a representation of the lasso estimator which shows how it achieves both $ell_1$ and $ell_2$ types of shrinkage simultaneously. Connections to other algorithms are also provided. The benefits of the Deterministic Bayesian Lasso algorithm are then illustrated on simulated and real data.
We develop an analytic model for the power spectra of polarized filamentary structures as a way to study the Galactic polarization foreground to the Cosmic Microwave Background. Our approach is akin to the cosmological halo-model framework, and repro duces the main features of the Planck 353 GHz power spectra. We model the foreground as randomly-oriented, three-dimensional, spheroidal filaments, accounting for their projection onto the sky. The main tunable parameters are the distribution of filament sizes, the filament physical aspect ratio, and the dispersion of the filament axis around the local magnetic field direction. The abundance and properties of filaments as a function of size determine the slopes of the foreground power spectra, as we show via scaling arguments. The filament aspect ratio determines the ratio of $B$-mode power to $E$-mode power, and specifically reproduces the Planck-observed dust ratio of one-half when the short axis is roughly one-fourth the length of the long axis. Filament misalignment to the local magnetic field determines the $TE$ cross-correlation, and to reproduce Planck measurements, we need a (three-dimensional) misalignment angle with a root mean squared dispersion of about 50 degrees. These parameters are not sensitive to the particular filament density profile. By artificially skewing the distribution of the misalignment angle, this model can reproduce the Planck-observed (and parity-violating) $TB$ correlation. The skewing of the misalignment angle necessary to explain $TB$ will cause a yet-unobserved, positive $EB$ dust correlation, a possible target for future experiments.
We briefly review our work about the polarized foreground contamination of the Cosmic Microwave Background maps. We start by summarizing the main properties of the polarized cosmological signal, resulting in electric (E) and magnetic (B) components o f the polarization tensor field on the sky. Then we describe our present understanding of sub-degree anisotropies from Galactic synchrotron and from extra-Galactic point sources. We discuss their contamination of the cosmological E and B modes.
Herding is a technique to sequentially generate deterministic samples from a probability distribution. In this work, we propose a continuous herded Gibbs sampler, that combines kernel herding on continuous densities with Gibbs sampling. Our algorithm allows for deterministically sampling from high-dimensional multivariate probability densities, without directly sampling from the joint density. Experiments with Gaussian mixture densities indicate that the L2 error decreases similarly to kernel herding, while the computation time is significantly lower, i.e., linear in the number of dimensions.
A core problem in cognitive science and machine learning is to understand how humans derive semantic representations from perceptual objects, such as color from an apple, pleasantness from a musical chord, or seriousness from a face. Markov Chain Mon te Carlo with People (MCMCP) is a prominent method for studying such representations, in which participants are presented with binary choice trials constructed such that the decisions follow a Markov Chain Monte Carlo acceptance rule. However, while MCMCP has strong asymptotic properties, its binary choice paradigm generates relatively little information per trial, and its local proposal function makes it slow to explore the parameter space and find the modes of the distribution. Here we therefore generalize MCMCP to a continuous-sampling paradigm, where in each iteration the participant uses a slider to continuously manipulate a single stimulus dimension to optimize a given criterion such as pleasantness. We formulate both methods from a utility-theory perspective, and show that the new method can be interpreted as Gibbs Sampling with People (GSP). Further, we introduce an aggregation parameter to the transition step, and show that this parameter can be manipulated to flexibly shift between Gibbs sampling and deterministic optimization. In an initial study, we show GSP clearly outperforming MCMCP; we then show that GSP provides novel and interpretable results in three other domains, namely musical chords, vocal emotions, and faces. We validate these results through large-scale perceptual rating experiments. The final experiments use GSP to navigate the latent space of a state-of-the-art image synthesis network (StyleGAN), a promising approach for applying GSP to high-dimensional perceptual spaces. We conclude by discussing future cognitive applications and ethical implications.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا