This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

An Entropy Law Governing Retail Shopper Decisions: Supplemental Material

Nicholas Sohre1 [email protected]    Alisdair Wallis2    Stephen J. Guy1 [email protected] University of Minnesota1, Computer Science & Engineering, Tesco PLC2

Appendix A Proofs

Lemma: Simulated Inversion rate of two items

Given two items having true geodesic distances bb and cc from the agent, their estimated distances b^\hat{b} and c^\hat{c} under our simulation model are

b^=b+wb,ϵb𝒩(0,αb)c^=c+wc,ϵc𝒩(0,αc)\begin{split}\hat{b}=b+w_{b},\epsilon_{b}\sim\mathcal{N}(0,\alpha*b)\\ \hat{c}=c+w_{c},\epsilon_{c}\sim\mathcal{N}(0,\alpha*c)\end{split} (1)

Which produce Gaussian random variables of the estimated distances B𝒩(b,αb)B\sim\mathcal{N}(b,\alpha*b) and C𝒩(c,αc)C\sim\mathcal{N}(c,\alpha*c). Suppose b<cb<c (that is, item bb is closer to the agent than item cc). Then, the probability of an inversion is the probability that the estimated distances swap in magnitude: C<BCB<0C<B\rightarrow C-B<0. Let Y=CBY=C-B be a new Gaussian random variable, then

Y𝒩(cb,(αc)2+(αb)2)=𝒩(cb,αc2+b2)\begin{split}Y\sim\mathcal{N}(c-b,\sqrt{(\alpha*c)^{2}+(\alpha*b)^{2}})\\ =\mathcal{N}(c-b,\alpha*\sqrt{c^{2}+b^{2}})\end{split} (2)

and the likelihood of inversion is P(Y<0)P(Y<0). For a given bb, cc, and α\alpha, this quantity can be computed analytically from the CDF of YY evaluated at 0:

P(Y<0)=12[1+erf(0(cb)α2(c2+b2))]P(Y<0)=\frac{1}{2}\left[1+\text{erf}\left(\frac{0-(c-b)}{\alpha*\sqrt{2(c^{2}+b^{2})}}\right)\right] (3)

Theorem: Simulated inversion chance increases with distance to the closer item

First, we note that 1+erf(x)1+\text{erf}(x) is a positive, increasing function of xx. Then, given equation 3, it suffices to show that for any bb (the distance to the closer item), the input to erf is increasing:

b[δδbbcα2(c2+b2)0], 0<bc\forall b\left[\frac{\delta}{\delta b}\,\frac{b-c}{\alpha*\sqrt{2(c^{2}+b^{2})}}\geq 0\right],\;0<b\leq c (4)

Proof: evaluating the partial derivative in equation 4 with respect to bb, we have

δδbbcα2(c2+b2)\displaystyle\frac{\delta}{\delta b}\,\frac{b-c}{\alpha*\sqrt{2(c^{2}+b^{2})}}
=1α2δδbbcc2+b2\displaystyle=\frac{1}{\alpha\sqrt{2}}\frac{\delta}{\delta b}\,\frac{b-c}{\sqrt{c^{2}+b^{2}}}
=1α2b2+c2(bc)bb2+c2(b2+c2)\displaystyle=\frac{1}{\alpha\sqrt{2}}\frac{\sqrt{b^{2}+c^{2}}-\frac{(b-c)b}{\sqrt{b^{2}+c^{2}}}}{(b^{2}+c^{2})}\;\; (quotient rule)
=1α2c(c+b)(b2+c2)32\displaystyle=\frac{1}{\alpha\sqrt{2}}\frac{c(c+b)}{(b^{2}+c^{2})^{\frac{3}{2}}}\;\; (simplify) (5)

Since equation 5 is positive whenever b,cb,c and α\alpha are positive, the derivative is positive with respect to bb and constrains equation 3 to be increasing with increasing bb.

Theorem: Simulated inversion chance decreases with increasing distance between items

Another important property for maintaining the relationship between difficulty and inversion chance is that the inversion chance must decrease with increasing W=cbW=c-b.
Proof: We wish to show that the derivative with respect to WW is always negative:

W=cb,[δδWbcα2(c2+b2)0], 0<bc\forall W=c-b,\left[\frac{\delta}{\delta W}\,\frac{b-c}{\alpha*\sqrt{2(c^{2}+b^{2})}}\leq 0\right],\;0<b\leq c (6)

