Resource saving via ensemble techniques for quantum neural networks

Quantum neural networks hold significant promise for numerous applications, particularly as they can be executed on the current generation of quantum hardware. However, due to limited qubits or hardware noise, conducting large-scale experiments often requires significant resources. Moreover, the output of the model is susceptible to corruption by quantum hardware noise. To address this issue, we propose the use of ensemble techniques, which involve constructing a single machine learning model based on multiple instances of quantum neural networks. In particular, we implement bagging and AdaBoost techniques, with different data loading configurations, and evaluate their performance on both synthetic and real-world classification and regression tasks. To assess the potential performance improvement under different environments, we conducted experiments on both simulated, noiseless software and IBM superconducting-based QPUs, suggesting these techniques can mitigate the quantum hardware noise. Additionally, we quantify the amount of resources saved using these ensemble techniques. Our findings indicate that these methods enable the construction of large, powerful models even on relatively small quantum devices.


Introduction
The emerging field of quantum machine learning [1] holds promise for enhancing the accuracy and speed of machine learning algorithms by utilizing quantum computing techniques.Although the potential of quantum machine learning is expected to be advantageous for certain classes of problems in chemistry, physics, material science, and pharmacology [2], its applicability to more conventional use cases remains uncertain [3].Notably, utilizable quantum machine learning algorithms generally need to be adapted to run on 'NISQ' devices [4], that are current noisy quantum computer, no error corrected and with modest number of qubits and circuit depth capabilities.In the quantum machine learning scenario, the quantum counterparts of classical neural networks, quantum neural networks [5], have emerged as the de facto standard model for solving supervised and unsupervised learning tasks in the quantum domain.
While quantum neural networks have generated much interest, they presently have some issues.The first is barren plateau [6] characterised by the exponentially-fast decay of the loss gradient's variance with increasing system size.This problem may be exacerbated by various factors, such as having overly-expressive quantum circuits [7].To address this issue, quantum neural networks need to be carefully designed [8] and to incorporate expressibility control techniques such as projection [9] and bandwidth control [10].The second problem, which is the one addressed in this work, concerns the amount of resources required to run quantum neural networks (the limited number of total qubits -currently up to over a hundred-and the low fidelity of operations on current quantum devices severely restrict the size of the quantum neural network in terms of input dimension and layers).
In order to address the latter issue, we propose employing of NISQappropriate implementation of ensemble learning [11], a widely used technique in classical machine learning for tuning the bias and variance of a specific machine learning mechanism via the construction of a stronger classifier using multiple weak components, such that the ensemble, as a whole, outperforms the best individual classifier.The effectiveness of ensemble systems has been extensively demonstrated empirically and theoretically [12], although there does not currently exist any overarching theoretical framework capable of e.g.covering the requirements of ensemble components diversity to guarantee its out-performance.We here seek to provide and quantify a motivation for employing classical ensemble techniques in relation to NISQ-based quantum neural networks, which we address via the following three arguments.
The first argument concerns the potential for the superior performance of an ensemble system composed of small quantum neural networks compared to a single larger quantum neural network.This notion is based on the rationale that while quantum neural networks are inherently powerful machine learning models, they exhibit intrinsic variance due to the nature of highly non-convex loss landscape, implying that different predictors will result from randomlyinitialised stochastic gradient descent training, in common with classical neural networks.Modern deep learning practice often deliberately overparametrises the network in order to render the loss more convex [13], with the asymptotic case of infinitely wide neural networks exhibiting a fully convex loss landscape, making it effectively a linear model [14].Although overparameterization in quantum neural networks has been studied theoretically [15][16][17] and has been shown to be beneficial to generalization performances within certain settings, the increase in resource requirements makes this approach almost completely impractical on NISQ devices.In the classical literature, however, it has been demonstrated that ensemble techniques can perform comparably to the largest (generally overparameterized) models with significantly fewer resources, especially in relation to overall model parameterization, c.f. for example [18,Figure 2].
The second argument pertains to the resource savings achievable by ensemble systems, particularly in terms of the number of qubits, gates, and training samples required.For example, the boosting ensemble technique involves progressive dividing of the training dataset into multiple, partially overlapping subsets on the basis of their respective impact on the performance of the cumulative ensemble classifier created by summing of the partial weak classifiers trained on previously-selected data subsets.This enables the ensemble quantum neural network to be constructed in parallel with individual quantum neural networks operating on datasets of reduced size.The random subspace technique, by contrast, trains each base predictor on a random subset of features, but also provides an advantage in terms of the overall number of qubits and gates required.Employing the random subspace technique in a quantum machine learning setting would parallel the various quantum circuit splitting techniques (c.f. for example [19]), and divide-and-conquer approaches, that have been utilized in the field of quantum chemistry [20] and quantum optimization [21].
Our third argument, which is specific to quantum computing, examines the potential of ensembles' noise-canceling ability.Previous works have demonstrated that ensembles can enhance the performance of several noisy machine-learning tasks (see [22]).Our investigation aims to determine whether and to what extent these techniques can reduce the impact of noise during the execution on a NISQ device at the applicative level.This approach differs from most current approaches, which aim to reduce noise at a lower level, as described in [23].
We here examine the impact of ensemble techniques based on bagging (bootstrap aggregation) and boosting ensembles in a quantum neural network setting across seven variant data loading schemes.Bagging techniques are selected for their applicability in high-variance settings, i.e. those exhibiting significant fluctuations in relation to differ initialisations and differ sample subselections; contrarily, boosting techniques are effective in relation to high-bias models, i.e. those which are relatively insensitive to data subsampling.
Our first objective is to quantify the amount of resources (in particular, the number of qubits, gates, parameters, and training samples) saved by the respective approaches.Secondly, we evaluate the performance using quantum neural networks as base predictors to solve a number of representative synthetic and real-world regression and classification tasks.Critically, the accuracy and loss performance of these approaches are assessed with respect to the number of layers of the quantum neural networks in a simulated environment.We thus obtain a layer-wise quantification of performance that addresses one of the fundamental questions in architecting deep neural systems, namely, how many layers of abstraction to incorporate?Note that this question is fundamentally different in a quantum setting compared to classical neural systems; in the latter, the possibility of multi-level feature learning exists, and thus the potential for indefinite performance improvement with neural layer depth [17].This contrasts with the quantum neural networks, in which an increase in the number of layers affects the expressibility of the ansatz and thus might introduce a barren plateau [7].
Finally, the noise-canceling capabilities of ensembles will be investigated by testing a synthetic linear regression task both on a simulated noisy environment mimicking IBM's real quantum device Lima and on IBM's superconductor-based quantum processing unit (QPU) Lagos.

