# How to estimate uncertainty of measurements of equivalent widths?

I'm measuring equivalent widths of absorption lines using a spectrum of a star. I make two or three measurements of each line by making reasonable gaussian fits of the line with IRAF's splot tool. Then I calculate the mean of the measurements, which serves as my final equivalent width estimate.

What is a good way of estimating the uncertainty of this measurement?

## My current method

I'm currently using half of the range for the uncertainty. For example, if I made two measurements 10 and 16 mA (milliangstrom), then the mean is 13 mA and uncertainty is 3 mA. This gives the estimate of equivalent width to be 13±3 mA. Do you see any problems with this method of estimating uncertainty?

What it does not do is evaluate the uncertainty in the EW caused by the quality or signal-to-noise ratio of the data itself. You might assess this using some rule-of-thumb formulae for a Gaussian line, e.g. $$\Delta EW \simeq 1.5 \frac{\sqrt{fp}}{{\rm SNR}},$$ (eqn 6 of Cayrel de Strobel 1988) where $$f$$ is the FWHM of the spectral line (in wavelength units), $$p$$ is the size of one pixel in wavelength units and SNR is the signal-to-noise ratio of the data in an average pixel. Or you could take a synthetic spectrum and add some artificial noise to it with the appropriate properties and measure the EW of several randomisations of the same spectrum, taking the standard deviation of your EW measurements to indicate the EW uncertainty for a particular level of signal-to-noise ratio.