pwSpCorr: Power spectra estimate via correlation

In rcosmo: Cosmic Microwave Background Data Analysis

Description

This function provides an angular power spectra estimate using the values of the sample correlations. The approach is based on Lawson-Hanson algorithm for non-negative least squares.

Usage

pwSpCorr(corcmb, lmax = 20 * length(corcmb$u))

Arguments

corcmb	An object of the class CMBCorrelation.
lmax	A number of angular power spectra components to estimate

Value

A data frame which first column is 1-d grid of l values from 0 to lmax. The second column is estimated angular power spectra components on this grid.

References

Formula (2.1) in Baran A., Terdik G. Power spectrum estimation of spherical random fields based on covariances. Annales Mathematicae et Informaticae 44 (2015) pp. 15–22.

Examples

## Download the map first
# downloadCMBMap(foreground = "smica", nside = 1024)
# df <- CMBDataFrame("CMB_map_smica1024.fits")
# Corrf <- corrCMB(df, max.dist = 0.1, num.bins = 30,
# sample.size=10000)
# pw <- pwSpCorr(Corrf)

rcosmo documentation built on Dec. 11, 2021, 9:29 a.m.