Contributions
Our contributions are the following: • We evaluate various ensemble schemes that incorporate bagging and boosting techniques into quantum neural networks, and quantify the benefits in terms of resource savings, including the number of qubits, gates, and training samples required for these approaches.• We apply our approach both to a simulated noisy environment and to the IBM Lagos superconductor-based quantum processing unit to investigate the potential advantages of bagging techniques in mitigating the effects of noise during the execution of quantum circuits on NISQ devices.• We conduct a layer-wise analysis of quantum neural network performance in the ensemble setting with a view to determining the implicit trade-off between ensemble advantage and layer-wise depth.

Related Works
The quest for quantum algorithms able to be executed on noisy small-scale quantum systems led to the concept of Variational Quantum Circuits (VQCs), i.e. quantum circuits based on a hybrid quantum-classical optimization framework [24,25].VQCs are currently believed to be promising candidates to harness the potential of QC and achieve a quantum advantage [26][27][28].VQCs rely on a hybrid quantum-classical scheme, where a parameterized quantum circuit is iteratively optimized with the help of a classical co-processor.This way, low-depth quantum circuits can be efficiently designed and implemented on the available NISQ devices; the noisy components of the quantum process are mitigated by the low number of quantum gates present in the VQCs.The basic structure of a VQC include a data encoding stage, where classical data are embedded into a complex Hilbert space as quantum states, a processing of such quantum states via an ansatz made of parameterized rotation gates and entangling gates, and finally a measurement of the circuit to retrieve the expected outcome.Many different circuit architectures and ansatzes have been proposed for VQCs [29][30][31][32], depending on the structure of the problem or on the underlying quantum hardware.VQCs demonstrated remarkable performances and a good resilience to noise in several optimization tasks and real-world applications.For example, researchers in [33] introduced a circuitcentric quantum classifier based on VQC that could effectively be implemented on a near-term quantum device.It correctly classified quantum encoded data and demonstrated to be robust against noise.Authors in [25] proposed a VQC that successfully approximated high-dimensional regression and classification functions with a limited number of qubits.
VQCs are incredibly well-suited for the realization of quantum neural networks with a constraint on the number of qubits [34].A quantum neural network is usually composed of a layered architecture able to encode input data into quantum states and perform heavy manipulations in a high-dimensional feature space.The encoding strategy and the choice of the circuit ansatz are critical for the achievement of superior performances over classical NNs: more complex data encoding with hard-to-simulate feature maps could lead to a concrete quantum advantage [35], but too expressive quantum circuits may exhibit flatter cost landscapes and result in untrainable models [7].An example of quantum neural network was given in [36], where a shallow NN was employed to perform classification and regression tasks using both simulators and real quantum devices.In [37], authors proposed a multi-layer Quantum Deep Neural Network (QDNN) with three variational layers for an image classification task.They managed to prove that QDNNs have more representation capacity with respect to classical deep NN.A hybrid Quantum-classical Recurrent Neural Network (QRNN) was presented in [38] to solve a time series prediction problem.The QRNN, composed of a quantum layer as well as two classical recurrent layers, demonstrated superior performances over the classical counterpart in terms of prediction error.However, quantum neural networks suffer from some non-negligible problems, which deeply affect their performances and limit their impact in the quantum ecosystem.Firstly, they are still subject to quantum noise, and it gets worse as the number of layers (i.e., the depth of the quantum circuit) increases [39,40].Secondly, barren plateaus phenomena may occur depending on the ansatz and the number of qubits chosen, reducing the trainability of such models [6,7,41].Finally, data encoding on NISQ devices continues to represent an obstacle when the number of features is considerable [34], making them hard to implement and train [38].
In classical ML, ensemble learning has been investigated for years to improve generalization and robustness over a single estimator [11,42].Ensembling is based on the so-called "wisdom of the crowd" principle, namely it combines the predictions of several base estimators with the same learning algorithm to build a single stronger model.Despite there are many different ensemble methods, the latter can be easily grouped into two different categories: bagging methods, which build and train several estimators independently and then compute an average of their predictions [43], and boosting methods, which in turn train the estimators sequentially so that the each one corrects the predictions of the prior models and output a weighted average of such predictions [44].Ensemble methods for NNs have also been extensively studied, yielding remarkable performances in both classification and regression tasks [45][46][47][48][49][50][51][52][53].Authors in [54] have shown that overparameterization renders an optimization problem easier to train.
In the quantum setting, the adoption of an ensemble strategy has received little consideration in the past few years, with very few approaches focusing on near-term quantum devices and VQC ensembles.In [55,56], the authors exploit the superposition principle to obtain an exponentially large ensemble wherein each instance is weighted according to its accuracy on the training dataset.However, they make use of a fault-tolerant approach rather than considering limited quantum resources.A similar approach is explored in [57], where authors create an ensemble of Quantum Binary Neural Networks (QBNNs) with reduced computational training cost without taking into consideration the amount of quantum resources necessary to build the circuit.An efficient strategy for bagging with quantum circuits is proposed in [58] instead.Very recently, [59] has proposed a distributed framework for ensemble learning on a variety of NISQ quantum devices, although it requires many NISQ devices to be actually implemented.A quantum Error-correcting output codes multiclass ensemble approach was proposed in [60].In [61], the authors investigated the performance enhancement of a majority-voting-based ensemble system in the quantum regime.Authors in [62] studied the role of ensemble techniques in the context of quantum reservoir computing.Finally, an analysis of robustness to hardware error as applied to quantum reinforcement learning, and presenting compatible results, is given in [63].
In this paper, we propose a classical ensemble learning approach to the outputs of several quantum neural networks in order to reduce the quantum resources for a given quantum model and provide superior performances in terms of error rate over single quantum neural network instances.To the best of our knowledge, no one has ever proposed such an ensemble framework for VQCs.We also compare both bagging and boosting strategy to provide an analysis on the most appropriate ensemble methods for quantum neural networks in a noiseless setting.An error analysis with respect to the number of layers of the quantum neural networks reveals that bagging models greatly outperform the baseline model with low number of layers, with remarkable performances as the number of layers increase; in fact, sufficiently complex bagging models allow to select better points on the bias-variance trade-off curve, such that one can maximise generalisability in a way not always possible with a single learner (especially ones with intrinsic constraints such as the QNN) [64,65].Finally, we apply our approach both to a simulated IBM Quantum Lima noisy backend and to the IBM Lagos superconductor-based QPU to investigate the potential advantages of bagging techniques in mitigating the effects of noise during the execution of quantum circuits on NISQ devices.

Background and Notation
We provide a brief introduction to the notation and concepts used in this work.The sets X and Y represent the set of features and targets, respectively.Typically, X is equal to R d , with d equal to the dimensionality in input, whereas Y is equal to R for regression tasks and Y is equal to {c 1 , ..., c k } for k-ary classification tasks.Sequences of elements are indexed in the apex with x (j) , where the i-th component is denoted as x i .The notation ϵ ∼ N (µ, σ 2 ) indicates that the value of ϵ is randomly sampled from a univariate normal distribution with mean µ and variance σ 2 .We use the function P to denote one when the predicate P is true and zero otherwise.

