Abstract
In this erratum, we retract the upper limits on the 21 cm power spectrum presented in the published article. The published article reported an upper limit on δ _{21} ^{2}(k ) of (22.4 mK)^{2} at z= 8.4 in the range 0.15 < k < 0.5h Mpc^{1}. This analysis underestimated the level of signal loss, or attenuation of the target cosmological 21 cm signal associated with the chosen power spectrum estimator, and also underestimated the statistical error on those estimates. A revised result, with a new analysis, is presented in M. K. Kolopanis et al. (2018, in preparation). Below, we briefly summarize the errors in the original analysis and how they are corrected. For an indepth analysis and discussion of the errors, we refer the reader to Cheng et al. (2018). Signal loss was expected in the original analysis because the covariance matrices, C, used to weight the unnormalized bandpower estimates, qa, in (Formula Presented) were empirically estimated from a timeaveraged finite ensemble of the data, x, such that (Formula Presented). While the true covariance C leads to an inherently unbiased lossless estimator of the power spectrum, using an empirically estimated Ĉ can lead to signal loss. Specifically, weighting data by an empirically estimated covariance carries the risk of overfitting and downweighting EoR fluctuations that are coupled to the data. In Cheng et al. (2018), it is shown that these couplings are especially strong in the fringerate filtered PAPER64 data set. The first and most impactful error relates to the method by which signal loss was estimated. To assess signal loss from the empirically estimated covariance matrix, different realizations of mock cosmological signals e of known amplitudes are added to the original data to form a new data vector, (Formula Presented). New covariance matrices, (Formula Presented) are used to estimate unnormalized bandpowers (Formula Presented), which can be written as (Formula Presented). The normalized power estimate can then be compared to the known injected power in e to estimate signal loss. The key error in the previous analysis was to assume that, since e was statistically independent of x, that the two middle crossterms in Equation (2) would average to zero in an ensemble. However, as shown in Cheng et al. (2018) and Switzer et al. (2015), these crossterms can contain significant negative power because Ĉ_{r} contains information that correlates the two vectors. Ignoring these crossterms leads to a significant underestimate of signal loss. As a result, we presented negligible signal loss in our original analysis, when in fact approximately 99.99% of the signal was removed (Cheng et al. 2018). Correcting for the actual signal loss is the biggest factor revising the upper limit on 21 D2 . The second mistake made in the original analysis was to underestimate the statistical errors in the reported power spectrum estimates. The original analysis used a bootstrap resampling technique on power spectral measurements over the baseline and time axes. However, fringerate filtering introduces significant correlations in the data along the time axis. As is discussed in Cheng et al. (2018), bootstrapping across correlated samples can result in a significant underestimate of the variation in the data if the number of resamplings is not equal to the number of independent samples in the data, as in the case of the original analysis. The error bars associated with this oversampling were underestimated by approximately a factor of 2 (in mK). The revised analysis in M. K. Kolopanis et al. (2018, in preparation) only applies bootstrap resampling across the baseline axis to avoid this problem. The mistake in estimating the statistical errors should have become apparent when comparing results to our theoretical thermal noise sensitivity. Unfortunately, a third miscalculation was made in estimating the thermal noise sensitivity. As detailed in Cheng et al. (2018), this miscalculation stemmed from numerous small mismatches between the idealized analysis pipeline used to estimate sensitivity and the actual analysis applied to the data. As a result, our estimated thermal noise sensitivity was approximately a factor of 3 low (in mK), leading to the mistaken impression that our error bars were consistent with the level of thermal noise. In summary, we retract the power spectrum results shown in Figures 18 and 20 in the published article. Results that relied on the original limits, including those presented in Figure 21, are retracted. Additionally, the companion paper to the original manuscript, Pober et al. (2015), used the original limits to place constraints on the spin temperature of the intergalactic medium (IGM) at z=8.4. Our revised limits do not place significant constraints on the IGM temperature, and the results of Figure 4 from Pober et al. (2015) should be disregarded. However, we note that their analysis would still be relevant should a future experiment place constraints on the 21 cm signal similar to those claimed in the published article. An updated analysis of this same data set is presented in M. K. Kolopanis et al. (2018, in preparation), where these revised results are put into context with measurements at other redshifts.
Original language  English (US) 

Article number  201 
Journal  Astrophysical Journal 
Volume  863 
Issue number  2 
DOIs 

State  Published  Aug 20 2018 
Fingerprint
ASJC Scopus subject areas
 Astronomy and Astrophysics
 Space and Planetary Science
