## A Hierarchical Bayesian Model for Size Recommendation in Fashion

Romain Guigourès romain.guigoures@zalando.de Zalando SE

Yuen King Ho yuen.king.ho@zalando.de Zalando SE

Evgenii Koriagin evgenii.koriagin@zalando.de Zalando SE

Abdul-Saboor Sheikh saboor.sheikh@zalando.de

Zalando Research

Urs Bergmann urs.bergmann@zalando.de Zalando Research

## ABSTRACT

We introduce a hierarchical Bayesian approach to tackle the challenging problem of size recommendation in e-commerce fashion. Our approach jointly models a size purchased by a customer, and its possible return event: 1. no return, 2. returned too small 3. returned too big. Those events are drawn following a multinomial distribution parameterized on the joint probability of each event, built following a hierarchy combining priors. Such a model allows us to incorporate extended domain expertise and article characteristics as prior knowledge, which in turn makes it possible for the underlying parameters to emerge thanks to sufficient data. Experiments are presented on real (anonymized) data from millions of customers along with a detailed discussion on the efficiency of such an approach within a large scale production system.

Reza Shirvany reza.shirvany@zalando.de Zalando SE

sizes themselves represent different physical measurements from one brand to another; 2. different ways of converting a local size system to another: in Europe, garment sizes are not standardized and brands don't always use the same conversion logic.

## KEYWORDS

Bayesian model; size recommendation; fashion; e-commerce

## ACMReference Format:

Romain Guigourès, Yuen King Ho, Evgenii Koriagin, Abdul-Saboor Sheikh, Urs Bergmann, and Reza Shirvany. 2018. A Hierarchical Bayesian Model for Size Recommendation in Fashion. In Proceedings of Twelfth ACM Conference on Recommender Systems (RecSys '18). ACM, New York, NY, USA, 5 pages. https://doi.org/10.1145/3240323.3240388

## 1 INTRODUCTION

Fashion is a way to express identity, moods and opinions. Customers also tend to use fashion to either emphasize certain parts of their body or hide others. In that context, size and fit have been shown to be among factors influencing the most the overall satisfaction [1]. Online customers have to buy before trying their clothes on. The sensory feedback phase about how the article fits via touch and visual cues is then delayed. Because of these uncertainties, a lot of consumers are still reluctant to engage in the purchase process.

To make matters worse, fashion articles including shoes and apparel have important sizing variations primarily due to: 1. different definitions of respective sizes from brands: the size systems used for specific categories are limited (e.g. S, M, L, etc.), however the

Permission to make digital or hard copies of part or all of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for third-party components of this work must be honored. For all other uses, contact the owner/author(s).

RecSys '18, October 2-7, 2018, Vancouver, BC, Canada

© 2018 Copyright held by the owner/author(s).

ACM ISBN 978-1-4503-5901-6/18/10.

https://doi.org/10.1145/3240323.3240388

A way to circumvent the confusion created by these variations is to use size tables which map physical body measurements to the article size system, requiring customers to have accurate measurements of their body. However, size tables themselves might suffer from a large variance, up to one inch within a single size. These differences stem from either different datasets used for size tables (e.g. German vs. UK population) or are due to vanity sizing, i.e deliberate size inconsistencies in brands targeting a specific focus group based on age, sportiness, etc. which represent major influences on the body measurements [2-4]. The combination of the above factors leaves the customers alone to face a highly challenging problem of determining the right size and fit during their purchase journey. In recent years, there has been a lot of interest in building recommender systems in fashion e-commerce with major focus on modeling style preferences based on customers past interactions, taste and affinities [5-7].ÂăHowever, few research work have been conducted to tackle the size recommendation problem.

The recommendation of size and fit has been recently studied in [8] where sparsity in purchased data is mentioned as a major issue, especially considering that articles have a limited stock. To minimize that problem, the authors propose to represent articles as a combination of brand, usage, size, and fit. A neural network is then trained to learn a latent vector describing each article defined as the combination of features mentioned before. Customer vector representation is obtained by aggregating over purchased articles and, finally, a gradient boosted classifier predicts the fit of an article to a customer.

Following a different approach, the authors of [9] propose a solution for determining if an article of a certain size would be fit, large, or small for a certain customer, using the purchase history. This is achieved by iteratively deducing the true sizes for customers and products, fitting a linear function based on the difference in sizes, and performing ordinal regression on the output of the function to get the loss. Extra features are simply included by addition to the linear function. To handle multiple persons behind a single account, hierarchical clustering is performed on each customer account before doing the above. An extension of that work has been very recently published proposing a Bayesian approach on a similar model [10]. Instead of learning the parameters in an iterative process, the updates are done with mean-field variational inference with Polya-Gamma augmentation. This method therefore

naturally benefits from the nice advantages of Bayesian modeling - the uncertainty outputs, and the use of priors. In this paper, we present two approaches: a baseline algorithm, launched on shoes in 2016 and on garments in 2017, which consists in inferring the size that a customer intends to buy and, independently, the article's sizing characteristics ; and a hierarchical Bayesian approach which aims at jointly modeling the purchases of one or multiple sizes of an article along with their possible return events: 1. no return (article is kept), 2. returned too small 3. returned too big.

In the size recommendation context, data sparsity is severe since it affects both articles and customers. To tackle this, we propose two design choices: a) building a hierarchy on top of parameters, exploiting prior knowledge on articles and customers; b) virtually treating the size as a continuous variable in the training phase.

