Getting the current sample variable value so it can be used to compute the observed inputs in a secondary model likelihood

    # Integrate the line area
    line_area_obs = pm.Deterministic('linearea',
                                     tt.sum(emission_obs[mask_loc]) / 
        (obs_wave[mask_loc[-1]] - obs_wave[mask_loc[0]]) / 1e15)

The dividing by 1e15 part.

1 Like