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(a)fas.harvard.edu
<http://1429236.signature1.mailinfo.com/confirm2.6/000C0108/010B024D/000D014
D/0111831384.jpg>
Show replies by date