## 2 METHODOLOGY

Customers experience in fashion e-commerce consists in selecting an article in a desired size, trying the article, forming an opinion on its size and returning or keeping it. To simulate the customers behavior towards sizing, we model the joint probability of a customer to pick a size and the resulting return status. Return status is described by three possible events: the customer keeps the article, the customers returns the article because it's too small and the customer returns the article because it is too big. We ignore the cases where the customer returns the article for any other reason.

## 2.1 Notation

Let us denote C the set of customers and A the set of articles. The size S i is a continuous random variable. The variable R indicates the return status described above. Orders O are defined by a customer, an article, the purchased size and the return status. Both approaches introduced in the paper model the joint probability p (S o , R o | C o , A ) o as detailed in the following.

## 2.2 Baseline Model

The baseline model makes a simplifying assumption that the size the customer chose and the return status are two independent events. Thus, the joint probability is defined as the product of the probability over sizes and the probability of return status:

$$p ( \mathcal { S }, \mathcal { R } \, | \, \mathcal { C }, \mathcal { A } ) = p ( \mathcal { S } \, | \, \mathcal { C }, \mathcal { A } ) p ( \mathcal { R } \, | \, \mathcal { C }, \mathcal { A } ) \quad \quad ( 1 ) \quad \text{be}$$

Probability over sizes. We assume that multiple persons can use a single account. The probability distribution over sizes is obtained by Gaussian Kernel Density Estimation. To avoid getting degenerate distributions when customers always purchase the same size, we set a lower limit on the variances. Let denote O j the set of n j orders and S j = { s i , i = 1 .. n j } the set of n j sizes purchased by the customer c j . The related probability density function is defined as:

$$p ( s \ | \ c _ { j } ) = \frac { 1 } { n _ { j } h _ { j } } \sum _ { i = 1 } ^ { n _ { j } } \phi \left ( \frac { s - s _ { i } } { h _ { j } } \right ) \quad \quad ( 2 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where ϕ is the normal density function and h j is the bandwidth parameter for that specific customer c j . The latter is obtained by minimizing the mean integrated squared error [11].

Probability over return status. Customers have different return behaviors. However, impact of customers on returns is assumed negligible compared to potential sizing issues of the article. The probability of each return status is consequently marginalized over customers: p (R | C , A) = p (R | A) . The probability over return status is the empirical distribution over the three possible events: article is kept, too big or too small. The case of one of the events being not observed in the training data may lead to a null probability in validation. To avoid that problem, we add one to the counts of each event. For an article a i , sold n i times, the probability of a return event r , observed n i , r times in the data, is defined as:

$$p ( r \ | \ a _ { i } ) = \frac { n _ { i, r } + 1 } { n _ { i } + 3 } \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ }$$