Models in quantum machine learning
We define the state of a quantum system as the density matrix ρ having unitary trace and belonging to the Hilbert space H ≡ C 2 n ×2 n where n is the number of qubits.The system starts in the state ρ 0 = |0⟩⟨0|.The evolution in a closed quantum system is described by a unitary transformation U = exp(−itH), t ∈ R, H Hermitian operator, and acts like ρ → U † ρU .The measurement of the system in its computational basis {Π i = |i⟩⟨i|} 2 n −1 i=0 applied to the system in the state ρ will give outcome i ∈ 0, 1, ..., 2 n − 1 with probability Tr[Π i ρΠ i ] after which the state collapses to ρ ′ = Π i ρΠ i / Tr[Π i ρΠ i ].The expectation value of a Hermitian operator associated to a physical observable, O = i λ i Π i , acting on the system on the state ρ, is given by the Born rule ⟨O⟩ = Tr [ρO].
Quantum computation can be described using a quantum circuit, a sequence of gates (i.e.elementary operations) acting on one or more qubits of the system terminating with the measurement operation over some or all of its qubits.The output of the measurement can be post-processed using a classical function."The set of gates available shall be universal ", i.e. the composition of such elementary operation allows the expression of any unitary transformation with arbitrary precision.An exemplar universal gate set is composed of parametric operators R (i) z ), and the operator CNOT (i,j) = exp(−i π 4 σ (i) x σ (j) x ).The gate I is the identity.The matrices are the Pauli matrices.The apex denotes explicitly the qubits in which the transformation acts.Quantum machine learning forms a broad family of algorithms, some of which require fault-tolerant quantum computation while others are ready to execute on current generation 'NISQ' (noisy) quantum devices.The family of NISQ-ready techniques of interest in this document is denoted variational quantum algorithms [24].These algorithms are based on the tuning of a cost function C(θ) dependent on a set of parameters θ ∈ [0, 2π] P and optimized classically (possibly via gradient descent-based techniques) to obtain the value θ * = arg min θ C(θ).Optimization through gradient-descent thus involves computation of the gradient of C.This can be done using finite difference methods or else the parameter-shift rule [66].The parameter-shift rule is particularly well-suited for NISQ devices as it can utilise a large step size relative to finite difference methods, making it less sensitive to noise in calculations.
In general, C(θ) is a function corresponding to a parametric quantum transformation U (θ) of a length polynomial in the number of qubits, the set of input states {ρ i }, and the set of observables {O k }.Specifically, a quantum neural network is a function in the form where ρ 0 is the initial state of the system, V (x) is a parametric quantum circuit depending on the input parameters x ∈ X , U (θ) is a parametric quantum circuit named an ansatz that depends on the trainable parameters θ ∈ [0, 2π) P , and O is an observable.Given the training dataset {(x (i) , y (i) )} M i=1 ∈ (X × Y) M , the cost function of a quantum neural network, being a supervised learning problem, is the empirical risk where ℓ : Y × Y → R is any convex loss function, e.g. the Mean Squared Error (MSE).
The quantum neural network constitutes a linear model in the Hilbert space of the quantum system as a consequence of the linearity of quantum dynamics.It behaves, in particular, as a kernel machine that employs the unitary V (x) as the feature map ρ → ρ x = V (x)ρV † (x), while the variational ansatz ρ → ρ θ = U (θ)ρU † (θ) adjusts the model weights.Note that although the model is linear in the Hilbert space of the quantum system, the measurement projection makes it nonlinear in the parameter space, enabling a set of rich dynamics; nevertheless, this is not the only way of introducing nonlinearity into the quantum model [67].Quantum neural networks can have a layer-wise structure, i.e., U (θ) = ℓ i=1 U i (θ i ), which provides it with further degrees of freedom for optimization (however, due to the lack of nonlinearity between the layers, the model does not possess the hierarchical feature learning capabilities of classical neural networks).
The selection of the ansatz is thus a crucial aspect in defining the quantum neural network, and it is required to adhere to certain classifier-friendly principles.Expressibility is one such, being the model's ability to approximate any quantum state in the Hilbert space.Although there are various ways to formalize expressibility, one of the most widely used definitions is based on the generation of the ensemble of states {ρ θ = U (θ)ρ 0 U † (θ) | θ ∈ Θ} and the standard ensemble of random states induced by the Haar measure over the corresponding unitary group.Expressible unitaries are those that make small the deviation between the former and the latter ensembles.However, expressible circuits are susceptible to the barren plateau problem, where the variance of the gradient decreases exponentially with the number of qubits, making parameter training infeasible.The varieties of ansatz and their expressibilities are presented in [68].Expressibility is tightly connected to the concept of controllability in quantum optimal control, and authors in [8] show that the asymptotic limit of the number of layers ℓ → ∞ in the expressible circuits are the controllable ones, i.e. those whose ansatz is underlied by a Lie algebra matching the space of skew-Hermitian matrices u(2 n ).

