code: def logpdf(self, value):
      res = np.sum(
          scipy.special.xlogy(self.alphas - 1, value)
          - scipy.special.gammaln(self.alphas),
          axis=-1,
      ) + scipy.special.gammaln(np.sum(self.alphas, axis=-1))
      return res
  
imports: import numpy
  import pytensor
  import pytensor.tensor
  import pytest
  import scipy
  import scipy.special
  from pytensor.graph.basic import equal_computations
  from pymc.distributions.continuous import Cauchy, ChiSquared
  from pymc.logprob.basic import conditional_logp, icdf, logcdf, logp
  from pymc.logprob.transforms import ArccoshTransform, ArcsinhTransform, ArctanhTransform, ChainedTransform, CoshTransform, ErfcTransform, ErfcxTransform, ErfTransform, ExpTransform, LocTransform, LogTransform, ScaleTransform, SinhTransform, TanhTransform, Transform
  from pymc.logprob.utils import ParameterValueError
  from pymc.testing import Rplusbig, Vector, assert_no_rvs
  from tests.distributions.test_transform import check_jacobian_det
label_desc: Use image node_id 3 for calling the DirichletScipyDist obj's underlying member method code with example usage: obj.logpdf(value) and returns: res