diff --git a/data-raw/US/create_MortalityTables_US_Annuities.R b/data-raw/US/create_MortalityTables_US_Annuities.R
index 96a1632b7df9304b596b66da9e4503671d1e2e0f..f5ffaef5bdedcc59e0440a2020c8702e23627a3c 100644
--- a/data-raw/US/create_MortalityTables_US_Annuities.R
+++ b/data-raw/US/create_MortalityTables_US_Annuities.R
@@ -723,6 +723,7 @@ USA2012IAR = array(
       ages = USA2012IAM.data$age,
       deathProbs = USA2012IAM.data$qx,
       improvement = USA2012IAM.data$G2x,
+      modification = function(q) { round(q, 6) },
       data = list(
         dim = list(table = "2012 IAR", sex = "m", collar = "Mortality", type = "Rententafel", data = "loaded", year = 2012, country = "USA")
       )
@@ -733,6 +734,7 @@ USA2012IAR = array(
       ages = USA2012IAM.data$age,
       deathProbs = USA2012IAM.data$qy,
       improvement = USA2012IAM.data$G2y,
+      modification = function(q) { round(q, 6) },
       data = list(
         dim = list(table = "2012 IAR", sex = "w", collar = "Mortality", type = "Rententafel", data = "loaded", year = 2012, country = "USA")
       )
@@ -742,8 +744,6 @@ USA2012IAR = array(
   dimnames = list(Sex = c("m", "f"))
 )
 
-
-
 save(USA2012IAM.basic, file = here::here("data", "USA2012IAM.basic.RData"))
 save(USA2012IAM, file = here::here("data", "USA2012IAM.RData"))
 save(USA2012IAR, file = here::here("data", "USA2012IAR.RData"))
diff --git a/data/USA2012IAR.RData b/data/USA2012IAR.RData
index 9290a8c602373abc348f81e5a519255025273db5..ccf41b196a4d7cec6e7017d7a7e465f17d62d929 100644
Binary files a/data/USA2012IAR.RData and b/data/USA2012IAR.RData differ
diff --git a/vignettes/using-the-mortalityTables-package.Rmd b/vignettes/using-the-mortalityTables-package.Rmd
index 4b17a1d7d01dd2ea000a7e6f2330ebb9a45e45b6..e50bee65a6e6c44be1cb93f0ed6cbc473c90b9b0 100644
--- a/vignettes/using-the-mortalityTables-package.Rmd
+++ b/vignettes/using-the-mortalityTables-package.Rmd
@@ -463,6 +463,24 @@ AVOe2005R.female.mod@name = "Modified table (lower bound of 3%)"
 plot(AVOe2005R.female, AVOe2005R.female.mod, title = "Original and modified table")
 ```
 
+## Rounding death probabilities
+
+Some tables like the US 2021 IAR table with projection scale G2 prescribe a certain
+rounding of death probabilities. In the 2012 IAR case, the $q_x$ are to be rounded 
+to six decimal digits after applying the projection scale to the basic 2012 mortalities.
+
+The `modification` slot of the `mortalityTable` class can be used for this case, too:
+
+```{r}
+mortalityTable.improvementFactors(
+  name = "USA 2012 IAR, male",
+  ages = USA2012IAR$m@ages,
+  deathProbs = USA2012IAR$m@deathProbs,
+  improvement = USA2012IAR$m@improvement,
+  modification = function(q) { round(q, 6) }
+)
+```
+
 # Creating mortality tables from data and modifying them using various helper functions
 
 The package MortalityTables not only provides the data structures and some