def filter_pheno_nan(x_species, y_pheno, pheno_name, verbose=False):
    """Filter X and Y to have same index and drop NANs
    """
    y = y_pheno[pheno_name]
    y.dropna(axis=0, how='all', inplace=True)
    X = x_species
    X = X.reindex(y.index)
    return X, y