# Estimation of Jointly Normally Distributed Demand for Cross-Selling Items in Inventory Systems with Lost Sales

Authors: Ren-Qian Zhang; Qi-Qi Wang; Yi-Ye Zhang; Hai-Tao Zheng; Jie Hu
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

Demand estimation is often confronted with incomplete information of censored demand because of lost sales. Many estimators have been proposed to deal with lost sales when estimating the parameters of demand distribution. This study introduces the cross-selling effect into estimations, where two items are cross-sold because of the positive externality in a newsvendor-type inventory system. We propose an approach to estimate the parameters of a jointly normally distributed demand for two cross-selling items based on an iterative framework considering lost sales. Computational results based on more than two million numerical examples show that our estimator achieves high precision. Compared with the point estimations without lost sales, all the relative errors of the estimations of demand expectation, standard deviation, and correlation coefficient are no larger than 2% on average if the sample size is no smaller than 800. In particular, for demand expectation, the error is smaller than 1% if the comprehensive censoring level is no larger than four standard deviations (implying a2σ-level of safety stock for each item), even if the sample size decreases to 50. This implies that the demand estimator should be competent in modern inventory systems that are rich in data.


---

## Body

## 1. Introduction

In newsvendor-type inventory models, it is frequently assumed that the demand distribution for items is known or can be estimated based on sales data [1, 2]. Estimating the demand directly based on sales data implies that the sales quantity must be exactly equal to the realistic demand amount. However, this is not the case in practice [3]. The order quantity of an item is finite; thus, the inventory is limited but demand is often uncertain. It is very difficult to meet an uncertain demand based on a limited inventory. As a result, stock-outs emerge. Studies show that more than 85% of customers choose to give up the products when the items are out of stock [4], and therefore, the sales amount is the censored information of demand. When sales data is considered, demand might be underestimated. To remedy this flaw, demand estimation methods for censored sales data have been proposed [5, 6].Another factor that should be considered in inventory management is the presence of externalities, including positive and negative ones [7]. As a type of positive externalities, cross-selling has not been considered in demand estimation. Cross-selling means that “a customer who has purchased a particular product may also be willing to purchase a related product” [8], and this implies that items are possibly purchased together owing to their unknown interior associations [9].This study considers the cross-selling effect between two items to establish a new method to estimate the jointly (bivariate) normally distributed demand for two cross-selling items based on demand information from sales data that is censored by inventory order quantities. In our study, the sales data of any one item are not only censored because of lost sales caused by its own stock-out, but also affected by another associated item because of cross-selling. The core task is to recover the original demand distributions from the incomplete data that is a mixture of censored demand and cross sales. The proposed estimation is entirely based on the sales data, which is easily implemented in a modern business scenario. To the best of our knowledge, this is the first time that the demand estimation for cross-selling items with unobservable lost sales is proposed. Our work will consolidate the inventory model with positive externalities, as in the models of Netessine and Zhang [7] and Zhang et al. [10], and it can be used to provide basic information for inventory decision making.The remainder of the paper is organized as follows. Section2 reviews the related work. Section 3 includes assumptions and notations. Section 4 presents the demand estimation method and provides estimators. Section 5 illustrates the method by a numerical example, and Section 6 conducts numerical experiments to evaluate the performance of the estimators based on numerical examples. The paper is concluded in Section 7.


## 2. Literature Review

Many studies have focused on single-variable demand estimation considering censored data based on parametric methods where the type of random distribution is given. For censored data of a normal demand, Fisher [11] developed estimators for the mean and standard deviation by maximizing the likelihood function, where two specifically predetermined tables are needed. Cohen [12] introduced new formulae to improve Fisher's method, based on which the special tables were excluded. Later, Gupta [13] investigated the problem and presented a better linear unbiased estimator. Nahmias [5] developed a parametric methodology for estimating the normal demand and validated the approach based on censored samples. The method was further extended to the case of a negative binomial distribution [6]. Berk et al. [14] estimated the demand with censored data of negative binomial, gamma, Poisson, and normal distributions by Bayesian updating, adopting an approximate posterior distribution that matches the first two moments of the exact posterior. Heeseab [15] provided a Bayesian method for demand estimation considering unobserved lost sales, which depends on the prior information to estimate the posterior distribution. Jain et al. [16] modeled the cumulative demand arrivals during a period as a stochastic process with an unknown parameter. The parameter is updated in a Bayesian fashion by the observation of sales and stock-out events, which is further simplified by introducing stock-out timing observation.Without a given type of demand distribution, nonparametric demand estimation methods are adopted. Beutel and Minner [17] proposed demand forecasting by solving a linear programming (LP) problem with fully observed demand, which does not need to know the random distribution. The method was further extended to deal with the censored observations of demand by Sachs and Minner [18]. They both assume that the demand is linearly dependent on multiple exogenous variables. For the nonlinear case, an intelligent forecasting algorithm based on artificial neural network (ANN) and conventional regression has been adopted. The method was applied for forecasting gasoline demand [19].The above estimators are focused on the case of a single stochastic variable and did not consider consumption externalities. For multiple variables, Dempster et al. [20] proposed the expectation-maximization (EM) algorithm for the case of incomplete data, which, in contrast to the previous approaches, is an iterative method. Dahiya and Korwar [21] provided a parameter estimation method for the likelihood equation in the bivariate normal case. They assumed that the missing data were from only one variable and no consumption externalities existed. Adamids and Loukas [22] considered the missing values and proposed bivariate Poisson estimations for a two-item inventory based on the EM method without externalities. Regarding the time-varying demand, multivariate auto regressive integrated moving average (ARIMA) models were used to forecast the demand for perishable goods [23], but consumption externalities were not considered either.Consumption externality means that the demand/consumption for one item affects the demand/consumption for other items, and there are two types of externality in inventory models: negative and positive [7]. Besanko et al. [24] presented a logit demand estimation based on the consumer choice model of substitutable items with negative externalities, where the endogenous price is the Nash equilibrium of wholesale and retail in a market. Anupindi et al. [25] and Smith and Agrawal [26] developed demand estimators considering stock-outs and multi-item substitutions, assuming Poisson arrivals of demand. Kök and Fisher [27] presented an EM-based method for estimating the parameters of substitution behavior and demand for products, which was used for optimizing the assortment planning of a retailer. Conlon and Mortimer [28] estimated the multinomial distribution demand under incomplete product availability by maximizing a likelihood function based on an EM algorithm. They assumed that the customer behavior can be described based on the discrete choice model of substitutable items. For multi-store multi-product substitution, Wan et al. [29] built demand estimation models considering that consumers may substitute a product for the stocked out item in the same store or switching to a neighboring store, where the parameters were estimated by the Markov chain Monte Carlo algorithm in a Bayesian manner.It is noticeable that the extant inventory models and demand estimations with externalities are mainly focused on substitutable items, where choice happens because of substitution, and the customer behavior can be modeled by the consumer choice theory. However, this is not the case for cross-selling scenarios, where positive externalities exist between complementary items and no choice happens even if some items are stocked out. For comparing our study with previous works, Table1 classifies the above-reviewed related literature into three groups.Table 1
Summary of demand-estimation related literature.
Factors considered by estimation method Literature (1) Single item with stock-outs Fisher [11], Cohen [12], Gupta [13], Agrawal and Smith [6], Berk et al. [14], Azadeh et al. [19], Heeseab [15], Beutel and Minner [17], Sachs and Minner [18], Jain et al. [16]  (2) Multiple items without externalities Dempster et al. [20], Dahiya and Korwar [21], Adamids and Loukas [22], Huber et al. [23]  (3) Multiple items with externalities of substitution Anupindi et al. [25], Besanko et al. [24], Smith and Agrawal [26], Kök and Fisher [27], Conlon and Mortimer [28], Wan et al. [29]The related literature in Table1 mainly includes three streams in the prior research: (1) estimation of the demand distribution of a single item considering censored data because of stock-out; (2) estimation of the jointly distributed demand for multiple items without considering externalities; and (3) estimation of the demand considering negative externalities because of substitution. We can observe that the extant research did not simultaneously consider the unobservable lost sales caused by stock-out and the positive externalities of cross-selling. Our study proposes a new method to estimate the jointly normally distributed demand for cross-selling items (because of positive externalities) based on censored sales data with unobservable lost sales.


## 3. Assumptions and Notations