Though, this method seems inelegant, it has a Bayesian grounding. Indeed, this is equivalent to taking the maximum a posteriori of a categorical distribution with a Dirichlet conjugate prior, which concentration parameter is equal to one, i.e the uniform distribution. In case of a cold start, i.e if a customer (resp. an article) is new, the marginal distribution over sizes of all the customers (resp. over return status of all articles) is used.

## 2.3 Hierarchical Bayesian Model

The baseline described above has a risk of specious parameter estimation and overfitting. Bayesian approaches conversely aim at providing a probability of the estimated parameters given a set of observed data, supporting the decision process when offering a recommendation to the customer. Therefore, using a Bayesian approach, we aim at modeling the joint probability of a size to be purchased and a return status to be observed without the simplifying hypothesis from the base approach. For each pair of customer and article, orders O are drawn following a categorical distribution.

$$\mathcal { O } \sim C a t ( p ( \mathcal { S }, \mathcal { R } \, | \, \mathcal { C }, \mathcal { A } ) )$$

Contrary to the baseline, both S and R are not assumed independent, instead the joint probability is factorized as:

$$p ( \mathcal { S }, \mathcal { R } \ | \ \mathcal { C }, \mathcal { A } ) = p ( \mathcal { S } \ | \ \mathcal { R }, \mathcal { C }, \mathcal { A } ) \times p ( \mathcal { R } \ | \ \mathcal { C }, \mathcal { A } ) \quad \ \ ( 5 )$$

It is worth noting that methods described in [8, 10] aim at modeling p (R | S , C A) , . Doing so requires discretizing the continuous variable S leading to an increase in the number of parameters to be inferred and making the model more susceptible to the sparsity of the data. The factorization we chose in Equation 5 allows us to model p (S | R , C A) , as a continuous distribution. This enables us to learn a smaller set of parameters specifying the distribution over all sizes, which helps to alleviate part of the sparsity problem.

Probability over return status. For the same reasons as explained in the Section 2.2, the probability of return status is marginalized over customers: p (R | C A) , = p (R | A) . Returns are assumed independent from one another, allowing us to model them using a categorical distribution. As the number of purchases might be low for some articles, a Dirichlet prior is used. The concentration parameter of the prior is based on the counts at the brand level and at the category level (e.g. dresses, t-shirts, sneakers, etc.). Let n K , n S , and n B denote the counts of kept articles, returned articles for being too small and too big at the article level respectively. In

a similar way, m K , m S and m B indicates the counts at the brand level and m ′ , m ′ and m fi B , at the category level.

K

S

$$p ( \mathcal { R } \, | \, \mathcal { A } ) \sim D i r i c h l e t ( \alpha )$$

with α = w · [ m K , m S , m B ] + w fi · [ m ′ K , m ′ S , m ′ B ] . The weights w and w ′ are learned under the assumption that they follow Beta distribution with a low first shape parameter and the second shape parameter equal to 1, in order to favor low weight values.

$$w \sim B e t a ( 0. 5, 1 ) \text{ and } w ^ { \prime } \sim B e t a ( 0. 1, 1 ) \text{ \quad \ \ } ( 7 ) \text{ \quad on}$$

Since the Dirichlet prior is the conjugate of the categorical distribution, the posterior probability can be analytically computed, easing the inference of the parameters of the model.

$$P ( \mathcal { R } = r \ | \ \mathcal { A }, \mathcal { O } ) = \frac { n _ { r } + \alpha _ { r } } { \sum _ { \substack { i \in \\ \{ K, S, B \} \\ \{ K, S, B \} } } n _ { i } + \sum _ { \substack { i \in \\ \{ K, S, B \} \\ \{ K, S, B \} \\ \} } } \quad \text{We} \\ P ( \mathcal { R } = r \ | \ \mathcal { A }, \mathcal { O } ) = \frac { \text{We} } { \text{so} } 1 \\ \text{Big} \\ \text{me} ^ { \prime } }$$