Ensemble techniques
The purpose of using ensemble systems is to improve the generalization performance through reducing the bias or variance of a decision system.Such a result is obtained by training several models and combining the outcomes according to a combination rule.A large body of literature on ensemble techniques exists; the reader is referred to [11] for a general overview.
The idea behind the ensemble system may be motivated by Condorcet's jury theorem [12]: a jury of m peers, each having probability p = 1 2 + ϵ, 0 < ϵ ≪ 1, of giving the correct answer, implies that the probability of the verdict given by majority voting to be correct is and quickly approaches 1 as m → ∞.The theorem, broadly interpreted, suggests that a combination of small, individually ineffective machine learning models h 1 , ..., h m (weak learners) can be combined to constitute a more powerful one, with arbitrarily good performance depending on the nature of data manifold and the base classifiers h ens (strong learner ).According to [11], three aspects Fig. 1: Taxonomy of the three aspects characterizing an ensemble system.characterize an ensemble system: a data selection strategy, the composition plus training strategies of the single model instances, and the combination rule of its output.Some of the possible choices are summarized in Figure 1.
The data selection strategy determines how the data should be distributed to the individual instances.If all instances are trained on the same dataset, their predictions will be highly correlated, resulting in similar output.The bootstrapping technique creates smaller, overlapping subsets by sampling with replacement from the dataset, which are then assigned to different instances.Alternatively, the pasting technique can be used for processing larger datasets by subsampling without replacement.Another approach is to divide the dataset by randomly assigning different sets of features with replacement, known as the random subspace technique (when the bootstrapping and random subspace techniques are combined, the result is the random patch technique).
There are numerous schemes for combining predictors, with bagging being the most straightforward and commonly used.Bagging, short for bootstrap aggregation, involves the creation of multiple homogeneous model instances trained on bootstrapped datasets.An instance of a bagging scheme is the random forest, which involves bagging decision trees trained on differing sample subsets (in some cases, random forests may favor a random patch data selection strategy over bagging).Another predictor combination scheme is boosting, which involves training a sequence of predictors via subsampling data according to the following strategy: an initial predictor is trained on a uniformly drawn subset of samples, while the i-th instance of the predictor is trained on a subset of elements that the previous ensemble classifier incorrectly predicted.The bagging ensemble trains the models in parallel over a subset of the dataset drawn uniformly; each prediction is then merged via an average function.The boosting ensemble trains the models sequentially, the first predictor draws the samples uniformly, and the subsequent models draw the elements from a probability distribution biased toward previously misclassified items.
The ensemble is itself the convex cumulative sum over predictors.Numerous variations of boosting exist, one of the most notable being AdaBoost [69].Contrary to vanilla boosting, AdaBoost employs an exponential loss such that the ensemble error function allows for the fact that it is only the sign of outcome that is significant.These two scheme are illustrated in Figure 2. The other major ensemble scheme is stacking in which a collection of heterogeneous classifiers trained on the same dataset are combined via an optimised meta-classifier.The combination rule merges the output of individual models h 1 , ..., h m .In classification tasks i.e.where the label output is discrete y ∈ C = {c 1 , ..., c k }, the most commonly used rule is majority voting.This is calculated as y ens = arg max c∈C m i=1 h i (x) = c .Where there exists prior knowledge regarding the performance of individual predictors, positive weights w i can be assigned, such that the output is a weighted majority vote.The ensemble prediction in this case will be y ens = arg max c∈C m i=1 w i h i (x) = c .Alternatively, the borda count method sorts labels in descending order by likelihood, with the ensemble prediction being the highest ranking sum.Nevertheless, averaging functions can also be utilised for ensemble classifiers for small ensemble sizes [70].For regression tasks where y ∈ R, common combination rules are (possibly weighted) mean, minimum, and maximum.

Discussion
Ensemble techniques, while well-established in the classical realm, have been largely overlooked in the quantum literature, leaving a number of open questions in this setting, such as whether bagging techniques, which reduce variance, can be deployed as effectively as boosting techniques, which reduce bias (both of which are also data-manifold and base-model dependent).It is also unclear as to the relative resource saving in terms of circuit size (number of qubits) and depth (number of gates), and also samples required for training, that can be obtained by using an ensemble of quantum neural networks instead of a single, large quantum network.Furthermore, it is not currently well understood the extent to which an ensemble system can mitigate hardware noise.Our experiments are designed to explore these questions.
To investigate the first two aspects, we conduct a suite of experiments within a simulation environment, employing seven distinct ensemble schemes with varying strategies for data selection, model training and decision combination applied to four synthetic and real-world datasets, encompassing both regression and classification tasks.Specifically, we analyze: a synthetic linear regression dataset, the Concrete Compressive Strength regression dataset, the Diabetes regression dataset, and the Wine classification dataset, which are widely used benchmarks for evaluating machine learning models.
Six of the proposed techniques are classified as bagging methods, employing bootstrapped data to generate the ensemble, while the seventh is a sequential boosting technique, namely AdaBoost.In particular, we implemented the AdaBoost.R2 version [71] for the regression tasks and the AdaBoost SAMME.R version [72] for the classification problem.The bagging ensembles are characterized by two parameters: the sample ratio r n ∈ [0, 1], which determines the percentage of training samples used for each base predictor (with replacement), and the feature ratio r f ∈ [0, 1], which indicates the percentage of features used for each predictor (without replacement).We test six bagging schemes by varying (r n , r f ) ∈ {0.2, 1.0} × {0.3, 0.5, 0.8}.For both the classification and regression tasks, the outputs of the base predictors are combined via averaging, as suggested in [70].In the case of the AdaBoost ensemble, the training set for each base predictor has the same size and dimensionality as the original training set.However, the samples are not uniformly drawn but are selected and weighted based on the probability of misclassification by previous classifiers composing the cumulative ensemble; single predictors are hence combined using a weighted average.Each ensemble system comprises 10 base predictors.The characteristics of these ensemble schemes are summarized in Table 1, where FM identifies the baseline quantum neural network model, whereas Bag r f r n represents a bagging model with r f percentage of the features and r n percentage of the samples.Our experiments aim to evaluate the performance of each of the ensemble frameworks in comparison to the baseline model, as well as to assess the overall resource saving, including the number of qubits and overall parametric requirements.
To investigate the impact of quantum hardware noise, we conduct two additional experiments.The first one is performed in a noisy simulated setting with PennyLane-Qiskit plugin, which allows to simulate of a noisy device by selecting a fake backend from IBM Quantum's suite.It mimics the behaviors of IBM Quantum systems using system snapshots.The latter contains information about the quantum system such as coupling map, basis gates, qubit properties (T1, T2, error rate, etc.).These elements are useful for incorporating a realistic noise model into our experiment and thus performing noisy simulations of the system.We selected IBM Quantum Lima as the quantum backend to mimic.It is a 5-qubit superconducting-based quantum computer, and its topology is depicted in Figure 3a.In this experiment, the baseline FM is compared to the Bagging models, with the aim to assess their robustness in the presence of quantum hardware noise and qubit coupling.The second noisy experiment was performed directly on the IBM Lagos QPU real hardware.Such a device is a 7-qubit superconducting-based quantum computer.The topology of Lagos is depicted in Figure 3b.Specifically, here we compare the performance of the baseline model FM with that of the Bag 0.8 0.2 configuration on the linear regression dataset.Our goal with these two experiments is to determine whether ensemble techniques can effectively mitigate quantum noise, and whether the difference in performance between single predictors and ensemble systems is more pronounced within a simulated, noise-free environment in comparison with noisy executions on quantum hardware.

Experimental setup
This section outlines experimental protocols used to evaluate the performance of the various ensemble approaches in terms of both the experimental structure and specific parameters/settings used to configure the algorithm and hardware.

Choice of quantum neural networks
We utilize a quantum neural network of the form f , which operates on n qubits, with n corresponding to the number of features in the classification/regression problem.For the feature map, we opted for the simple parametric transformation y (x i ).This choice was motivated by the findings in [73], suggesting that more complex feature maps can lead to unfavorable generalization properties, the incorporation of which may thus unnecessarily bias our findings.(In [74], various feature maps are compared).
The ansatz is implemented with the parametric transformations structured layer-wise with, for ℓ the number of layers, a total of 3ℓn parameters.It is thus defined as: The role of CNOT gates is the introduction of entanglement in the system, which would otherwise be efficiently classical simulable.We select as the observable z , which operates on a single qubit.Local observables like this one are less susceptible to the barren plateau problem than global ones, for example, z (as noted in [41]).Note that the measurement on the first qubit depends nontrivially on the overall state due to the entanglement introduced by the chosen unitary.The quantum neural network described in our investigation is pictured in Figure 4.

