## Operator on Operator Regression in Quantum Probability

Suprio Bhar

IIT Kanpur

Department of Mathematics and Statistics

Kanpur 208016, India

Email: suprio@iitk.ac.in

Subhra Sankar Dhar

IIT Kanpur

Department of Mathematics and Statistics

Kanpur 208106, India

Email: subhra@iitk.ac.in

Soumalya Joardar

IISER Kolkata

Department of Mathematics and Statistics

Mohanpur 741246, India

Email : Soumalya@iiserkol.ac.in

August 2, 2024

## Abstarct

This article introduces operator on operator regression in quantum probability. Here in the regression model, the response and the independent variables are certain operator valued observables, and they are linearly associated with unknown scalar coefficient (denoted by β ), and the error is a random operator. In the course of this study, we propose a quantum version of a class of estimators (denoted by M estimator) of β , and the large sample behaviour of those quantum version of the estimators are derived, given the fact that the true model is also linear and the samples are observed eigenvalue pairs of the operator valued observables.

## 1 Introduction

Regression analysis is one of the most well-known toolkit in statistical modelling. More specifically, regression analysis studies the 'relationship' between the dependent/response variable and one or more independent/explanatory variables. This 'relationship' can be linear, nonlinear or even non-parametric as well (see, e.g., [4] and a few relevant references therein), and using various statistical techniques, one may estimate the unknown parameters involved in the 'relationship'. Recently, the concepts of regression have been extended to functional valued random elements also; for instance, the interested reader may look at [3] and a few relevant references therein regarding function on function regression. However, to the best of our knowledge, regression analysis is an area unexplored in the setting of quantum probability. This work aims to investigate various statistical issues involved in operator on operator regression in quantum probability. We first provide a literature review and then brief the research problem without much technical details in the following.

In the context of other statistical concepts studied in the language of quantum probability theory, recently, [11] extended the concepts of sufficient statistic and Rao-Blackwellization in the framework of quantum probability, and a few years before, [7] overviewed and reviewed the problem of statistical sufficiency in the setting of quantum statistics. They provided sufficient and necessary conditions for sufficiency and established many results from traditional statistical theory in the language of noncommutative geometry, and at the same time, [6] derived the noncommutative version of well-known factorization theorem in traditional Statistics (see, e.g, [10]).

Let us now consider a pair of operator valued observables ( X,Y ) (see, e.g., [5] for details on observables) with observed n pairs of eigen values ( λ , µ 1 1 ) , . . . , ( λ , µ n n ). Observe that here in the setting of quantum probability, we are observing the eigen values of ( X,Y ), unlike the observations on the random element ( X,Y ) as in traditional Statistics/classical Probability. Suppose that unobserved operators Y and X are linearly associated in true sense with scalar valued coefficient (see (2.1)), and in the sense of regression modelling, the operator valued observables are assumed to be linear with scalar valued coefficient having some operator valued

error (see (2.2)). In this work, based on ( λ , µ 1 1 ) , . . . , ( λ , µ n n ), we estimate (denoted by ˆ β n as an estimator) the unknown parameter β involved in the model described in (2.2). In the course of this study, we establish the connection between the concepts in quantum probability and classical probability (see Section 2.1), and in view of this relation, we investigate the large sample properties ˆ β n (see Section 3).

The rest of the article is structured as follows. Section 2 describes the preliminaries of the model in terms of quantum probability, and for this problem, a key connection between the quantum probability and classical probability is investigated in Section 2.1, and the problem is reformulated in terms of classical probability in Section 2.2. The large sample properties of the proposed estimator is studied in Section 3, and Section 4 contains a few relevant concluding remarks. Finally, necessary technical details are provided in the Appendix in Section 5.

## 2 Preliminaries : Model

Let H be a real separable Hilbert space with inner-product 〈· , ·〉 H and norm ‖ · ‖ H . If T and S are two bounded linear operators on H related by T = βS for some non-zero real scalar β , then given any eigen-value λ of S , we have an eigen-value µ = βλ of T with the same eigen vector of unit norm. Again, given any eigen-value µ of T , we have an eigen-value λ = µ β of S with the same eigen vector of unit norm. In light of this fact, one may expect that the linear regression model associating two operator valued observables can be represented by their corresponding eigen values. The research problem is foramlly described in the following.

