def model(data, tau_like, mu_base, tau_base, alpha, K = 20):
    
    mdl_infty_weights = dp.sticks(alpha, K)
    
    mdl_gaussian = gaussian_mixture(data, mdl_infty_weights['pi'], 
                     tau_like, mu_base, tau_base, K)
    
    return [mdl_infty_weights, mdl_gaussian]
def model(data, tau_like, mu_base, tau_base, alpha, K):    
    mdl_fixed_weights = fixed_mixture_weights(alpha, K)    
    mdl_gaussian = gaussian_mixture(data, mdl_fixed_weights['pi'], 
                     tau_like, mu_base, tau_base, K)
    
    return [mdl_fixed_weights, mdl_gaussian]