one more debug line

Author: jjacobson95

coderbuild/utils/calc_pdx_metrics.py

@@ -320,14 +320,19 @@ def lmm(time, volume, treatment, drug_name):
     print("[DBG][LMM] volume<=0:", (data['volume']<=0).sum(), " nonfinite log_volume:", (~np.isfinite(np.log(data['volume']))).sum())
     print("[DBG][LMM] groups:", data['model_id'].nunique())
 
+    print("[DBG][LMM] DATA:", data)
     fit = model.fit()
 
+    print("[DBG][LMM] fit:", fit.params)
+    
     # Get the coefficient for the interaction term 'time:exp_type'
     #interaction_term = 'time:exp_type'
 #    if interaction_term in fit.params:
 #    time_coef_value = fit.params['time']
     #print(fit.params)
     i_coef_value = fit.params['time:exp_type[T.'+drug_name+']']
+    
+    print("[DBG][LMM] i_coef_value:", i_coef_value)
     #i_coef_value = fit.params['time:exp_type['+drug_name+']']
    # else:
    #     coef_value = None  # Handle the case when the interaction term is not present