Training of the model
To train models, we utilize a standard state-of-the-art gradient descent-based algorithm, ADAM.The MSE was selected as the loss function and error metric to evaluate the performances of the models in the regression tasks, as it is a standard error metric in supervised learning.MSE was selected as the loss function to train the networks because it is more sensitive to larger errors.Categorical Cross Entropy (CCE) was used as the loss function for the classification task instead, while Accuracy score was employed as the error metric to assess the goodness of the classification.Given the output f of the model, the computation of its gradient ∇f , which is required to calculate the gradient of the loss function, is accomplished using the parameter-shift rule Fig. 4: Quantum Neural Network used to classify the linear regression dataset, having 5 qubits and ℓ = 1 layers.The rotational gates parameterized by the feature x i form the feature map, while those parameterized via the θs form the ansatz.
[66], since the commonly-used finite difference method ∇f (x; θ) ≈ (f (x; θ) − f (x; θ + ϵ))/ϵ is highly susceptible to hardware noise.The optimization hyperparameters used are the learning rate, set to 0.1, and the number of training epochs, which were selected through empirical investigation.Specifically, we carry out 150 training epochs to obtain the simulated noise-free results, while for noisy simulations and real QPU-based results, we perform just 100 and 10 epochs respectively, due to technological constraints on current hardware.

Datasets
We assess the performance of our approach using both synthetic and realworld datasets, across both regression and classification problems.The linear regression dataset is artificially generated with parametric control over the number of samples n, the dimensionality d, and the noise variance σ.It is procedurally generated by randomly sampling a weight vector w uniformly over [−1, 1] d such that the training set {(x (i) , y (i) )} n i=1 is constructed with x (i) uniformly sampled from [−1, 1] d , y (i) = w • x (i) + ϵ (i) , and ϵ (i) sampled from a normal distribution with zero mean and variance σ.In our case we have n = 250 (jointly the training and testing datasets), d = 5 and σ = 0.1.The other datasets involved in the experiments are the Concrete Compressive Strength dataset, the Diabetes dataset, and the Wine dataset.The first of these is a multivariate regression problem calculating the strength of the material based on its age and ingredients.The second is a multivariate regression problem correlating the biological and lifestyle characteristic of patients to their insulin levels.The third one is a multivariate, three-class classification problem investigating the geographic origin of wine samples from their chemical characteristics.All are freely available and open source.Table 2 summarizes the characteristics of these datasets; more details can be found in [75][76][77]

Simulated Noiseless Experiments
Initially, we evaluate our method in a simulated environment, one free of noise, such that the output estimation is infinitely precise.This differs significantly from execution on a NISQ quantum processing unit, which introduces various types of hardware error (such as decoherence and infidelity of operations) as well as sampling error caused via the measurement operation.We examine the performance of both the baseline models and ensemble systems in a scenario where the number of layers (i.e.quantum neural network depth) is gradually increased.To establish robustness to random initialization of parameters (that is, susceptibility to local minima effects), each simulation is repeated ten times.

Experiment I
The first experiment seeks to perform linear regression on a synthetic noisy 5-dimensional dataset.The function generating the targets is as follows: y = w•x+ϵ, where x ∈ (−1, 1) 5 ⊆ R 5 , w ∈ R 5 is randomly generated from a uniform distribution having as support the range −1 to 1, and ϵ is a Gaussian noise of mean zero and standard deviation 0.1.The total number of samples composing this synthetic dataset is 250.Each experimental data point instantiates a layer number, a number of bagged features, and a percentage of training data points available to the ensemble.The results of the first experiment are indicated in Figure 5.Both FM and AdaBoost achieve the lowest MSE on test of about 0.021 at 10 layers, reaching a performance plateau at 5 layers.The bagging models utilising 80% of the features are able to reach satisfactory results with 10 layers, which are only 0.03 -0.05 points higher than the error obtained by the best performing models.In general, it appears that quantum bagging models with a high number of features are able to generalize well on unseen data in this setting, even with only 20% of the training samples (unsurprisingly, the performance of bagging models with only 20% of training samples are worse than those of the counterparts using 100% of the training samples).Nevertheless, they still achieve remarkable results and show impressive generalization capabilities, confirming the effectiveness of bagged quantum models in generalizing well with relatively little training data [78].
It is also notable that all of the bagging models have a lower MSE test error as compared to FM and AdaBoost when the number of layers is low.In particular, with just 1 layer, all of the bagging models outperform FM and AdaBoost.However, as the number of layers increases, the performances of bagging models begin to plateau more rapidly than FM and Adaboost which, in contrast, continue their trend of decreasing error with increasing circuit depth.This is due to bagging models reaching their performance saturation point more quickly as the number of layers increases, so that additional layers do not significantly improve their predictive power.The primary cause for this phenomenon is the impossibility of bagging models in utilizing all the available information from the features of the samples; an increment in the number of bagging estimators may alleviate such phenomenon.On the other hand, FM and AdaBoost have the advantage of fully leveraging all the features of the samples for the prediction.As a result, they gain benefits from higher circuit depth, enabling them to progressively diminish their error rates.
Finally, the decreasing error trend seen in the more complex bagging models as well as the FM and AdaBoost models is not visible in relation bagging with 30% of the features.We conjecture that since this bagging configuration utilises only 1 qubit, it cannot appropriately model the evolution of the quantum state with respect to the input.Hence, despite leveraging 10 different submodels of 1 qubit (i.e., one feature) each, the performance of bagging models with 30% of the features cannot improve as the number of layers increases (adding more layers in this case translates in performing rotations on the single qubit only, without the possibility of further CNOTs or other entangling gate operations).This result hence highlights the importance of entanglement in quantum neural network models as a means of improving performance.In order to provide a comprehensive evaluation of the models' performances, we hereafter report the results obtained by varying not only the characteristics of the single base learners (i.e. the number of layers, the number of samples r n , and the number of features r f ), but also the ensemble size (i.e., the number of base learners that are averaged).As illustrated in Table 3, the ensemble works as expected, since the test error decreases while increasing the ensemble size.

