Power spectrum estimators for large CMB datasets