Cite this
Erratum : PAPER64 constraints on reionization: The 21 cm power spectrum at z = 8.4 (Astrophysical Journal (2015) 809 (61) DOI: 10.1088/0004637X/809/1/61). / Ali, Zaki S.; Parsons, Aaron R.; Zheng, Haoxuan; Pober, Jonathan C.; Liu, Adrian; Aguirre, James E.; Bradley, Richard F.; Bernardi, Gianni; Carilli, Chris L.; Cheng, Carina; Deboer, David R.; Dexter, Matthew R.; Grobbelaar, Jasper; Horrell, Jasper; Jacobs, Daniel; Klima, Pat; Macmahon, David H.E.; Maree, Matthys; Moore, David F.; Razavi, Nima; Stefan, Irina I.; Walbrugh, William P.; Walker, Andre.
In: Astrophysical Journal, Vol. 863, No. 2, 201, 20.08.2018.Research output: Contribution to journal › Comment/debate
}
TY  JOUR
T1  Erratum
T2  PAPER64 constraints on reionization: The 21 cm power spectrum at z = 8.4 (Astrophysical Journal (2015) 809 (61) DOI: 10.1088/0004637X/809/1/61)
AU  Ali, Zaki S.
AU  Parsons, Aaron R.
AU  Zheng, Haoxuan
AU  Pober, Jonathan C.
AU  Liu, Adrian
AU  Aguirre, James E.
AU  Bradley, Richard F.
AU  Bernardi, Gianni
AU  Carilli, Chris L.
AU  Cheng, Carina
AU  Deboer, David R.
AU  Dexter, Matthew R.
AU  Grobbelaar, Jasper
AU  Horrell, Jasper
AU  Jacobs, Daniel
AU  Klima, Pat
AU  Macmahon, David H.E.
AU  Maree, Matthys
AU  Moore, David F.
AU  Razavi, Nima
AU  Stefan, Irina I.
AU  Walbrugh, William P.
AU  Walker, Andre
PY  2018/8/20
Y1  2018/8/20
N2  In this erratum, we retract the upper limits on the 21 cm power spectrum presented in the published article. The published article reported an upper limit on δ 21 2(k ) of (22.4 mK)2 at z= 8.4 in the range 0.15 < k < 0.5h Mpc1. This analysis underestimated the level of signal loss, or attenuation of the target cosmological 21 cm signal associated with the chosen power spectrum estimator, and also underestimated the statistical error on those estimates. A revised result, with a new analysis, is presented in M. K. Kolopanis et al. (2018, in preparation). Below, we briefly summarize the errors in the original analysis and how they are corrected. For an indepth analysis and discussion of the errors, we refer the reader to Cheng et al. (2018). Signal loss was expected in the original analysis because the covariance matrices, C, used to weight the unnormalized bandpower estimates, qa, in (Formula Presented) were empirically estimated from a timeaveraged finite ensemble of the data, x, such that (Formula Presented). While the true covariance C leads to an inherently unbiased lossless estimator of the power spectrum, using an empirically estimated Ĉ can lead to signal loss. Specifically, weighting data by an empirically estimated covariance carries the risk of overfitting and downweighting EoR fluctuations that are coupled to the data. In Cheng et al. (2018), it is shown that these couplings are especially strong in the fringerate filtered PAPER64 data set. The first and most impactful error relates to the method by which signal loss was estimated. To assess signal loss from the empirically estimated covariance matrix, different realizations of mock cosmological signals e of known amplitudes are added to the original data to form a new data vector, (Formula Presented). New covariance matrices, (Formula Presented) are used to estimate unnormalized bandpowers (Formula Presented), which can be written as (Formula Presented). The normalized power estimate can then be compared to the known injected power in e to estimate signal loss. The key error in the previous analysis was to assume that, since e was statistically independent of x, that the two middle crossterms in Equation (2) would average to zero in an ensemble. However, as shown in Cheng et al. (2018) and Switzer et al. (2015), these crossterms can contain significant negative power because Ĉr contains information that correlates the two vectors. Ignoring these crossterms leads to a significant underestimate of signal loss. As a result, we presented negligible signal loss in our original analysis, when in fact approximately 99.99% of the signal was removed (Cheng et al. 2018). Correcting for the actual signal loss is the biggest factor revising the upper limit on 21 D2 . The second mistake made in the original analysis was to underestimate the statistical errors in the reported power spectrum estimates. The original analysis used a bootstrap resampling technique on power spectral measurements over the baseline and time axes. However, fringerate filtering introduces significant correlations in the data along the time axis. As is discussed in Cheng et al. (2018), bootstrapping across correlated samples can result in a significant underestimate of the variation in the data if the number of resamplings is not equal to the number of independent samples in the data, as in the case of the original analysis. The error bars associated with this oversampling were underestimated by approximately a factor of 2 (in mK). The revised analysis in M. K. Kolopanis et al. (2018, in preparation) only applies bootstrap resampling across the baseline axis to avoid this problem. The mistake in estimating the statistical errors should have become apparent when comparing results to our theoretical thermal noise sensitivity. Unfortunately, a third miscalculation was made in estimating the thermal noise sensitivity. As detailed in Cheng et al. (2018), this miscalculation stemmed from numerous small mismatches between the idealized analysis pipeline used to estimate sensitivity and the actual analysis applied to the data. As a result, our estimated thermal noise sensitivity was approximately a factor of 3 low (in mK), leading to the mistaken impression that our error bars were consistent with the level of thermal noise. In summary, we retract the power spectrum results shown in Figures 18 and 20 in the published article. Results that relied on the original limits, including those presented in Figure 21, are retracted. Additionally, the companion paper to the original manuscript, Pober et al. (2015), used the original limits to place constraints on the spin temperature of the intergalactic medium (IGM) at z=8.4. Our revised limits do not place significant constraints on the IGM temperature, and the results of Figure 4 from Pober et al. (2015) should be disregarded. However, we note that their analysis would still be relevant should a future experiment place constraints on the 21 cm signal similar to those claimed in the published article. An updated analysis of this same data set is presented in M. K. Kolopanis et al. (2018, in preparation), where these revised results are put into context with measurements at other redshifts.
AB  In this erratum, we retract the upper limits on the 21 cm power spectrum presented in the published article. The published article reported an upper limit on δ 21 2(k ) of (22.4 mK)2 at z= 8.4 in the range 0.15 < k < 0.5h Mpc1. This analysis underestimated the level of signal loss, or attenuation of the target cosmological 21 cm signal associated with the chosen power spectrum estimator, and also underestimated the statistical error on those estimates. A revised result, with a new analysis, is presented in M. K. Kolopanis et al. (2018, in preparation). Below, we briefly summarize the errors in the original analysis and how they are corrected. For an indepth analysis and discussion of the errors, we refer the reader to Cheng et al. (2018). Signal loss was expected in the original analysis because the covariance matrices, C, used to weight the unnormalized bandpower estimates, qa, in (Formula Presented) were empirically estimated from a timeaveraged finite ensemble of the data, x, such that (Formula Presented). While the true covariance C leads to an inherently unbiased lossless estimator of the power spectrum, using an empirically estimated Ĉ can lead to signal loss. Specifically, weighting data by an empirically estimated covariance carries the risk of overfitting and downweighting EoR fluctuations that are coupled to the data. In Cheng et al. (2018), it is shown that these couplings are especially strong in the fringerate filtered PAPER64 data set. The first and most impactful error relates to the method by which signal loss was estimated. To assess signal loss from the empirically estimated covariance matrix, different realizations of mock cosmological signals e of known amplitudes are added to the original data to form a new data vector, (Formula Presented). New covariance matrices, (Formula Presented) are used to estimate unnormalized bandpowers (Formula Presented), which can be written as (Formula Presented). The normalized power estimate can then be compared to the known injected power in e to estimate signal loss. The key error in the previous analysis was to assume that, since e was statistically independent of x, that the two middle crossterms in Equation (2) would average to zero in an ensemble. However, as shown in Cheng et al. (2018) and Switzer et al. (2015), these crossterms can contain significant negative power because Ĉr contains information that correlates the two vectors. Ignoring these crossterms leads to a significant underestimate of signal loss. As a result, we presented negligible signal loss in our original analysis, when in fact approximately 99.99% of the signal was removed (Cheng et al. 2018). Correcting for the actual signal loss is the biggest factor revising the upper limit on 21 D2 . The second mistake made in the original analysis was to underestimate the statistical errors in the reported power spectrum estimates. The original analysis used a bootstrap resampling technique on power spectral measurements over the baseline and time axes. However, fringerate filtering introduces significant correlations in the data along the time axis. As is discussed in Cheng et al. (2018), bootstrapping across correlated samples can result in a significant underestimate of the variation in the data if the number of resamplings is not equal to the number of independent samples in the data, as in the case of the original analysis. The error bars associated with this oversampling were underestimated by approximately a factor of 2 (in mK). The revised analysis in M. K. Kolopanis et al. (2018, in preparation) only applies bootstrap resampling across the baseline axis to avoid this problem. The mistake in estimating the statistical errors should have become apparent when comparing results to our theoretical thermal noise sensitivity. Unfortunately, a third miscalculation was made in estimating the thermal noise sensitivity. As detailed in Cheng et al. (2018), this miscalculation stemmed from numerous small mismatches between the idealized analysis pipeline used to estimate sensitivity and the actual analysis applied to the data. As a result, our estimated thermal noise sensitivity was approximately a factor of 3 low (in mK), leading to the mistaken impression that our error bars were consistent with the level of thermal noise. In summary, we retract the power spectrum results shown in Figures 18 and 20 in the published article. Results that relied on the original limits, including those presented in Figure 21, are retracted. Additionally, the companion paper to the original manuscript, Pober et al. (2015), used the original limits to place constraints on the spin temperature of the intergalactic medium (IGM) at z=8.4. Our revised limits do not place significant constraints on the IGM temperature, and the results of Figure 4 from Pober et al. (2015) should be disregarded. However, we note that their analysis would still be relevant should a future experiment place constraints on the 21 cm signal similar to those claimed in the published article. An updated analysis of this same data set is presented in M. K. Kolopanis et al. (2018, in preparation), where these revised results are put into context with measurements at other redshifts.
UR  http://www.scopus.com/inward/record.url?scp=85052436705&partnerID=8YFLogxK
UR  http://www.scopus.com/inward/citedby.url?scp=85052436705&partnerID=8YFLogxK
U2  10.3847/15384357/aad7b4
DO  10.3847/15384357/aad7b4
M3  Comment/debate
AN  SCOPUS:85052436705
VL  863
JO  Astrophysical Journal
JF  Astrophysical Journal
SN  0004637X
IS  2
M1  201
ER 