Experiment II
The second experiment seeks to assess the performance of the respective ensemble techniques on the Concrete Compressive Strength dataset, which consists Table 3: MSE of bagging ensembles by varying the number of bagging estimators in Experiment I.Only even number of layers are reported for brevity. in 1030 samples of 8 features.The target value to predict in this regression case is hence the concrete compressive strength, measured in Megapascal (MPa), a highly nonlinear function of age and composition of the material.
The results of the regression experiment are in line with the findings of Experiment I, and are reported in Figure 6.FM, AdaBoost and the two bagging models applied in relation to 80% of features achieve comparable results at 10 layers, with the Bag.0.8 1.0 configuration obtaining the lowest MSE, followed by Bag.0.8 0.2, FM and finally by AdaBoost.Also in this case, the differential between bagging models with 20% of samples and with 100% of samples is marginal, confirming the effectiveness of bagging quantum models in relation to reduced training dataset size.In contrast with Experiment I, bagging models having 30% of available features now have 2 qubits, and therefore demonstrate a relative improvement in test error when l = 2.However, their expressive power soon saturates and their error curves plateau.
In general, the generalization capability of bagging models decreases monotonically with the number of layers, in contrast to FM and AdaBoost.In fact, they exhibit episodes of increased test errors when utilising 5 (and up to 7) layers, while bagging appears to be able to evade this outcome.Such an increase in the test error of FM and AdaBoost may be due to the complex landscape of the cost function, which in some cases could be challenging to optimize and lead to bad generalization.
All of the bagging models analyzed still outperform FM and AdaBoost at a low number of layers, suggesting that they may be the right choice for implementation on NISQ devices, or else when there is any necessity of implementing low-depth quantum circuits.As in the first experiment, it is also of interest to note that all the bagging models with l = 1 here have very similar MSE values, while their performances vary as the number of layers increases.This may indicate that the MSE value reached at l = 1 is the optimal for that family of bagging models, given their expressibility.Moreover, a sharp decrease in MSE beyond the first layers would appear to be a common pattern, both with respect to the ensembles and the FM model.For example, at l ≥ 3, the MSE of FM and AdaBoost dramatically decrease, while bagging models with 50% of the features exhibit this trend between l = 1 and l = 2. (A future analysis of this topic might seek to exploit this characteristic in order to predict a priori how many layers one would need to attain an error level within a given bound).
As for Experiment I, we report here the MSE of the bagging models with respect to the variation of the ensemble size (Table 4).It confirms the trend for which as the ensemble size increases, the MSE on test decreases.

Experiment III
The dataset used in Experiment III is the reference Diabetes dataset from Scikit-learn, consisting of 10 numerical features, including age, sex, body mass index, blood serum measurements, and also a target variable, a quantitative measure of disease progression one year after baseline.The dataset is composed of 442 instances and is often used for non-trivial regression analysis in ML.
Figure 7 illustrates the results of this experiment.The performance of the quantum models is notably different from those of the previous two experiments.It may be seen that the best-performing models are the bagging models containing 80% of the features for almost any number of layers, while FM and AdaBoost achieve satisfactory results up to 6 layers, at which point their MSE Table 4: MSE of bagging ensembles by varying the number of bagging estimators in Experiment II.Only even number of layers are reported for brevity.begins to increase.At l = 10, every model has stabilized, however.Bag.0.8 1.0 and Bag.0.8 0.2 have an MSE of respectively 8.8% and 6.1% lower than that of FM.AdaBoost has an MSE comparable to the error of Bag.0.3 1.0, being only 0.9% higher than FM.Bagging models with 50% of the features have Fig. 6: Evolution of MSE with respect to the number of quantum neural network layers in Experiment II.surprisingly good results, better than those of FM and very close to bagging models with 80% of the features.
As in Experiment I and II, a very sharp MSE reduction between l = 1 and l = 3 is evident for all of the models.Less complex models like bagging with 30% and 50% of the features immediately reach a plateau, while the error curves for bagging with 80% of the features, FM and AdaBoost evolves as the number of parameters increases.Considering layer numbers between l = 6 and l = 8, it is possible that the capacity of FM and AdaBoost models saturates as the number of model parameters increases, and thus they perform poorly on both training and test data; they struggle in navigating the cost function and escaping from local minima during training.In particular, their learning is constrained to such an extent that they almost reach the same performance level of the simplest bagging models with 30% of the features.The latter show no indication of bad generalization however, in common with bagging models having 50% of the features.Bagging with 80% of the features shows light signs of bad generalization and learning saturation when l > 6, but still achieve the best results from among all of the tested algorithms.The robustness of bagging models to learning saturation with respect to AdaBoost and FM arises from their ability to reduce variance via averaging of decorrelated error across the predictions of each submodel.By contrast, when the number of layers is high, AdaBoost and FM utilise a model that is too complex and expressive for the underlying task, leading to learning saturation, bad generalization performances and difficulties in approximating the underlying function.In concordance with Experiment II, this result confirm the effectiveness of bagging models in improving the predictive performance of QNN models, especially in cases where the optimization of the cost function becomes challenging, in line with the classical counterpart.
In addition, this experiment also highlights more markedly the discrepancy between the error level of bagging models with the same number of features but a distinct number of training samples.The difference between the MSE of the bagging model with 30% and 20% of samples and that with 100% of samples is now far more apparent, suggesting that when the variance of the dataset is very high, even bagging models require a sufficient threshold of training samples to perform well in the NISQ setting.
Regarding the performances of the bagging models with different ensemble sizes, the results are displayed in Table 5.The findings further support the previously observed trend that as the size of the ensemble increases, the MSE decreases on the test dataset.

Experiment IV
For the classification task in Experiment IV, we used the reference UCI Wine dataset.It is a multi-class classification dataset corresponding to the results of Table 5: MSE of bagging ensembles by varying the number of bagging estimators in Experiment III.Only even number of layers are reported for brevity.a chemical analysis of wines grown within a specific region of Italy.It consists of 13 numerical features representing various chemical properties, such as alcohol, malic acid, and ash content, and a target variable indicating the class of the wine.The dataset has 178 samples and is a common baseline ML benchmark for low-parametric complexity classifiers.
Results from Experiment IV are reported in Figure 8.Although they cannot be directly compared to the previous results due to the intrinsically different nature of the problem, there are few comparative insights that can be gained from the respective plot of Accuracy curves.First, all the models except bagging with 30% of the features achieve the same accuracy score of 97.2% using 10 layers.The performances of Bag.0.3 0.2 and Bag.0.3 1.0 are still relatively strong, however, having an accuracy score of 94.2% and 96.9% respectively.Given the very low complexity of these two models, this is a striking result.
A further notable aspect of the Accuracy curves is that all ensemble models converge with far fewer layers than FM.In particular, they require 3 layers in order to reach a performance plateau on average, after which they saturate and the accuracy score reaches saturation as well.By contrast, FM struggles to achieve a comparable accuracy score, only achieving accuracy greater than 90% when l ≥ 7.This means that the ensemble models are able to learn and capture the complex relationships between the input features far more efficiently than FM, which requires a much deeper architecture to attain comparable results.This observation is particularly relevant when considering the implementation of these models on NISQ devices, where the number of qubits and the coherence time are severely limited.
As expected, bagging models with 100% of the samples obtain almost everywhere (especially with few layers), a higher accuracy score than their counterparts with 20% of the features given the same number of layers.This suggests that using more training samples can improve the performance of ensemble models provided that the number of layers is low, as it allows them to better capture the underlying patterns of class discriminability in the data.
Finally, Table 6 shows the performance of the bagging models with varying ensemble sizes.The outcomes confirm the previously observed pattern that increasing the ensemble size leads to an increment in the Accuracy calculated on the test set.This indicates that the use of larger ensembles can contribute to improving the accuracy and robustness of the quantum predictive models.

