better fix for #58 (commit 7390235 did not take care of an arbitrary …

…number of these special case)
ycroissant · Jul 2, 2024 · 03aa426 · 03aa426
1 parent 23fbd4d
commit 03aa426
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/R/tool_vcovG.R b/R/tool_vcovG.R
@@ -503,7 +503,7 @@ vcovG.plm <- function(x, type = c("HC0", "sss", "HC1", "HC2", "HC3", "HC4"),
         # drop groups in Ti that are now empty (group had 1 observation before first-differencing, hence 0 after)
         # and adjust n0 due to same reason
         Ti <- Ti[!drop]
-        n0 <- n0 - 1
+        n0 <- n0 - sum(drop)
       }
       t0 <- t0 - 1
 
@@ -1033,7 +1033,7 @@ vcovBK.plm <- function(x, type = c("HC0", "HC1", "HC2", "HC3", "HC4"),
         # drop groups in Ti that are now empty (group had 1 observation before first-differencing, hence 0 after)
         # and adjust n0 due to same reason
         Ti <- Ti[!drop]
-        n0 <- n0 - 1
+        n0 <- n0 - sum(drop)
       }
       t0 <- t0 - 1
     }