diff --git a/NEWS.md b/NEWS.md index 9249961..682224e 100644 --- a/NEWS.md +++ b/NEWS.md @@ -1,6 +1,6 @@ # POMADE 0.1.0.9999 -* Updated version of power calculation of the CE-RVE model in accordance with corrigendum +* Corrected formula for expectation of tau-squared estimator in the CE-RVE model, in accordance with corrigendum. # POMADE 0.1.0 diff --git a/R/power_MADE.R b/R/power_MADE.R index bdc9e84..cc3632e 100644 --- a/R/power_MADE.R +++ b/R/power_MADE.R @@ -469,7 +469,7 @@ power_MADE_single <- a <- sum(1/sigma2j) # Equation 17 in Vembye, Pustejovsky, & Pigott (2023, corrigendum) - tau2_e <- tau^2 + omega^2 * (1 - 1/a^2 * sum(1 / (kj * sigma2j^2))) / (1 - 1/a^2 * sum(1 / sigma2j^2)) + tau2_e <- tau^2 + omega^2 * (1 - sum(1 / (kj * sigma2j^2)) / a^2) / (1 - sum(1 / sigma2j^2) / a^2) # _dd = dotdot wj_dd <- 1 / (sigma2j + tau2_e)