Let X and Y be two compact self-adjoint operator (on H ) valued observables related by

$$Y = \beta _ { 0 } X,$$

where β 0 ∈ R is unknown. We now consider the operator on operator regression model

$$Y = \beta X + \epsilon,$$

where /epsilon1 is a random operator valued error, and β ∈ R is an unknown constant. In this work, our objective is to propose an estimator (denoted by ˆ ) of the unknown β n β based on the observed

eigenvalues of ( X,Y ). The technical assumptions are stated in Section 3, where the large sample properties of ˆ β n are studied.

Suppose that ( λ , µ 1 1 ) , . . . , ( λ , µ n n ) are n pairs of observed eigen-values of ( X,Y ). Motivated by our earlier observation on deterministic operators T and S , we assume that each pair ( λ , µ j j ) appear with a common eigen-vector v j of unit norm. Here, as mentioned earlier, using the observed ( λ , µ j j ) n j =1 , we would like to estimate the unknown parameter β . The estimation procedure of β is described in the following.

In order to carry our the estimation of β , we first reformulate the problem, from originally stated in terms of X and Y in the quantum probability setup to their eigen-value pairs ( λ, µ ) in the classical probability setup. Next, using the classical probability model thus obtained the problem is reduced to a standard linear regression problem and we can then use classical statistical tools in order to estimate β . Section 2.1 studies the eigen decomposition of operator valued observables X and Y , which is the key concept in reformulating the problem in the setup of classical probability.

## 2.1 Distribution of eigen decomposition of X and Y : From quantum to classical probability

Assume that X and Y are compact self-adjoint operators on H with discrete spectrum given by σ X ( ) and σ Y ( ), respectively. As per our discussion above, the eigenspaces corresponding to each feasible eigen-value pairs are the same. Consequently, we have the following spectral decompositions (see [8]):

$$\vdots & X = \sum _ { \alpha \in \sigma ( X ) } \alpha P _ { \alpha }, \ \ Y = \sum _ { \alpha \in \sigma ( Y ) } \alpha P _ { \alpha }, \\ \text{able and } P \ \text{ denotes the projection onto the ei}$$

/negationslash where the sum is countable, and P α denotes the projection onto the eigen-space of any eigenvalue α . Here, the eigenvalues are all real and the eigen-spaces corresponding to 0 = α ∈ σ X ( ) are finite-dimensional. Now, given any state φ (unit trace), we have the expectations in the setting of quantum probability as

$$E X = t r ( \phi X ) = \sum _ { \alpha \in \sigma ( X ) } \alpha \, t r ( \phi P _ { \alpha } ), \ \ E Y = t r ( \phi Y ) = \sum _ { \alpha \in \sigma ( Y ) } \alpha \, t r ( \phi P _ { \alpha } ).$$

Note that tr ( φP α ) ≥ 0 , ∀ α ∈ σ X ( ). Moreover, ∑ α ∈ σ X ( ) tr ( φP α ) = tr ( ρ ) = 1. A similar statement is true when X is replaced by Y .

Remark 2.1 Observe that one can think of λ j 's drawn independently from a discrete probability distribution supported on σ X ( ) and with the probability mass function

$$f _ { \lambda } ( \alpha ) = t r ( \phi P _ { \alpha } ) 1 _ { ( \alpha \in \sigma ( X ) ) }.$$

Similarly, µ j 's are drawn independently from a discrete probability distribution supported on σ Y ( ) and with the probability mass function

$$f _ { \mu } ( \alpha ) = t r ( \phi P _ { \alpha } ) 1 _ { ( \alpha \in \sigma ( Y ) ) }.$$

This fact enables us to reformulate the problem in the set up of traditional statistics.

## 2.2 A reformulation : estimation of β

Given n pairs ( λ , µ 1 1 ) , . . . , ( λ , µ n n ) of observed eigen-values of ( X,Y ), assume that v , v 1 2 , · · · , v n be the corresponding eigen-vectors of unit norm. Then,

