I present you with a function that solves a problem that has bugged me for
many years. I think the problem may be general enough to at least consider
adding this function, or a revamped version of it, to the 'stats' package,
with the other multinomial functions reside.
I'm using R to export
Ugh, I made a typo at the very heart of my message:
when I preprocess each line in R as p-a/sum(a), occasionally a line will
sum to 0.999, 1.002, or the like
should be
when I preprocess each line in R as p-round(a/sum(a),3) occasionally a
line will sum to 0.999, 1.002, or the like
Also,
