diff --git a/NEWS.md b/NEWS.md
index 93002dee8..5c1af8ee2 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -5,6 +5,7 @@
 ## New Features
 * rbmi now uses the `lifecycle` package 
 * Added en-GB spell-check and a corresponding test to the package
+* `ancova_single()` now return residual variance.
 
 ## Documentation
 * Documented the return value (`@return`) of all exported functions and marked
diff --git a/R/ancova.R b/R/ancova.R
index 6c0b3605c..6ff548bbd 100644
--- a/R/ancova.R
+++ b/R/ancova.R
@@ -226,6 +226,11 @@ ancova_single <- function(
     lsm1 <- do.call(lsmeans, args)
 
     x <- list(
+        var = list(
+            est = summary(mod)$sigma^2,
+            se = summary(mod)$sigma^2 / sqrt(df.residual(mod) / 2),
+            df = df.residual(mod)
+        ),
         trt = list(
             est = coef(mod)[[group]],
             se = sqrt(vcov(mod)[group, group]),
diff --git a/tests/testthat/_snaps/pool.md b/tests/testthat/_snaps/pool.md
index 7557ddb71..531adbeac 100644
--- a/tests/testthat/_snaps/pool.md
+++ b/tests/testthat/_snaps/pool.md
@@ -13,12 +13,15 @@
         =========================================================
             parameter      est    se     lci    uci      pval    
         ---------------------------------------------------------
+           var_visit_1    0.153  0.052  0.234  0.131  < 0.000001 
            trt_visit_1    0.078  0.024  0.058  0.12    0.07057   
          lsm_ref_visit_1  0.031  0.005  0.031  0.034  < 0.000001 
          lsm_alt_visit_1  0.071  0.025  0.081  0.101  < 0.000001 
+           var_visit_2    0.333  0.076  0.329  0.406  < 0.000001 
            trt_visit_2    0.104  0.025  0.14   0.088   0.00057   
          lsm_ref_visit_2  0.064  0.016  0.052  0.088  < 0.000001 
          lsm_alt_visit_2  0.068  0.009  0.073  0.068  < 0.000001 
+           var_visit_3    0.668  0.21   0.515  1.006  < 0.000001 
            trt_visit_3    0.107  0.023  0.094  0.136   0.00093   
          lsm_ref_visit_3  0.096  0.021  0.097  0.115  < 0.000001 
          lsm_alt_visit_3  0.087  0.021  0.074  0.117  < 0.000001 
diff --git a/tests/testthat/_snaps/print.md b/tests/testthat/_snaps/print.md
index 96c3e676f..1a4fd1237 100644
--- a/tests/testthat/_snaps/print.md
+++ b/tests/testthat/_snaps/print.md
@@ -226,12 +226,15 @@
       Analysis Function: ancova
       Delta Applied: FALSE
       Analysis Estimates:
+          var_visit_1
           trt_visit_1
           lsm_ref_visit_1
           lsm_alt_visit_1
+          var_visit_2
           trt_visit_2
           lsm_ref_visit_2
           lsm_alt_visit_2
+          var_visit_3
           trt_visit_3
           lsm_ref_visit_3
           lsm_alt_visit_3
@@ -297,9 +300,11 @@
       Analysis Function: rbmi::ancova
       Delta Applied: TRUE
       Analysis Estimates:
+          var_visit_1
           trt_visit_1
           lsm_ref_visit_1
           lsm_alt_visit_1
+          var_visit_3
           trt_visit_3
           lsm_ref_visit_3
           lsm_alt_visit_3
@@ -357,12 +362,15 @@
       Analysis Function: ancova
       Delta Applied: FALSE
       Analysis Estimates:
+          var_visit_1
           trt_visit_1
           lsm_ref_visit_1
           lsm_alt_visit_1
+          var_visit_2
           trt_visit_2
           lsm_ref_visit_2
           lsm_alt_visit_2
+          var_visit_3
           trt_visit_3
           lsm_ref_visit_3
           lsm_alt_visit_3
@@ -419,12 +427,15 @@
       Analysis Function: ancova
       Delta Applied: FALSE
       Analysis Estimates:
+          var_visit_1
           trt_visit_1
           lsm_ref_visit_1
           lsm_alt_visit_1
+          var_visit_2
           trt_visit_2
           lsm_ref_visit_2
           lsm_alt_visit_2
+          var_visit_3
           trt_visit_3
           lsm_ref_visit_3
           lsm_alt_visit_3
diff --git a/tests/testthat/test-ancova.R b/tests/testthat/test-ancova.R
index 5226131f3..3286025e0 100644
--- a/tests/testthat/test-ancova.R
+++ b/tests/testthat/test-ancova.R
@@ -192,9 +192,11 @@ test_that("ancova", {
     expect_equal(
         names(result_actual),
         c(
+            "var_visit 1",
             "trt_visit 1",
             "lsm_ref_visit 1",
             "lsm_alt_visit 1",
+            "var_visit 2",
             "trt_visit 2",
             "lsm_ref_visit 2",
             "lsm_alt_visit 2"