tests/moPrep-test.R

library(Amelia)

x <- rnorm(100)
s <- x + rnorm(100, 0, 0.1)
vv <- rep(0.1^2/var(s), 100)

df <- data.frame(x, s)
mop <- moPrep(df,s ~ s,error.proportion = vv)
a.out <- amelia(mop)
IQSS/Amelia documentation built on April 12, 2024, 11:36 a.m.