Thank you so much! Really helpful.  ... View more

Thank you all so much for your input and time. I think I had a better idea to handle my questions.   ... View more

Hi Paige , Thank you for your explanation.  I get your meaning, but since this is my first time to do this, could you please show me some examples, on how to "think a function"--you don't need to give me any detailed description, only some study case, or links you happen to know. Thank you!    Best,  Heather ... View more

Yes, that's what I want to know, but you r right, doesn't make sense......thanks.  ... View more

Hi,  I read the blog. It sample random variable using ranfun and Bernoulli.  For my question, I have a correlation and covariance matrix calculated from a real dataset (proc corr).  and I plan to generate random variables using this matrix.  As Dr. Wicklin explain in his blog: https://blogs.sas.com/content/iml/2011/01/12/sampling-from-the-multivariate-normal-distribution.html.     After generate the standard normal variables, these are not the real data.  I am not sure whether I can use these variable fit the logistic or OSL model directly? For example binary variable, such as gender, I need to "transfer back" to 1/0,  and variable age I can't use negative value -- but how can I do that?       Thanks.     ... View more

Hi I have not read this link  yet. Let me read it first. thanks a lot!  ... View more

I see.. How about "call streaminit(321)"? Does it use seed 321 for 1000 iterations?    ... View more

Sure.  I need to run 1000 times, however with the current code, they all have the same seed = 12345.  I want to create a macro variable which change every time I pull the sample.  One idea is to set the seed macro variable as the time we create the samples.  for example for the first iteration: if time is 10am 1 minute 1 second   then %let seed = 100101; second iteration seed become   100111 if each simulation takes 10 seconds.  I remember SAS has a function of current time...system....I forget. It create a macro variable of the current time.  Thanks    %macro uniform(datain =, dataout =, size = ); /*pulling a random sample from original data*/ data &dataout. (drop=obsleft sampsize); sampsize= &size.; obsleft=totobs; do while(sampsize>0); pickit+1; if ranuni(12345)<sampsize/obsleft then do; set &datain. point=pickit nobs=totobs; output; sampsize=sampsize-1; end; obsleft=obsleft-1; end; stop; run; %mend uniform; ... View more

Chris, Thank you for your guidance. I will follow your advice for next question.  I'd solved this question this time.  ... View more