jointfit function

Fitting Multivariate Double Hierarchical Generalized Linear Models using h-likelihood Approach