Hi,
I think you should direct this question to the Amelia email list since
it's about Amelia rather than Zelig.
Kosuke
On Tue, 27 May 2008, Donald Braman wrote:
I'm attempting to use AmeliaII in the following
way:
vars_to_impute <- gundata[,c("progun", "egalitarianism",
"individualism",
"crfear", "victim", "female", "RACE",
"income", "URBANKID", "URBANNOW",
"RELIGION", "iss", "democrat",
"conservative")]
imputed <- amelia(data=vars_to_impute, p2s=2,
noms=c("RACE","RELIGION","URBANKID","URBANNOW"),
outname="imputation",
ords=c("democrat", "conservative")
The problem that I run into is that in about half of the imputation attempts
fail due to non-invertible covariance matrices. I'm curious if there is any
way to deal with this aside from removing the most highly covariant
variables? For example, given that it produces imputed data about 1/5 of
the time, would it be acceptable to use successful imputations? E.g., can I
just set m=100 and use as many imputations as I like from the resulting set
of ~25 imputations?
If I do need to remove the covariate variables, do you know of a simple way
to check for that among a given set of variables?
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)lists.gking.harvard.edu
[un]subscribe Options:
http://lists.gking.harvard.edu/?info=zelig
Zelig program information:
http://gking.harvard.edu/zelig/