$$X v _ { j } = \lambda _ { j } v _ { j } \text{ for all } j = 1, \dots, n,$$

where for j = 1 , . . . , n , || v j || H = 1, and similarly, we have

$$Y v _ { j } = \mu _ { j } v _ { j } \text{ for all } j = 1, \dots, n.$$

Next, it follows from (2.2), (2.3) and (2.4) that for all j = 1 , . . . , n , Y v j = βXv j + /epsilon1v j , which yields

$$\mu _ { j } v _ { j } = \beta \lambda _ { j } v _ { j } + \epsilon v _ { j },$$

and consequently,

$$\mu _ { j } = \langle v _ { j }, \mu _ { j } v _ { j } \rangle _ { \mathcal { H } } = \langle v _ { j }, \beta \lambda _ { j } v _ { j } + \epsilon v _ { j } \rangle _ { \mathcal { H } } = \beta \lambda _ { j } + \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } \,, \, \forall j = 1, \dots, n.$$

We now estimate β , involved in (2.6), based on ( λ , µ j j ) n j =1 in various ways, by assuming conditions on the error terms 〈 v , /epsilon1v j j 〉 H , j = 1 , . . . , n . Now, recall Remark 2.1 that λ j 's are i.i.d. random variables following a discrete probability distribution supported on σ X ( ) and with the probability mass function

$$f _ { \lambda } ( \alpha ) = t r ( \phi P _ { \alpha } ) 1 _ { ( \alpha \in \sigma ( X ) ) },$$

and µ j 's are i.i.d. random variables following a discrete probability distribution supported on σ Y ( ) and with the probability mass function

$$f _ { \mu } ( \alpha ) = t r ( \phi P _ { \alpha } ) 1 _ { ( \alpha \in \sigma ( Y ) ) }.$$

For the data ( λ , µ j j ) ( j = 1 , . . . , n ), the M -estimator of β (see (2.6)), which is denoted by ˆ , is defined as β n

$$\L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L \L
 \hat { \beta } _ { n } \ & = \ a r g \min _ { \beta } \sum _ { j = 1 } ^ { n } \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } ) \\ & = \ a r g \min _ { \beta } \sum _ { j = 1 } ^ { n } \rho ( \mu _ { j } - \beta \lambda _ { j } ) \ ( \uplus \text{sing} \ ( 2. 6 ) ), \quad \ ( 2. 7 ) \\ \text{s a certain convex function, and few more assumptions will be stated in the sub-}$$

where ρ . ( ) is a certain convex function, and few more assumptions will be stated in the subsequent section. Observe that for ρ x ( ) = x 2 , ˆ β n coincides with the well-known least squares estimator (see, e.g., [12], pp. 42-43) based on the data ( λ , µ j j ) n j =1 , and when ρ x ( ) = | x | , ˆ β n coincides with the well-known least absolute deviation estimator (see, e.g., [12], pp. 4243) based on the data ( λ , µ j j ) n j =1 . Besides, it also includes Huber's estimator with ρ x ( ) = x 2 1 ( | x |≤ c ) / 2 + ( c x | | -c / 2 2)1 ( | x &gt;c | ) , where c &gt; 0 (see, e.g., [12], p. 43), L q regression estimate with ρ x ( ) = | x | q , where q ∈ [1 , 2] (see, e.g., [9]) and well-known α -th regression quantiles with ρ x ( ) = | x | +(2 α -1) x , where α ∈ (0 , 1) (see, e.g., [12], p. 43, equation (5.5)). This is one of the reasons to consider the technique of M -estimation, which includes many well-known estimators of the unknown regression parameter.

## 3 Large sample properties of ˆ βn

In order to study the large sample properties of ˆ β n defined in (2.7), one needs to assume the following technical conditions.

(A1) ρ . ( ) is a strictly convex function.

(A2) The error random variables 〈 v , /epsilon1v j j 〉 H ( j = 1 , . . . , n ) are i.i.d., and the common distribution of 〈 v , /epsilon1v j j 〉 H , which is denoted by F (depends on joint distribution of ( λ, µ )), is such that F ( D ) = 0. Here D is the set of discontinuity points of ρ ′ ( ), . where ρ ′ ( ) denotes the derivative . of ρ . ( ).

