14 Apr
2015
14 Apr
'15
9:59 a.m.
Dear List Members,
I experience an error when running tscsPlot on an Amelia output with the message:
Error in output$theta[, , ceiling(i/drawsperimp)] :
subscript out of bounds
The error is due to the ceiling function, which for some settings returns wrong results
due to floating point arithmetic. In my case m = 6 and draws = 100 resulting in
ceiling(100 / ( (1/6) * 100)) = 7 != 6
so index is out of bounds.
I attached a patch which only uses integer arithmetic and is thus not susceptible to
floating point arithmetic errors.
Hope that helps!
Kind regards,
David