Confidence Intervals



next up previous contents
Next: Parameter Recovery Up: Simulations Previous: Simulations

Confidence Intervals

 

The validity of the confidence intervals returned by the Bayesian approach in ELFIT2D can be tested by running ELFIT2D on a number of noise realizations of the same galaxy rotation curve spectrum. A simulation of a typical SIS/LORAL3 spectrum (t = 23600sec, N = 8 e) of a z = 0.35 galaxy with V sin = 130 km/s, disk scale length of 2.0h kpc, [OII] ratio of 1.0, total [OII] flux of 1053 DU and background level of 0 DU was created. Random noise was added to this galaxy spectrum using the IRAF task MKNOISE in the ARTDATA package. Fifty galaxy spectra were created using a different seed for the noise each time. ELFIT2D was run with the same initial conditions on all fifty spectra. The results are given in Table gif and shown in Figure gif.

Columns 2-5 in Table gif show that the measured median values are scattered in an unbiased way around the true parameter values. Column 8 gives for each parameter the fraction of spectra in which the true parameter value was within the 68 confidence interval derived by ELFIT2D from the posterior probability distribution . These confidence intervals include the true value between 0.56 and 0.72 of the time according to the values. The expected value of 0.68 for corresponds to 34 spectra out of 50, and the = 5.8 Poisson error can account for the fluctuations seen in . Therefore, the confidence intervals given by ELFIT2D are valid. Figure gif shows for three noise realizations. It is interesting to note that the skewness of the distributions exhibits the same behavior as in the suspicious coin experiment (see Figure gif).

 

 

 



Luc Simard
Mon Sep 2 12:37:40 PDT 1996