Probability over sizes conditionally on the return status. For a given customer and article, the probability distribution of the customer buying a size is a mixture of Gaussians. Since the number of users of an account is unknown, we decide to use an infinite mixture model. However, we assume that the number of distinct persons using a single account is low. That's why we opt for a Dirichlet process with a truncation level fixed to four. In order to ease the inference, we use a truncated stick-breaking process [12].

$$\L t e n i \text{ inference, we use a truncate} \, \text{suck-breaking process} \, [ \text{12} ]. \\ p ( s \, | \, r, c, a ) = \sum _ { i = 1 } ^ { 4 } \pi _ { i } \phi ( s \, | \, \mu _ { i }, \sigma _ { i } ^ { 2 } ) \\ \pi _ { i } = b _ { i } \prod _ { j = 1 } ^ { i - 1 } ( 1 - b _ { j } ) \\ b _ { i } \sim B e t a ( 1, \alpha ) \, \text{for} \, i = 1.. 3 \, \text{and} \, b _ { 4 } = 1 \\ \text{namer} \, \pi \, \text{is the mixin} \sigma \, \text{nronortion} \, \text{that can be interpreted}$$

The parameter π is the mixing proportion, that can be interpreted as the probability of person i using the account of the customer c . The shape parameter α of the Beta distribution is acting as a concentration parameter of the Dirichlet process: the case α = 1 is equivalent to the uniform distribution, thus favoring the scenario of multiple persons sharing a single account; conversely, the case α → 0 produces high density around 1 which favors the scenario of a single person placing all the orders. In the context of size recommendation, α is fixed at 0 5. .

The function ϕ in equation 9 is the normal probability density function over sizes for the person i using the account of the customer c , buying an article a , resulting in a return status r . The parameter µ is a combination of three parameters µ = µ C + µ A + η R : ·

- the average size of a person µ C ,

$$\mu _ { C } \sim \mathcal { N } ( \mu _ { 0 }, \sigma _ { 0 } ^ { 2 } ) \quad \quad ( 1 0 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where hyperparameters µ 0 and σ 2 0 depend on the category, the gender of the article and the size system;

- · the average offset of the article µ A ,

$$\mu _ { A } \sim \mathcal { N } ( 0, 1 ) \quad \quad \quad ( 1 1 )$$

where assumption is made that most articles have an accurate size, i.e. an offset of zero;

- · a shifting parameter η = { η K , η S , η B } for each return status: resp. article is kept, returned too small and returned too big,

$$\eta _ { K } = 0 \ ; \, \eta _ { S } \sim \mathcal { N } ( - 1, 1 ) \ ; \, \eta _ { B } \sim \mathcal { N } ( 1, 1 ) \quad \quad ( 1 2 )$$

shifting parameter is fixed at 0 for the case the article is kept, while it's sampled using a Gaussian distribution centered to 1 (resp. -1) when the customer has returned the article because it is too big (resp. too small).

We assume the parameter σ in equation 9 depends on the customer only. It is sampled following an Inverse Gamma distribution, with the shape parameter γ 1 and the scale parameter γ 2.

$$\sigma _ { C } ^ { 2 } \sim \Gamma ^ { - 1 } ( \gamma _ { 1 }, \gamma _ { 2 } )$$

We fix the parameters of the distribution to γ 1 = 1 and γ 2 = 2, so that the mode of the Inverse Gamma distribution is equal to 1. Figure 1 represents the graphical model of the approach.

Inference. Monte-Carlo Markov Chains are popular sampling methods for Bayesian inference. But those approaches are often slow to converge. Variational inference methods run faster than sampling based methods but also introduce an approximation bias that may lead to a bad estimation of the parameters. However, those approaches have demonstrated reasonable performances on Dirichlet processes [13] and are well suited for problems involving large amount of data. The inference is consequently done using mean-field approximation.

Figure 1: Graphical model of the Bayesian approach

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00825v1-figure-93f75140-2d58-43e6-a1ad-9470e03dc004.jpg)

