topic Simulating Poisson data in Statistical Procedures

Simulating Poisson data

CHELS — Fri, 17 Jan 2025 19:32:32 GMT

Hi everybody,

I am trying to simulate Poisson data but not sure if I am doing it correctly.

I know that in a given population, an event occurs at a rate of 1.2 per year.

If I would like to simulate 1000 set of data each containing 50 patients, is this how the code should be written?

Thanks!

data Poisson;
call streaminit(4321);
lambda = 0.8;

do SimID = 1 to 1000;
do i = 1 to 50;
x = rand("Poisson", lambda);
output;
end;
end;
run;

Re: Simulating Poisson data

quickbluefish — Fri, 17 Jan 2025 21:15:11 GMT

I think you want to set lambda to 1.2 in your case:

data test;
do i=1 to 10000;
  x=rand('poisson',1.2);
  output;
end;
run;

proc univariate data=test noprint;
histogram x / midpoints=0 to 7 by 1;
inset mean median max / position=NE;
run;

Re: Simulating Poisson data

Ksharp — Sat, 18 Jan 2025 08:31:51 GMT

https://blogs.sas.com/content/iml/2024/08/07/simulate-poisson-regression.html