From 988b96578348bcc6de4fa65478f9910ccad993dd Mon Sep 17 00:00:00 2001
From: Merlise Clyde <clyde@stat.duke.edu>
Date: Thu, 30 Nov 2023 21:54:10 -0500
Subject: [PATCH] initialize dispersion

---
 src/family.c  | 12 ++++++------
 src/glm_fit.c |  2 +-
 2 files changed, 7 insertions(+), 7 deletions(-)

diff --git a/src/family.c b/src/family.c
index f0500859..d48d0811 100644
--- a/src/family.c
+++ b/src/family.c
@@ -124,11 +124,11 @@ double poisson_dispersion(double *resid,  double *weights, int n, int rank) {
 /* Gamma */
 
 double gamma_loglik(double *Y, double*mu, double *wts, double devb, int n) {
-  int i;
+  int i=0;
   double ll = 0.0;
-  double disp, ntemp = 0.0;
+  double disp= 0.0, ntemp = 0.0;
+  
   
-
   for (i = 0; i < n; i++) {
     ntemp += wts[i];
   }
@@ -144,7 +144,7 @@ double gamma_loglik(double *Y, double*mu, double *wts, double devb, int n) {
 
 void gamma_variance(double *mu, double *var, int n) {
   
-  int i;
+  int i=0;
   
   for (i = 0; i<n; i++) {
     var[i] = pow (mu[i],2.0);
@@ -154,8 +154,8 @@ void gamma_variance(double *mu, double *var, int n) {
 
 void gamma_dev_resids(double *ry, double *rmu, double *rwt, double *rans, int n)
 {
-  int i;
-  double mui, yi, wti;
+  int i=0;
+  double mui=1.0, yi=0.0, wti=0.0;
   
   for (i = 0; i < n; i++) {
     mui = rmu[i];
diff --git a/src/glm_fit.c b/src/glm_fit.c
index d93ce150..d2661042 100644
--- a/src/glm_fit.c
+++ b/src/glm_fit.c
@@ -52,7 +52,7 @@ SEXP glm_bas(SEXP RX, SEXP RY, glmstptr *glmfamily, SEXP Roffset, SEXP Rweights,
 		*variance=REAL(Rvariance);
 
 	double  one = 1.0,  tol, devold, devnew;
-	double disp;
+	double disp= 1.0;
 	
 	int   i, j, l, rank=1, *pivot=INTEGER(Rpivot), conv=0;