Resource efficiency of quantum neural network ensembles
Besides performance, resource efficiency is a key argument for the utilization of quantum neural network ensembles.Efficiency can be measured by various metrics: for example, number of qubits, gates, parameters, and training samples required to achieve comparable performance.
To determine the potential savings in the number of qubits we here deploy the random subspace technique (also known as attribute bagging or attribute bootstrap aggregation).Our experiments (cf Figure 9) suggest a potential saving of 20% to 80% of the total qubit budget via this approach.However, such a saving is made at the cost of the ensemble was a whole having the potential for less rich class-discrimination behaviour, dependent on both the sampling required to achieve full feature coverage and the nature of the underlying data manifold.A positive consequence of reducing the number of qubits, though, is that each quantum circuit will have fewer gates and parameters, resulting in improved noise robustness on real hardware (i.e less decoherence, higher overall fidelity), as well as faster gradient calculation (individual gradient calculations require P + 1 quantum circuit evaluations for P parameters).This allows for Fig. 8: Evolution of Accuracy score with respect to quantum neural network depth in Experiment IV. a saving of the parameter budget of up to 75% in the indicated experimental regime, while the saving on gates corresponds proportionately (cf Figure 4).Savings for each dataset and ensemble technique are as depicted in Figure 9.

Simulated Noisy Experiments
For the noisy simulated experiment, we compare the performance of the FM baseline with all the Bagging techniques on the same synthetic linear regression dataset used in Experiment I, in order to assess the ensemble's noise mitigation properties.Unfortunately, further experiments on the other datasets are currently unfeasible due to technological constraints, since noisy simulations with Qiskit's backend on large datasets are very time-consuming above 5-6 qubits.Moreover, such constraints allow to simulate up to 3 layers only; simulation of more complex models would demand a significant amount of runtime and computational resources that are currently unfeasible.Nevertheless, they are enough to demonstrate the behavior and goodness of the ensemble techniques in the presence of quantum hardware noise, as well as the evolution of the models' performance as the number of layers increases.Each experiment is repeated 5 times to ensure statistical validity.From Figure 10a, it is evident that FM is highly affected by noise, while sufficiently complex Bagging techniques are more resilient to such errors and manage to achieve better results, even outperforming FM in the case of Bagging with 80% of the features.This discrepancy in performance is expected to accentuate as the number of layers (i.e., the depth) of the circuit or the number of ensemble estimators increase.These results confirm the goodness of the Bagging approach in mitigating the effect of noise for QNNs if compared to the outcomes of the same experiment in an ideal setting, as illustrated in Figure 10b for the best-performing techniques (FM and Bagging with 80% of the features).

Experiments executed on superconducting-based QPU
For the real-hardware evaluation, we compare the performance of the baseline quantum neural network with the Bag 0.8 0.2 ensemble on the same synthetic linear regression dataset used in Experiment I. We selected the Bag 0.8 0.2 model as representative ensemble technique for its outstanding performance in the simulated experiments despite the low number of training samples.To ensure statistical validity, we repeat each experiment 10 times.However, due to technological constraints on real quantum hardware, we analyze only the linear dataset with a quantum neural network having a single layer.Figure 11 presents the real-world experimental findings, which indicate that the bagging ensemble reduces the expected MSE by one-third and the expected variance by half when executed on quantum hardware, compared to the baseline model.Such results demonstrate that the noise-canceling capabilities of ensemble technique can be effectively exploited to work on NISQ devices in realistic settings.Additionally, the performance of the ten bagging models varied significantly, underlining the need to reinitialise the ensemble multiple times and validate it against a suitable validation dataset to ensure that the best model is selected.

Conclusion
We propose the use of ensemble techniques for practical implementation of quantum machine learning models on NISQ hardware.In particular, we justify the application of these techniques based on their capacity for significant reduction in resource usage, including in respect to the overall qubit, parameter, and gate budget, which is achieved via the random subspace (attribute bagging) technique.This resource-saving is especially crucial for noisy hardware, which is typically limited to a small number of qubits, being vulnerable to decoherence, noise, and operational errors.Consequently, the contribution of ensemble techniques may be seen as a form of quantum noise reduction.
To establish this, we evaluated and compared various configurations of bagging and boosting ensemble techniques on synthetic and real-world datasets, tested in a simulated noise-free environment, in a simulated noisy setting and on a superconducting-based QPU by IBM, and subtending a range of layer depths.
Our experimental findings showed that bagging ensembles can effectively train quantum neural network instances using fewer features and qubits, which leads to ensemble models with superior performance compared to the baseline model.Reducing the number of features in bagging models of quantum neural networks directly translates into a reduction in the number of qubits, that is a desirable characteristics for practical quantum applications.Ensembles of quantum neural network can also help addressing some of the toughest challenges associated with noise and decoherence in NISQ devices, as well as to mitigate barren plateau effects.These can be key considerations in the development of quantum machine learning models, particularly when working with limited resources on modern quantum systems.
Moreover, bagging models were found to be extremely robust to learning saturation, being able to effectively capture the underlying patterns in the data with high generalization ability.This makes them better suited for tasks where generalization is important, such as in real-world applications.As in the classical case, the increase of the ensemble size in case of bagging also helps to achieve better performances for QNNs.However, it is important to notice that the effectiveness of bagging quantum models diminishes with a decrement in the number of features, which suggests that complex bagging models are still needed to obtain satisfactory results.Using only a subset of the features can reduce the computational complexity of the model and help in the optimization process, but it may also result in a loss of information and a decrease in performance.On the contrary, the number of training samples do not seem to have a deep impact on bagging quantum models, hence this bagging strategy may be used when executing quantum neural network instances on real hardware in order to deal with long waiting queues and job scheduling issues.In this regard, having a low number of training data leads to faster training procedures and quantum resource savings.The training of ensembles can also be done in parallel on multiple QPUs in a distributed learning fashion.Therefore, it is important to strike a balance between model complexity and performance to achieve the best possible outcomes.
Additionally, the fact that the bagging models outperform FM and AdaBoost at low number of layers suggests that the former models are better suited for low-depth quantum circuits, which have limited capacity and are prone to noise and errors.For quantum machine learning tasks with NISQ devices, using bagging models with a low number of layers may be a good strategy to achieve good generalization performance while minimizing the impact of noise and errors in the circuit.
Overall, our results suggest that ensembles of quantum neural network models can be a promising avenue for the development of practical quantum machine learning applications on NISQ devices, both from a performance and resource usage perspective.A careful evaluation of the trade-offs between model complexity, performance, quantum resources available and explainability may be necessary to make an informed decision.
In a future work, we plan to further investigate the relationship between ensembles and quantum noise, which is a key consideration when developing quantum neural network models.In addition, it would be relevant to discuss how our proposed approach compares to classical ensembles for real-world applications; it would contribute to a more comprehensive and insightful understanding of the potential advantages and limitations of our approach.To sum up, our findings could potentially contribute to the development of more efficient and accurate quantum machine learning algorithms, which could have significant implications for real-world applications.to result in lower MSE or higher Accuracy, although more tests should be conducted to further study such a phenomenon and obtain statistically significant results with more samples.Table B4: Spearman correlation coefficient between Accuracy and ensemble size in Experiment IV (p-value ≤ 0.05).Zero values indicate that the alternative hypothesis is rejected, i.e. the coefficient is 0 and hence there is no correlation.This outcome was observed exclusively with the more complex ensemble, and can be attributed to its ability to attain favorable Accuracy values even when employing a relatively small number of estimators.Only even number of layers are reported for brevity.

