From: Steven on
Hi All,

Can anyone tell me the right way to use binomial function in proc
model? The following did not work for me.


data s;
_name_ = "z";
z = 1; y = .5; /* first row */
output;
_name_ = "y";
z = .5; y = 1; /* Second row */
output;
run;

data all ;
y = 0 ;
z = 0 ;
output ;
run ;


/* Simulate bivariate with binormial vs poisson distributions */ proc
model data=all ;
y = 0;
errormodel y ~poisson( 20 );
z = 0;
errormodel z ~binomial(0.15 );
solve y z / random = 1000 sdata = s out =
sim_poisson_normal(where=(_rep_>0)) ; run;

Thanks!


Thanks!

Steven