Hi Matt---

 

I uncovered a couple small glitches in AMELIA…

 

(1)  If I identify a nominal variable with exactly 10 unique values, AMELIA gives the following message:

 

"Warning in amcheck(x = x, m = m, idvars = numopts$idvars, priors = priors, :

 The number of catagories in one of the variables marked nominal has greater than 10 categories. Check nominal specification"

 

I suspect that AMELIA is counting "NA" as one of the unique values.

 

This can be corrected by rewriting the statement that counts the number of unique values, as follows:

 

     categories <- length( unique(IVs) [ !is.na (unique (IVs) ) ] )

 

I'm not sure whether AMELIA only checks the number of unique values for the variables passed in the "bounds" input file, or if it also checks the number of unique values passed in the "noms" and "ords" input files… If it checks all of these, then there may be more than one statement that needs to be modified.

 

(2) Also, the error message cited above has a misspelled word ("catagories")…

 

 

Wayne A. Thornton

 

Harvard University

 

thornton@fas.harvard.edu