code: def Y(self,value):
          
          self._N,self._P = value.shape
          self._Y = value
          self.clear_cache('Ystar1','Ystar','Yhat','LRLdiag_Yhat',
                           'beta_grad','Xstar_beta_grad','Zstar','DLZ')
docs: set phenotype