Should change the probability in rand().   %let NumGuesses = 10; data method_3; set death; z=count/5;   L=int(z); if mod(count, 5) ne 0 then do; do i = 1 to &NumGuesses; if mod(count, 5)=1 then x = L + RAND('BERNoulli',0.2); if mod(count, 5)=2 then x = L + RAND('BERNoulli',0.4); if mod(count, 5)=3 then x = L + RAND('BERNoulli',0.6); if mod(count, 5)=4 then x = L + RAND('BERNoulli',0.8); y=5*x; output; end; end; if mod(count, 5)= 0 then do; y=count; output; end; run;   ... View more

Thanks Rick for your hint. Here is my code. %let NumGuesses = 10; data method_3; set death; z=count/5; L=int(z); if mod(count, 5) ne 0 then do; do i = 1 to &NumGuesses; if mod(count, 5)=1 then x = L + RAND('BERNoulli',0.8); if mod(count, 5)=2 then x = L + RAND('BERNoulli',0.6); if mod(count, 5)=3 then x = L + RAND('BERNoulli',0.4); if mod(count, 5)=4 then x = L + RAND('BERNoulli',0.2); y=5*x; output; end; end; if mod(count, 5)= 0 then do; y=count; output; end; run; Any suggestions? Thank you very much! ... View more

It works. Not sure if that's what I exactly want. Thanks a lot! ... View more

Thank you very much Rick! The probability proportion is what I need. I twisted your code posted at Predict a random integer as the below. %let NumGuesses = 10; data Guess(keep = x PredRound diffRound PredRand diffRand Frac count k Frac); set table; call streaminit(12345); array prob[4] _temporary_ (0.2, 0.4, 0.6, 0.8); /* P(X=i) */ k = int(count/5)*5; Frac = (count - k)/10; /* distance from E(X) to x. divided by 10 to fit the Bernoulli */ do i = 1 to &NumGuesses; PredRound = round(count); /* guess the nearest integer */ PredRand = k + rand("Bern", Frac)*5; /* random guesses between k and k+1, weighted by Frac */ /* The guesses are made. Now generate a new instance of X and compute residual difference */ x = rand("Table", of prob[*]); diffRound = x - PredRound; /* rounding estimate */ diffRand = x - PredRand; /* unbiased estimate */ output; end; ; run;   Seem the PredRand is generated from the Bernoulli distribution; the x is generated from random probability. So I should use x as the random rounding number. Am I right?   I am not sure how to revise x = rand("Table", of prob[*]) for my case. Can you proviea bit more instruction? Thank you so much!!!!     ... View more