common random number when satisfying certain conditions

9 visualizaciones (últimos 30 días)
Hi all,
I'm trying to generate a matrix of random numbers that has to satisfy the restriction that when certain conditions are satisfied, the random numbers must be the same. Consider the following example with two individuals (id takes the values 1 and 2) and up to six options (defined by ch)
id=[1 1 1 1 1 1 2 2 2 2 2 2]';
ch=[1 2 3 1 2 3 3 4 5 3 5 6]';
I would like to generate random numbers such that whenever id==i and ch==j, the numbers are the same (and the same is true for all values of i and j).
I was able to do this in the following way
RN=zeros(12,5);
U=unique(id);J=unique(ch);
RNo=randn(kron(max(size(U)),max(size(J))),5); %for each combination I need 5 random numbers. Note that this generates random numbers that will never be used because there are some combinations of id and ch that do not appear in the data.
DJ=repmat(J,max(U),1); % Define a vector with all possible alternatives
DI=kron(U,ones(max(J),1)); % Define a vector with id repeated as many times as possible alternatives
for k=1:max(size(id))
ii=id(k,1);jj=ch(k,1);
for s=1:5,
RN(id==ii & ch==jj,s)=RNo(DI==ii & DJ==jj,s);
end
end
Now RN has the random numbers and these are the same whenever a combination of id and ch is repeated. However, doing this at a larger scale (my data consists of 500,000 observations that combine 900 values for id with up to 25 values for ch), with even one random number (instead of 5), takes forever.
Is there any way to generate the matrix with random numbers that does not require all this manipulation? I have tried to use indices but I haven't figured out a clever way to do it.
Thanks,
Fernando
  2 comentarios
Walter Roberson
Walter Roberson el 20 de Mzo. de 2013
50000 observations -- do you mean you want 50000 random numbers generated? 900 * 25 < 50000, so is the task to verify that the observations meet the condition if id&ch being consistent ?
Fernando
Fernando el 20 de Mzo. de 2013
Hi Walter,
I have 900 individuals and up to 25 alternatives. Individuals are observed over a lot of periods, meaning that each id may be observed with the same ch several times. This generates the 500,000 observations.
I need to generate the 900*25 random numbers and then assign these numbers accordingly to id and ch, generating a vector of 500,000x1 of repeated random numbers. The numbers will be repeated every time a combination of id and ch is repeated. My problem is how to assign the 900*25 random numbers to the data, according to the combination of id and ch.
Thanks,

Iniciar sesión para comentar.

Respuesta aceptada

Walter Roberson
Walter Roberson el 20 de Mzo. de 2013
Editada: Walter Roberson el 20 de Mzo. de 2013
[uid, uida, uidc] = unique(id);
[uch, ucha, uchc] = unique(ch);
numid = length(uid);
numch = length(uch);
randmat = rand(numid, numch);
randobs = randmat( sub2ind([numid, numch], uidc, uchc) );
Note: the above will fail if the id and ch vectors are not the same length.
  2 comentarios
Fernando
Fernando el 20 de Mzo. de 2013
Awesome, thanks!
Walter Roberson
Walter Roberson el 20 de Mzo. de 2013
Sorry, I missed accessing randmat() at those locations; I have edited that into place.

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Creating and Concatenating Matrices en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by