code: def test_join_mixed_ndim_supp():
      base1_rv = pt.random.normal(size=3, name="base1")
      base2_rv = pt.random.dirichlet(np.ones(3), name="base2")
      y_rv = pt.concatenate((base1_rv, base2_rv), axis=0)
  
      y_vv = y_rv.clone()
      with pytest.raises(
          ValueError,
          match="Joined logps have different number of dimensions",
      ):
          logp(y_rv, y_vv)
  
imports: import numpy
  import pytensor
  import pytest
  from pytensor import tensor
  from pytensor.graph import RewriteDatabaseQuery
  from pytensor.graph.rewriting.basic import in2out
  from pytensor.graph.rewriting.utils import rewrite_graph
  from pytensor.tensor.basic import Alloc
  from scipy import stats
  from pymc.logprob.basic import conditional_logp, logp
  from pymc.logprob.rewriting import logprob_rewrites_db
  from pymc.logprob.tensor import naive_bcast_rv_lift
  from pymc.testing import assert_no_rvs
label_desc: Use image node_id 9 for calling a global function with example usage: test_join_mixed_ndim_supp() without return types