code: def residuals(self,Y):
                   
  
          return (Y-np.dot(self._create_B(Y),self._create_Z(Y)))
docs: Creates the model residuals
  
          Parameters
          ----------
          Y : np.array
              The dependent variables Y
  
          Returns
          ----------
          The model residuals