## 2.4 Providing a Size Recommendation

The set of sizes of an article is a finite set and as a consequence the probability density function needs to be discretized. For a customer c and an article a with a set of k sizes S = { s i , i = 1 .. k } , the probability over sizes is discretized as follows:

$$p ( s, r ) = p ( r ) \frac { \int _ { s - \frac { 1 } { 2 } \epsilon } ^ { s + \frac { 1 } { 2 } \epsilon } f ( x ) d x } { \int _ { s _ { 1 } - \frac { 1 } { 2 } \epsilon } ^ { s _ { k } + \frac { 1 } { 2 } \epsilon } f ( x ) d x }$$

where f is the probability over sizes, marginal for the baseline and conditionally to the return status for the Bayesian model ; /uni03F5 is equal to the step between two sizes. To provide a size recommendation, we choose the size having the highest probability to be kept by the customer.

## 3 EXPERIMENTS

For this experiment, anonymized purchase data is collected for adult shoes. The data consists of 14.5 million purchases, 3 million distinct customers and 73,000 distinct articles. As the data has a strong temporal component, cross validation is performed under the following conditions [14]: 1. validation data occurs later than training data, 2. a period of three weeks - corresponding to the time to collect most returns from customers - is ignored between train set and validation set and 3. validation sets must not overlap. Both models are trained and cross-validated on the same data. We reported the average logarithm of the joint probability p s ( , r ) over all observations in the validation set in Table 1. Higher numbers show better performances of the model.

Table 1: Average log joint probability

|                         | Baseline        | Bayesian        |
|-------------------------|-----------------|-----------------|
| incl. unknown customers | -2.85 (Âś 0.15) | -2.35 (Âś 0.11) |
| excl. unknown customers | -3.32 (Âś 0.26) | -1.83 (Âś 0.19) |

Table 1 compares likelihood for both approaches, including and excluding customers not observed in the training data. Bayesian approach shows better results in both cases, where results are the best when excluding unknown customers. Conversely, the baseline performs better when unknown customers are included. This is mainly due to the fact that the baseline overfits by putting high probability density on the events observed in the training set, and very low density on unseen events.

In the context of size recommendation, two indicators play a key role in the decision process: coverage and accuracy. The coverage is the percentage of purchases for which the algorithm is confident making a decision. The accuracy is the number of correctly predicted sizes and return status, over all predictions. Figure 2 (Top) shows the accuracy versus the coverage for several values of a sliding threshold on the joint probability for both models. On Figure 2 (Bottom), the results are presented for the Bayesian model where we also include a threshold on the posterior probability of parameters.

Figure 2 (Top) shows that accuracy decreases as coverage increases, when changing the threshold on the joint probability. The performances of the Baseline and the Bayesian model are similar, when the decision is based on the value of the joint probability. In Figure 2 (Bottom), by putting a threshold on the posterior probability of parameters, we prevent the model from recommending article sizes to customers, where the parameters are poorly estimated. Performances are slightly better, however only 13% of the purchases can be covered.

It is worth mentioning that both approaches need to filter out a lot of purchases before starting to show reasonable accuracy levels. From the computational complexity point of view, inference of the Bayesian approach is more costly compared to the baseline model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00825v1-figure-48914bdb-7466-4b97-9cad-c6c31a9c4b97.jpg)

Accuracy

Figure 2: Accuracy versus coverage. (Top) Baseline and Bayesian model. (Bottom) Bayesian model with and without threshold on the posterior probability of the parameters.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00825v1-figure-2412e6e0-b9d0-460f-bc75-8da117eba641.jpg)

Results are encouraging and demonstrate the complexity of the size recommender topic, motivating a deeper research work in the field.

## 4 CONCLUSION

