# LocalReadoutMitigator

*class *`qiskit.result.LocalReadoutMitigator(assignment_matrices=None, qubits=None, backend=None)`

Bases: `BaseReadoutMitigator`

1-qubit tensor product readout error mitigator.

Mitigates `expectation_value()`

and `quasi_probabilities()`

. The mitigator should either be calibrated using qiskit experiments, or calculated directly from the backend properties. This mitigation method should be used in case the readout errors of the qubits are assumed to be uncorrelated. For *N* qubits there are *N* mitigation matrices, each of size $2 x 2$ and the mitigation complexity is $O(2^N)$, so it is more efficient than the `CorrelatedReadoutMitigator`

class.

Initialize a LocalReadoutMitigator

**Parameters**

**assignment_matrices**(*List**[**ndarray**] | None*) – Optional, list of single-qubit readout error assignment matrices.**qubits**(*Iterable**[**int**] | None*) – Optional, the measured physical qubits for mitigation.**backend**– Optional, backend name.

**Raises**

**QiskitError** – matrices sizes do not agree with number of qubits

## Attributes

### qubits

The device qubits for this mitigator

### settings

Return settings.

## Methods

### assignment_matrix

`assignment_matrix(qubits=None)`

Return the measurement assignment matrix for specified qubits.

The assignment matrix is the stochastic matrix $A$ which assigns a noisy measurement probability distribution to an ideal input measurement distribution: $P(i|j) = \langle i | A | j \rangle$.

**Parameters**

**qubits** (*List**[**int**] | None*) – Optional, qubits being measured for operator expval.

**Returns**

the assignment matrix A.

**Return type**

np.ndarray

### expectation_value

`expectation_value(data, diagonal=None, qubits=None, clbits=None, shots=None)`

Compute the mitigated expectation value of a diagonal observable.

This computes the mitigated estimator of $\langle O \rangle = \text{Tr}[\rho. O]$ of a diagonal observable $O = \sum_{x\in\{0, 1\}^n} O(x)|x\rangle\!\langle x|$.

**Parameters**

**data**(*Counts*) – Counts object**diagonal**(*Callable**|**dict**|**str**|**ndarray**| None*) – Optional, the vector of diagonal values for summing the expectation value. If`None`

the default value is $[1, -1]^\otimes n$.**qubits**(*Iterable**[**int**] | None*) – Optional, the measured physical qubits the count bitstrings correspond to. If None qubits are assumed to be $[0, ..., n-1]$.**clbits**(*List**[**int**] | None*) – Optional, if not None marginalize counts to the specified bits.**shots**(*int**| None*) – the number of shots.

**Returns**

the expectation value and an upper bound of the standard deviation.

**Return type**

**Additional Information:**

The diagonal observable $O$ is input using the `diagonal`

kwarg as a list or Numpy array $[O(0), ..., O(2^n -1)]$. If no diagonal is specified the diagonal of the Pauli operator :math`O = mbox{diag}(Z^{otimes n}) = [1, -1]^{otimes n}` is used. The `clbits`

kwarg is used to marginalize the input counts dictionary over the specified bit-values, and the `qubits`

kwarg is used to specify which physical qubits these bit-values correspond to as `circuit.measure(qubits, clbits)`

.

### mitigation_matrix

`mitigation_matrix(qubits=None)`

Return the measurement mitigation matrix for the specified qubits.

The mitigation matrix $A^{-1}$ is defined as the inverse of the `assignment_matrix()`

$A$.

**Parameters**

**qubits** (*List**[**int**] |* *int* *| None*) – Optional, qubits being measured for operator expval. if a single int is given, it is assumed to be the index of the qubit in self._qubits

**Returns**

the measurement error mitigation matrix $A^{-1}$.

**Return type**

np.ndarray

### quasi_probabilities

`quasi_probabilities(data, qubits=None, clbits=None, shots=None)`

Compute mitigated quasi probabilities value.

**Parameters**

**data**(*Counts*) – counts object**qubits**(*List**[**int**] | None*) – qubits the count bitstrings correspond to.**clbits**(*List**[**int**] | None*) – Optional, marginalize counts to just these bits.**shots**(*int**| None*) – Optional, the total number of shots, if None shots will be calculated as the sum of all counts.

**Returns**

**A dictionary containing pairs of [output, mean] where “output”**

is the key in the dictionaries, which is the length-N bitstring of a measured standard basis state, and “mean” is the mean of non-zero quasi-probability estimates.

**Return type**

**Raises**

**QiskitError** – if qubit and clbit kwargs are not valid.

### stddev_upper_bound

`stddev_upper_bound(shots, qubits=None)`

Return an upper bound on standard deviation of expval estimator.

**Parameters**

**shots**(*int*) – Number of shots used for expectation value measurement.**qubits**(*List**[**int**] | None*) – qubits being measured for operator expval.

**Returns**

the standard deviation upper bound.

**Return type**