This work is focused on inventory systems that are controlled based on the single-stage newsvendor strategy for multiple periods, but each period is independent. It means that neither the surplus product nor backordering is carried over from the former period to the next one. We focus on the cross-selling effect between two items, and thus the two-item inventory system is considered. Assumption1 summarizes the above considerations.Assumption 1.
The two-item inventory system is operated based on the single-stage newsvendor policy for each product, and there are no carryovers (surplus or backordered items) between adjacent periods.Consider a two-item inventory system, where itemi and item j affect each other in sales. In other words, besides the original “naked” demands for the two items (denoted by Di and Dj, respectively), the two items can also be demanded accompanying the sales of the counterpart item. For example, if a customer buys one unit of item i, she/he will possibly buy a certain quantity of the item j because of the cross-selling effect with item i, and vice versa. Let μi,σi and μj,σj denote the mean and standard deviation of the demand for item i and item j, respectively, and let ρ=cov⁡(Di,Dj)/σjσi represent the correlation coefficient. It is noted that externalities do not imply statistical correlations between the stochastic distributions of demand.It could be considered that the cross-selling demand caused by the original “naked” demand might lead to the secondary “new” additional demand, which will make the problem too complex to solve. However, in practice, vendors mainly pay close attention to the cross-sold items caused by the original “naked” demand. The secondary “new” additional part is trivial and therefore is overlooked. Following this consideration, the following Assumption2 holds.Assumption 2.
Cross-selling demand is just caused by the original “naked” demand, and the cross-selling demand will not lead to a secondary “new” additional demand.To quantitatively measure the cross-selling effect, we adopt the concept of association rules where cross-selling coefficients can be determined by data mining based on a transaction database [9, 30]. The cross-selling coefficient is defined as the amount of the cross-selling demand of one item caused by one unit another item's sale in average and vice versa. For clarity, we introduce Assumption 3.Assumption 3.
Cross-selling coefficients are known and can be determined by mining of association rules, based on which the demand for an item has a linearly proportional increase by cross-selling.Suppose that the sale of one unit of itemi (of the “naked” demand) leads to the sale of rij units of item j, and the sale of one unit of item j (of the “naked” demand) leads to the sale of rji units of item i, in average. From Assumption 3, we can see that if there are no stock-outs, the total demand for the two items is Xi=Di+rji∙Dj and Xj=Dj+rij∙Di, respectively.Moreover, regarding the ordering policy of the inventory system, the following Assumption4 is relevant in practice.Assumption 4.
For each item, the order quantity is known and is constant for all periods.LetQi and Qj represent the order quantities of items i and j, respectively. Denote xit,xjt the corresponding observed demands for item i and item j in period t, which are the realized sales that can be derived from transaction records. Evidently, the realized sales are bounded by order quantities, so that we have xit=min⁡(Qi,Xit)≤Qi and xjt=min⁡(Qj,Xjt)≤Qj.Note that ifxit=Qi or xjt=Qj, it could be considered that some sales are lost and the recorded sales amount is not equal to the actual demand. In this case, the sampled sales data may be censored as in Table 2.Table 2
Sampled sales data being censored.
Periodt 1 2 3 4 5 6 7 8 9 10 … t … T  Sales of itemi x i 1 x i 2 Q i Q i x i 5 x i 6 x i 7 x i 8 Q i x i 10 … x i t … x i T Sales of itemj Q j x j 2 x j 3 Q j x j 5 x j 6 Q j Q j x j 9 Q j … x j t … x j TIt is noticeable that even if the original “naked” demand is smaller than the order quantity, this does not mean that there are no lost sales. Because of cross-selling, the sale of one item is influenced by another item's sales and order quantity. Suppose that the realized “naked” demand of itemsi and j in period 3 (see Table 2) is ξi3 and ξj3, respectively. Even ξi3≤Qi, it is possible that xi3=ξi3+rji∙min⁡Qj,ξj3>Qi, resulting in lost sales of item i because of the additional demand caused by cross-selling. As mentioned above, the cross-selling coefficients rji and rij have been determined by association rule mining, so we need to estimate the distribution parameters: μi, μj, σi, σj, and ρ, of the original demand Di and Dj. We summarize all notations in Table 3.Table 3
Parameters and notations.
Parameters/notations Meaning i , j Indices of items D i , D j Unknown original “naked” demand for itemsi and j ξ i t , ξ j t The realized original “naked” demand of itemsi and j x i t , x j t Observed demand, i.e., the realized quantity of sales of itemsi and j μ i , μ j Means of demand (demand expectations) for itemsi and j to be estimated σ i , σ j Standard deviation of demand for itemsi and j to be estimated ρ Correlation coefficient of the original demand to be estimated Q i, Qj Order quantities of itemsi and j r i j , r j i Cross-selling coefficients x - i , x - j Means of the sampled sales data of itemsi and j μ ~ i , μ ~ j Estimators ofμi and μj σ ~ i , σ ~ j Estimators ofσi and σj ρ ~ Estimator ofρ z / z Quantile of standard normal distribution T Total number of sampled data (sample size)


## 4. Estimation Methodology

### 4.1. Estimators for One Item with Another Demand Given

Sorting increasingly all the sampled sales data, we will have two ordered data series as(1)Itemi:xi1,xi2,…,xiu,xi,u+1,…,xiT,(2)Itemj:xj1,xj2,…,xjv,xj,v+1,…,xjT,where xiu<Qi;xi,u+1=xi,u+2=…=xiT=Qi and xjv<Qj;xj,v+1=xj,v+2=…=xjT=Qj. u and v are the total number of xit and xjt that meet xit<Qi and xjt<Qj, respectively.Define(3)x-u=1u∑τ=1uxiτ,(4)x-v=1v∑τ=1vxjτ,(5)Su2=1u-1∑τ=1uxiτ-x-u2,(6)Sv2=1v-1∑τ=1vxjτ-x-v2,(7)Pi=uT,(8)Pj=vT,where x-u and x-v represent the mean values of the observed demand that are smaller than the order quantities, and Su2,Sv2 and Pi,Pj are the corresponding sample variances and proportions, respectively. Regarding the distribution estimation of the original demand, we have the following results:Proposition 5.
Given demand expectationμj and standard deviation σj of item j, and with a known correlation coefficient ρ, the value of μi and σi of item i can be estimated by Eq. (9) and Eq. (10), respectively.(9)σ~i=ρ2-1rji2σj2+Su21-ziϕzi/Pi-ϕzi2/Pi2-ρrjiσj,(10)μ~i=x-u+ϕziσ~i2+rji2σj2+2ρrjiσ~iσjPi-rjiμj,where ϕ(∙) and Φ(∙) are the density function and cumulative probability function of the standard normal distribution, respectively, and zi=Φ-1(Pi) is the quantile of normal distribution corresponding to probability Pi.Proof.
See AppendixAProposition 6.
Givenμi, σi, and ρ, the value of μj and σj can be estimated by Eq. (11) and Eq. (12), respectively. (11)σ~j=ρ2-1rij2σi2+Sv21-zjϕzj/Pj-ϕzj2/Pj2-ρrijσi,(12)μ~j=x-v+ϕzjσ~j2+rij2σi2+2ρrijσiσ~jPj-rijμi,where zj=Φ-1(Pj).Proof.
It is straightforward from Proposition5, because of the symmetry of the two items.Evidently, the second terms at the right-hand sides of Eq. (9) and Eq. (11), -ρrjiσj and -ρrijσi, are both negative values for calculating the standard deviations of items i and j. Not only that, the first terms under the radical sign, (ρ2-1)rji2σj2 and (ρ2-1)rij2σi2, are also negative. This implies that if we do not consider the cross-selling effect (the negative terms will not emerge), the standard deviation of the original demand might be overestimated. Similarly, the third terms at the right-hand sides of Eq. (10) and Eq. (12) also show that the cross-selling might lead to overestimation of the demand expectation. Thus, for a cross-selling inventory system, the association between items must be considered to mitigate the overestimation of the original demand. In other words, the cross-selling will increase either the observable demand amount or the observable demand uncertainty.Besides, if we do not consider lost sales, the deviation and the mean of the original demand will be estimated based onSu2, Sv2, x-u, and x-v. In this case, even if there is no cross-selling, the estimator will inevitably give underestimated outputs because of incomplete information from the censored sales data. This implies that lost sales decrease both the quantity and uncertainty of the effective demand.


### 4.2. Estimators of the Correlation Coefficient

To estimate the correlation coefficientρ, define set q={(xiτ,xjτ)|xiτ<Qiandxjτ<Qj} representing all the sample points whose sales amounts are smaller than the order quantities. Evidently, for (xiτ,xjτ)∈q, the sales data are not censored so that we have(13)xiτ=ξiτ+λjiξjτxjτ=ξjτ+λijξiτ,where ξiτ and ξjτ are the realized “naked” demand for items i and j, respectively. Solving the above bivariate linear equation group yields(14)ξiτ=xiτ-λjixjτ1-λijλjiξjτ=xjτ-λijxiτ1-λijλji.As the correlation is independent of the position where sample points are distributed, we will use {ξiτ,ξjτ}τ∈q to estimate the correlation coefficient. Let |q| be the cardinality of set q, and it follows that the estimator of the correlation coefficient is(15)ρ~=1q-1∑τ∈qξiτ-ξ-iqSiqξjτ-ξ-jqSjq,where ξ-iq and ξ-jq are the means of {ξiτ}τ∈q and {ξjτ}τ∈q, respectively, and Siq=(1/(|q|-1))∑τ∈qξiτ-ξ-iq2 and Sjq=(1/(|q|-1))∑τ∈qξjτ-ξ-jq2 are the corresponding sample standard deviations.


### 4.3. Iterative Estimation for Two Items