(A3) E ρ [ ′ ( 〈 v , /epsilon1v 1 1 〉 H + )] = c ac + ( o | c | ) as | c | → 0. Here a and c are two some constants.

(A4) g c ( ) := E ρ [ ′ ( 〈 v , /epsilon1v 1 1 〉 H + ) c -ρ ′ ( 〈 v , /epsilon1v 1 1 〉 H )] exists for some c ∈ ( -δ, δ ), where δ &gt; 0 is a some constant. Moreover, g . ( ) is a continuous function at 0.

$$( \mathbf A 5 ) \, E [ \{ \rho ^ { \prime } ( \langle v _ { 1 }, \epsilon v _ { 1 } \rangle _ { \mathcal { H } } ) \} ^ { 2 } ] \coloneqq D > 0.$$

$$\left ( A 6 \right ) S _ { n } \coloneqq \sum _ { j = 1 } ^ { n } \lambda _ { j } ^ { 2 } > 0 \text{ for some } n \geq n _ { 0 } \text{ and } \\ \lambda ^ { 2 }$$

$$d _ { n } ^ { 2 } \coloneqq \max _ { 1 \leq j \leq n } \frac { \lambda _ { j } ^ { 2 } } { S _ { n } } \to 0 \ a s \ n \to \infty.$$

