Conditional Random number generation

Question

wang syr 2017-3-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/327656-conditional-random-number-generation

评论： Bruno Luong 2023-8-12

采纳的回答： Roger Stafford

Hello there, For example; If I want to generate 5 random integer numbers with a sum of 20, how can I do that?

" ... example = ceil(10*rand(100, 5)) ... "

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Rik 2020-12-20

Why did you edit away your question? It is stored on Google cache anyway, so it's easy to restore.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Roger Stafford 2017-3-4

4
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/327656-conditional-random-number-generation#answer_257296

在 MATLAB Online 中打开

   function R = randfixedsumint(m,n,S);
   % This generates an m by n array R.  Each row will sum to S, and
   % all elements are all non-negative integers.  The probabilities
   % of each possible set of row elements are all equal.
   % RAS - Mar. 4, 2017
   if ceil(m)~=m|ceil(n)~=n|ceil(S)~=S|m<1|n<1|S<0
    error('Improper arguments')
   else
    P = ones(S+1,n);
    for in = n-1:-1:1
     P(:,in) = cumsum(P(:,in+1));
    end
    R = zeros(m,n);
    for im = 1:m
     s = S;
     for in = 1:n
      R(im,in) = sum(P(s+1,in)*rand<=P(1:s,in));
      s = s-R(im,in);
     end
    end
   end
   return

6 个评论
显示 4更早的评论隐藏 4更早的评论

Yu Takahashi 2021-2-9

编辑：Walter Roberson 2021-2-10

Wondering whether it is possible to specify the max and min of the devided value? i.e., something like what you kindly provided in the randfixedsum function, thanks!

Ref

https://viewer.mathworks.com/?viewer=plain_code&url=https%3A%2F%2Fwww.mathworks.com%2Fmatlabcentral%2Fmlc-downloads%2Fdownloads%2Fsubmissions%2F9700%2Fversions%2F1%2Fcontents%2Frandfixedsum.m&embed=web

Bruno Luong 2023-8-12

@Walter Roberson "Wondering whether it is possible to specify the max and min of the devided value?"

https://www.mathworks.com/matlabcentral/fileexchange/133762-ricb

请先登录，再进行评论。

Answer 2

Walter Roberson 2017-3-2

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/327656-conditional-random-number-generation#answer_256941

See https://www.mathworks.com/matlabcentral/fileexchange/9700-random-vectors-with-fixed-sum

9 个评论
显示 7更早的评论隐藏 7更早的评论

Walter Roberson 2017-3-2

在 MATLAB Online 中打开

I just thought of a completely different method for the case where you have a fixed sum of a fixed number of positive integers:

random_posint_partition = @(of_what, num_rows, number_of_partitions) accumarray( [repmat((1:num_rows).',of_what, 1), randi(number_of_partitions, num_rows*of_what, 1)], 1);

Example run:

>> random_posint_partition(20,7,5)
ans =
     4     3     4     5     4
     5     3     3     5     4
     3     6     1     6     4
     3     6     4     5     2
     5     3     6     2     4
     4     6     8     1     1
     5     1     6     4     4
>> sum(ans,2)
ans =
    20
    20
    20
    20
    20
    20
    20

That is, the first parameter gives the number which is to be partitioned, the second parameter gives the number of different times you want the generation to be done, and the third parameter gives the number of pieces to partition into.

... I thought it would be a lot harder to come up with a fair method!

Oops, I just discovered that this can generate entries with a count of 0. :( For example,

6 5 5 4 0

The adjusted version is

random_posint_partition = @(of_what, num_rows, number_of_partitions) 1 + accumarray( [repmat((1:num_rows).',of_what-number_of_partitions, 1), randi(number_of_partitions, num_rows*(of_what-number_of_partitions), 1)], 1);

Roger Stafford 2017-3-4

编辑：Roger Stafford 2017-3-4

@Walter. The distribution of values one gets with this method you have described depends very heavily upon the assumption that is made about the underlying statistics involved. By my calculations, using your method, the probability of getting the result [4,4,4,4,4] is 305,540,235,000 times as likely as that of getting the result [20,0,0,0,0]. (I’m assuming your original method, which would allow zero integers.) It is equivalent to tossing 20 pebbles into 5 cans at random and recording the sum of pebbles in each can. In other words, the distribution is heavily weighted in favor of roughly equal values among the five integers.

However, suppose you assume instead that you are in a five dimensional “hyper-cube” of integer space, each integer of which can vary from 0 to 20, and that each integer-valued point in this cube has an equal a priori probability, namely, 1/(21^5). Then suppose, as a conditional probability, that we restrict ourselves to the subset in which the five integers must have 20 as their sum. Such a conditional probability distribution will give [4,4,4,4,4] and [20,0,0,0,0] an equal conditional probability, as opposed to the above enormous difference. This latter kind of conditional probability is analogous to the kind of conditional probability assumption that is the basis of my ‘randfixedsum’ algorithm #9700 for continuous distributions. If we were to use an integer adjustment to the results of randfixedsum, this latter is roughly the kind of integer distribution we would obtain.

Walter Roberson 2017-3-4

Ah. I don't think I know how to implement your suggestion, though, at least not without generating all of the possible choices that sum to 20 and then picking one at random.

John D'Errico's https://www.mathworks.com/matlabcentral/fileexchange/12009-partitions-of-an-integer can calculate all of the possible partitions; a question is whether we can avoid having to take that step.

Walter Roberson 2017-3-4

https://en.wikipedia.org/wiki/Partition_(number_theory)#Restricted_part_size_or_number_of_parts talks about restricted partitioning briefly, and ties it to change making problems, which does indeed sound equivalent to the approach I was taking. Those are in turn tied to knapsack problems.

请先登录，再进行评论。

Answer 3

Bruno Luong 2020-8-10

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/327656-conditional-random-number-generation#answer_477376

在 MATLAB Online 中打开

m = 5;
n = 3;
s = 10;

This will generate uniform distribution with sum criteria

% generate non-negative integer random (m x n) array row-sum to s 
[~,r] = maxk(rand(m,s+n-1),n-1,2);
z = zeros(m,1);
r = diff([z, sort(r,2), (s+n)+z],1,2)-1;

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Bimal Ghimire 2020-10-4

https://www.mathworks.com/matlabcentral/answers/327656-conditional-random-number-generation#answer_257296

While generating conditional random numbers, how can we generate random numbers that has a limit of some maximum value and have certain specified sum value?

请先登录，再进行评论。