R Under development (unstable) (2025-04-26 r88181 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(msme)
Loading required package: MASS
Loading required package: lattice
> 
> # library(msme, lib.loc="lib")
> 
> data(medpar)
> 
> med.nb.g <- ml_glm3(los ~ hmo + white,
+                    family = "gNegBinomial", 
+                    link = "log",
+                    group = medpar$provnum, 
+                    data = medpar)
> 
> summary(med.nb.g)

Call:
ml_glm3(formula = los ~ hmo + white, data = medpar, family = "gNegBinomial", 
    link = "log", group = medpar$provnum)

Deviance Residuals:
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
 -51.14   51.14   51.14   37.66   51.14   51.14 

Coefficients (all in linear predictor):
            Estimate     SE      Z        p    LCL    UCL
(Intercept)   0.9293 0.0758 12.258 1.53e-34  0.781 1.0779
hmo          -0.0701 0.0515 -1.361    0.173 -0.171 0.0309
white        -0.0270 0.0653 -0.413     0.68 -0.155 0.1011

Null deviance: NA  on  1494 d.f.
Residual deviance: 3909331  on  1492 d.f.
AIC:  8992.173

Number of optimizer iterations:  96 

> 
> 
> proc.time()
   user  system elapsed 
   1.18    0.25    1.37