Fig. 2 :
Fig.2: Comparison between bagging (left) and 'vanilla' boosting (right) techniques.The bagging ensemble trains the models in parallel over a subset of the dataset drawn uniformly; each prediction is then merged via an average function.The boosting ensemble trains the models sequentially, the first predictor draws the samples uniformly, and the subsequent models draw the elements from a probability distribution biased toward previously misclassified items.

Fig. 5 :
Fig. 5: Evolution of MSE with respect to the number of quantum neural network layers in Experiment I.Each experimental data point instantiates a layer number, a number of bagged features and a percentage of training data points available to the ensemble.

Fig. 7 :
Fig. 7: Evolution of MSE with respect to the number of quantum neural network layers in Experiment III.

Fig. 10 :
Fig. 10: Comparison of average performance of the baseline model and the Bagging ensemble models on simulated IBM Lima quantum hardware.(10a) shows the difference in terms of MSE over 5 executions with all the bagging techniques.(10b) shows the performance of FM and bagging models with 80% of the features both in an ideal and simulated noisy setting.

Fig. 11 :
Fig. 11: Comparison of average performance of the baseline model and the Bag 0.8 0.2 ensemble technique on IBM quantum hardware.(11a) shows the difference in terms of MSE over 10 executions.(11b) shows the performance of the bagging model with respect to its estimators.

Fig. A1 :
Fig. A1: Comparison of the performance of the baseline model and ensemble systems on the Linear Regression dataset.It exhibits the MSE and standard deviation, with a semi-transparent area, of the ensemble schemes in comparison to the baseline models.The top-left image shows ensembles with Random Subspace at 30% of the features, top-right shows ensembles with Random Subspace at 50%, bottom-left displays ensembles with Random Subspace at 80%, and bottom-right illustrates AdaBoost.

Fig. A2 :
Fig. A2: Comparison of the performance of the baseline model and ensemble systems on the Concrete Compressive Strength dataset.It exhibits the MSE and standard deviation, with a semi-transparent area, of the ensemble schemes in comparison to the baseline models.The top-left image shows ensembles with Random Subspace at 30% of the features, top-right shows ensembles with Random Subspace at 50%, bottom-left displays ensembles with Random Subspace at 80%, and bottom-right illustrates AdaBoost.

Fig. A3 :
Fig. A3: Comparison of the performance of the baseline model and ensemble systems on the Diabetes dataset.It exhibits the MSE and standard deviation, with a semi-transparent area, of the ensemble schemes in comparison to the baseline models.The top-left image shows ensembles with Random Subspace at 30% of the features, top-right shows ensembles with Random Subspace at 50%, bottom-left displays ensembles with Random Subspace at 80%, and bottomright illustrates AdaBoost.

Fig. A4 :
Fig. A4: Comparison of the performance of the baseline model and ensemble systems on the Wine dataset.It exhibits the average accuracy and standard deviation, with a semi-transparent area, of the ensemble schemes in comparison to the baseline models.The top-left image shows ensembles with Random Subspace at 30% of the features, top-right shows ensembles with Random Subspace at 50%, bottom-left displays ensembles with Random Subspace at 80%, and bottom-right illustrates AdaBoost.

Table 1 :
Characteristics of the baseline benchmark model (0) and ensemble systems (I to VII).The ensemble system is identified by its broad data loading method (BST for Boosting and RSBS for Random Subspace), predictor composition & training type (Ensemble), number of base predictors (#BP), composition rule (Rule, with Avg representing the average function and W.Avg representing weighted average).

Table 2 :
. Every dataset is divided into 80% train samples and 20% test samples.Moreover, in a data preprocessing phase, raw data were scaled in the range [−1, 1] to best suit the output of the quantum neural networks; the scaler was fitted using training data only.No other preprocessing technique, i.e.PCA, has been applied.Characteristics of the datasets analyzed.UCI stands for the open source UCI Repository available at archive.ics.uci.edu.Scikit-Learn is an opensource software library for Python3.The number of features does not include the target.
(3,4)ls Our implementation 1 is written in Python3, and utilizes Pennylane as a framework to define and simulate quantum circuits, with the Pennylane-Qiskit plugin used to execute circuits on IBM Quantum devices via the Qiskit software stack.To improve simulation times, we employed the JAX linear algebra framework as the simulation backend.By using JAX, the quantum circuit can be just-in-time compiled to an intermediate representation called XLA, which can significantly speed up simulation times (by up to a factor of 10).Our simulations were run on a commercial computer with an AMD Ryzen 7 5800X (8-core CPU with a frequency of 3.80 GHz) and 64 GB of RAM.The first experiment on the noise canceling properties of ensemble systems was simulated with Qiskit's FakeLimaV2 backend, which mimics the behaviour of IBM Quantum's Lima QPU.It consists of 5 qubits arranged in the topology {(0, 1);(1,2);(1,3);(3,4)}.The single-gate and CNOT fidelities of this QPU were 4.79e −4 and 1.07e −2 , respectively.The second experiment on the noise canceling properties was conducted on the ibm lagos quan-

Table 6 :
Accuracy of bagging ensembles by varying the number of bagging estimators in Experiment IV.Only even number of layers are reported for brevity.

Table B1 :
Spearman correlation coefficient between MSE and ensemble size in Experiment I (p-value ≤ 0.05).Only even number of layers are reported for brevity.

Table B2 :
Spearman correlation coefficient between MSE and ensemble size in Experiment II (p-value ≤ 0.05).Only even number of layers are reported for brevity.

Table B3 :
Spearman correlation coefficient between MSE and ensemble size in Experiment III (p-value ≤ 0.05).Only even number of layers are reported for brevity.