Noting that (b2+c2)=W2+2cb=W2+2b(b+W)(b^{2}+c^{2})=W^{2}+2cb=W^{2}+2b(b+W), we can substitute into equation 6 and get

δδWWα2(W2+2b(b+W))\frac{\delta}{\delta W}\,\frac{-W}{\alpha*\sqrt{2(W^{2}+2b(b+W))}} (7)

While we cannot write the derivative in terms of only WW, we can treat bb as a positive constant and take the partial with respect to WW. If the result is negative for any value of bb, then the derivative with respect to WW is negative regardless of bb and the property is satisfied:

δδW[Wα2(W2+2b(b+W))]\displaystyle\frac{\delta}{\delta W}\,\left[-\frac{W}{\alpha*\sqrt{2(W^{2}+2b(b+W))}}\right]
=1α2δδW[WW2+2b(b+W)]\displaystyle=\frac{1}{\alpha\sqrt{2}}\frac{\delta}{\delta W}\,\left[-\frac{W}{\sqrt{W^{2}+2b(b+W)}}\right]
=1α2[W2+2b(b+W)W(W+b)W2+2b(b+W)W2+2b(b+W)]\displaystyle=\frac{1}{\alpha\sqrt{2}}\left[-\frac{\sqrt{W^{2}+2b(b+W)}-\frac{W(W+b)}{\sqrt{W^{2}+2b(b+W)}}}{W^{2}+2b(b+W)}\right]\;\; (quotient rule)
=1α2[W(W+b)(W2+2b(b+W))32W2+2b(b+W)(W2+2b(b+W))32]\displaystyle=\frac{1}{\alpha\sqrt{2}}\left[\frac{W(W+b)}{(W^{2}+2b(b+W))^{\frac{3}{2}}}-\frac{W^{2}+2b(b+W)}{(W^{2}+2b(b+W))^{\frac{3}{2}}}\right]\;\; (simplify)
=1α2b(2b+W)(W2+2b(b+W))32\displaystyle=\frac{1}{\alpha\sqrt{2}}\frac{-b(2b+W)}{(W^{2}+2b(b+W))^{\frac{3}{2}}}\;\; (simplify) (8)

Since WW and bb are both positive non-zero quantities, the resulting derivative is always negative as desired.

Theorem: Simulation inversion chance increases monotonically with difficulty

To show this property, it is sufficient to show that difficulty also increases monotonically with decreasing WW and increasing bb, since both these two values fully specify both the inversion rate (given an α\alpha) and the difficulty.
Proof: First, we note that WW and bb are sufficient to fully describe difficulty:

difficulty=log2(AW)=log2(b+WW)\text{{difficulty}}=log_{2}\left(\frac{A}{W}\right)=log_{2}\left(\frac{b+W}{W}\right) (9)

Then we can construct the partial derivatives with respect to both WW and bb for difficulty and see that they are always positive and negative respectively:

δδWAW=δδWb+WW=bW2\frac{\delta}{\delta W}\frac{A}{W}=\frac{\delta}{\delta W}\frac{b+W}{W}=\frac{-b}{W^{2}} (10)
δδbAW=δδbb+WW=b+WW2\frac{\delta}{\delta b}\frac{A}{W}=\frac{\delta}{\delta b}\frac{b+W}{W}=\frac{b+W}{W^{2}} (11)

Thus, both inversion rate and difficulty monotonically increase with increasing bb and decrease with increasing WW. Since WW and bb are both sufficient to fully specify both inversion rate and difficulty, we cannot make one smaller or larger (by adjusting WW or bb) without having the same effect on the other. Therefore, both these quantities will have a monotonic relationship with each other.

Appendix B Simulation Details

Input: itemsToRetrieve;
Output: itemOrder;
itemOrder = [];
alpha = 0.30;
while While itemsToRetrieve.length << 1 do
       distances = [];
       for i in 0 : itemsToRetrieve.length-1  do
             trueDistance = getGeodesicDistance(itemsToRetrieve[i]);
             noisyDistance = trueDistance + sampleNormal(0, alpha * trueDistance);
             distances.push(noisyDistance);
            
       end for
      itemsByEstimatedDist = sort(itemsToRetrieve, by = distances);
       itemOrder.push(itemsByEstimatedDist[0]);
       itemsToRetrieve.remove(itemsByEstimatedDist[0]);
      
end while
itemOrder.push(itemsToRetrieve[0]);
return itemOrder;
Algorithm 1 Basket Simulation