Dear All,
I have looked at this problem with *Zelig* previously, but still cannot
find how to combine the five datasets produced by *AmeliaView*. In
our data it is critical to set bounds for factors, otherwise the imputed
data produces out of range values for standard blood tests. How exactly do
you create the *output object* that is needed to combine the five
datasets that are produced by AmeliaView?
Although page 27 of the *Amelia* user's manual gives the instructions
for setting bounds without using the *AmeliaView* GUI, the syntax seems
cumbersome "* bds <- matrix(c(3, 30, 40), nrow = 1, ncol = 3) "* if this
must be done for each variable. The alternative presented, to allow
unbounded imputation that might produce negative values, that would
then be set to "0" would produce unacceptable results.... such patients
would be dead.
I very much appreciate your efforts with this program. I and my
colleagues at the Wistar Institute of the University of Pennsylvania have
successfully used the older version of Amelia and Zelig for a
manuscript on pediatric treatment response for HIV that has been accepted
pending revisions by PLoS- One, so I hope to be able to continue to use
the program. I even tried to load the older, more user friendly version
of *Zelig* , but it would not work with the latest version of *R*. Any
additional guidance as to how the bounded datasets produced by *AmeliaView*
can be turned into an output object and then combined would be much
appreciated. Perhaps a new updated version of this excellent package
(named "*Hoffa*" ?) might simplify this process.
Skip Barbour
>
--
*There is nothing so fatal to character as half finished tasks.*
David Lloyd George