def statistics(<EOL>graph, ps, spanning_cluster=True, model='<STR_LIT>', alpha=alpha_1sigma, runs=<NUM_LIT><EOL>):
my_microcanonical_averages = microcanonical_averages(<EOL>graph=graph, runs=runs, spanning_cluster=spanning_cluster, model=model,<EOL>alpha=alpha<EOL>)<EOL>my_microcanonical_averages_arrays = microcanonical_averages_arrays(<EOL>my_microcanonical_averages<EOL>)<EOL>return canonical_averages(ps, my_microcanonical_averages_arrays)<EOL>
Helper function to compute percolation statistics

See Also
--------

canonical_averages

microcanonical_averages

sample_states