exprob: Threshold exceedance probability

Description Usage Arguments Value Examples

Description

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

Usage

1
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
## 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)

frycast/rcosmo documentation built on June 4, 2019, 5:13 p.m.