pwSpCorr: Power spectra estimate via correlation

Description Usage Arguments Value References Examples

View source: R/Statistics.R

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

1
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

1
2
3
4
5
6
## 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 Feb. 14, 2020, 5:06 p.m.