Fixing the HI Width with Planck Data
By Elijah Bernstein-Cooper, August 12, 2015, 0 comments.

Table of Contents



HI and CO spectra

Yesterday’s post incorrectly showed the CO spectra for each cloud. Below are the corrected CO spectra included within the region outlined for each cloud.


Perseus

Taurus

California

Figure 1

Median HI spectra with model fit in purple, and the HI velocity range used as the gray shaded region. The velocity widths are consistent with what was done in Imara et al. (2012).


Binning Errors

We should make sure the pixels used in the MLE calculation have reasonable errors and follow a linear trend. In previous posts I was incorrectly displaying the errors in this plot. After

When binning the image, we will quote the mean of the pixels within the bin as the binned value, and the standard deviation of the pixels as the uncertainty of the binned value.

In the case where the values of the pixels vary greatly compared to the error of each pixel, the uncertainty of the bin will be dominated by the spread in the pixel values. However in the case where the values of pixels are constant compared to the error of each pixel, the uncertainty of the binned value will be dominated by the errors of each pixel.

We are interested in two different quantities:

  • Standard error, : This is the error on our estimate of the mean. This depends on the sample size. where is the standard deviation for each element. If the variances are the same .

  • Standard deviation, : the spread in the population. This is an intrinsic property of the data which will not change with sample size.

Since these two properties are independent of one another, we can add their uncertainties quadratically to get the total uncertainty of the sample, :

Lee+12 IRIS , residual masking


Figure 2

Left: masked vs. N(HI), right: N(H) vs. N(HI) for Perseus region.