numpyro.render_model render unexpected graph with handler.do

Similar to #1710.
The reproduced code is as follows.

import numpyro
import numpyro.distributions as dist

from numpyro.handlers import do

def model(y=None):
    alpha = numpyro.sample("alpha", dist.Normal(0., 1.))
    beta = numpyro.sample("beta", dist.Normal(alpha, 1.))
    gamma = numpyro.sample("gamma", dist.HalfNormal(beta))

    return gamma
    
intervention = {
    "beta": 0.,
}
model_intervention = do(model, intervention)

numpyro.render_model(model_intervention)

This code render the graph below.

But expected is the edge between alpha and beta is cut and the edge between beta and gamma exists.

Start by running the reproduced Python snippet with numpyro.render_model and handlers.do, then compare its graph with the expected image and the related behavior in #1710. Done means the intervention graph cuts the alpha-to-beta edge while retaining the beta-to-gamma edge.