In this paper, a network model has been proposed to control dengue disease transmission considering host-vector dynamics in n patches. The control of mosquitoes is performed by SIT. In SIT, the male insects are sterilized in the laboratory and released into the environment to control the number of offsprings. The basic reproduction number has been computed. The existence and stability of various states have been discussed. The bifurcation diagram has been plotted to show the existence and stability regions of disease-free and endemic states for an isolated patch. The critical level of sterile male mosquitoes has been obtained for the control of disease. The basic reproduction number for n patch network model has been computed. It is evident from numerical simulations that SIT control in one patch may control the disease in the network having two/three patches with suitable coupling among them.
1.
Introduction
As pointed out in Claessens et al. (2014), global financial crises underline the importance of innovative modelling approaches to financial markets. The quantum mechanics approach suggests an alternative way to describe the unpredictable stock market behaviour (see e.g., Baaquie (2009)).
This paper is motivated by the quantum mechanics approach to the actuarial modelling and risk analysis as initiated in Tamturk and Utev (2018) and developed in Lefèvre et al. (2018), Tamturk and Utev (2019). We modify and adapt their methods to derive and analyze a quantum-type financial modelling. Specifically, the Schrödinger, Heisenberg, Feynman, Dirac approaches in quantum mechanics (see, e.g., Griffiths and Schroeter (2018); Parthasarathy (2012); Plenio (2002)) will be applied to the well-known option pricing models of Black and Scholes (1973) and Cox et al. (1979).
Moreover, motivated by the data analysis of the quantum reserve process proposed in Lefèvre et al. (2018), we model the financial data as eigenvalues of certain 1 or 2-step observable operators. These data then are analyzed to identify price jumps using supervised machine learning tools such as k-fold cross-validation techniques (see, e.g., Bishop (2006); Hastie et al. (2009); Wittek (2014)).
Option pricing models using quantum techniques discussed, for example, in Baaquie(2004, 2014); Bouchaud and Potters (2003); Haven (2002) are often based on the Schrödinger wave function with Hamiltonian operator H and are mainly oriented to continuous-time markets. For discrete-time markets as considered here, following Chen(2001, 2004), , we choose the discrete-time formalism and analyze the quantum version of the Cox-Ross-Rubinstein binomial model. Then, we establish the limit of the spectral measures providing the convergence to the geometric Brownian motion model. We also identify the limit of the N-step non-self adjoint bond market as a planar Brownian motion.
The paper is organized as follows. Section 2 deals with heterogeneous quantum binomial markets. In Section 3, two convergence properties to continuous-time quantum markets are obtained. In Section 4, discrete and continuous-time quantum mechanics techniques are applied to the problem of option pricing. Section 5 is devoted to an analysis of stock market data.
2.
Quantum modelling in finance
The motivation to our models come from the quantum mechanics approach to insurance proposed in Tamturk and Utev (2018, 2019), the insurance claim data analysis via quantum tools introduced in Lefèvre et al. (2018) and the non-traditional financial modelling initiated in Ma and Utev (2012) and developed in Karadeniz and Utev(2015, 2018). The quantum modelling approach is partly inspired by Baaquie (2004); Chen(2001, 2004), and Parthasarathy (2012).
2.1. Share price operators
We begin by outlining some standard arguments for a quantum-type modelling (see e.g., Lefèvre et al. (2018); Parthasarathy (2012)). An observable is a linear operator (matrix) on a certain Hilbert space. The quantum product of two independent observables A,B is implemented by the tensor product of the observables A⊗B. So, ln(A⊗B) acts as a quantum sum of two independent observables; in particular, A⊗A is the quantum product of two independent identical observables.
The 1-step quantum geometric random walk is defined as a 2×2 matrix A with eigenvalues eu,ed. Thus, A⊗N (N≥0) models the N-step geometric random walk.
2.1.1. Quantum binomial model
The quantum type modification of the classical Cox-Ross-Rubinstein model is originated in Chen(2001, 2004). The dynamics per period is defined by two moves: eu (the share price goes up) and ed (the share price goes down) with d<0<u.
The quantum binomial model over N periods is then represented by the share price operator
where the main 1-step observable H=A is a 2×2 self-adjoint (hermitian) matrix with eigenvalues eu,ed and representation
where U is a 2×2 unitary matrix.
In the sequel, the quantum binomial model discussed will be heterogeneous with share price operator
2.1.2. Quantum actuarial-type model
The motivation to this relatively new financial model comes from Lefèvre et al. (2018) quantum mechanics approach to non-life insurance and the Lamplighter group approximation to the financial modelling suggested in Ma and Utev (2012). Based on the approach to the financial data analysis developed in Ma and Utev (2012); Karadeniz and Utev(2015, 2018) (see also references therein), we treat the data as having big jumps, say eu, small jumps, say ed or no jump. Moreover, the financial stock is observed at fixed times Δk but the number of jumps occurred during the time period ((k−1)Δ,kΔ] is not observed. For simplicity and following Lefèvre et al. (2018), we then assume there are at most two jumps per period.
In this circumstance, the main observable operator is given, similarly to the insurance case Lefèvre et al. (2018), by
where In is a n×n identity matrix, A is the 2×2 matrix representing the 1-step geometric move operator (2) with eigenvalues eu,ed, and P0,P1,P2 are 3×3 matrix projection operators corresponding to the 0,1,2 claim occurrence operators and defined by
where Di+1|3 is a 3×3 diagonal matrix which has all its elements equal to 0 except the (i+1,i+1)-th with value 1, and W is a 3×3 unitary matrix.
The share price operator over N periods is then defined as
Furthermore, when constructing the density operator, we consider the following two cases.
Maxwell-Boltzmann statistics. In this case, the jump sizes are i.i.d. with a two-point distribution. More precisely, there are 0,1,2 jumps with probabilities δ0,δ1,δ2 given via a Poisson process, and each jump size has two possible values ed,eu with probabilities q,p (see also later in Section 4.1).
Bose-Einstein statistics. In this case, the claim sizes are dependent but remain independent of the claim occurrences.
2.1.3. Quantum trinomial model
This model makes a bridge between the traditional binomial model and the actuarial type model. In this case, the dynamics per period is defined by three moves: no change, down and up. The trinomial type financial modelling is a well established topic in finance (see e.g., Boyle (1986); Tian (1993); Leisen and Reimer (1996)).
Now, the 1-step observable operator H2=B is the 3×3 self-adjoint matrix with eigenvalues eu,1,ed corresponding to these three moves and is given by
where U is now a 3×3 unitary matrix.
Then, the share price operator over N periods is defined by
Remarks. Notice that the heterogeneous versions of the trinomial model and the actuarial type model are available with the representation similar to (3). However, they are not treated in this paper because the main purpose for considering these two models is the non-standard data analysis (to be presented in Section 5). Although the data analysis of time dependent models is a fascinating topic, it is out of the scope of this paper.
2.2. Quantum binomial price
As mentioned in above, the quantum binomial model (1) is originated in Chen(2001, 2004). Our presentation is somewhat different and based on the algebraic tensor product properties. In addition, since heterogeneity is an important topic in option pricing (see e.g., Benninga and Mayshar (2000)), we treat a slightly more general time dependent or heterogeneous quantum binomial model (3).
A quantum state ρ is defined as a positive self-adjoint operator with trace tr(ρ)=1. We recall the following properties.
Lemma 2.1. Let A,B,C,D be self-adjoint operators, U a unitary matrix, f a function of observables and ρ a quantum state. Then,
For simplicity, we choose a quantum state ρ as a tensor product, i.e.
where each ρi is a self-adjoint non-negative 2×2 matrix such that tr(ρi)=1. From Lemma 2.1, we see that ρ is a proper quantum state.
The risk-neutral world of the quantum Black-Scholes model consists of self-adjoint non-negative 2×2 matrices ρi that satisfy
where ri is the risk-free interest rate for the period i.
Let us define
Note that the ~ρi are also quantum states. In addition, ˜ρi has non-negative diagonal elements q(i)u,q(i)d and it can be shown to have the representation
where
The transformed operator ˜ρ is then obtained from Lemma 2.1 (third property) as
after using (12), and ˜ρ is again a quantum state. For the classical probability case, that is when all matrices are commutative, this transform is the form of the change of measure.
Moreover, the quantum no arbitrage condition (11) is satisfied for the transformed density, i.e.
2.3. Quantum binomial option pricing
Consider a payoff function f and a discount factor
From (3), the price of the general option without arbitrage OP(f(HSN)) for the N-step quantum binomial model is then defined by
Applying properties given in Lemma 2.1, we then obtain
after using (15). Since ˜ρ is a tensor product and f(S0D1⊗...⊗DN) is a diagonal matrix, we deduce the option pricing formula (20) below.
Theorem 2.2. For the heterogeneous quantum binomial model,
where the index σ denotes any feasible path yσ of the form
which occurs with the probability qσ=qσ1…qσN where qσi∈{qui,qdi}.
In particular, for the homogeneous case (the quantum Cox-Ross-Rubinstein model) where for all i, ri=r,ui=u,di=d with qui=qu,qdi=qd, the formula (20) reduces to
Non-self adjoint quantum binomial market. Assume that Hi are invertible 2×2 matrices with two different eigenvalues ed and eu but no more self-adjoint, in general. In this case the transformed density ˜ρi are no more proper states, in general. However, they still have same diagonal elements qui and qdi and so the transformed matrix ˜ρ defined in (15) again has same diagonal elements as in the self-adjoint case, but is not a proper state, in general. By inspecting the proof, we see that Theorem 2.2 holds true in this case as well.
3.
Convergence to continuous-time markets
In this section, we consider the homogeneous quantum binomial model, and we discuss two examples on the convergence to continuous-time markets, namely the Black-Scholes model and the planar Brownian motion.
3.1. Convergence to the Black-Scholes model
Theorem 3.1. Let μN be the measure of the eigenvalues of H⊗N with respect to the quantum state ρ⊗N. Suppose that r=λ/N and u=−d=σN−1/2. Then, as N→∞,
where μ is a lognormal distribution (i.e. μ(x)=P(ea+σZ≤x) for suitable constants a,σ and Z a standard normal variable).
Proof. We begin with the representation via the spectral measure μN. Observe that
where μN is the measure of the eigenvalues λσ of H⊗N with respect to the quantum probability tr(ρ⊗NH⊗N). Since λσ=eσ1…eσN with σi∈{ui,di}, we get
where qσi∈{qu,qd} respectively.
Let us move on to the weak convergence desired. We can write that
where the Yi are i.i.d. variables with P(Yi=eu)=qu,P(Yi=ed)=qd. Thus,
with TN≡lnY1+…+lnYN. As N→∞, we obtain from the central limit theorem that
for certain constants a,b and Z a standard normal. This gives the limit result (23). ⋄
Remarks. (ⅰ) The limiting measure μ corresponds to that of the geometric Brownian motion St=S0e(ρ−σ2/2)t+σBt when S0=1, t=1 and ρ−σ2/2=a.
(ⅱ) It would be interesting to compare the technique with the semi-classical approximation, such as expanding the action around the classical path (see (48) in Contreras et al. (2010)). Another interesting question is to analyse connection with arbitrage as discussed in Haven (2002) and Contreras et al. (2010).
(ⅲ) The Cauchy transform is an alternative approach to deal with quantum probabilities (see Mudakkar and Utev (2013)). In particular,the convergence of spectral measures is reduced to the convergence of the Cauchy transforms. We recall that the Cauchy transform for the measure μ is defined by
for all z∈C∖R={z∈C:ℑz≠0} and open intervals (a,b) with μ({a,b})=0. The goal is then to show that Sμn(z)→Sμ(z). However,the common approach of the moment expansion does not work in this case since
3.2. Convergence to the planar Brownian motion
In this part,by treating the relatively simple bond market,we show that non-commutative markets provide a richer class of financial models (see also Haven (2002); Contreras et al. (2010); Herscovich (2016)).
Bond market. Assume that returns are non-risky,that is the outcomes are equal. In the quantum setup,we suppose that matrix H has a single eigenvalue,eu say.
Self-adjoint quantum bond market. Let ρ⊗N denote a quantum state defined as before. Assume that in addition operator H is self-adjoint. Then H=euI2 and
In particular,any option claim f(HSN) is commutative with the state ρ⊗N which implies that the self-adjoint quantum bond market is commutative,that is equivalent to the simple classical probability financial market with two non-risky assets. The no-arbitrage condition becomes global to give eu=tr(ρH)=1+r since tr(ρ)=1 which restates that under no-arbitrage the non-risky returns are equal. Moreover,the option price for f(HSN) is provided by
Non self-adjoint quantum bond market. Now,assume that H is no longer self-adjoint,for example a product of two noncommutative self adjoint observables. Thus,non self-adjoint quantum bond market is noncommutative,in general. In this case,we show that the limit is sensitive to the density state ρ and the representation of H. For simplicity,we assume that the basic observable H and the state ρ are defined by
i.e.,Hu is now a 3×3 Jordan matrix. The bond price process is then defined by
To satisfy the no-arbitrage condition,we ask that
Theorem 3.2. Suppose that r=λ/N,u=a/N and δ=−Δ/N with Δ≥0. Then,for any positive integer k,as N→∞,
regardless of a, where i is the imaginary unit and BΔ is a Brownian motion at time Δ.
Proof. Notice that
From the definition of HSN and ρ, we then get
after using the assumptions made on r,u,δ. This provides the limit result (36). ⋄
Remarks. (ⅰ) This representation is useful in computing the European call option OP([HSN−K)+]) via the Fourier techniques combined with the Monte Carlo approximation. For the distribution of this process, we view R2 as the complex plane and the planar Brownian motion ~Bt=(B1(t),B2(t)) is then interpreted as a complex-valued Brownian motion.
(ⅱ) Although observables are traditionally considered to be self-adjoint, the non self-adjoint data is modelled by considering an observable such as H=AB, i.e. the product of two self-adjoint matrices A,B where A and B may represent the non-trade and trade time changes, respectively.
(ⅲ) Similarly to the convergence to the Black-Scholes model, the limit does not depend on the shift parameter a. However, now the limit depends on the mysterious characteristic Δ.
4.
Quantum mechanics in finance
In this section, we will apply some methods of quantum mechanics to finance. Discrete and continuous-time markets must be treated separately because of different stochastic behaviors. To simplify the presentation, we assume that the interest rate is 0 and that the risky processes for share prices are martingales.
4.1. Discrete-time quantum approach
The use of Dirac-Feynman quantum mechanics techniques for insurance risk modelling was initiated in Tamturk and Utev (2018) and then developed in Lefèvre et al. (2018); Tamturk and Utev (2019). We adapt this approach to the problem of option pricing in finance, in particular for the pricing of path-dependent options.
In the Dirac formalism, bra-ket notation is a standard way of describing quantum states. Consider a class of n×n matrices treated as C∗ algebra. A column vector x corresponds to a ket-vector |x>. An associated bra-vector <x| is a row vector defined as its Hermitian conjugate. The usual inner product is denoted by <x|y>, while the outer product |x><y| is the operator/matrix defined by
In the Feynman path integral methods, the transition probability P(xj→xj+1) is computed as the propagator <xj|Aj+1|xj+1> when Aj is a Markovian operator. Thus, the typical path is written as |x0>→|x1>→…→|xn>, and its probability is given by
The main ingredient is then the path calculation formula that calculates the probability P(x0→xn) via the sum of the probabilities on all the appropriate paths, i.e.
It remains to define a suitable propagator for discrete time. For simplicity, we take the operators Aj all equal to A and the time intervals Δti all equal to Δt. In a similar way to e.g. Baaquie (2004); Tamturk and Utev (2018), we assume that the operator A is defined via an Hamiltonian operator H such as A=e−ΔtH, where −H is a Markovian generator called Markovian Hamiltonian. Thus, P(xi→xi+1)=<xi|e−ΔtH|xi+1> which can be computed applying the Fourier transform to the momentum space (e.g., Griffiths and Schroeter (2018); Tamturk and Utev (2018)). Specifically, let |p> be a basis in that space, and write <x|p>=eipx and <p|x>=e−ipx. Then, we get
where {|α>,Kα} is the set of eigenvectors and eigenvalues of the Hamiltonian operator H (i.e., H|α>=Kα|α>). Therefore, we deduce from (41), (42) that
For a more detailed overview of the theory, we refer the reader to the books by Feynman and Hibbs (2010); Griffiths and Schroeter (2018); Parthasarathy (2012) and Plenio (2002).
Option pricing formula. Consider a claim of the form C=f(S0,S1,...,SN), using the notation of Section 2. From (43), we obtain for the corresponding option price
The discrete-time approach followed to derive (43) can then be easily applied to the current formula (44).
A modified Cox-Ross-Rubinstein model. Consider a discrete-time market in which, during each i-th time interval Δt, the share price Si can
(.) have 1 jump giving Si+1=Sieu,Sied with probabilities p,q=1−p, or
(.) have 2 jumps giving Si+1=Sieu+d,Sie2u,Sie2d with probabilities 2pq,p2,q2, or
(.) remain the same giving Si+1=Si,
where d,u are integers with d<0<u. Furthermore, the possible jumps arrive according to a Poisson process of parameter λ so that
(.) δ1≡P[N(Δt)=1]=e−λΔt(λΔt),
(.) δ2≡P[N(Δt)=2]=e−λΔt(λΔt)2/2,
(.) δ0≡P[N(Δt)=0]≈1−δ1−δ2.
Define xi=ln(Si). The transition probabilities P(Si→Si+1) are equivalent to P(xi→xi+1). Set Δt=1, say. We can now apply (42) with A=e−H where the set {|α>,Kα} is found by solving the Schrödinger equation e−H|α>=e−Kα|α> in which
Notice that the Hamiltonian H is not Markovian, but (42) is still applicable to the discrete times kΔt. Moreover, by construction, the martingale probabilities p=qu and q=qd for the Cox-Ross-Rubinstein model (without interest) yield a martingale in the present situation too since we have
We have illustrated numerically the option pricing results obtained for the model. The tables and figures are however too large to be included here.
4.2. Continuous-time quantum approach
This short part is mostly a review of the application of quantum mechanics approach to the continuous time markets and closely follows Baaquie (2004).
The continuous time formalism is based on the Fourier transform of tempered distributions on the basis |p> in the momentum space. Consider a risky asset price St that evolves in function of an Hamiltonian operator H. First, the method is applied to compute the pricing kernel
Then, the option price at time t for the claim Q≡Q(ST), T>t, given St=x is defined by
Black-Scholes model. In this classical approach, the stock price St follows a geometric Brownian motion, i.e. St=S0e˜Bt where ˜Bt=μt+σBt (μ is the drift, σ the volatility and Bt a standard Brownian motion).
We apply the quantum mechanics approach. Motivated by Baaquie (2004), we work with ˜Bt rather than with St. The corresponding Hamiltonian for the Brownian motion is Hf=−(1/2)f″ (computed, for example, via the Itô formula). The Brownian motion kernel is then defined from (47) by the normal density function
To obtain the option price (48), we set x = S_0e^u and Q(x) = Q(S_0e^u)\equiv g(u) . Then, the Feynman-Kac formula yields
Via path integrals. To find the pricing kernel p(x, \tau; x') = < x|e^{-\tau H}|x' > for \tau = T-t , an alternative method consist in using path integral methods. For this, we discretize the time in N intervals of length \Delta and consider the x_i = x(t_i) where t_i = i \Delta . We then proceed as in the situation in discrete-time. The pricing kernel for (x, x') = (x_0, x_N) becomes
Applications to non-life insurance. Consider Feynman's modification \tilde{H} of the Brownian motion Hamiltonian H by adding the potential V , i.e. for the Hamiltonian \tilde{H} = H+V . For example, choose a ruin level B and take V(x) = +\infty for x < B . Then, the path calculation formula (43) allows us to compute ruin probabilities when B = 0 (Tamturk and Utev (2018, 2019)) and exotic options with barriers when B > 0 .
Finally, let us mention that numerically, the binomial model formula (20) and the path integral approach (51) were found to give results close to the Black-Scholes formula, even for relatively small values of N (of order 40 ).
5.
Analysis of stock market data
We are going to analyse a generated set of financial data by choosing two different quantum models described in Section 2.1. First, we consider the quantum actuarial-type model (see Section 2.1.3 and in 4.1) with the N -step observable defined in (6). Then, we consider the quantum trinomial model (see Section 2.1.2) with the N -step observable defined in (8). Note that matching 1 -step for the actuarial case to two steps for the trinomial model is natural, since in the actuarial case, we choose at most two jumps.
5.1. Methodology for data analysis
We apply supervised machine learning methods such as developed e.g. in the books by Bishop (2006); Hastie et al. (2009) and Wittek (2014).
Overall approach. The dataset is supposed to come from a non-ordered class of randomly perturbed observables. Each data is the observation of an eigenvalue \lambda of the observable perturbed by i.i.d. error terms. The observable is the 1 -step operator for the actuarial-type model and the 2 -step oparator for the trinomial model. Matching 1 -step in the actuarial case to 2 -step in the trinomial case is natural since there are at most two jumps in the actuarial case considered.
First, the data are classified in classes G_\lambda with respect to the eigenvalues \lambda . Then, the probabilities p_\lambda are estimated by maximum likelihood using Maxwell-Boltzmann or Bose-Einstein statistics. Finally, the \lambda are estimated via the weighted L_1 -norm risk error function.
Let us explain in more detail for the actuarial-type model, for example.
Step 1. An initial (u = u_0, u = d_0) is chosen randomly.
Step 2. For a given (u, d) , the data is classified and labeled against the eigenvalues of the observable using a nearest neighbor algorithm. This leads to the classes G_{\lambda} .
Step 3. For the same (u, d) , the estimates \hat{p} and \hat{q} are obtained by maximizing the likelihood function L(p, q).
Step 4. The (u, d) is updated by minimizing a weighted L_1 -norm risk error function F(u, d) \equiv F(\lambda) defined by
Step 5. The loop of steps 2 to 4 is repeated until the relative error becomes smaller than a selected difference M , i.e. when
k- fold cross-validation. To reduce the risk of error, we use a k- fold cross-validation strategy. The dataset is randomly divided into k subsets of equal size. One of the subsets is chosen as the training set and the others as test sets. The process is repeated k times, each subset constituting a training element. At each iteration, steps 1 to 5 above are applied to the training data and the results obtained are then checked in the test data. Finally, the estimates used are an average of those obtained on the k iterations.
Numerical example. As a simple illustration, we will consider the following dataset
5.2. Data analysis via the actuarial-type model
From the assumptions of the model, the 1 -step observable H has the eigenvalues
with probabilities respectively given by
For the probabilities p_{d+u}, q_{2d}, p_{2u} , we consider two possible statistics often used in the analysis of quantum observables.
Maxwell-Boltzmann independence. This case yields the binomial model since
Bose-Einstein dependence. In this case, corresponds to probabilities
As pointed out in Lefèvre et al. (2018), both statistics admit a formal construction, via the proper choice of the density operator \rho for the number of occurrences and the density projection operators \rho_\lambda , such that
Likelihood functions. Denote by \#x the number of x observed in the data set. For the Maxwell-Boltzmann statistics, the likelihood is defined by the probabilities p , q , \delta_0 , \delta_1 and \delta_2 such that
For the Bose-Einstein statistics, the likelihood function is modified as
Risk functions. From (55), the scaled share price spectrum of H_{S_1} is given by S_0\{\lambda\} = \{S_0, S_0e^{u}, S_0 e^{d}, S_0e^{u+d}, S_0e^{2d}, S_0e^{2u}\} . Thus, using the risk function (52), we get for the Maxwell-Boltzmann case
and for the Bose-Einstein case,
Numerical illustration. The data (54) is treated as the eigenvalues of the 1 -step observable H defined in (6) for N = 1 . Let us assume that the Poisson process rate is 1 and the length of time is \Delta t = 1 .
We estimate the values d, u and the probabilities q, p by applying the algorithm of Section 5.1 and using the Maxwell-Boltzmann statistics. Choose, for instance, S_0 = 70 , (u_0, d_0) = (0.8, -0.1) and M = 0.000001 . The results obtained from (60), (62) are presented in Table 1.
We then also apply a 4 -foldcross-validation procedure. Let V_i , 1\leq i \leq 4 , before randomly chosen subsets of V such that \cup_{i = 1}^{4} V_i = V . For each iteration, \{V\setminus V_i\} and V_i are treated as the training set and the test set, respectively. First, the maximum likelihood estimation and the risk error computation are executed for the training data. Then, the obtained estimates are implemented in the test data. This gives the results of Table 2.
We observe that the risk errors are small for the test data but are larger for the training data. Thus, the total risk errors are significant. Similar numerical calculations have also been performed under the Bose-Einstein assumption.
5.3. Data analysis via the trinomial model
The spectrum of H^{\otimes N}_2 is given by
Thus, for the case N = 2 , the set of observables is exactly the same as (55). However, the associated probabilities differ from (56) and are equal to
Likelihood and risk functions. Using the Maxwell-Boltzmann statistics, the likelihood is defined by the probabilities p_1 , p_2 , p_3 as
Note that for the Cox-Ross-Rubinstein model, p_2 = 0 so that the likelihood is simplified to L(p_1, p_3) = (2p_1p_3) ^{\#d+u}(p_3^2) ^{\#2d}(p_1^2) ^{\#2u} .
The corresponding risk function is then defined by
Numerical illustration. We process the data (54) again with S_0 = 70 , (u_0, d_0) = (0.8, -0.1) and M = 0.000001 . The results obtained using the algorithm of Section 5.1 with the functions (66), (67) are given in Table 3.
Then, we apply a 4 -fold cross-validation procedure as previously done. The obtained results are shown in Table 4. Note that this method allows to reduce the risk errors.
The observable operators of the quantum actuarial-type and trinomial models give us different results as expected. Which model to choose? One possible approach might be to consider a mixture of quantum models via the mixture of Hamiltonians (see, e.g., Wittek (2014)).
6.
Conclusion
Several quantum type financial models are constructed that benefit from the physical interpretation of the unpredictable stock market behaviour and associated dependences. The models provide a general physical type framework for pricing of derivatives and a possibility to construct quantum trading strategies. Moreover, it is revealed that certain quantum type models are applied both in actuarial and financial sciences.
Acknowledgements
C. Lefèvre received support from the Chair DIALog sponsored by CNP Assurances. S. Utev benefited from a Scientific Mission granted by the Belgian FNRS. We thank the referees for the careful reading of the paper which greatly helped to improve the presentation.
Conflict of interest
All authors declare no conflicts of interest in this paper.