Fix crash in AR when point estimate<0 and CI requested

droodman · Nov 10, 2023 · 7943a02 · 7943a02
1 parent 11e4cdd
commit 7943a02
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/interface.jl b/src/interface.jl
@@ -379,7 +379,7 @@ Function to perform wild-bootstrap-based hypothesis test
 * `clustid::AbstractVecOrMat{<:Integer}`: data vector/matrix of error and bootstrapping cluster identifiers; see notes 
 * `nbootclustvar::Integer=size(clustid,2)`: number of bootstrap-clustering variables
 * `nerrclustvar::Integer=nbootclustvar`: number of error-clustering variables
-* `issorted:Bool=false`: time-saving flag: data matrices are already sort by column types 2, then 3, then 1 (see notes)
+* `issorted:Bool=false`: time-saving flag: data matrices are already sorted by column types 2, then 3, then 1 (see notes)
 * `hetrobust::Bool=true`: true unless errors are treated as iid
 * `feid::AbstractVector{<:Integer}`: data vector for one-way fixed effect group identifier
 * `fedfadj::Integer`: degrees of freedom that fixed effects (if any) consume; defaults to number of FEs

diff --git a/src/plot-CI.jl b/src/plot-CI.jl
@@ -52,7 +52,7 @@ function plot!(o::StrBootTest{T}) where T
   if o.arubin
 		# p = o.dist[1] * o.multiplier
 		# p = ccdf(Chisq{T}(T(o.dof)), o.sqrt ? p^2 : p)
-		halfwidth = abs.(o.confpeak * sqrtNaN((o.dist[1] * o.multiplier) / o.dof) / Phi)  # abs.(o.confpeak) * quantile(Normal{T}(zero(T),one(T)), p/2) / Phi
+		halfwidth = abs.(o.confpeak * sqrtNaN((abs(o.dist[1]) * o.multiplier) / o.dof) / Phi)  # abs.(o.confpeak) * quantile(Normal{T}(zero(T),one(T)), p/2) / Phi
   else
 		halfwidth = T(-1.5) * Phi .* sqrtNaN.(diag(getV(o)))
 		any(isnan.(halfwidth)) && return