Propositions5 and 6 give the demand estimation of one of the two cross-selling items when the other item's demand and correlation coefficient are given. However, the other item's original demand is unknown before accomplishing the estimation, and thus a circular deadlock is created. Here, we present an iterative procedure to obtain the estimation of the demand.Specifically, we first estimate the correlation coefficientρ by Eq. (15) and obtain the estimator ρ~. Then, we estimate the two items' demand without considering cross-selling based on the sampled data. This step is easy to realize based on Eqs. (9)–(12) by setting rij=rji=0, which yields (16)σ~i0=PiSuPi2-Piziϕzi-ϕzi2,μ~i0=x-u+ϕziσ~iPi,(17)σ~j0=PjSvPj2-Pjzjϕzj-ϕzj2,μ~j0=x-v+ϕzjσ~jPj.Subsequently, substituting (μ~i0,μ~j0,σ~i0,σ~j0) and ρ~ into Eqs. (9)–(12), we obtain the first iterative value of the estimators as(18)σ~i1=ρ~2-1rji2σ~j02+Su21-ziϕzi/Pi-ϕzi2/Pi2-ρ~rjiσ~j0,(19)μ~i1=x-u+ϕziσ~i12+rji2σ~j02+2ρ~rjiσ~i1σ~j0Pi-rjiμ~j0,(20)σ~j1=ρ~2-1rij2σ~i12+Sv21-zjϕzj/Pj-ϕzj2/Pj2-ρ~rijσ~i1.(21)μ~j1=x-v+ϕzjσ~j12+rij2σ~i12+2ρ~rijσ~i1σ~j1Pj-rijμ~i1.Again, we substitute (μ~i1,μ~j1,σ~i1,σ~j1) and ρ~ into Eqs. (9)–(12), yielding (μ~i2,μ~j2,σ~i2,σ~j2). Generally, for step k≥1, we have(22)σ~ik=ρ~2-1rji2σ~j,k-12+Su21-ziϕzi/Pi-ϕzi2/Pi2-ρ~rjiσ~j,k-1,(23)μ~ik=x-u+ϕziσ~ik2+rji2σ~j,k-12+2ρ~rjiσ~ikσ~j,k-1Pi-rjiμ~j,k-1,(24)σ~jk=ρ~2-1rij2σ~ik2+Sv21-zjϕzj/Pj-ϕzj2/Pj2-ρ~rijσ~ik,(25)μ~jk=x-v+ϕzjσ~jk2+rij2σ~ik2+2ρ~rijσ~ikσ~jkPj-rijμ~ik.The above procedures are repeated until a certain converging criterion is satisfied. We summarize the above procedure as follows.Step 1.
Calculateρ~, the estimated value of correlation coefficient from Eq. (15).Step 2.
Determine initial values of demand expectation and standard deviation for itemsi,j from Eqs. (16) and (17).Step 3.
For thek-th iteration (k≥1), estimate the demand expectation and standard deviation for items i and j from Eqs. (22)–(25).Step 4.
When the relative changes ofμi, μj, σi, and σj between two adjacent iterations are smaller than δ for q consecutive iterations, stop; otherwise, repeat Step 3.δ is a small positive real number and q is a positive integer. This study sets δ=10-6 and q=10. Regarding the convergence of the above procedure, we have Propositions 7 and 8.Proposition 7.
For positive correlated demand, i.e.,ρ≥0, the sequence σ~ik,σ~jk defined by Eqs. (22) and (24) are convergent.Proof.
See AppendixB.Proposition 8.
If the sequenceσ~ik,σ~jk are convergent, the sequence μ~ik,μ~jk defined by Eqs. (23) and (25) is convergent.Proof.
See AppendixB.It is very difficult to theoretically prove the convergence in case of negative correlated demand(ρ<0) because of the mathematical intractability. But the numerical experiment in Section 6 shows that, for each of more than two million numerical examples, the iterative procedures can quickly converge to a stable output even ρ<0. We conclude that, from a computational point of view, the iterative procedures based on our estimators are stable and convergent.


## 4.1. Estimators for One Item with Another Demand Given

Sorting increasingly all the sampled sales data, we will have two ordered data series as(1)Itemi:xi1,xi2,…,xiu,xi,u+1,…,xiT,(2)Itemj:xj1,xj2,…,xjv,xj,v+1,…,xjT,where xiu<Qi;xi,u+1=xi,u+2=…=xiT=Qi and xjv<Qj;xj,v+1=xj,v+2=…=xjT=Qj. u and v are the total number of xit and xjt that meet xit<Qi and xjt<Qj, respectively.Define(3)x-u=1u∑τ=1uxiτ,(4)x-v=1v∑τ=1vxjτ,(5)Su2=1u-1∑τ=1uxiτ-x-u2,(6)Sv2=1v-1∑τ=1vxjτ-x-v2,(7)Pi=uT,(8)Pj=vT,where x-u and x-v represent the mean values of the observed demand that are smaller than the order quantities, and Su2,Sv2 and Pi,Pj are the corresponding sample variances and proportions, respectively. Regarding the distribution estimation of the original demand, we have the following results:Proposition 5.
Given demand expectationμj and standard deviation σj of item j, and with a known correlation coefficient ρ, the value of μi and σi of item i can be estimated by Eq. (9) and Eq. (10), respectively.(9)σ~i=ρ2-1rji2σj2+Su21-ziϕzi/Pi-ϕzi2/Pi2-ρrjiσj,(10)μ~i=x-u+ϕziσ~i2+rji2σj2+2ρrjiσ~iσjPi-rjiμj,where ϕ(∙) and Φ(∙) are the density function and cumulative probability function of the standard normal distribution, respectively, and zi=Φ-1(Pi) is the quantile of normal distribution corresponding to probability Pi.Proof.
See AppendixAProposition 6.
Givenμi, σi, and ρ, the value of μj and σj can be estimated by Eq. (11) and Eq. (12), respectively. (11)σ~j=ρ2-1rij2σi2+Sv21-zjϕzj/Pj-ϕzj2/Pj2-ρrijσi,(12)μ~j=x-v+ϕzjσ~j2+rij2σi2+2ρrijσiσ~jPj-rijμi,where zj=Φ-1(Pj).Proof.
It is straightforward from Proposition5, because of the symmetry of the two items.Evidently, the second terms at the right-hand sides of Eq. (9) and Eq. (11), -ρrjiσj and -ρrijσi, are both negative values for calculating the standard deviations of items i and j. Not only that, the first terms under the radical sign, (ρ2-1)rji2σj2 and (ρ2-1)rij2σi2, are also negative. This implies that if we do not consider the cross-selling effect (the negative terms will not emerge), the standard deviation of the original demand might be overestimated. Similarly, the third terms at the right-hand sides of Eq. (10) and Eq. (12) also show that the cross-selling might lead to overestimation of the demand expectation. Thus, for a cross-selling inventory system, the association between items must be considered to mitigate the overestimation of the original demand. In other words, the cross-selling will increase either the observable demand amount or the observable demand uncertainty.Besides, if we do not consider lost sales, the deviation and the mean of the original demand will be estimated based onSu2, Sv2, x-u, and x-v. In this case, even if there is no cross-selling, the estimator will inevitably give underestimated outputs because of incomplete information from the censored sales data. This implies that lost sales decrease both the quantity and uncertainty of the effective demand.


## 4.2. Estimators of the Correlation Coefficient

To estimate the correlation coefficientρ, define set q={(xiτ,xjτ)|xiτ<Qiandxjτ<Qj} representing all the sample points whose sales amounts are smaller than the order quantities. Evidently, for (xiτ,xjτ)∈q, the sales data are not censored so that we have(13)xiτ=ξiτ+λjiξjτxjτ=ξjτ+λijξiτ,where ξiτ and ξjτ are the realized “naked” demand for items i and j, respectively. Solving the above bivariate linear equation group yields(14)ξiτ=xiτ-λjixjτ1-λijλjiξjτ=xjτ-λijxiτ1-λijλji.As the correlation is independent of the position where sample points are distributed, we will use {ξiτ,ξjτ}τ∈q to estimate the correlation coefficient. Let |q| be the cardinality of set q, and it follows that the estimator of the correlation coefficient is(15)ρ~=1q-1∑τ∈qξiτ-ξ-iqSiqξjτ-ξ-jqSjq,where ξ-iq and ξ-jq are the means of {ξiτ}τ∈q and {ξjτ}τ∈q, respectively, and Siq=(1/(|q|-1))∑τ∈qξiτ-ξ-iq2 and Sjq=(1/(|q|-1))∑τ∈qξjτ-ξ-jq2 are the corresponding sample standard deviations.


## 4.3. Iterative Estimation for Two Items

