Abortion-related deaths with no abortion indicator in the data

[jarathomas]

In probbase

Pr(abortion-related death | i310o) = 0.00001

note: i310o is a parent item to the abortion indicators i334o and i335o

Potential fix:

Pr(abortion-related deaths | i310o = 1) = 0

where i310o is

"Please confirm: When she died, she was NEITHER pregnant NOR had recently been pregnant NOR had recently delivered when she died - is that right?"

R code would be something like

sci.new <- probbaseV5
sci.new[1, 3] <- "Customized with abortion fix"
# update sci to make P(b0901 | i310) = 0
#                                 P(b0902 | i310) = 0
sci.new[which(sci.new[,1] == "i310o"), "b_0901"] <- "N"
sci.new[which(sci.new[,1] == "i310o"), "b_0902"] <- "N"
fit <- codeVA(..., sci = sci.new)