I have defined set W and parameter Coef(W) as follows:
Code: Select all
Set W /1*10/;
Parameter Coef(W);
Code: Select all
[1 1 0 1 1 0 1 1 0 1]
or
[0 1 1 1 0 0 1 1 1 1]
or
...
Kind regards,
Amirhossein
Code: Select all
Set W /1*10/;
Parameter Coef(W);
Code: Select all
[1 1 0 1 1 0 1 1 0 1]
or
[0 1 1 1 0 0 1 1 1 1]
or
...
Code: Select all
*reset the seed for the random number generator.
execseed = 1 + gmillisec(jnow);
Set W /1*10/
select(W);
alias(W,Wp);
select(W)=NO;
Parameter Coef(W)
p
A /70/
;
loop(W$(ord(w) le floor(A*card(W)/100)),
p = uniformInt(1,card(W));
while( sum(Wp$(ord(Wp) eq p),1) eq sum(Wp$(ord(Wp) eq p and select(Wp)),1),
p = uniformInt(1,card(W));
);
select(Wp)$(ord(Wp) eq p) = yes;
Coef(Wp)$(ord(Wp) eq p) = 1;
);
Coef(Wp)$(not select(Wp)) = 0;
display Coef;
Code: Select all
Set W /1*10/;
Parameter Coef(W);
scalar a / 0.7 /;
* set the seed in some random way, based on the time
execseed = gmillisec(jnow);
coef(w) = [uniform(0,1) < a];
display coef;
Code: Select all
set i /1*10/;
parameter c(i);
scalar k 'number of 1s in c' / 4 /;
c(i)$[ord(i) <= k] = 1;
* set the seed in some random way, based on the time
execseed = gmillisec(jnow);
option shuffle = c;
display c;