Propositions5 and 6 give the demand estimation of one of the two cross-selling items when the other item's demand and correlation coefficient are given. However, the other item's original demand is unknown before accomplishing the estimation, and thus a circular deadlock is created. Here, we present an iterative procedure to obtain the estimation of the demand.Specifically, we first estimate the correlation coefficientρ by Eq. (15) and obtain the estimator ρ~. Then, we estimate the two items' demand without considering cross-selling based on the sampled data. This step is easy to realize based on Eqs. (9)–(12) by setting rij=rji=0, which yields (16)σ~i0=PiSuPi2-Piziϕzi-ϕzi2,μ~i0=x-u+ϕziσ~iPi,(17)σ~j0=PjSvPj2-Pjzjϕzj-ϕzj2,μ~j0=x-v+ϕzjσ~jPj.Subsequently, substituting (μ~i0,μ~j0,σ~i0,σ~j0) and ρ~ into Eqs. (9)–(12), we obtain the first iterative value of the estimators as(18)σ~i1=ρ~2-1rji2σ~j02+Su21-ziϕzi/Pi-ϕzi2/Pi2-ρ~rjiσ~j0,(19)μ~i1=x-u+ϕziσ~i12+rji2σ~j02+2ρ~rjiσ~i1σ~j0Pi-rjiμ~j0,(20)σ~j1=ρ~2-1rij2σ~i12+Sv21-zjϕzj/Pj-ϕzj2/Pj2-ρ~rijσ~i1.(21)μ~j1=x-v+ϕzjσ~j12+rij2σ~i12+2ρ~rijσ~i1σ~j1Pj-rijμ~i1.Again, we substitute (μ~i1,μ~j1,σ~i1,σ~j1) and ρ~ into Eqs. (9)–(12), yielding (μ~i2,μ~j2,σ~i2,σ~j2). Generally, for step k≥1, we have(22)σ~ik=ρ~2-1rji2σ~j,k-12+Su21-ziϕzi/Pi-ϕzi2/Pi2-ρ~rjiσ~j,k-1,(23)μ~ik=x-u+ϕziσ~ik2+rji2σ~j,k-12+2ρ~rjiσ~ikσ~j,k-1Pi-rjiμ~j,k-1,(24)σ~jk=ρ~2-1rij2σ~ik2+Sv21-zjϕzj/Pj-ϕzj2/Pj2-ρ~rijσ~ik,(25)μ~jk=x-v+ϕzjσ~jk2+rij2σ~ik2+2ρ~rijσ~ikσ~jkPj-rijμ~ik.The above procedures are repeated until a certain converging criterion is satisfied. We summarize the above procedure as follows.Step 1.
Calculateρ~, the estimated value of correlation coefficient from Eq. (15).Step 2.
Determine initial values of demand expectation and standard deviation for itemsi,j from Eqs. (16) and (17).Step 3.
For thek-th iteration (k≥1), estimate the demand expectation and standard deviation for items i and j from Eqs. (22)–(25).Step 4.
When the relative changes ofμi, μj, σi, and σj between two adjacent iterations are smaller than δ for q consecutive iterations, stop; otherwise, repeat Step 3.δ is a small positive real number and q is a positive integer. This study sets δ=10-6 and q=10. Regarding the convergence of the above procedure, we have Propositions 7 and 8.Proposition 7.
For positive correlated demand, i.e.,ρ≥0, the sequence σ~ik,σ~jk defined by Eqs. (22) and (24) are convergent.Proof.
See AppendixB.Proposition 8.
If the sequenceσ~ik,σ~jk are convergent, the sequence μ~ik,μ~jk defined by Eqs. (23) and (25) is convergent.Proof.
See AppendixB.It is very difficult to theoretically prove the convergence in case of negative correlated demand(ρ<0) because of the mathematical intractability. But the numerical experiment in Section 6 shows that, for each of more than two million numerical examples, the iterative procedures can quickly converge to a stable output even ρ<0. We conclude that, from a computational point of view, the iterative procedures based on our estimators are stable and convergent.


## 5. An Example for Illustration

### 5.1. The Sample Data

Suppose that the original “naked” demand for two cross-selling items is jointly distributed as(Di,Dj)~N(μi,σi2;μj,σj2;ρ), where μi=300,σi=10;μj=400,σj=20;ρ=0.6. The two items are cross-sold with cross-selling coefficients rij=0.2 and rji=0.3. The order quantities are Qi=430,Qj=470 in all periods.If the original “naked” demand is realized as summarized in Table4, the total demand without lost sales for the two items should be as summarized in Table 5, considering that Xi=Di+rji∙Dj and Xj=Dj+rij∙Di.Table 4
Original realized “naked” demand.
t 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 ξ i t 313 311 289 303 292 296 290 272 293 286 305 294 308 315 283 ξ j t 395 418 391 411 420 377 388 363 369 356 412 422 412 425 367  t 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30  ξ i t 309 290 299 299 299 291 295 284 286 292 301 309 300 295 303 ξ j t 404 366 390 396 375 418 397 388 375 418 382 422 392 386 387  t 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45  ξ i t 299 307 317 290 301 289 306 306 282 305 301 287 286 290 321 ξ j t 387 428 416 378 405 410 437 402 374 401 406 392 406 380 421  t 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60  ξ i t 294 269 285 314 302 306 305 309 295 288 321 321 305 291 297 ξ j t 381 350 380 419 420 403 396 382 395 386 401 438 428 369 393Table 5
Total demand (without lost sales) for two items.
t 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15  x i t 431.5 436.4 406.3 426.3 418.0 409.1 406.4 380.9 403.7 392.8 428.6 420.6 431.6 442.5 393.1 x j t 457.6 480.2 448.8 471.6 478.4 436.2 446.0 417.4 427.6 413.2 473 480.8 473.6 488.0 423.6  t 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30  x i t 430.2 399.8 416.0 417.8 411.5 416.4 414.1 400.4 398.5 417.4 415.6 435.6 417.6 410.8 419.1 x j t 465.8 424.0 449.8 455.8 434.8 476.2 456.0 444.8 432.2 476.4 442.2 483.8 452.0 445.0 447.6  t 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45  x i t 415.1 435.4 441.8 403.4 422.5 412.0 437.1 426.6 394.2 425.3 422.8 404.6 407.8 404.0 447.3 x j t 446.8 489.4 479.4 436.0 465.2 467.8 498.2 463.2 430.4 462.0 466.2 449.4 463.2 438.0 485.2  t 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60  x i t 408.3 374.0 399.0 439.7 428.0 426.9 423.8 423.6 413.5 403.8 441.3 452.4 433.4 401.7 414.9 x j t 439.8 403.8 437.0 481.8 480.4 464.2 457.0 443.8 454.0 443.6 465.2 502.2 489.0 427.2 452.4Because of the limitation of order quantities, the observed demand (sales amount) will be censored as summarized in Table6 according to xit=min⁡(Qi,Xit) and xjt=min⁡(Qj,Xjt).Table 6
Observed demand censored by order quantities.
t 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15  x i t 430 430 406.3 426.3 418.0 409.1 406.4 380.9 403.7 392.8 428.6 420.6 430 430 393.1 x j t 457.6 470 448.8 470 470 436.2 446.0 417.4 427.6 413.2 470 470 470 470 423.6  t 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30  x i t 430 399.8 416.0 417.8 411.5 416.4 414.1 400.4 398.5 417.4 415.6 406 430 419 424 x j t 465.8 424.0 449.8 455.8 434.8 470 456.0 444.8 432.2 470 452 442.2 452.0 445.0 447.6  t 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45  x i t 415.1 430 430 403.4 422.5 412.0 430 426.6 394.2 425.3 422.8 404.6 407.8 404.0 430 x j t 446.8 470 470 436.0 465.2 467.8 470 463.2 430.4 462.0 466.2 449.4 463.2 438.0 470  t 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60  x i t 408.3 374.0 399.0 430 428.0 426.9 423.8 423.6 413.5 403.8 430 430 430 401.7 414.9 x j t 439.8 403.8 437.0 470 470 464.2 457.0 443.8 454.0 443.6 465.2 470 470 427.2 452.4We will estimateμi,σi,μj,σj, and ρ based on the observed sales amount in Table 6.


### 5.2. The Results

In total, there are 60 sampled data, among which 39 points meetxit<Qi and xjt<Qj, and thus the coefficient ρ is estimated from Eq. (15) as ρ~=0.5913. Moreover, there are 14 data points of item i censored and 18 data points of item j censored so that u=46 and v=42. Thus, Pi=u/T=46/60=0.7667 and Pj=v/T=42/60=0.7000. We, therefore, have zi=Φ-1Pi=Φ-10.6833=0.7279 and zj=Φ-1Pj=Φ-10.6000=0.5244. Subsequently, we determine the initial value of (μ~i,σ~i,μ~j,σ~j) from Eqs. (16) and (17), yielding(26)μ~i0=305.03,σ~i0=12.00;μ~j0=374.51,σ~j0=20.88.Starting the iterative procedure, we will obtain the estimated results when the iteration converges, which yields(27)μ~i=298.41,σ~i=12.01;μ~j=396.76,σ~j=20.86.Figure 1 presents the converging process of the output by the estimators on each round of iteration.Figure 1
Estimated demand expectations and standard deviations.In Figure1, the lines with triangles represent thepoint estimations obtained based on the sampling original “naked” demand ξit,ξjt in Table 4. It is easy to find that the point estimations are(28)μ~ξi=1T∑ξit=298.18,σ~ξi=1T-1∑ξit-ξ-i2=11.26,μ~ξj=1T∑ξjt=396.77,σ~ξj=1T-1∑ξjt-ξ-j2=20.41.Evidently, the results of our method are very close to the corresponding point estimations. The relative error is smaller than 1%. It shows that the performance of our estimators is as good as that of thepoint estimators in terms of the data sampled from the original demand. Similarly, we estimate ρ based on ξit,ξjt and obtain ρ~ξ=0.6557. Our estimator gives the estimation ρ~=0.5913, which is also close to ρ~ξ.It is noticeable that we in fact cannot directly observe the realized original demandξit,ξjt in practice. However, for the numerical experiments, we can use the point estimator based on the ξit,ξjt data as a benchmark to evaluate the accuracy of the new estimators proposed herein.


## 5.1. The Sample Data

