hi all,
i have a series of probability vector that i'd like to feed into
multinomial to get an array of vector outcomes back. for example,
given:
p = array([[ 0.9 , 0.05, 0.05],
[ 0.05, 0.05, 0.9 ]])
i'd like to call multinomial like this:
multinomial(1, p)
to get a vector of multinomial samplers, each using the nth list in
'p'. something like:
array([[1, 0, 0], [0, 0 1]]) in this case. is this possible? it seems
like 'multinomial' takes only a one dimensional array. i could write
this as a "for" loop of course but i prefer a vectorized version since
speed is crucial for me here. the 'multinomial' function from scipy
does not support this, it seems. is there a way to do this in pygsl?
thanks very much.
------------------------------------------------------------------------------
Come build with us! The BlackBerry(R) Developer Conference in SF, CA
is the only developer event you need to attend this year. Jumpstart your
developing skills, take BlackBerry mobile applications to market and stay
ahead of the curve. Join us from November 9 - 12, 2009. Register now!
http://p.sf.net/sfu/devconference
_______________________________________________
pygsl-discuss mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/pygsl-discuss