Remark 3.1 Condition (A1) is required to have unique minima of the equation described in (2.7) , and it is a common assumption across the literature of M-estimation (see, e.g., [12]) in traditional Statistics. Condition (A2) indicates certain smoothness of the function ρ along with the fact that the error terms (in the forms of inner product) are i.i.d. random variables. Conditions (A3) and (A4) imply the existence of the derivative of the criterion function of M estimation in some neighbourhood of 0 , which is necessary to solve minimization problem involved in M estimation. Next, Condition (A5) imply that the M -estimator ˆ β (see (2.7) ), after appropriate normalization, will weakly converge to a non-degenerate random variable (see the statement of Theorem 3.2 below, and finally, Condition (A6) is required to have the asymptotic normality of ˆ β n , and from statistical point of view, one can interpret as the variation explained by any eigen value will not be a dominating factor.

Now, we state the consistency and the asymptotic normality results associated with ˆ . β n

Theorem 3.1 Under (A1)-(A6), ˆ β n p → β as n → ∞ . Here p → denotes the convergence in probability, and ˆ β n and β are the same as defined in (2.7) and (2.6) , respectively.

Theorem 3.2 Under (A1)-(A6),

$$T _ { n } ^ { - 1 / 2 } K _ { n } ( \hat { \beta } _ { n } - \beta ) \stackrel { w } { \rightarrow } Z,$$

where Z is a random variable associated with standard normal distribution, and w → denotes the convergence in distribution. Here T n = DS n and K n = aS n , where D a , and S n are the same as defined in (A5), (A3) and (A6), respectively.

Remark 3.2 In the set up of quantum probability, the observations ( λ , µ j j ) n j =1 are n eigen value pairs of ( X,Y ) unlike the classical probability, where the observations are ( X , Y j j ) n j =1 obtained from the distribution of ( X,Y ) . Despite this issue, Theorems 3.1 and 3.2 assert that ˆ β n is a consistent estimator of β , and after appropriate normalization, ( ˆ β n -β ) converges weakly to a standard normal random variable.

## 4 Concluding Remarks

In this work, observe that the true model (see (2.1)) is also linear with scalar valued coefficient, i.e., in other words, we work on specified model. However, the study will be entirely different if the model becomes mis-specified, which is of great interest in many statistical problem (see, e.g., [1] and a few relevant references therein). Investigating this issue may be of interest to future research. Besides, we have assumed our operators to be compact and self-adjoint. The classical analog would be to assume that the real random variables are discrete. It will be of interest when the involved operators are just self-adjoint. Then the spectrum may be continuous and in this case, our technique would not suffice. Moreover, in this work as the parameter β is scalar, the operators involved here are apriori simultaneously disgonalizable and therefore they can be measured simultaneously.

In the course of this study, it is assumed that the error random variables 〈 v , /epsilon1v j j 〉 H ( j = 1 , . . . , n ) involved in (2.6) are i.i.d., and the assertions in the main results (i.e., Theorems 3.1 and 3.2) rely on this assumption. However, this assumption can be weakened by following the approach considered in [13]. In this work, we did not dig in this issue as establishing the large sample properties of ˆ β n is not the main theme of the work, and therefore, we kept the assumption simpler so that the readers can appreciate the notion on quantum probability conveniently for operator on operator regression.

Acknowledgement: Subhra Sankar Dhar gratefully acknowledges his core research grant (CRG/2022/001489), Government of India. Suprio Bhar gratefully acknowledges the Matrics grant MTR/2021/000517 from the Science and Engineering Research Board (Department of Science &amp; Technology, Government of India). Soumalya Joardar acknowledges support from the SERB MATRICS grant (MTR 2022/000515).

## References

- [1] Bagchi, P. and Dhar, S. S. (2024). Characterization of the least squares estimator: misspecified multivariate isotonic regression model with dependent errors. Theory Probab. Math. Statist. , (110):143-158.
- [2] Bai, Z. D., Rao, C. R., and Wu, Y. (1992). M-estimation of multivariate linear regression parameters under a convex discrepancy function. Statistica Sinica , 2(1):237-254.
- [3] Dette, H. and Tang, J. (2024). Statistical inference for function-on-function linear regression. Bernoulli , 30(1):304-331.
- [4] Dhar, S. S., Jha, P., and Rakshit, P. (2022). The trimmed mean in non-parametric regression function estimation. Theory Probab. Math. Statist. , (107):133-158.
- [5] Dirac, P. A. M. (1945). On the analogy between classical and quantum mechanics. Rev. Modern Phys. , 17:195-199.

- [6] Jenˇ cov´, A. and Petz, D. (2006a). a Sufficiency in quantum statistical inference. Comm. Math. Phys. , 263(1):259-276.
- [7] Jenˇ cov´, A. and Petz, D. (2006b). a Sufficiency in quantum statistical inference. A survey with examples. Infin. Dimens. Anal. Quantum Probab. Relat. Top. , 9(3):331-351.
- [8] Knapp, A. W. (2005). Advanced real analysis . Cornerstones. Birkh¨user a Boston, Inc., Boston, MA. Along with a companion volume Basic real analysis .
- [9] Lai, P. Y. and Lee, S. M. S. (2005). An overview of asymptotic properties of lp regression under general classes of error distributions. Journal of the American Statistical Association , 100(470):446-458.
- [10] Shao, J. (2003). Mathematical statistics . Springer Texts in Statistics. Springer-Verlag, New York, second edition.
- [11] Sinha, K. B. (2022). Sufficient statistic and Rao-Blackwell theorem in quantum probability. Infin. Dimens. Anal. Quantum Probab. Relat. Top. , 25(4):Paper No. 2240005, 16.
- [12] van der Vaart, A. W. (1998). Asymptotic statistics , volume 3 of Cambridge Series in Statistical and Probabilistic Mathematics . Cambridge University Press, Cambridge.
- [13] Wu, W. B. (2007). M-estimation of linear models with dependent errors. The Annals of Statistics , 35(2):495 - 521.

## 5 Appendix : Technical Details

Proof of Theorem 3.1: The arguments in the proof are parallel to the proof of Theorem 2.2 in [2]. Without loss of generality, we assume that β = 0. It is now enough to show that

$$P [ | \hat { \beta } _ { n } | \geq a _ { n } ] \rightarrow 0$$

as n →∞ , where { a n } n ≥ 1 is such that a n →∞ as n →∞ . Observe that it follows from (3.11) in [2] that there exists a sequence a ′ n such that a ′ n →∞ as n →∞ and a ′ n ≤ a n for all n ∈ N .

Moreover, in view of it, we also have

$$\text{Moreover, in view of it, we also have} \\ \sup _ { | \beta | \leq a ^ { \prime } _ { n } } & \left | \sum _ { j = 1 } ^ { n } [ \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } - \beta \lambda _ { j } ) - \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } ) + \beta \lambda _ { j } \rho ^ { \prime } ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } ) ] - \frac { \beta ^ { 2 } K _ { n } } { 2 } \right | \stackrel { p } { \to } 0 \quad \ \ ( 5. 1 ) \\ \text{as } n \to \infty, \, \text{where } K _ { n } \text{ is the same as defined in the statement of Theorem 3.2.} \\ \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{-- } \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--}" \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--}\text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--} \text{--}$$

