Re: [amelia] Error in tscsPlot due to floating point representation

On Apr 14, 2015, at 15:59, David K. &lt;code(a)rotinaj.org&gt; wrote: 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 <Amelia_tscsPlot_fpbug.patch>-- Amelia mailing list served by HUIT [Un]Subscribe/View Archive: http://lists.gking.harvard.edu/?info=amelia More info about Amelia: http://gking.harvard.edu/amelia Amelia mailing list Amelia(a)lists.gking.harvard.edu To unsubscribe from this list or get other information: https://lists.gking.harvard.edu/mailman/listinfo/amelia