# NOT RUN {
## 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)
# }
Run the code above in your browser using DataCamp Workspace