Further, when | β | = a ′ n for some n ∈ N , then observe that

$$\frac { \beta ^ { 2 } K _ { n } } { 2 } \geq \frac { a _ { n } ^ { ^ { \prime } 2 } \max _ { 1 \leq j \leq n } \lambda _ { j } } { 2 },$$

which follows from the definition of of K n defined in the statement of Theorem 3.2 and S n defined in (A6) before Remark 3.1. Moreover, observe that using (A4) on (5.1), we have

$$\sum _ { j = 1 } ^ { n } \lambda _ { j } \rho ^ { \prime } ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } ) = O _ { p } ( 1 ),$$

and hence,

Therefore, using (5.1), (5.2), (5.3) and (5.4), we have

$$\sup _ { | \beta | \leq a _ { n } ^ { \prime } } \left | \beta \sum _ { j = 1 } ^ { n } \lambda _ { j } \rho ^ { \prime } ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } ) \right | = O _ { p } ( a _ { n } ^ { \prime } ). \\, \, \text{using} \, ( 5. 1 ), ( 5. 2 ), ( 5. 3 ) \text{ and } ( 5. 4 ), \, \text{we have}$$

$$I \text{ nerefore, using } ( 0. 1 ), ( 0. 2 ), ( 0. 3 ) \text{ and } ( 0. 4 ), \text{ we have } \\ P \left ( \inf _ { | \beta | = d _ { n } ^ { \prime } } \sum _ { j = 1 } ^ { n } \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } - \beta \lambda _ { j } ) - \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } ) \leq 0 \right ) \to 0 \\ n \to \cap \text{ and next } \underset { n \infty } { \text{ in} } ( A 1 ) \text{ we have }$$

as n →∞ , and next, using (A1), we have as n →∞ .

$$n & \to \infty, \, \text{and next, using } ( A 1 ), \, \text{we have} \\ & P \left ( \inf _ { | \beta | \geq d _ { n } ^ { \prime } } \sum _ { j = 1 } ^ { n } \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } - \beta \lambda _ { j } ) - \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } ) \leq 0 \right ) \to 0 \\ n & \to \infty.$$

Finally, by the definition of ˆ β n in (2.7), we have

$$P ( | \hat { \beta } _ { n } | \geq a _ { n } ^ { ^ { \prime } } ) \to 0$$

as n →∞ . As a ′ n ≤ a n for all n ∈ N , the proof follows.

/square

Proof of Theorem 3.2: Without loss of generality, we here consider β = 0. First of all, using (A6) and in view of Lindeberg Feller CLT (see [12]), we have

$$\overset { \sim } { \sum _ { j = 1 } ^ { n } } \lambda _ { j } \rho ^ { \prime } ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } ) \coloneqq K _ { n } \widetilde { \beta } _ { n } \overset { w } { \rightarrow } Z, \\ 1 1$$

as n → ∞ , where Z is a random variable associated with standard normal distribution, K n is the same as defined in the statement of Theorem 3.2, and ' w → ' denotes weak convergence. Hence, in order to prove this theorem, it is enough to show that

$$\hat { \beta } _ { n } - \widetilde { \beta } _ { n } \xrightarrow { p } 0 & & ( 5. 9 ) \\ \intertext { are the same as defined in ( 5. 8 ) \text{ and (2.7), respectively.  Observe }$$

$$& \to \infty, \, \text{where} \, \rho _ { n } \, \text{and} \, \rho _ { n } \, \text{are} \, \text{same as unclein} \, \{ 0. 0 } \, \text{and} \, ( 2. i ), \, \text{respectively} \, \text{,} \, \text{overse} \\ & \, \text{using} \, ( 5. 1 ) \, \text{and} \, \text{the definition of} \, \widetilde { \beta } _ { n } \, \text{in} \, ( 5. 8 ), \, \text{we have} \\ & \quad \sum _ { j = 1 } ^ { n } [ \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } - \widetilde { \beta } _ { n } \lambda _ { j } ) - \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } ) ] + \frac { \beta \widetilde { \beta } _ { n } K _ { n } } { 2 } \, \stackrel { p } { \to } 0 \\ & \quad \to \infty, \, \text{and for some sequence} \, \{ b _ { n } \} _ { n \geq 1 } \, \text{such that} \, b _ { n } \to \infty \, \text{as} \, n \to \infty \, \text{and for some arbitrary}$$

