Takes a matrix in which rows sum to one. For each row, performs a single multinomial draw from amongst the columns, weighted by their values in that column

samp_from_mat(M)

Arguments

M

a matrix whose rows are reals, each summing to one

Value

a vector length = nrow(M) of indices, with each element being the column that was chosen in that row's sampling