exprob: Threshold exceedance probability

In rcosmo: Cosmic Microwave Background Data Analysis

Description

This function returns an estimated exceedance probability for the specified CMBDataFrame column intensities, threshold level alpha and CMBWindow region.

Usage

exprob(cmbdf, win, alpha, intensities = "I")

Arguments

cmbdf	A CMBDataFrame.
win	A CMBWindow
alpha	A numeric threshold level.
intensities	The name of the column in cmbdf that contains the measured values.

Value

Estimated threshold exceedance probability

Examples

## Download the map first
# downloadCMBMap(foreground = "smica", nside = 1024)
# df <- CMBDataFrame("CMB_map_smica1024.fits")
# cmbdf <- sampleCMB(df, sample.size = 1000)

# intensities <- "I"
# alpha <- mean(cmbdf[,intensities, drop = TRUE])
# alpha

# win1 <- CMBWindow(theta = c(0,pi/2,pi/2), phi = c(0,0,pi/2))
# exprob(cmbdf, win1, alpha)

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