code: def markov_blanket(y, mean, scale, shape, skewness):
          
          return ss.t.logpdf(x=y, df=shape, loc=mean, scale=scale)
docs: Markov blanket for the Poisson distribution
  
          Parameters
          ----------
          y : np.ndarray
              univariate time series
  
          mean : np.ndarray
              array of location parameters for the t distribution
  
          scale : float
              scale parameter for the t distribution
  
          shape : float
              tail thickness parameter for the t distribution
  
          skewness : float
              skewness parameter for the t distribution
  
          Returns
          ----------
          - Markov blanket of the t family