contingency_table = pd.crosstab(filtered_df['pms'], filtered_df['neuroticismvariable'])
chi2, p_value, dof, expected = chi2_contingency(contingency_table)
print('Likelihood ratio:', chi2)