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"
Al
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 da
