I'm given these two equations:
Set the scalar pi to a value after you invoke IML.
pi = constant('pi');
IML has build-in function, you can use it .
proc iml;
x=j(10000,1,.);
call randseed(12345678);
call randgen(x, 'CAUCHY');
create want var{x};
append;
close;
quit;
You asked, "how do I account for the [-infinity, infinity] part?" Here is a hint: What is the range of the expression
uniform*pi - pi/2 ?
Ask yourself how the TAN function behaves on this interval.
For more about the inverse CDF method of simulating data, see "The inverse CDF method for simulating from a distribution."
Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.
Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.
Learn how to run multiple linear regression models with and without interactions, presented by SAS user Alex Chaplin.
Find more tutorials on the SAS Users YouTube channel.