tests/bounds.R

In IQSS/Amelia: A Program for Missing Data

library(Amelia)

data(freetrade)

bds <- matrix(c(3, 30, 32), nrow = 1, ncol = 3)

set.seed(12345)
a.out.bds <- amelia(freetrade, ts = "year", cs = "country", bounds = bds,
                    max.resample = 10, p2s = 0)

out <- range(a.out.bds$imputations$imp1[is.na(freetrade[,3]),3])



if (out[1] < 30) {
  stop("lower bounds not working")
}

if (out[2] > 32) {
  stop("upper bounds not working")
}