The readout.py module allows you to estimate the measurement confusion matrix.
 get_flipped_program(program) For symmetrization, generate a program where X gates are added before measurement. estimate_confusion_matrix(qc, qubit, num_shots) Estimate the readout confusion matrix for a given qubit. estimate_joint_confusion_in_set(qc, qubits, …) Measures the joint readout confusion matrix for all groups of size group_size among the qubits. estimate_joint_reset_confusion(qc, qubits, …) Measures a reset ‘confusion matrix’ for all groups of size joint_group_size among the qubits. marginalize_confusion_matrix(…) Marginalize a confusion matrix to get a confusion matrix on only the marginal_subset.