/*
**   count4.e
**
**   negative binomial regression
**
*/

library count;
#include count.ext;
countset;

dep = { wars };
ind = { age, party, unem };
dataset = "count";

output file = count4.out reset;

{ b,vc,ll } = negbin(dataset,dep,ind,0);

call countCLPrt(b,MAXTlimits(b,vc),ll);

output off;