Suppose that the original “naked” demand for two cross-selling items is jointly distributed as(Di,Dj)~N(μi,σi2;μj,σj2;ρ), where μi=300,σi=10;μj=400,σj=20;ρ=0.6. The two items are cross-sold with cross-selling coefficients rij=0.2 and rji=0.3. The order quantities are Qi=430,Qj=470 in all periods.If the original “naked” demand is realized as summarized in Table4, the total demand without lost sales for the two items should be as summarized in Table 5, considering that Xi=Di+rji∙Dj and Xj=Dj+rij∙Di.Table 4
Original realized “naked” demand.
t 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 ξ i t 313 311 289 303 292 296 290 272 293 286 305 294 308 315 283 ξ j t 395 418 391 411 420 377 388 363 369 356 412 422 412 425 367  t 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30  ξ i t 309 290 299 299 299 291 295 284 286 292 301 309 300 295 303 ξ j t 404 366 390 396 375 418 397 388 375 418 382 422 392 386 387  t 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45  ξ i t 299 307 317 290 301 289 306 306 282 305 301 287 286 290 321 ξ j t 387 428 416 378 405 410 437 402 374 401 406 392 406 380 421  t 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60  ξ i t 294 269 285 314 302 306 305 309 295 288 321 321 305 291 297 ξ j t 381 350 380 419 420 403 396 382 395 386 401 438 428 369 393Table 5
Total demand (without lost sales) for two items.
t 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15  x i t 431.5 436.4 406.3 426.3 418.0 409.1 406.4 380.9 403.7 392.8 428.6 420.6 431.6 442.5 393.1 x j t 457.6 480.2 448.8 471.6 478.4 436.2 446.0 417.4 427.6 413.2 473 480.8 473.6 488.0 423.6  t 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30  x i t 430.2 399.8 416.0 417.8 411.5 416.4 414.1 400.4 398.5 417.4 415.6 435.6 417.6 410.8 419.1 x j t 465.8 424.0 449.8 455.8 434.8 476.2 456.0 444.8 432.2 476.4 442.2 483.8 452.0 445.0 447.6  t 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45  x i t 415.1 435.4 441.8 403.4 422.5 412.0 437.1 426.6 394.2 425.3 422.8 404.6 407.8 404.0 447.3 x j t 446.8 489.4 479.4 436.0 465.2 467.8 498.2 463.2 430.4 462.0 466.2 449.4 463.2 438.0 485.2  t 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60  x i t 408.3 374.0 399.0 439.7 428.0 426.9 423.8 423.6 413.5 403.8 441.3 452.4 433.4 401.7 414.9 x j t 439.8 403.8 437.0 481.8 480.4 464.2 457.0 443.8 454.0 443.6 465.2 502.2 489.0 427.2 452.4Because of the limitation of order quantities, the observed demand (sales amount) will be censored as summarized in Table6 according to xit=min⁡(Qi,Xit) and xjt=min⁡(Qj,Xjt).Table 6
Observed demand censored by order quantities.
t 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15  x i t 430 430 406.3 426.3 418.0 409.1 406.4 380.9 403.7 392.8 428.6 420.6 430 430 393.1 x j t 457.6 470 448.8 470 470 436.2 446.0 417.4 427.6 413.2 470 470 470 470 423.6  t 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30  x i t 430 399.8 416.0 417.8 411.5 416.4 414.1 400.4 398.5 417.4 415.6 406 430 419 424 x j t 465.8 424.0 449.8 455.8 434.8 470 456.0 444.8 432.2 470 452 442.2 452.0 445.0 447.6  t 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45  x i t 415.1 430 430 403.4 422.5 412.0 430 426.6 394.2 425.3 422.8 404.6 407.8 404.0 430 x j t 446.8 470 470 436.0 465.2 467.8 470 463.2 430.4 462.0 466.2 449.4 463.2 438.0 470  t 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60  x i t 408.3 374.0 399.0 430 428.0 426.9 423.8 423.6 413.5 403.8 430 430 430 401.7 414.9 x j t 439.8 403.8 437.0 470 470 464.2 457.0 443.8 454.0 443.6 465.2 470 470 427.2 452.4We will estimateμi,σi,μj,σj, and ρ based on the observed sales amount in Table 6.


## 5.2. The Results

In total, there are 60 sampled data, among which 39 points meetxit<Qi and xjt<Qj, and thus the coefficient ρ is estimated from Eq. (15) as ρ~=0.5913. Moreover, there are 14 data points of item i censored and 18 data points of item j censored so that u=46 and v=42. Thus, Pi=u/T=46/60=0.7667 and Pj=v/T=42/60=0.7000. We, therefore, have zi=Φ-1Pi=Φ-10.6833=0.7279 and zj=Φ-1Pj=Φ-10.6000=0.5244. Subsequently, we determine the initial value of (μ~i,σ~i,μ~j,σ~j) from Eqs. (16) and (17), yielding(26)μ~i0=305.03,σ~i0=12.00;μ~j0=374.51,σ~j0=20.88.Starting the iterative procedure, we will obtain the estimated results when the iteration converges, which yields(27)μ~i=298.41,σ~i=12.01;μ~j=396.76,σ~j=20.86.Figure 1 presents the converging process of the output by the estimators on each round of iteration.Figure 1
Estimated demand expectations and standard deviations.In Figure1, the lines with triangles represent thepoint estimations obtained based on the sampling original “naked” demand ξit,ξjt in Table 4. It is easy to find that the point estimations are(28)μ~ξi=1T∑ξit=298.18,σ~ξi=1T-1∑ξit-ξ-i2=11.26,μ~ξj=1T∑ξjt=396.77,σ~ξj=1T-1∑ξjt-ξ-j2=20.41.Evidently, the results of our method are very close to the corresponding point estimations. The relative error is smaller than 1%. It shows that the performance of our estimators is as good as that of thepoint estimators in terms of the data sampled from the original demand. Similarly, we estimate ρ based on ξit,ξjt and obtain ρ~ξ=0.6557. Our estimator gives the estimation ρ~=0.5913, which is also close to ρ~ξ.It is noticeable that we in fact cannot directly observe the realized original demandξit,ξjt in practice. However, for the numerical experiments, we can use the point estimator based on the ξit,ξjt data as a benchmark to evaluate the accuracy of the new estimators proposed herein.


## 6. Numerical Experiments

We conduct numerical experiments with two aims: (1) to evaluate the estimation accuracy of the proposed method by comparing our estimators with the point estimation of the sampled original demandξit,ξjt; (2) to determine in a sensitivity analysis how the parameters impact the estimation.


### 6.1. Scheme for Generating Numerical Examples

The numerical experiment will be carried out based on artificially generated data. To investigate as many cases as possible, we will rearrange the values of the parameters from low to high, including the expectation, coefficient of variation (CV), and correlation coefficient of demand distribution, along with the cross-selling coefficient, censoring levels (order quantities), and sample size. Specifically, we select values of parameters for generating numerical examples, as summarized in Table7.Table 7
Selected values for parameters.
Parameter Notation Values to be selected Expectation of itemi’s original demand μ i 100, 300, 500 Expectation of itemj’s original demand μ j 100, 300, 500 Coefficient of variation of itemi C V i = σ i / μ i 0.05, 0.15, 0.25 Coefficient of variation of itemj C V j = σ j / μ j 0.05, 0.15, 0.25 Cross-selling coefficient from itemi to j r i j 0.1, 0.3, 0.5, 0.7, 0.9 Cross-selling coefficient from itemj to i r j i 0.1, 0.3, 0.5, 0.7, 0.9 Correlation coefficient of demand ρ −0.9, −0.7, −0.5, −0.3, −0.1, 0.1, 0.3, 0.5, 0.7, 0.9 Order quantities (censoring level) ofi Q i μ X i, μXi+σXi, μXi+2σXi, μXi+3σXi Order quantities (censoring level) ofj Q j μ X j, μXj+σXj, μXj+2σXj, μXj+3σXj Sample size T 50, 100, 300, 500, 700, 900, 1000Note thatμXi,μXj and σXi,σXj in Table 7 represent the expectation and standard deviation of total demand Xi=Di+rji∙Dj and Xj=Dj+rij∙Di, respectively. According to Table 7, in total, we will test 32×32×52×10×42×7=2,268,000 numerical examples.


### 6.2. Computational Results

#### 6.2.1. Relative Errors and Censoring Levels

As mentioned in Section5.2, we will evaluate the performance of our estimation method by comparing the results with thepoint estimation of the data ξit,ξjt that is sampled from the original demand. We define the relative error of the estimated expectation of demand as(29)eμi=μ~i-μ~ξiμ~ξi×100%,eμj=μ~j-μ~ξjμ~ξj×100%,where μ~ξi,μ~ξj are the corresponding point estimations. Similarly, the relative errors of the standard deviation and the correlation coefficient are(30)eσi=σ~i-σ~ξiσ~ξi×100%,eσj=σ~j-σ~ξjσ~ξj×100%,eρ=ρ~-ρ~ξρ~ξ×100%.Moreover, the censoring level is discriminated by the deviation of order quantities over the corresponding demand expectation, i.e., the multiples ofσXi and σXj. To be concise, we adopt the sum of the two multiples to comprehensively discriminate different censoring levels. The results will be reported considering the seven comprehensive censoring levels as those in the first column of Table 8.Table 8
Comprehensive censoring levels.
Comprehensive censoring levels Multiples ofσXi and σXj 0 σ 0 σ X i ,  0 σ X j 1 σ 0 σ X i , 1 σ X j  o r  1 σ X i , 0 σ X j 2 σ 0 σ X i , 2 σ X j  o r  1 σ X i , 1 σ X j  o r  2 σ X i , 0 σ X j 3 σ 0 σ X i , 3 σ X j  o r  1 σ X i , 2 σ X j  o r  2 σ X i , 1 σ X j  o r  3 σ X i , 0 σ X j 4 σ 1 σ X i , 3 σ X j  o r  2 σ X i , 2 σ X j  o r  3 σ X i , 1 σ X j 5 σ 2 σ X i , 3 σ X j  o r  3 σ X i , 2 σ X j 6 σ 3 σ X i , 3 σ X j


