Given a matrix x below, I would like to have a resultant matrix y having 1 for each element of x<0 else 0: proc iml; x = {1 2, 1 -3, 2 4}; IF x<0 THEN y=1; ELSE y=0; print y; Resultant matrix should be: y= {0 0, 0 1, 0 0} The result however is a scalar y=0   ... View more

You are right thanks. I'll have tofigure out a way to get around this. ... View more

/*NORMINV of Random Binomial Probability of Default*/ %let N = 10; %let M= 3; proc iml; call randseed(123); x = j(&N,&M); cdf=j(&N,&M); q=j(&N,&M); p = {.100 .200 .500}; p_mat = repeat(p, &N); stdev_mat=p_mat#(1-p_mat); N_mat = repeat(&N, &N, &M); call randgen(x, "Binomial", p_mat, N_mat); x=x/&N;/*Random probability of default*/ q = quantile("Normal", x, p_mat, stdev_mat); /*Quantile of the probability from normal distribution (Distance to default)*/ mean=mean(x); max=x[<>,]; min=x[><,]; print p,mean,max,min,x,q; Quit; Hi Rick_SAS,   I have refined the code with your help and now I am passing on a matrix of probabilities to the function using the formula you suggested. However, I am getting an error on Quantile function again i.e. ERROR: (execution) Invalid argument to function.  Just to give you some background: I am trying to generate random variables using the mean probability of default ( matrix x) and then want to calculate distance to default by finding the quantile ( matrix q) assuming normal distribution. Thanks for your help in advance.   ... View more

Thanks a million. It worked.     ... View more

Thanks for your advice. I realize that IML is different but quite convenient too if one is good at it. ... View more

Thanks. I realized the mistake 🙂 ... View more

Thanks for your detailed explanation on parameters as vectors/matrices.   Although the RANDGEN function should work on a vector of parameters, it doesn't give me the intended results. Even after generating a full matrix of parameter p in the code below by repeating the row vector (effectively a N*M matrix), the results clearly suggest that p is being taken as a scalar i.e. only the first element of p=0.1 from the matrix is used. Do you have any suggestions? %let N = 100; %let M= 3; proc iml; call randseed(123); x = j(&N,&M); y = j(&N,&M); p= j(1, &M);/*specify row vector*/ p= {.1 .2 .5}; /*Probability*/ p_mat = repeat(p, &N); /*Repeat probability parameter for j(&N,&M) matrix*/ call randgen(x, "Binomial",p_mat,&N ); x=x/&N; mean=mean(x); max=max(x); print p,mean,max; Quit;   RESULT:   p 0.1 0.2 0.5 mean 0.1014 0.1001 0.1075 max 0.2 ... View more

Thank you all for your quick responses. I didn't know we had such a cooperative SAS  community.   Creating matrices of parameters (mu_mat and stdev_mat) is quite intuitive and seems to be working fine. However, it seems that I am having some problem with RANDGEN too. Although I am specifying the parameters as vectors they seem to be taken as scalars. Specifically, RANDGEN seems to be taking only the first parameter (.01) in the vector below and applies it to whole matrix.   mu= j(1, &M); mu= {.01 .1 .6};   It would also help if you could elaborate a bit more on what do you mean by 'row vector (apply to cols), column vector (apply to rows)'. I was goin gthrough the doc and it says: "if the parameters contain m elements, the jth column of the result matrix consists of random values drawn from the distribution with parameters param1[j], param2[j], and param3[j]".   Thanks in advance.     ... View more