diff --git a/EPCsAndCarbon/carbonCosts.Rmd b/EPCsAndCarbon/carbonCosts.Rmd
index e4d958104d237cb9ece360d912cab5107921f41c..e2f804b7b5965816bd7b9e28c2b2978eb5982def 100644
--- a/EPCsAndCarbon/carbonCosts.Rmd
+++ b/EPCsAndCarbon/carbonCosts.Rmd
@@ -172,18 +172,20 @@ sotonEpcMSOA_DT <- sotonPostcodesDT[, .(nEPCs = .N,
                                     ] 
 ```
 
-```{r, compareEstimates}
+Join the estimates together at MSOA level for comparison. There are `r uniqueN(sotonElecDT$MSOACode)` MSOAs in Southampton.
 
+```{r, compareEpcEstimates}
+# 32 LSOAs in Soton
 # add deprivation
 setkey(sotonElecDT, MSOACode)
 setkey(deprivationDT,MSOACode)
-setkey(sotonMSOA_DT,MSOACode)
 sotonMSOA_DT <- deprivationDT[sotonElecDT]
 sotonMSOA_DT <- sotonDep_DT[sotonMSOA_DT]
 
 
 t <- sotonMSOA_DT[, .(nHouseholds2011 = sum(totalHouseholds),
                       nElecMeters2018 = sum(METERS),
+                      nEPCs = sum(),
                       sumMWh = sum(KWH)/1000)]
 
 kableExtra::kable(t, caption = "Comparison of different estimates of the number of dwellings")