A portion of the disclosure of this patent document contains material which is subject to copyright protection. The copyright owner has no objection to the facsimile reproduction by anyone of the patent document or the patent disclosure, as it appears in the Patent and Trademark Office patent file or records, but otherwise reserves all copyright rights whatsoever.
Embodiments of the invention are generally related to recommender systems and collaborative filtering, and are particularly related to logistic matrix factorization of implicit feedback data, with application to media environments or streaming services.
Users of modern e-commerce services are often presented with a myriad of product choices. To enhance the user experience, recommender systems can be used in conjunction with such services to provide personalized recommendations. A recommender system can be content-based, such that it examines features directly associated with users and products, for example a user's age or the release date of a music album; or can be based on collaborative filtering, such that it examines users' past behaviors or feedback to predict how users might act in the future. However, explicit feedback data, such as user-provided ratings, can be difficult to obtain. Additionally, there is increased interest in the use of implicit feedback data, such as click-throughs, which can be collected faster and with greater scale, and without requiring users to provide an explicit indication of sentiment.
In accordance with an embodiment, described herein is a system and method for logistic matrix factorization of implicit feedback data, with application to media environments or streaming services. While users interact with an environment or service, for example a music streaming service, usage data reflecting implicit feedback can be collected in an observation matrix. A logistic function can be used to determine latent factors that indicate whether particular users are likely to prefer particular items. Exemplary use cases include providing personalized recommendations, such as personalized music recommendations, or generating playlists of popular artists.
As described above, recommender systems can be used in conjunction with e-commerce services to provide personalized recommendations. However, explicit feedback data, such as user-provided ratings, can be difficult to obtain. Additionally, there is increased interest in the use of implicit feedback data, such as click-throughs, which can be collected faster and with greater scale, and without requiring users to provide an explicit indication of sentiment.
In accordance with an embodiment, described herein is a system and method for logistic matrix factorization of implicit feedback data, with application to media environments or streaming services. While users interact with an environment or service, for example a music streaming service, usage data reflecting implicit feedback can be collected in an observation matrix. A logistic function can be used to determine latent factors that indicate whether particular users are likely to prefer particular items. Exemplary use cases include providing personalized recommendations, such as personalized music recommendations, or generating playlists of popular artists.
In accordance with an embodiment, each entry in the observation matrix represents the number of times that a particular user has interacted with a particular item. For example, if U=(u1, . . . , un) is a group of n users, and I=(i1, . . . , im) is a group of m items, then an observation matrix R=(rui)n×m can be determined wherein each entry rui∈R≥0 represents the number of times that user u has interacted with item i.
In accordance with an embodiment, the observation matrix can be factorized by lower-dimensional matrices Xn×f and Ym×f wherein f is a number of latent factors, the rows of X are latent factor vectors that represent a user's taste, and the columns of YT are latent factor vectors that represent an item's implicit characteristics. The user and item vectors can be initially configured with random noise, and an alternating gradient ascent procedure performed, including alternating between solving for the vectors while modifying the observation matrix according to a gradient. During each iteration, the system can solve for item vectors using solved user vectors, and solve for user vectors using solved item vectors. After completing a number of iterations, latent factors can be determined for users and items. The latent factors can then be used to determine recommendations, for example one or more items for a user with which they have not yet interacted, i.e., for which rui=0, but which they are likely to enjoy based on the latent factors.
In accordance with an embodiment, the process can be adapted for parallel processing and for scale, including sharding the observation matrix and vectors into partitions, working locally with each partition using a data collection and aggregation processor, such as an Apache Hadoop environment, and then joining or otherwise combining the results of each partition.
Media Environments
In accordance with an embodiment, a media environment or streaming service which enables users to select a media device to stream, play, or otherwise access a particular media content, can utilize logistic matrix factorization of implicit feedback data to enhance the user experience.
In accordance with an embodiment, the media device and media server can include physical device or computer hardware resources 104, 144, such as one or more processors (CPU), physical memory, network components, or other types of hardware resources. The media device can optionally include a user interface 106, which is adapted to display media options and to determine a user interaction or input, for example to download, stream or otherwise provide access to a corresponding particular media content item. The media device can further include a media application 108, together with an in-memory media content buffer 110, and media play logic 122, which controls the playback of media content received from the media server, for playing either at a requesting client device, or at a controlled client device. A connected media environment firmware or logic 120 enables the media device to participate within the media environment.
In accordance with an embodiment, the media server can be provided at a computer system, including an operating system 146 or other processing environment which supports execution of a media server application 150 that can be used, for example, to stream music, video, or other forms of media content to a media device. The media server can provide a service which enables a media device to communicate with and to receive media content from the media server, in the form of media-access requests. A media application interface 148 can receive requests from media devices, or from other systems, to retrieve media content from the media server. Databases 162 can include a context database that stores data associated with the presentation of media content by a particular media device; a media content database that stores media content such as music, songs, videos, movies, or other types of media content; and an advertisement database that stores advertisement or sponsor-directed content. A media streaming logic 151 can retrieve or otherwise provide access to media content items, in response to requests from media devices or other systems, and populate a media content buffer 181, at a media delivery module or stream service 153, with streams 182, 184, 186 of corresponding media content data, which can then be returned to the requesting device or controlled device.
In accordance with an embodiment, media devices and media servers can communicate with one another using a network, for example the Internet 190. A user 192 can interact 194 with the media application at a media device, and issue requests to access media content. For example, the user's selection of a particular media option can be communicated 196 to the media server, via its media application interface. In response, the media server can populate 204 its media content buffer with corresponding items of media or sponsor-directed content 206, for example as one or more streams of media content and/or advertisement or other sponsor-directed content. The media server can then communicate 208 the selected media content to the user's media device, or to the controlled device as appropriate, where optionally it can be buffered (206′) in the media content buffer, for subsequent playing at the device.
In accordance with an embodiment, a recommender logic 210 can provide recommendations to a user based on logistic matrix factorization of implicit feedback, as described in further detail below.
In accordance with an embodiment, while users 250, 252 use a media server within a regional media server environment, for example by interacting 254, 256 with a media application playlist or a search function to retrieve, play, stream, or otherwise request media content items 261, 262, 263, 264, usage data 266, 268 can be collected and provided to the data collection and aggregation processor. An observation matrix 282 can be determined, and processed by an implicit matrix factorization logic 284, optionally including sharding 285 and MapReduce 286 functionality as appropriate, for use in determining latent factors 288. The latent factors can be used to configure a recommender logic with recommendations 290, which can be provided or streamed 292 to users as appropriate, including for example one or more items for a user with which they have not yet interacted, but which they are likely to enjoy based on the latent factors.
The Spotify Connect environment is generally described herein as an illustrative example of a media environment or streaming service which can utilize logistic matrix factorization of implicit feedback data. In accordance with other embodiments the techniques described herein can also be used with other types of media environment or streaming service.
Collaborative Filtering
As described above, recommender systems can be used in conjunction with e-commerce services to provide personalized recommendations. However, explicit feedback data, such as user-provided ratings, can be difficult to obtain. Additionally, there is increased interest in the use of implicit feedback data, such as click-throughs, which can be collected faster and with greater scale, and without requiring users to provide an explicit indication of sentiment.
Collaborative filtering techniques can be used to examine users' past behaviors or feedback, including implicit feedback, to predict how users might act in the future. A context-aware model takes into account the context with which users provide implicit feedback, such as the time of day, day of the week, or location. A probabilistic model considers the probability of a user choosing to interact with an item, for example by modelling the distribution of an explicit ratings matrix, given user and item latent factor vectors, as being distributed according to a Gaussian or other distribution. Models may adopt a neighborhood-based approach, which generally focuses on the relationships between users or between items, and attempts to provide recommendations based on ratings of neighboring items by a same user. Alternatively, a latent factor approach attempts to characterize users and items based on multiple factors inferred from the usage patterns.
Matrix factorization models are a popular choice for recommender systems that include latent factor models, due to their simplicity, superior performance to neighborhood based methods, and ability to incorporate additional information.
In accordance with an embodiment, a user's past behavior can be analyzed in order to predict how they will act in the future. Implicit feedback can be provided in many forms, for example click-throughs, page views, or media streaming counts, but in all cases it can be assumed that there is a set of non-negative feedback values associated with each pair of users and items. If U=(u1, . . . , un) is a group of n users, and I=(i1, . . . , im) is a group of m items, then an observation matrix R=(rui)n×m can be determined wherein each entry rui∈R≥0 represents the number of times that user u has interacted with item i.
The rui values are not required to be integers, but instead can be any non-negative real values, to support contextual or temporal weighting of observations. For example, in the case of a music streaming service such as Spotify, explicit clicks by a user may be weighted by a higher weight than streams without explicit clicks. As another example, more recent user streams may be weighted higher than older streams, since a user's taste may change slightly over time.
In most cases the observation matrix R is a very sparse matrix, since most users only interact with a small number of items in I. For any entries rui where user u does not interact with item i the system can place 0's. Using this approach, a value of 0 does not necessarily imply that the user does not prefer the item, but could simply imply that the user does not know about the item. The system can then be used, for example, to find the top recommended items for each user for each item that they have not yet interacted with (i.e., where rui=0).
In accordance with an embodiment, the observation matrix can be factorized by lower-dimensional matrices Xn×f and Ym×f wherein f is a number of latent factors, the rows of X are latent factor vectors that represent a user's taste, and the columns of YT are latent factor vectors that represent an item's implicit characteristics.
In accordance with an embodiment, this relationship can be modeled by a logistic function. If lu,i denotes the event that user u has chosen to interact with item i (i.e., user u prefers item i), then the probability of this event occurring can be distributed according to a logistic function parameterized by the sum of the inner product of user and item latent factor vectors, and user and item biases:
In accordance with an embodiment, the βi and βj terms represent user and item biases which account for variation in behavior across both users and items. Some users may have a tendency to interact with a diverse assortment of items in I, while others may only interact with a small subset. Similarly, some items may be very popular and have a high expectation of being interacted with across a broad audience, while other items may be less popular and only apply to a niche group. The bias terms are latent factors associated with each user u∈U and item i∈I that offset the behavior and popularity biases. Given this formulation, the non-zero entries of the observation matrix rui≠0 represent positive observations, and the zero entries rui=0 represent negative observations. Additionally, a measure of confidence in the entries of R can be defined, where α is a tuning parameter, as:
c=αrui
Each nonzero element rui≠0 serves as c=αrui positive observations, and each zero element rui=0 serves as a single negative observation. Increasing α places more weight on the non-zero entries, while decreasing α places more weight on the zero entries. Choosing α to balance the positive and negative observations generally yields better results. In accordance with other embodiments, other confidence functions can replace c. To remove the power user bias that comes from a dataset where a small minority of users contribute the majority of the weight, a log scaling function can be used, such as:
c=1+α log(1+rui/ϵ)
By making the assumption that all entries of R are independent, the system can derive the likelihood of observations R given the parameters X, Y, and β as:
In accordance with an embodiment, zero-mean spherical Gaussian priors can be placed on the user and item latent factor vectors to regularize the model and avoid over-fitting to the training data:
Taking the log of the posterior, and replacing constant terms with a scaling parameter λ, the following can be determined:
Then, the process can attempt to learn the X, Y, and β that maximize the log posterior. A local maximum of the objective defined below can be found by performing an alternating gradient ascent procedure:
In accordance with an embodiment, in each iteration the system first fixes the user vectors X and biases β, and takes a step towards the gradient of the item vectors Y and biases β. Next, the system fixes the item vectors Y and biases β, and takes a step towards the gradient of the user vectors X and biases β. The partial derivatives for the user vectors and biases can be given by:
Each iteration is linear in the number of users |U| and items |I|, which for larger domains this can become a limitation. In situations where linear computation is not possible, then fewer negative samples (rui=0) can be sampled, together with a decreased α in response. This provides enough approximation that near-optimal vectors can be solved with much less computation time. Additionally, the number of iterations required for convergence can be greatly reduced by choosing the gradient step sizes adaptively. For example, if xtu denotes the value of xu at iteration t, and gtxu denotes the gradient of xu at iteration t, then at iteration t the following update to xu can be performed:
Each iteration of the alternating gradient descent procedure includes computing the gradient for all latent factor vectors, and then taking a step towards the positive direction of the gradient. Each of these gradients includes sum of a set of functions that each depend on a single user and item.
In accordance with an embodiment, the process can be adapted for parallel processing and for scale, including sharding the observation matrix and vectors into partitions, working locally with each partition using a data collection and aggregation processor, such as an Apache Hadoop environment, and then joining or otherwise combining the results of each partition.
In accordance with an embodiment, in the map phase, all observations rui, user vectors xu, and item vectors yi from the same block are provided to the same mapper. For each pair of users and items (u and i), the following is computed in parallel:
In accordance with an embodiment, in the reduce phase, the system can key off u (or off i if performing an item iteration), such that each vui and bui that map to the same user u (or to the same item i if performing an item iteration) are sent to the same reducer. It follows that ∂/∂xu=Σivui, and ∂/∂βu=Σibui and so these summations can be efficiently aggregated in parallel in the reduce phase. Finally, once the partial derivatives ∂/∂xu and ∂/∂βu have been computed, the system can update xu and βu.
At step 304, an observation matrix R is factorized by lower-dimensional matrices Xn×f and Ym׃, wherein f is a number of latent factors, and wherein the rows of X are latent factor vectors that represent a user's taste, and wherein the columns of YT are latent factor vectors that represent an item's implicit characteristics.
At step 306, for a particular user, recommended items are determined, for example items with which the user has not yet interacted, for which rui=0, but which they are likely to enjoy based on the latent factors.
At step 308, the method can optionally include parallelizing or scaling model training using, for example, sharding and a data collection and aggregation processor (e.g., an Apache Hadoop environment).
Example Results
The described approach was evaluated using a dataset of user listening behavior associated with a music streaming service and including usage data for |X|=50,000 active users, and the top |I|=10,000 most popular artists.
For each user, the number of times they listened to each artist in I was tracked over a period of time. A “listen” was defined as any continuous stream of an artist's song for more than 30 seconds, which removes the bias that results from users sampling artists that they end up not liking. For each user-item pair in U×I a tuple of the form (<user>, <artist>, <listens>) was created wherein <user>∈U, <artist>∈I, and rui=<listens>. In total, the observation matrix R contained 86 million non-zero entries and 414 million zero entries. Choosing α to balance the positive and negative observations yielded the best results.
A difficulty with evaluating collaborative filtering models with implicit feedback datasets is that the data only encodes positive feedback. If a user has not listened to a particular artist this does not necessarily mean that they dislike the artist, but could instead mean that they have never heard of the artist. However, if a user has listened to an artist a large number of times then this is a strong positive signal that the user likes that artist. Due to the lack of negative feedback, a recall-based evaluation metric Mean Percentage Ranking (MPR) can be used to evaluate a user's satisfaction with an ordered list of recommended items.
A uniform random subset of 10% of the entries in R was selected as a test set (Rtest) having entries rtui, and the system trained on the remaining 90%. Removing random entries in R, instead of random listens, avoids the bias that users will tend to revisit the same artists. By removing full entries, the dataset consists of user-item pairs for which the training set has no listening data, which more closely reflects the use case of predicting top new items for a user, rather than those top items the user has already interacted with.
For each user u∈U a ranked list of the items in I was generated and sorted by preference. A rankui denoted the percentile ranking of item i for user u, wherein a rankui=0% signified that i is predicted as the highest recommended item for u; while a rankui=100% signified that i is predicted as the lowest recommended item for u. The percentile ranking is then evenly distributed among the remaining items in the list by steps of 100%/|I|. The MPR expected percentile ranking of a user listening to the artists in the test set can be defined as:
Lower values of MPR are more desirable as they indicate that the user listened to artists higher in their predicted lists. Conversely, higher values of MPR indicate that users listened to artists lower in their predicted lists. Randomly produced lists would have an expected MPR of 50%.
The model was evaluated using the MPR evaluation metric for different numbers of latent factors ranging from 5 to 100. An Implicit Matrix Factorization (IMF) model was also implemented for the same set of varying latent factors. Both models were run until the change in MPR was minimal, amounting to around 10-20 iterations for IMF and 30-40 iterations for logistic matrix factorization, depending on the number of latent factors. A popularity baseline was also evaluated in which all ranked lists were based on the globally most popular artists over the training period, and which achieved an MPR of 14.9%, already an improvement over a purely random model.
Both logistic matrix factorization and IMF eventually converged to a similar MPR at around 100 latent factors (6.27 for IMF and 5.99 for logistic matrix factorization). The improved performance of logistic matrix factorization comes from its ability to outperform IMF with fewer number of latent factors. For example, using just 10 latent factors IMF achieved an MPR of 16.9% which is just above the popularity baseline; while logistic matrix factorization achieved an MPR of 8.065%. The ability to perform well using less latent factors is a desirable characteristic since, as the number of latent factors increases, more training data is needed to effectively learn near optimal vectors—in situations where there is insufficient amounts of training data, increasing the number of latent factors may not yield performance improvements. Additionally, some recommendation techniques, such as nearest-neighbor search require maintaining a data structure such as a k-d tree in memory, and the size of these data structures is proportional to the number of items as well as the dimensionality of the latent factors.
In addition to the quantitative evaluation described above, a qualitative analysis was performed. One of the benefits to using latent factor models is that it places all items in I into a lower-dimensional space where similarities can be computed between items. A cosine similarity can be used to measure the similarity between item vectors in a latent factor space, for example given two vectors a and b, the cosine similarity between the vectors is defined as:
The top related items for a given item a can be defined as the top items b∈I ranked by simab. Both logistic matrix factorization and IMF were trained using the same listening history dataset and 60 latent factors, and the top related items were computed for several artists in I using the learned item vectors. Qualitatively, the logistic matrix factorization was able to successfully capture a music sentiment, for use in choosing complimentary music and artists.
Example Software Implementation
In accordance with an embodiment, the present invention may be conveniently implemented using one or more conventional general purpose or specialized digital computer, computing device, machine, or microprocessor, including one or more processors, memory and/or computer readable storage media programmed according to the teachings of the present disclosure. Appropriate software coding can readily be prepared by skilled programmers based on the teachings of the present disclosure, as will be apparent to those skilled in the software art.
The software code illustrated in
In some embodiments, the present invention includes a computer program product which is a non-transitory storage medium or computer readable medium (media) having instructions stored thereon/in which can be used to program a computer to perform any of the processes of the present invention. Examples of the storage medium can include, but is not limited to, any type of disk including floppy disks, optical discs, DVD, CD-ROMs, microdrive, and magneto-optical disks, ROMs, RAMs, EPROMs, EEPROMs, DRAMs, VRAMs, flash memory devices, magnetic or optical cards, nanosystems (including molecular memory ICs), or any type of media or device suitable for storing instructions and/or data.
The foregoing description of embodiments of the present invention has been provided for the purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the precise forms disclosed. Many modifications and variations will be apparent to the practitioner skilled in the art.
The embodiments were chosen and described in order to best explain the principles of the invention and its practical application, thereby enabling others skilled in the art to understand the invention for various embodiments and with various modifications that are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the following claims and their equivalents.
This application claims the benefit of priority to U.S. Provisional Patent Application No. 61/947,262, titled “SYSTEM AND METHOD FOR LOGISTIC MATRIX FACTORIZATION FOR IMPLICIT FEEDBACK DATA”, filed Mar. 3, 2014; which is herein incorporated by reference.
Number | Name | Date | Kind |
---|---|---|---|
5872850 | Klein et al. | Feb 1999 | A |
7398274 | Ittycheriah et al. | Jul 2008 | B2 |
7627605 | Lamere et al. | Dec 2009 | B1 |
7797446 | Heller et al. | Sep 2010 | B2 |
7849092 | Slaney et al. | Dec 2010 | B2 |
8306976 | Handman et al. | Nov 2012 | B2 |
8370351 | Kalasapur et al. | Feb 2013 | B2 |
8407226 | Das | Mar 2013 | B1 |
8473368 | Gronow et al. | Jun 2013 | B2 |
8874574 | Purdy | Oct 2014 | B2 |
9110955 | Bernhardsson | Aug 2015 | B1 |
20020019858 | Kaiser et al. | Feb 2002 | A1 |
20060015193 | Kato | Jan 2006 | A1 |
20090034750 | Ayoub et al. | Feb 2009 | A1 |
20090037355 | Brave et al. | Feb 2009 | A1 |
20110153663 | Koren | Jun 2011 | A1 |
20110295843 | Ingrassia et al. | Dec 2011 | A1 |
20120030020 | Bucak | Feb 2012 | A1 |
20120059788 | Sekino | Mar 2012 | A1 |
20120078830 | Bodor | Mar 2012 | A1 |
20120303713 | Harbick et al. | Nov 2012 | A1 |
20130159243 | Wei et al. | Jun 2013 | A1 |
20130262469 | Whitman | Oct 2013 | A1 |
20130311163 | Somekh et al. | Nov 2013 | A1 |
20130332842 | Bernhardsson et al. | Dec 2013 | A1 |
20140129500 | Nice et al. | May 2014 | A1 |
20140214751 | Damera-Venkata | Jul 2014 | A1 |
20140280181 | Rodger et al. | Sep 2014 | A1 |
20150248618 | Johnson | Sep 2015 | A1 |
20160328409 | Ogle et al. | Nov 2016 | A1 |
20160335266 | Ogle et al. | Nov 2016 | A1 |
Number | Date | Country |
---|---|---|
2312516 | Apr 2011 | EP |
2012107762 | Aug 2012 | WO |
Entry |
---|
Wikipedia, “Collaborative filtering; https://en.wikipedia.org/wiki/Collaborative_filtering”, captured on Jul. 3, 2015; retrieved from http://web.archive.org/web/20150703051148/https://en.wikipedia.org/wiki/Collaborative_filtering on Feb. 6, 2019, 8 pages. |
Dato, “GraphLab Factorization Recommender; https://dato.com/products/create/docs/generated/graphlab.recommender.factorization_recommender.FactorizationRecommender.html” captured on Sep. 30, 2015; retrieved from “https://web.archive.org/web/20150508070324/http://dato.com/products/create/docs/generated/graphlab.recommender.factorization_recommender.FactorizationRecommender.html”, 3 pages. |
The Echo Nest, “Acoustic Attributes Overview; http://developer.echonest.com:80/acoustic-attributes.html”, captured on Oct. 27, 2015; retrieved from “https://web.archive.org/web/20151027062454/http://developer.echonest.com:80/acoustic-attributes.html” on Feb. 6, 2019, 2 pages. |
Google, “Word2vec; https://code.google.com/p/word2vec/”, captured on Jan. 22, 2016; retrieved from “https://web.archive.org/web/20160122041211/https://code.google.com/p/word2vec/” on Feb. 6, 2019, 5 pages. |
Hu, et al., “Collaborative Filtering for Implicit Feedback Datasets”, Data Mining 2008, ICDM Eighth IEEE International Conference, Dec. 15-19, 2008, pp. 263-272,10 pages. |
International Search Report in connection with PCT Application No. PCT/US2016/042821, dated Aug. 31, 2016, 3 pages. |
International Search Report and Written Opinion of the International Search Authority in connection with PCT Application No. PCT/US2016/042830, dated Sep. 6, 2016, 13 pages. |
Johnson, Christopher C., “Logistic Matrix Factorization for Implicit Feedback Data”, Conference on Neural Information Processing Systems 2014, 9 pages. |
Koren, et al., “Matrix Factorization Techniques for Recommender Systems”, Computer, Issue 8, IEEE Computer Society, Aug. 2009, pp. 42-49, 8 pages. |
O'Mahony, et al., “Detecting Noise in Recommender System Databases”, Proceedings of the 11th International Conference on Intelligent User Interfaces, Sydney, Australia, Jan. 29-Feb. 1, 2006, pp. 109-115, 7 pages. |
Rendle, Steffen, “Factorization Machines”, ICDM '10, Proceedings of the 2010 IEEE International Conference on Data Mining, Dec. 13, 2010, pp. 995-1000, 6 pages. |
Scikit-learn, “Generalized Linear Models—scikit-learn 0.17 documentation; http://scikit-learn.org:80/stable/modules/linear_model.html”, captured on Feb. 1, 2016; retrieved from “https://web.archive.org/web/20160201220706/http://scikit-learn.org:80/stable/modules/linear_model.html” on Feb. 6, 2019, 26 pages. |
United States Patent and Trademark Office, Office Action dated Dec. 22, 2016, in connection with U.S. Appl. No. 15/212,961, 18 pages. |
United States Patent and Trademark Office, Office Action dated Jun. 15, 2017, in connection with U.S. Appl. No. 15/212,961, 25 pages. |
United States Patent and Trademark Office, Office Action dated Sep. 26, 2018, in connection with U.S. Appl. No. 15/212,961, 13 pages. |
Number | Date | Country | |
---|---|---|---|
20150248618 A1 | Sep 2015 | US |
Number | Date | Country | |
---|---|---|---|
61947262 | Mar 2014 | US |