#### 6.2.2. Results

For the seven censoring levels, we report the estimation errors ofμi, μj, σi, and σj, under different sample sizes, correlation coefficients, cross-selling coefficients, and coefficients of variation, respectively.(1) Errors vs. Sample Size ( T ).Figure 2 presents the average error under different numbers of samples, i.e., T=50,100,300,500,700,900,1000.Figure 2
Variations in estimation errors with sample size.Figure 3
Variation in estimation errors with correlation coefficient.We observe from Figure2 that the demand expectation estimators yield very good estimation for the demand expectations of both items. Even for 50 sample points, the error is a maximum of 2.5%. With the standard deviation estimators, the error can reach more than 20% when the sample size is too small. However, if the sample size is larger than 800, the estimator performs much better, with an error no larger than 5%. We consider that for many modern business applications of inventory management, the data is no longer an obstruction. In particular, with the development of data acquisition and storage technologies, a large amount of data can be obtained from POS/MIS/MES, and big data systems [31]. The accuracy of our estimators can be further improved based on practical large sample data.Moreover, it is expected that the estimation error will become smaller when censoring is alleviated. The results show that if the comprehensive censoring level reaches4σ (implying that the safety stock of the two items is in average 2σ-level with a fill-rate of 95.45%), the error of the demand expectation estimator decreases from about 2.0% of 0σ to one no larger than 1.0%, and even to an error no larger than 0.5% with a 6σ comprehensive censoring level (corresponding to a 3σ-level safety stock with a fill-rate of 99.74%). We also see that, with high fill-rate, the standard deviation error decreases from larger than 15% of small sample data to smaller than 2% of large sample data.(  2) Errors vs. Correlation Coefficient ( ρ ).For the value of demand expectation, the estimation error slightly increases with the correlation coefficient, except for the censoring level of 0σ. However, almost all estimation errors of the mean demand are no larger than 1.4%, and they even decrease to 0.2% when censoring is sufficiently alleviated with a comprehensive censoring level no smaller than 4σ (corresponding to 2σ-level fill rate for each item in average), irrespective of the correlation coefficient.It can also be observed that, under a higher safety stock level (with comprehensive censoring level no smaller than4σ), the relative error of the standard deviation estimators is basically decreasing with the correlation coefficient. However, when censoring is aggravated (with censoring level no larger than 4σ), for weak correlation cases, the estimation error of the standard deviation reaches the maximal value, but it decreases when the correlation becomes stronger (irrespective of being positive or negative). The explanation is that, for a larger value of ρ, the estimator could capture more impact of the cross-selling effect (see Eqs. (9) and (11)), which might be beneficial to estimate the standard deviation of the demand.(  3) Errors vs. Cross-Selling Coefficient ( r i j , r j i ).The cross-selling is of dual-direction, so that for each one of the seven censoring levels we must consider two impact factors: rij and rji, which together will lead to a lengthy data report consisting of 28 two-dimensional tables. To be concise, we present a comprehensive index, r=(rij+rji)/2, to discriminate different levels of the cross-selling effect, as summarized in Table 9.Table 9
Comprehensive cross-selling coefficients.
Comprehensive cross-selling coefficientr=(rij+rji)/2 Values of(rij,rji) 0.1 ( 0.1,0.1 ) 0.2 0.1,0.3  o r  ( 0.3,0.1 ) 0.3 ( 0.1,0.5 )  o r  ( 0.5,0.1 ) 0.4 0.1,0.7  o r  0.3,0.5  o r  0.5,0.3  o r  ( 0.7,0.1 ) 0.5 ( 0.1,0.9 )  o r  ( 0.3,0.7 )  o r  ( 0.7,0.3 )  o r  ( 0.9,0.1 ) 0.6 ( 0.3,0.9 )  o r  ( 0.5,0.7 )  o r  ( 0.7,0.5 )  o r  ( 0.9,0.3 ) 0.7 ( 0.5,0.9 )  o r  ( 0.9,0.5 ) 0.8 ( 0.7,0.9 )  o r  ( 0.9,0.7 ) 0.9 ( 0.9,0.9 )For the different comprehensive cross-selling coefficients in Table9, the average relative errors of μi, μj, σi, and σj are presented in Figure 4.Figure 4
Variation in estimation errors with cross-selling effect.The estimation errors of the demand expectation increase with the cross-selling coefficient, but they do not exceed 4.50% for all censoring levels. It means that the estimation forμi and μj has a wide range of applicability for handling different levels of cross-selling. For the standard deviation, the error tends to increase with the cross-selling coefficient, but it is still smaller than 5% for censoring levels of 5σ and 6σ in most cases. In particular, for 6σ-level censoring, the relative error is not really increasing with the cross-selling effect, and it even slightly decreases with larger cross-selling coefficients.(  4) Errors vs. Coefficient of Variation C V i , C V j. The two coefficients of variation, CVi and CVj, affect the estimation. Similarly, to be concise, we introduce a comprehensive coefficient of variation, cv=(CVi+CVj)/2, to discriminate different variation levels, as summarized in Table 10.Table 10
Comprehensive coefficient of variation.
Comprehensive coefficient of variationcv=(CVi+CVj)/2 Values of(CVi,CVj) 0.05 ( 0.05,0.05 ) 0.1 0.05,0.15  o r  ( 0.15,0.05 ) 0.15 0.05,0.25  o r  0.15,0.15  o r  ( 0.25,0.05 ) 0.2 ( 0.15,0.25 )  o r  ( 0.25,0.15 ) 0.25 ( 0.25,0.25 )For the different comprehensive coefficients of variation in Table10, the average relative errors of μi, μj, σi, and σj are presented in Figure 5.Figure 5
Estimation errors with the coefficient of variation (CV).Figure 6
Estimation errors of correlation coefficient.Figure5 presents that the estimation errors of μi and μj increase with the coefficient of variation, but even for a 20% variation both are still smaller than 1.5%. However, regarding the standard deviations, σi and σj, the increase in order quantity for alleviating censoring is still an important means to obtain a good estimation. Only when the comprehensive censoring level reaches more than5σ, the relative error becomes smaller than 5%.(4) Estimation Errors of Correlation Coefficient ( ρ ). We can see that the sample size is still a positive factor for the estimators of correlation coefficient, and it is evident that the more slightly the data is censored, the better the estimator will perform. We also observe that the error is smaller in negative correlations than in positive correlations on the overall trend. However, the estimation of ρ still reaches a high precision provided that the fill rate of the inventory is high enough. Remarkably, for the 6σ comprehensive censoring level, the error is no larger than 1.0% on average.Combining all the above results, it is revealed that for inventory systems with at least2σ-level safety stock (corresponding to 4σ comprehensive censoring level), the estimation generates high-precision outcomes for larger sample data. It is worth noting that 2σ-level safety stock, i.e., 95% fill rate, is a basic service level for daily operations of a store. Hence, our approach could be of high practical relevance.


## 6.1. Scheme for Generating Numerical Examples

The numerical experiment will be carried out based on artificially generated data. To investigate as many cases as possible, we will rearrange the values of the parameters from low to high, including the expectation, coefficient of variation (CV), and correlation coefficient of demand distribution, along with the cross-selling coefficient, censoring levels (order quantities), and sample size. Specifically, we select values of parameters for generating numerical examples, as summarized in Table7.Table 7
Selected values for parameters.
Parameter Notation Values to be selected Expectation of itemi’s original demand μ i 100, 300, 500 Expectation of itemj’s original demand μ j 100, 300, 500 Coefficient of variation of itemi C V i = σ i / μ i 0.05, 0.15, 0.25 Coefficient of variation of itemj C V j = σ j / μ j 0.05, 0.15, 0.25 Cross-selling coefficient from itemi to j r i j 0.1, 0.3, 0.5, 0.7, 0.9 Cross-selling coefficient from itemj to i r j i 0.1, 0.3, 0.5, 0.7, 0.9 Correlation coefficient of demand ρ −0.9, −0.7, −0.5, −0.3, −0.1, 0.1, 0.3, 0.5, 0.7, 0.9 Order quantities (censoring level) ofi Q i μ X i, μXi+σXi, μXi+2σXi, μXi+3σXi Order quantities (censoring level) ofj Q j μ X j, μXj+σXj, μXj+2σXj, μXj+3σXj Sample size T 50, 100, 300, 500, 700, 900, 1000Note thatμXi,μXj and σXi,σXj in Table 7 represent the expectation and standard deviation of total demand Xi=Di+rji∙Dj and Xj=Dj+rij∙Di, respectively. According to Table 7, in total, we will test 32×32×52×10×42×7=2,268,000 numerical examples.


## 6.2. Computational Results

### 6.2.1. Relative Errors and Censoring Levels

