[R] Error message using mi() in mi package

[Andrew Miles]

On Jul 6, 2010, at 2:15 PM, Erik Iverson wrote:

>
>
>>>
>>> This looks suspiciously like a syntax problem.
>>> I would get my text editor to search for 'c14ordered'
>>> in the code. You might have missed some punctuation.
>>>
>>> -Peter Ehlers
>> A good thought.  I checked my own code (the stuff coding the data  
>> and using the mi package) and, for good measure, the code of the  
>> mi.polr function in the mi package, but the string "c14ordered"  
>> never appears.
>
> And what does
>
> >traceback()
>
> tell you?


Here is the original error message:

Error while imputing variable: a8 , model: mi.polr
  Error in eval(expr, envir, enclos) : could not find function  
"c14ordered"

I confess I'm not good at reading traceback() output, but here is what  
it says:

18: eval(expr, envir, enclos)
17: eval(predvars, data, env)
16: model.frame.default(formula = form, data = data,  
drop.unused.levels = FALSE)
15: model.frame(formula = form, data = data, drop.unused.levels = FALSE)
14: model.frame(formula = form, data = data, drop.unused.levels = FALSE)
13: eval(expr, envir, enclos)
12: eval(expr, p)
11: eval.parent(m)
10: bayespolr(formula = form, data = data, start = 0, method =  
c("logistic"),
         drop.unused.levels = FALSE, n.iter = n.iter)
9: mi.polr(formula = "a8 ~ rural_now + a3 + a7 . . . [more arguments  
here including the rest of the variables in the data set and a list of  
all the actual data], start = NULL,
        n.iter = 100, draw.from.beta = TRUE)
8: do.call(model.type, args = c(list(formula = info[[CurrentVar]] 
$imp.formula,
        data = dat, start = if (!is.null(start.val[[i]][[jj]])) {
            start.val[[i]][[jj]]
        } else {
            NULL
        }), info[[CurrentVar]]$params))
7: eval(expr, envir, enclos)
6: eval(substitute(expr), data, enclos = parent.frame())
5: with.default(data = dat, do.call(model.type, args = c(list(formula  
= info[[CurrentVar]]$imp.formula,
        data = dat, start = if (!is.null(start.val[[i]][[jj]])) {
            start.val[[i]][[jj]]
        } else {
            NULL
        }), info[[CurrentVar]]$params)))
4: with(data = dat, do.call(model.type, args = c(list(formula =  
info[[CurrentVar]]$imp.formula,
        data = dat, start = if (!is.null(start.val[[i]][[jj]])) {
            start.val[[i]][[jj]]
        } else {
            NULL
        }), info[[CurrentVar]]$params)))
3: .local(object, ...)
2: mi(imp.data, info = info2, n.iter = 10)
1: mi(imp.data, info = info2, n.iter = 10)