def initial_mean_index_tuples(n_mixtures, factors):
    """Index tuples for initial_mean.

    Args:
        n_mixtures (int): Number of elements in the mixture distribution of the factors.
        factors (list): The latent factors of the model

    Returns:
        ind_tups (list)

    """
    ind_tups = []
    for emf in range(n_mixtures):
        for factor in factors:
            ind_tups.append(("initial_states", 0, f"mixture_{emf}", factor))
    return ind_tups