A hierarchical Bayesian approach was proposed to tackle the challenging problem of size recommendation in e-commerce fashion. The size purchased by a customer and its possible return events were jointly modeled thanks to a Bayesian approach. Experimental results were presented on real (anonymized) data from millions of customers along with a detailed discussion and comparison with a baseline approach with simplified hypothesis. It was shown that the Bayesian approach outperforms the baseline approach while providing better theoretical framework for gaining deeper understanding of the predictions. Future work consists in exploring different approaches to learn the joint probability, making use of additional article related data from fashion industry, along with deeper dives in the segmentation of customers and articles.

A Hierarchical Bayesian Model for Size Recommendation in Fashion

## REFERENCES

- [1] Gina Pisut and Lenda Jo Connell. Fit preferences of female consumers in the usa. Journal of Fashion Marketing and Management: An International Journal , 11(3):366-379, 2007.
- [2] Darko Ujević, Lajos Szirovicza, and Isak Karabegović. Anthropometry and the comparison of garment size systems in some european countries. Collegium antropologicum , 29(1):71-78, 2005.
- [3] Su-Jeong Hwang Shin and Cynthia L Istook. The importance of understanding the shape of diverse ethnic female consumers for developing jeans sizing systems. International Journal of Consumer Studies , 31(2):135-143, 2007.
- [4] Marie-Eve Faust and Serge Carrier. Designing Apparel for Consumers: The Impact of Body Shape and Size . Woodhead Publishing, 2014.
- [5] Yang Hu, Xi Yi, and Larry S Davis. Collaborative fashion recommendation: A functional tensor factorization approach. In Proceedings of the 23rd ACM international conference on Multimedia , pages 129-138. ACM, 2015.
- [6] Sagar Arora and Deepak Warrier. Decoding fashion contexts using word embeddings. In KDD Workshop on Machine learning meets fashion , 2016.
- [7] Christian Bracher, Sebastian Heinz, and Roland Vollgraf. Fashion dna: Merging content and sales data for recommendation and article mapping. arXiv preprint arXiv:1609.02489 , 2016.
- [8] GMohammedAbdulla and Sumit Borar. Size recommendation system for fashion e-commerce.
- [9] Vivek Sembium, Rajeev Rastogi, Atul Saroop, and Srujana Merugu. Recommending product sizes to customers. In Proceedings of the Eleventh ACM Conference on Recommender Systems , pages 243-250. ACM, 2017.
- [10] Vivek Sembium, Rajeev Rastogi, Lavanya Tekumalla, and Atul Saroop. Bayesian models for product size recommendations. In Proceedings of the 2018 World Wide Web Conference , WWW'18, pages 679-687, 2018.
- [11] Bernard W Silverman. Density estimation for statistics and data analysis . Routledge, 2018.
- [12] Jayaram Sethuraman. A constructive definition of dirichlet priors. Statistica sinica , pages 639-650, 1994.
- [13] David M Blei, Michael I Jordan, et al. Variational inference for dirichlet process mixtures. Bayesian analysis , 1(1):121-143, 2006.
- [14] Sylvain Arlot, Alain Celisse, et al. A survey of cross-validation procedures for model selection. Statistics surveys , 4:40-79, 2010.
Romain Guigourès Yuen King Ho Evgenii Koriagin Abdul-Saboor Sheikh Urs Bergmann Reza Shirvany
A Hierarchical Bayesian Model for Size Recommendation in Fashion
We introduce a hierarchical Bayesian approach to tackle the challenging
problem of size recommendation in e-commerce fashion. Our approach jointly
models a size purchased by a customer, and its possible return event: 1. no
return, 2. returned too small 3. returned too big. Those events are drawn
following a multinomial distribution parameterized on the joint probability of
each event, built following a hierarchy combining priors. Such a model allows
us to incorporate extended domain expertise and article characteristics as
prior knowledge, which in turn makes it possible for the underlying parameters
to emerge thanks to sufficient data. Experiments are presented on real
(anonymized) data from millions of customers along with a detailed discussion
on the efficiency of such an approach within a large scale production system.