As mentioned in Section5.2, we will evaluate the performance of our estimation method by comparing the results with thepoint estimation of the data ξit,ξjt that is sampled from the original demand. We define the relative error of the estimated expectation of demand as(29)eμi=μ~i-μ~ξiμ~ξi×100%,eμj=μ~j-μ~ξjμ~ξj×100%,where μ~ξi,μ~ξj are the corresponding point estimations. Similarly, the relative errors of the standard deviation and the correlation coefficient are(30)eσi=σ~i-σ~ξiσ~ξi×100%,eσj=σ~j-σ~ξjσ~ξj×100%,eρ=ρ~-ρ~ξρ~ξ×100%.Moreover, the censoring level is discriminated by the deviation of order quantities over the corresponding demand expectation, i.e., the multiples ofσXi and σXj. To be concise, we adopt the sum of the two multiples to comprehensively discriminate different censoring levels. The results will be reported considering the seven comprehensive censoring levels as those in the first column of Table 8.Table 8
Comprehensive censoring levels.
Comprehensive censoring levels Multiples ofσXi and σXj 0 σ 0 σ X i ,  0 σ X j 1 σ 0 σ X i , 1 σ X j  o r  1 σ X i , 0 σ X j 2 σ 0 σ X i , 2 σ X j  o r  1 σ X i , 1 σ X j  o r  2 σ X i , 0 σ X j 3 σ 0 σ X i , 3 σ X j  o r  1 σ X i , 2 σ X j  o r  2 σ X i , 1 σ X j  o r  3 σ X i , 0 σ X j 4 σ 1 σ X i , 3 σ X j  o r  2 σ X i , 2 σ X j  o r  3 σ X i , 1 σ X j 5 σ 2 σ X i , 3 σ X j  o r  3 σ X i , 2 σ X j 6 σ 3 σ X i , 3 σ X j


### 6.2.2. Results

For the seven censoring levels, we report the estimation errors ofμi, μj, σi, and σj, under different sample sizes, correlation coefficients, cross-selling coefficients, and coefficients of variation, respectively.(1) Errors vs. Sample Size ( T ).Figure 2 presents the average error under different numbers of samples, i.e., T=50,100,300,500,700,900,1000.Figure 2
Variations in estimation errors with sample size.Figure 3
Variation in estimation errors with correlation coefficient.We observe from Figure2 that the demand expectation estimators yield very good estimation for the demand expectations of both items. Even for 50 sample points, the error is a maximum of 2.5%. With the standard deviation estimators, the error can reach more than 20% when the sample size is too small. However, if the sample size is larger than 800, the estimator performs much better, with an error no larger than 5%. We consider that for many modern business applications of inventory management, the data is no longer an obstruction. In particular, with the development of data acquisition and storage technologies, a large amount of data can be obtained from POS/MIS/MES, and big data systems [31]. The accuracy of our estimators can be further improved based on practical large sample data.Moreover, it is expected that the estimation error will become smaller when censoring is alleviated. The results show that if the comprehensive censoring level reaches4σ (implying that the safety stock of the two items is in average 2σ-level with a fill-rate of 95.45%), the error of the demand expectation estimator decreases from about 2.0% of 0σ to one no larger than 1.0%, and even to an error no larger than 0.5% with a 6σ comprehensive censoring level (corresponding to a 3σ-level safety stock with a fill-rate of 99.74%). We also see that, with high fill-rate, the standard deviation error decreases from larger than 15% of small sample data to smaller than 2% of large sample data.(  2) Errors vs. Correlation Coefficient ( ρ ).For the value of demand expectation, the estimation error slightly increases with the correlation coefficient, except for the censoring level of 0σ. However, almost all estimation errors of the mean demand are no larger than 1.4%, and they even decrease to 0.2% when censoring is sufficiently alleviated with a comprehensive censoring level no smaller than 4σ (corresponding to 2σ-level fill rate for each item in average), irrespective of the correlation coefficient.It can also be observed that, under a higher safety stock level (with comprehensive censoring level no smaller than4σ), the relative error of the standard deviation estimators is basically decreasing with the correlation coefficient. However, when censoring is aggravated (with censoring level no larger than 4σ), for weak correlation cases, the estimation error of the standard deviation reaches the maximal value, but it decreases when the correlation becomes stronger (irrespective of being positive or negative). The explanation is that, for a larger value of ρ, the estimator could capture more impact of the cross-selling effect (see Eqs. (9) and (11)), which might be beneficial to estimate the standard deviation of the demand.(  3) Errors vs. Cross-Selling Coefficient ( r i j , r j i ).The cross-selling is of dual-direction, so that for each one of the seven censoring levels we must consider two impact factors: rij and rji, which together will lead to a lengthy data report consisting of 28 two-dimensional tables. To be concise, we present a comprehensive index, r=(rij+rji)/2, to discriminate different levels of the cross-selling effect, as summarized in Table 9.Table 9
Comprehensive cross-selling coefficients.
Comprehensive cross-selling coefficientr=(rij+rji)/2 Values of(rij,rji) 0.1 ( 0.1,0.1 ) 0.2 0.1,0.3  o r  ( 0.3,0.1 ) 0.3 ( 0.1,0.5 )  o r  ( 0.5,0.1 ) 0.4 0.1,0.7  o r  0.3,0.5  o r  0.5,0.3  o r  ( 0.7,0.1 ) 0.5 ( 0.1,0.9 )  o r  ( 0.3,0.7 )  o r  ( 0.7,0.3 )  o r  ( 0.9,0.1 ) 0.6 ( 0.3,0.9 )  o r  ( 0.5,0.7 )  o r  ( 0.7,0.5 )  o r  ( 0.9,0.3 ) 0.7 ( 0.5,0.9 )  o r  ( 0.9,0.5 ) 0.8 ( 0.7,0.9 )  o r  ( 0.9,0.7 ) 0.9 ( 0.9,0.9 )For the different comprehensive cross-selling coefficients in Table9, the average relative errors of μi, μj, σi, and σj are presented in Figure 4.Figure 4
Variation in estimation errors with cross-selling effect.The estimation errors of the demand expectation increase with the cross-selling coefficient, but they do not exceed 4.50% for all censoring levels. It means that the estimation forμi and μj has a wide range of applicability for handling different levels of cross-selling. For the standard deviation, the error tends to increase with the cross-selling coefficient, but it is still smaller than 5% for censoring levels of 5σ and 6σ in most cases. In particular, for 6σ-level censoring, the relative error is not really increasing with the cross-selling effect, and it even slightly decreases with larger cross-selling coefficients.(  4) Errors vs. Coefficient of Variation C V i , C V j. The two coefficients of variation, CVi and CVj, affect the estimation. Similarly, to be concise, we introduce a comprehensive coefficient of variation, cv=(CVi+CVj)/2, to discriminate different variation levels, as summarized in Table 10.Table 10
Comprehensive coefficient of variation.
Comprehensive coefficient of variationcv=(CVi+CVj)/2 Values of(CVi,CVj) 0.05 ( 0.05,0.05 ) 0.1 0.05,0.15  o r  ( 0.15,0.05 ) 0.15 0.05,0.25  o r  0.15,0.15  o r  ( 0.25,0.05 ) 0.2 ( 0.15,0.25 )  o r  ( 0.25,0.15 ) 0.25 ( 0.25,0.25 )For the different comprehensive coefficients of variation in Table10, the average relative errors of μi, μj, σi, and σj are presented in Figure 5.Figure 5
Estimation errors with the coefficient of variation (CV).Figure 6
Estimation errors of correlation coefficient.Figure5 presents that the estimation errors of μi and μj increase with the coefficient of variation, but even for a 20% variation both are still smaller than 1.5%. However, regarding the standard deviations, σi and σj, the increase in order quantity for alleviating censoring is still an important means to obtain a good estimation. Only when the comprehensive censoring level reaches more than5σ, the relative error becomes smaller than 5%.(4) Estimation Errors of Correlation Coefficient ( ρ ). We can see that the sample size is still a positive factor for the estimators of correlation coefficient, and it is evident that the more slightly the data is censored, the better the estimator will perform. We also observe that the error is smaller in negative correlations than in positive correlations on the overall trend. However, the estimation of ρ still reaches a high precision provided that the fill rate of the inventory is high enough. Remarkably, for the 6σ comprehensive censoring level, the error is no larger than 1.0% on average.Combining all the above results, it is revealed that for inventory systems with at least2σ-level safety stock (corresponding to 4σ comprehensive censoring level), the estimation generates high-precision outcomes for larger sample data. It is worth noting that 2σ-level safety stock, i.e., 95% fill rate, is a basic service level for daily operations of a store. Hence, our approach could be of high practical relevance.


## 6.2.1. Relative Errors and Censoring Levels

As mentioned in Section5.2, we will evaluate the performance of our estimation method by comparing the results with thepoint estimation of the data ξit,ξjt that is sampled from the original demand. We define the relative error of the estimated expectation of demand as(29)eμi=μ~i-μ~ξiμ~ξi×100%,eμj=μ~j-μ~ξjμ~ξj×100%,where μ~ξi,μ~ξj are the corresponding point estimations. Similarly, the relative errors of the standard deviation and the correlation coefficient are(30)eσi=σ~i-σ~ξiσ~ξi×100%,eσj=σ~j-σ~ξjσ~ξj×100%,eρ=ρ~-ρ~ξρ~ξ×100%.Moreover, the censoring level is discriminated by the deviation of order quantities over the corresponding demand expectation, i.e., the multiples ofσXi and σXj. To be concise, we adopt the sum of the two multiples to comprehensively discriminate different censoring levels. The results will be reported considering the seven comprehensive censoring levels as those in the first column of Table 8.Table 8
Comprehensive censoring levels.
Comprehensive censoring levels Multiples ofσXi and σXj 0 σ 0 σ X i ,  0 σ X j 1 σ 0 σ X i , 1 σ X j  o r  1 σ X i , 0 σ X j 2 σ 0 σ X i , 2 σ X j  o r  1 σ X i , 1 σ X j  o r  2 σ X i , 0 σ X j 3 σ 0 σ X i , 3 σ X j  o r  1 σ X i , 2 σ X j  o r  2 σ X i , 1 σ X j  o r  3 σ X i , 0 σ X j 4 σ 1 σ X i , 3 σ X j  o r  2 σ X i , 2 σ X j  o r  3 σ X i , 1 σ X j 5 σ 2 σ X i , 3 σ X j  o r  3 σ X i , 2 σ X j 6 σ 3 σ X i , 3 σ X j


