/*
**    count10.e
**
**    bootstrap Poisson regression
**
**    data in a matrix in memory.
*/

library count,pgraph;
#include count.ext;
#include pgraph.ext;
countset;
graphset;

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

output file = count10.out reset;

_cn_ZeroTruncate = 0;
_cn_Inference = "BOOT";
__title = " Bootstrapped Poission Regression ";
__output = 5;

{b,cov,f} = poisson(data,1,2|3|4);

call countprt(b,cov,f);

call MAXDensity(_max_BootFname,0);

output off;