as n →∞ , where ˜ β n and ˆ β n are the same as defined in (5.8) and (2.7), respectively. Observe that using (5.1) and the definition of β n in (5.8), we have as n →∞ , and for some sequence { b n } n ≥ 1 such that b n →∞ as n →∞ and for some arbitrary constant δ &gt; 0, we have

$$\text{constant } \dots > 0, \, \text{we have} \\ \sup _ { | \beta | \leq b _ { n } + \delta } \left | \sum _ { j = 1 } ^ { n } [ \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } - \beta \lambda _ { j } ) - \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } ) ] + \beta \widetilde { \beta } _ { n } K _ { n } - \frac { \beta ^ { 2 } K _ { n } } { 2 } \right | \rightarrow 0 \quad \ \ ( 5. 1 1 ) \\ \text{as } n \to \infty. \, \text{Now, it follows from } ( 5. 1 0 ) \, \text{and } ( 5. 1 1 ), \, \text{we have} \\ \quad \ \. \quad \ \.$$

$$& \text{as} \, n \to \infty. \, \text{Now, it follows from} \, \{ 0. 1 0 \} \text{ and} \, \{ 0. 1 1 \}, \, \text{we have} \\ & \sup _ { | \beta - \widetilde { \beta } _ { n } | = \delta } \left | \sum _ { j = 1 } ^ { n } [ \rho ( \langle v _ { j }, \epsilon v _ { j } ) _ { \mathcal { H } } - \beta \lambda _ { j } ) - \rho ( \langle v _ { j }, \epsilon v _ { j } \rangle _ { \mathcal { H } } - \widetilde { \beta } _ { n } \lambda _ { j } ) ] - \frac { ( \beta - \widetilde { \beta } _ { n } ) ^ { 2 } K _ { n } } { 2 } \right | \, \to 0 \quad \ \ ( 5. 1 2 ) \\ & \text{as} \, n \to \infty. \, \text{Moreover, observe that when} \, | \beta - \widetilde { \beta } _ { n } | = \delta, \, \text{we have}$$

as n →∞ . Moreover, observe that when | β -β n | = , we have δ

Therefore, using (5.12) and (5.13), we have

$$\text{er, observe that when } | \beta - \widetilde { \beta } _ { n } | = \delta, \, \text{we have} \\ ( \beta - \widetilde { \beta } _ { n } ) ^ { 2 } K _ { n } \geq \left ( \min _ { 1 \leq j \leq n } \lambda _ { j } \right ) \delta ^ { 2 }. \\. \text{12} \right ) \text{and} \left ( 5. 1 3 \right ), \, \text{we have}$$

$$P ( | \beta - \widetilde { \beta } _ { n } | > \delta ) \to 0 & & ( 5. 1 4 ) \\ \text{is an arbitrary constant. It completes the proof.} & & \Box$$

as n →∞ , where δ &gt; 0 is an arbitrary constant. It completes the proof. /square
Suprio Bhar Subhra Sankar Dhar Soumalya Joardar
Operator on Operator Regression in Quantum Probability
This article introduces operator on operator regression in quantum
probability. Here in the regression model, the response and the independent
variables are certain operator valued observables, and they are linearly
associated with unknown scalar coefficient (denoted by $\beta$), and the error
is a random operator. In the course of this study, we propose a quantum version
of a class of estimators (denoted by $M$ estimator) of $\beta$, and the large
sample behaviour of those quantum version of the estimators are derived, given
the fact that the true model is also linear and the samples are observed
eigenvalue pairs of the operator valued observables.