## 6.2.2. Results

For the seven censoring levels, we report the estimation errors ofμi, μj, σi, and σj, under different sample sizes, correlation coefficients, cross-selling coefficients, and coefficients of variation, respectively.(1) Errors vs. Sample Size ( T ).Figure 2 presents the average error under different numbers of samples, i.e., T=50,100,300,500,700,900,1000.Figure 2
Variations in estimation errors with sample size.Figure 3
Variation in estimation errors with correlation coefficient.We observe from Figure2 that the demand expectation estimators yield very good estimation for the demand expectations of both items. Even for 50 sample points, the error is a maximum of 2.5%. With the standard deviation estimators, the error can reach more than 20% when the sample size is too small. However, if the sample size is larger than 800, the estimator performs much better, with an error no larger than 5%. We consider that for many modern business applications of inventory management, the data is no longer an obstruction. In particular, with the development of data acquisition and storage technologies, a large amount of data can be obtained from POS/MIS/MES, and big data systems [31]. The accuracy of our estimators can be further improved based on practical large sample data.Moreover, it is expected that the estimation error will become smaller when censoring is alleviated. The results show that if the comprehensive censoring level reaches4σ (implying that the safety stock of the two items is in average 2σ-level with a fill-rate of 95.45%), the error of the demand expectation estimator decreases from about 2.0% of 0σ to one no larger than 1.0%, and even to an error no larger than 0.5% with a 6σ comprehensive censoring level (corresponding to a 3σ-level safety stock with a fill-rate of 99.74%). We also see that, with high fill-rate, the standard deviation error decreases from larger than 15% of small sample data to smaller than 2% of large sample data.(  2) Errors vs. Correlation Coefficient ( ρ ).For the value of demand expectation, the estimation error slightly increases with the correlation coefficient, except for the censoring level of 0σ. However, almost all estimation errors of the mean demand are no larger than 1.4%, and they even decrease to 0.2% when censoring is sufficiently alleviated with a comprehensive censoring level no smaller than 4σ (corresponding to 2σ-level fill rate for each item in average), irrespective of the correlation coefficient.It can also be observed that, under a higher safety stock level (with comprehensive censoring level no smaller than4σ), the relative error of the standard deviation estimators is basically decreasing with the correlation coefficient. However, when censoring is aggravated (with censoring level no larger than 4σ), for weak correlation cases, the estimation error of the standard deviation reaches the maximal value, but it decreases when the correlation becomes stronger (irrespective of being positive or negative). The explanation is that, for a larger value of ρ, the estimator could capture more impact of the cross-selling effect (see Eqs. (9) and (11)), which might be beneficial to estimate the standard deviation of the demand.(  3) Errors vs. Cross-Selling Coefficient ( r i j , r j i ).The cross-selling is of dual-direction, so that for each one of the seven censoring levels we must consider two impact factors: rij and rji, which together will lead to a lengthy data report consisting of 28 two-dimensional tables. To be concise, we present a comprehensive index, r=(rij+rji)/2, to discriminate different levels of the cross-selling effect, as summarized in Table 9.Table 9
Comprehensive cross-selling coefficients.
Comprehensive cross-selling coefficientr=(rij+rji)/2 Values of(rij,rji) 0.1 ( 0.1,0.1 ) 0.2 0.1,0.3  o r  ( 0.3,0.1 ) 0.3 ( 0.1,0.5 )  o r  ( 0.5,0.1 ) 0.4 0.1,0.7  o r  0.3,0.5  o r  0.5,0.3  o r  ( 0.7,0.1 ) 0.5 ( 0.1,0.9 )  o r  ( 0.3,0.7 )  o r  ( 0.7,0.3 )  o r  ( 0.9,0.1 ) 0.6 ( 0.3,0.9 )  o r  ( 0.5,0.7 )  o r  ( 0.7,0.5 )  o r  ( 0.9,0.3 ) 0.7 ( 0.5,0.9 )  o r  ( 0.9,0.5 ) 0.8 ( 0.7,0.9 )  o r  ( 0.9,0.7 ) 0.9 ( 0.9,0.9 )For the different comprehensive cross-selling coefficients in Table9, the average relative errors of μi, μj, σi, and σj are presented in Figure 4.Figure 4
Variation in estimation errors with cross-selling effect.The estimation errors of the demand expectation increase with the cross-selling coefficient, but they do not exceed 4.50% for all censoring levels. It means that the estimation forμi and μj has a wide range of applicability for handling different levels of cross-selling. For the standard deviation, the error tends to increase with the cross-selling coefficient, but it is still smaller than 5% for censoring levels of 5σ and 6σ in most cases. In particular, for 6σ-level censoring, the relative error is not really increasing with the cross-selling effect, and it even slightly decreases with larger cross-selling coefficients.(  4) Errors vs. Coefficient of Variation C V i , C V j. The two coefficients of variation, CVi and CVj, affect the estimation. Similarly, to be concise, we introduce a comprehensive coefficient of variation, cv=(CVi+CVj)/2, to discriminate different variation levels, as summarized in Table 10.Table 10
Comprehensive coefficient of variation.
Comprehensive coefficient of variationcv=(CVi+CVj)/2 Values of(CVi,CVj) 0.05 ( 0.05,0.05 ) 0.1 0.05,0.15  o r  ( 0.15,0.05 ) 0.15 0.05,0.25  o r  0.15,0.15  o r  ( 0.25,0.05 ) 0.2 ( 0.15,0.25 )  o r  ( 0.25,0.15 ) 0.25 ( 0.25,0.25 )For the different comprehensive coefficients of variation in Table10, the average relative errors of μi, μj, σi, and σj are presented in Figure 5.Figure 5
Estimation errors with the coefficient of variation (CV).Figure 6
Estimation errors of correlation coefficient.Figure5 presents that the estimation errors of μi and μj increase with the coefficient of variation, but even for a 20% variation both are still smaller than 1.5%. However, regarding the standard deviations, σi and σj, the increase in order quantity for alleviating censoring is still an important means to obtain a good estimation. Only when the comprehensive censoring level reaches more than5σ, the relative error becomes smaller than 5%.(4) Estimation Errors of Correlation Coefficient ( ρ ). We can see that the sample size is still a positive factor for the estimators of correlation coefficient, and it is evident that the more slightly the data is censored, the better the estimator will perform. We also observe that the error is smaller in negative correlations than in positive correlations on the overall trend. However, the estimation of ρ still reaches a high precision provided that the fill rate of the inventory is high enough. Remarkably, for the 6σ comprehensive censoring level, the error is no larger than 1.0% on average.Combining all the above results, it is revealed that for inventory systems with at least2σ-level safety stock (corresponding to 4σ comprehensive censoring level), the estimation generates high-precision outcomes for larger sample data. It is worth noting that 2σ-level safety stock, i.e., 95% fill rate, is a basic service level for daily operations of a store. Hence, our approach could be of high practical relevance.


## 7. Conclusions

Demand distribution is key information for inventory decision making. Many estimators estimate the parameters of demand distribution for a single item based on the observed demand and considering unobservable lost sales. Methods focused on the demand estimation of multiple substitutable items have also been proposed. However, few studies consider the demand estimation for cross-selling items in inventory systems with lost sales, which leave a research gap in the extant literature.This study extends the problem to the case of two cross-selling items. In the new scenario, the demand/sales of the two items are intertwined, and the lost sales may not be observed. We proposed a demand estimation method based on an iterative framework, which can deal with the cross-selling effect and lost sales. Numerical computations show that our estimators perform very well if the sample size is large enough. The methods will be competent in a modern inventory system that has substantial amounts of data and pursues a high-level service fill rate. The contribution of this study is summarized as follows:(1) We propose an approach to estimate the jointly normally distributed demands for cross-selling items in inventory systems with lost sales. The approach can estimate the parameters of the original demand distribution from incomplete information that is a mixture of censored demand and cross sales.(2) A closed-form estimator for an individual item's demand is developed, so that the parameters, including the means, deviations, and correlation coefficient of the cross-selling items' demands, can be efficiently calculated based on an iterative framework.(3) The calculation of the proposed approach mainly utilizes the (censored) sales data, which should be easily implemented to provide demand information for inventory decisions in practical applications.Our work can be further extended considering the following research lines. First, the method has not considered the case of multiple items. Second, we can further investigate different forms of demand distribution, for instance, Poisson distribution and negative binomial distribution. Third, the current method assumes that the expectation and standard deviation of the demand distribution are stable, but in many cases, the demand varies over time. Adapting the method to deal with a time-varying demand for cross-selling items will lead to more extensive researches.



---

*Source: 7219326-2019-07-16.xml*