// Test case 1
clear all 
use accident3
mypoisson5 accidents cvalue kids traffic 

assert `"`e(cmd)'"'        == `"mypoisson5"'
assert `"`e(predict)'"'    == `"mypoisson5_p"'
assert `"`e(properties)'"' == `"b V"'

assert         e(N)    == 505
assert         e(rank) == 4

qui {
mat T_b = J(1,4,0)
mat T_b[1,1] = -.6558870690222316
mat T_b[1,2] = -1.009016972473914
mat T_b[1,3] =  .1467114650851019
mat T_b[1,4] =  .5743541247442324
}
matrix C_b = e(b)
assert mreldif( C_b , T_b ) < 1e-12
_assert_streq `"`: rowfullnames C_b'"' `"y1"'
_assert_streq `"`: colfullnames C_b'"' `"cvalue kids traffic _cons"'
mat drop C_b T_b

qui {
mat T_V = J(4,4,0)
mat T_V[1,1] =  .0049911902167341
mat T_V[1,2] =  .0002953642487161
mat T_V[1,3] = -.0000506909358346
mat T_V[1,4] = -.0089523155601508
mat T_V[2,1] =  .0002953642487161
mat T_V[2,2] =  .0065280055261688
mat T_V[2,3] =  .0002050149836939
mat T_V[2,4] = -.0054776138886792
mat T_V[3,1] = -.0000506909358346
mat T_V[3,2] =  .0002050149836939
mat T_V[3,3] =  .0009844631577381
mat T_V[3,4] = -.0075052131640854
mat T_V[4,1] = -.0089523155601508
mat T_V[4,2] = -.0054776138886792
mat T_V[4,3] = -.0075052131640854
mat T_V[4,4] =  .0806284655627814
}
matrix C_V = e(V)
assert mreldif( C_V , T_V ) < 1e-12
_assert_streq `"`: rowfullnames C_V'"' `"cvalue kids traffic _cons"'
_assert_streq `"`: colfullnames C_V'"' `"cvalue kids traffic _cons"'
mat drop C_V T_V