/*
**    count7.e
**
**    Poisson regression
**
**    data in a matrix in memory.
*/

library count;
#include count.ext;
countset;

rndseed 131;
n = 100;
l = exp(rndn(n,1));
data = 1+rndp(n,1,l)~rndu(n,3);

output file = count7.out reset;

__title = " EXAMPLE WITHOUT TRUNCATION ";
call countprt(poisson(data,1,2|3|4));

_cn_ZeroTruncate = 0;
__title = " EXAMPLE WITH TRUNCATION-AT-ZERO ";
call countprt(poisson(data,1,2|3|4));


output off;