Unveiling Hidden Implicit Similarities For Cross-Domain Recommendation

This article has been accepted for publication in a future issue of this journal, but has not been
fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TKDE.2019.2923904, IEEE
Transactions on Knowledge and Data Engineering
1
Unveiling Hidden Implicit Similarities for

Cross-Domain Recommendation
Quan Do, Wei Liu, Member, IEEE, Jin Fan, Dacheng Tao, Fellow, IEEE
Abstract—E-commerce businesses are increasingly dependent on recommendation systems to introduce personalized services and
products to targeted customers. Providing effective recommendations requires sufficient knowledge about user preferences and
product (item) characteristics. Given the current abundance of available data across domains, achieving a thorough understanding of
the relationship between users and items can bring in more collaborative filtering power and lead to a higher recommendation accuracy.
However, how to effectively utilize different types of knowledge obtained across domains is still a challenging problem. In this paper, we
propose to discover both explicit and implicit similarities from latent factors across domains based on matrix tri-factorization. In our
research, common factors in a shared dimension (users or items) of two coupled matrices are discovered, while at the same time,
domain-specific factors of the shared dimension are also preserved. We will show that such preservation of both common and
domain-specific factors are significantly beneficial to cross-domain recommendations. Moreover, on the non-shared dimension, we
propose to use the middle matrix of the tri-factorization to match the unique factors, and align the matched unique factors to transfer
cross-domain implicit similarities and thus further improve the recommendation. This research is the first that proposes the transfer of
knowledge across the non-shared (non-coupled) dimensions. Validated on real-world datasets, our approach outperforms existing
algorithms by more than two times in term of recommendation accuracy. These empirical results illustrate the potential of utilizing both
explicit and implicit similarities for making across-domain recommendations.
Index Terms—Cross-Domain Learning; Recommendation System; Matrix Factorization.
1 I NTRODUCTION require adequate ratings to analyze the latent similarities of

the users and the items [4]. This thorough understanding is
Product or content providers usually offer a wide variety
a foundation for achieving reliable recommendation results.
of products. On the one hand, this massive product selec-
The requirement of having sufficient ratings leads to two
tion meets a wide range of different consumer needs and
problems. The first is that newly established businesses
tastes. On the other hand, browsing over a long product
need time to collect ratings. The second one is the question
list is not a user-friendly task for any consumer to choose
raised after they have enough data: once they can achieve a
products of her preferences. Effective matching between
reliable recommendation, will their performance be further
product properties and consumer interests enhances user
improved if they have additional knowledge from external
satisfaction and eventually leads to more sale. Realizing that
data sources? These two problems are motivations for us to
recommendation systems can be a competitive advantage,
conduct this research.
lots of providers have been employing them to analyze
The above problems are addressable by joint analy-
users’ past behaviors to provide personalized product rec-
sis of information coming from different sources [5], [6],
ommendations [1].
[7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18].
Recommendation systems have gained their importance
Recent revolutions on the Internet have opened a great
and popularity among product providers. Two key tech-
opportunity for this joint analysis by making many public
niques are widely chosen for developing recommendation
and highly correlated datasets available. Finding a closely
systems: content-based approach [2] and collaborative filter-
related dataset from another domain is easier than ever. For
ing (CF) based approach [3]. The former focuses on informa-
example, some users buy mobile phones on Amazon, and
tion of users or items for making recommendations whereas
some other do so at Walmart; they may provide feedback on
the latter is based on latent similarities of the user interests
how much they like them. These ratings from the different
and those of the item characteristics for predicting items
users (i.e., Amazon and Walmart users) for the same set of
particular users would be interested. This paper focuses on
items can be formed two rating matrices coupled in their
CF based approach.
item dimensions. Thus, they can be jointly used to better
Even though CF based methods have shown their ef- understand user preferences and item characteristics. How
fectiveness in providing accurate recommendations, they to joint analysis of this rich information across domains to
improve recommendation performance is still a question
Q. Do and W. Liu are with Advanced Analytics Institute, University of
Technology Sydney, Australia. e-mail: {Quan.Do, Wei.Liu}@uts.edu.au. that needs to be better addressed.
J. Fan is with the College of Computer Science and Technology, Hangzhou Coupled datasets across domains possess explicit sim-
Dianzi University, China. e-mail: FanJin@hdu.edu.cn. ilarities. In case of the aforementioned rating matrices on
D. Tao is with the UBTECH Sydney Artificial Intelligence Center and Amazon and Walmart, they both contain user preferences
the School of Information Technologies, the Faculty of Engineering and
Information Technologies, the University of Sydney, Australia. e-mail: for the same set of items. Thus, explicit features in the
Dacheng.Tao@Sydney.edu.au. identical item dimension are conventionally used to do
1041-4347 (c) 2019 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/TKDE.2019.2923904, IEEE
2
High income Low crime

LGA rate LGA
Low income High crime

LGA rate LGA
(a) Percentage of high income family per LGA (b) Number of “break and enter dwelling” incidents per LGA
Fig. 1: (Best viewed in color) Implicit correlations between income and crime rate. Local government areas (LGA) in New South Wales (NSW)
state with more high-income families have a lower “break and enter dwelling” incidents. Data in a) is from Australian Bureau of Statistics and b)
is from NSW Bureau of Crime Statistics and Research.
coupled learning between datasets [19], [20], [21]. Besides factor for both datasets, is that cross-domain datasets have
these explicit similarities, we hypothesize that cross-domain their unique knowledge of their domains. Thus, HISF as-
datasets also have other implicit correlations in their re- sumes the coupled factors not to be equal, but they contain
maining dimensions. Our intuition comes from consumer common parts, which are shared between datasets across
segmentation which is a popular concept in marketing. domains, and domain-specific parts, which are unique for
Consumers can be grouped by their behaviors, for example, each domain.
“tech leaders” is a group of consumers with a strong desire Also, our another contribution is with implicit similar-
to own new smartphones with the latest technologies or ities. The fact that non-coupled dimensions in the afore-
another user group who only changes a new phone when mentioned example contain non-overlapping users prevents
the current one is out of order. Although users on Amazon a direct knowledge sharing in their non-coupled factors.
and Walmart are different, users sharing similar interests However, their latent behaviors are correlated and should
may share similar behaviors on the same set of products. be shared. These latent behaviors can be captured in low-
We think that this type of correlation also implicitly exists rank factors by matrix tri-factorization. As factorization is
in many other scenarios in the real world, for example, equivalent to spectral clustering [25], different users with
local suburbs with more high-income families may lead to a similar preferences are grouped in non-coupled user factors.
lower crime rate as shown in Fig. 1, or users with interest on Developing on this concept, we hypothesize that latent
action books may also like related movie genres (e.g., action clusters in these non-coupled factors may have a close
films). Thus, using both explicit and implicit similarities relationship. Therefore, we align correlated clusters in non-
properly will have a great potential to be applied to many coupled factors to be as close as possible. This idea matches
applications. the fundamental concept of CF in the sense that similar user
Different approaches have been proposed to perform a groups who rate similarly will continue to do so.
joint analysis of multiple datasets [22]. However, all of the In short, our key contributions are:
existing algorithms use explicit similarities as a bridge to 1) Sharing common latent variables while preserving
collaborate among datasets. The popular Collective Matrix unique ones in the coupled factors (Sect. 4.1): We propose
Factorization (CMF) [19] jointly analyzes datasets by as- coupled factors to include both common and their own
suming them to have an identical low-rank factor in their unique parts (e.g., common V(0) and domain-specific V(1) ,
coupled dimension. In this case, the shared identical factor V(2) in Figure 2.) This model better captures the true explicit
captures the explicit similarities across domains. Li et al. characteristics of cross-domain datasets.
[23] suggest correlated datasets to share explicit hidden 2) Aligning implicit similarities in non-coupled factors
rating patterns. The similarities between rating patterns are (Sect. 4.2): We introduce a method to leverage implicit
then transfered from one to another dataset. Gao et al. [24] similarities. HISF is the first factorization algorithm that
extend Li’s idea to include unique patterns in each dataset. utilizes both explicit and implicit similarities across domains
Weike et al. [4] regularizes factors of the target user-by- for recommendation accuracy improvement.
item matrix with those, called principal coordinates, from 3) Introducing an algorithm that optimizes all factors
user profile and item information matrices. Although these (Sect. 4.3): All factors are optimized by an algorithm follow-
explicit similarities showed their effectiveness in improving ing alternating least squared approach (Algorithm 1). We
recommendation, there are still rich implicit features that test it with real-world datasets, and the empirical results
were not used but have great potential to further improve suggest our proposed HISF the best choice for joint analysis
the recommendation. of datasets across domains (Sect. 6).
Motivated by this literature gap, we propose a Hid- 4) Utilizing knowledge from potentially unlimited
den Implicit Similarities Factorization model (HISF) as the number of sources (Sect. 5): We propose a generalized
first algorithm to utilize both explicit and implicit simi- model utilizing both similarities from multiple datasets.
larities across datasets. Our idea for explicit similarities, Our experiments demonstrate its advantages of leveraging
differed from CMF which assumes an identical coupled correlations from more than two datasets, suggesting the
3
item
 
user of X(1) specific parts
domain 1
[ ]
T
V(0)
V(1)
common parts
X(1) U(1) S(1) (explicitly shared)
user clusters alignment
(implicitly shared)
item
user of  
domain 2 X(2) specific parts
[ ]
T
V(0)
X(2) U(2) S(2)
V(2)
Fig. 2: (Best viewed in color) Our Hidden Implicit Similarities Factorization model for two rating matrices X(1) and X(2) . They are coupled in
their item dimension. Our proposed joint analysis decomposes X(1) into U(1) , S(1) , common V(0) and its specific V(1) . At the same time, X(2)
is factorized into U(2) , S(2) , common V(0) and its specific V(2) . Note that our model also matches clusters of users with similar interests in
non-coupled factors U(1) and U(2) by aligning correlated clusters (captured in their columns) as close as possible.
idea can be extended to multi-sources (Sect. 6). 2.2 Our notation

A rating matrix from n users for m items is denoted by a
boldface capital, e.g., X. Each row of the matrix is a vector of
2 BACKGROUND AND OUR NOTATION a user’s ratings for all items while each column is a vector of
ratings from all users for a specific item. We denote vectors
This section provides a background on factorization and a with boldface lowercases, i.e., u. A boldface capital and
reference for the notation we use. lowercase with indices in their subscript are respectively
used for an entry of a matrix and a vector. Table 1 lists all
other symbols we thoroughly use in this paper.
2.1 Rating matrix
Recommendation systems are based on past user prefer-
2.3 Matrix Tri-Factorization
ences such as feedbacks users provide to the systems. For
example, users of Facebook click ”thumb-up” to like a post. Matrix tri-factorization decomposes an input matrix into
For another example, Amazon’s users rate products from 1 three low-rank matrices, called factor matrices or simply
to 5 stars on Amazon website. These kinds of rating are factors. When the given matrix is complete, decomposing
often formed a matrix with users in one dimension and it into factors can be done by Singular Vector Decompo-
items in the other. Its entries are ratings that users provided. sition (SVD). However, where it is incomplete, computing
Rating matrices are usually sparse as the users only rate a its exact decomposition is an intractable task [26]. Thus, a
few items. more efficient and feasible approach is to approximate the
incomplete X of n users by m item as a matrix multiplication
TABLE 1: Symbols and their descriptions
of U ∈ Rn×r , S ∈ Rr×r and VT ∈ Rr×m :
X ≈ U × S × VT
Symbol Description
(i) where r is the rank of the factorization, UT × U = I and
X Rating matrix from i-th dataset
U(i) The first dimension factor of X(i)
VT × V = I.
V(0) Common parts of the coupled factors In this case, U is the user factor, V is the item factor
V(i) Domain-specific parts of the coupled factor of X(i) and S is the weight between U and V. These factors can be
S(i) Weighting factor of X(i) found by solving the optimization of the following:
AT Transpose of A 2
A+ minL = X − U × S × VT

Moore-Penrose pseudo inverse of A (1)
I The identity matrix
kAk Frobenius norm
n, m, p Dimension length
c Number of common clusters in coupled factors 2.4 Joint analysis of different datasets
r Rank of decomposition
ΩX Number of observations in X
Joint analysis of two matrices coupled in one dimension can
∂
Partial derivative with respect to x be done by minimizing a coupled loss function [19]:
∂x
L Loss function 2 2
L = X(1) − U(1) × S(1) × VT + X(2) − U(2) × S(2) × VT

λ Regularization parameter
× Multiplication
where V is the common coupled factor for both datasets.
4
3 P ROBLEM DEFINITION 4.1 Explicit similarities - Sharing common and preserv-

ing domain-specific item latent variables in V
We first define our problem formally here and then discuss
our solution in the next section. Suppose X(1) and X(2) are Even though X(1) and X(2) contain the same m items,
two correlated rating matrices from two different domains. it is unreasonable to force them to share the same low-
Also, suppose that they have an explicitly coupled relation- rank item factors (i.e., V(1) equals V(2) ). The fact that both
ship in their second dimension. For instance, X(1) is a rating matrices capture ratings of the same items indicates that
matrix from n users for m items and X(2) is another one they are strongly correlated. Nevertheless, they also have
from p users for the same m items. In this case, they have their own unique features of their domain. For this reason,
an identical item dimension. This assumption is reasonable we include both common and unique parts in item factors to
as we can find many cross-domain datasets that possess this better capture correlations of the coupled dimension among
coupled relationship. different datasets. Thus, the coupled loss function becomes:
Given X(1) and X(2) are coupled in the item dimension, (0) T 2 (0) T 2
(1) (1) V + X(2) − U(2) S(2) V
(1)
one can combine them (on their coupled dimension) to form X − U S
L= (1)
+ λθ
V V(2)
a big matrix and decompose it. By doing so, the factor in
the coupled dimension can be updated by both datasets. where V(0) ∈ Rm∗c is the common part, V(1) and V(2) ∈
As a result, this combination is the same as CMF [19]. This Rm∗(r−c) are domain-specific parts, r is rank of the de-
is probably not a proper approach especially when X(1) composition, c is the number of common rows such that
and X(2) are from two different domains. Their coupled 1 ≤ c ≤ r and θ is a squared L2 regularization term. These
factors have some explicit similarities, yet they also possess common V(0) and domain-specific V(1) , V(2) parts in cou-
their own domain’s characteristics. Thus, assuming their pled factors are illustrated in Figure 2.
coupled factors to be identical potentially loses the unique
yet critical knowledge of their own domains. Moreover, 4.2 Implicit similarities - Aligning related user latent
X(1) and X(2) also have implicit correlations in their non- clusters across domains
coupled user dimension. These similarities, if used properly, Besides the explicit similarities in common parts of coupled
may provide rich insights for better understanding their factors (as described above), the non-coupled dimension
underlying structure. To this end, two questions are raised: of datasets from different sources also possess a strong
how to better use these explicit similarities? And how we relationship. Even though users in X(1) and X(2) are dif-
can exploit the implicit similarities to fully understand ferent users, they may have similar behaviors which can
correlations between X(1) and X(2) ? be grouped by latent factors. For instance, different users
As for the explicit similarities, our key hypothesis is that with similar preferences on sci-fi movies can be grouped
two cross-domain datasets also have their unique patterns. as sci-fi fans. This intuition inspires us to hypothesize that
Our idea is to come up with a way to integrate these non-coupled user factors share closely related latent user
domain-specific patterns together with the common parts in clusters. Thus, we first find cluster matchings between user
the coupled factors. We formulate the coupled item factors clusters in U(1) and those in U(2) . Based on these matchings,
not to be identical, but to contain both common and domain- centroids of their correlated clusters are then regularized to
specific parts. Figure 2 illustrates an example of common be as close as possible. Figure 2 demonstrates a case of user
V(0) and domain-specific V(1) , V(2) in the coupled factors. clusters alignment between U(1) , U(2) .
Besides, our another key hypothesis is with implicit We use a few examples to explain how challenging it is to
similarities. Since factorization is equivalent to spectral clus- match user clusters across domains and how clustering can
tering [27], we propose non-coupled user factors to share be achieved with matrix factorization. We then introduce
closely related latent user clusters. Even though users in our solution that first finds related user clusters across
X(1) and X(2) are different users, they are clustered by their domains and then aligns them together.
preferences. Two implicitly related user clusters in the two
domains potentially provide some insights to understand 4.2.1 Factorization as a clustering method
their behaviors. We thus align them to be as close as possible
Bauckhage [27] proved that matrix factorization is equiv-
so the knowledge between user clusters can be exploited.
alent to k-means clustering when each row of the factor
Figure 2 demonstrates a case of user clusters alignment
contains (k − 1) 0 elements and only one element that is
between U(1) , U(2) .
1. This condition does not apply in the setting of this paper,
but we can assign 1 to the element with the maximum value
4 O UR H IDDEN I MPLICIT S IMILARITIES FACTOR - of each row and 0 to the rest to help guide the clustering.
We now make an example where X(1) contains ratings
IZATION M ODEL (HISF) of two groups of users with similar preferences as shown
In this section, we introduce how we exploit explicit and im- in Fig. 3: one group giving brown circle ratings and another
plicit similarities among different datasets. We first explain one giving green triangle ratings. When X(1) is factorized,
a case where X(1) and X(2) are coupled in their second these two user groups are captured in columns of the user
dimension. An extension to more matrices will be then dis- factor, i.e., one column of the user factor contains users with
cussed in the next section. We make two key improvements: brown circles and another one captures users with green
integrating both common and domain-specific parts in the triangles. This example illustrates that factorization can
item low-rank factors and aligning similar user clusters in cluster users based on their preferences. However, it does
the user factors. not provide any information on whether the first column is
5
item
   
1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 U(1) U(1) U(1) U(1)
user 1  
T T
X(1) U(1) S(1) V(1) U’(1) S’(1) V’(1)
a) b)
Fig. 3: (Best viewed in color) Matrix factorization as a clustering U(2) U(2) U(2) U(2)
method. Suppose X(1) contains ratings of users for items and there
are two user groups represented by their rating similarities (brown
circles and green triangles). When X(1) is factorized with matrix tri- (a) (b) (c) (d)
factorization, these two user groups are captured in columns of the
user factor. Two possible cases can happen: a) users with brown circles Fig. 5: (Best viewed in color) User clusters’ alignment between X(1) in
are in the first column and those with green triangles are in the second
Fig. 3 and X(2) in Fig. 4. There are four possible cases: in a) and d) the
column of U(1) ; or b) users with brown circles are in the second column
user cluster in the first column of U(1) matches the user cluster in the
and those with green triangles are in the first column of U0 (1) . first column of U(2) and the user cluster in the second column of U(1)
matches the user cluster in the second column of U(2) ; in b) and c) the
item user cluster in the first column of U(1) matches the user cluster in the
 1
 1 1 1 1
 1
 1 1 1 second column of U(2) and the user cluster in the second column of
1
1 1 1 1 1 1 1 1
U(1) matches the user cluster in the first column of U(2) . These show
user 2   the challenge to determine how clusters in U(1) and U(2) are aligned
at each iteration.
T T
X(2) U(2) S(2) V(2) U’(2) S’(2) V’(2)
Fig. 4: (Best viewed in color) Suppose X(2) is obtained from another

common and different columns. These columns of V define
domain where different users rate the same set of items as that of X(1) .
X(2) also has two user groups having similar preferences with those in new coordinates of low dimensional (linear transformation)
X(1) : circle and triangle preferences. column-wise (item side) information of X. The common
columns of V(1) and V(2) thus provide new coordinates
of common item information of X(1) and X(2) ; each of new
the group of users with brown circles or green triangles as coordinates (each common columns of V(1) and V(2) ) de-
both U(1) ×S(1) ×V(1) T and U0 (1) ×S0(1) ×V0 (1) T are equally good fines a cluster of items (linear transformation of the original
solutions. item side of X) having the same characteristics (e.g., sci-fi
movies or comedy movies). The rows of S(1) and S(2) are
4.2.2 A challenge to align similar clusters across domains the values (weights) of different user clusters (captured in
In case there is another dataset from a different domain X(2) columns of U(1) and U(2) respectively) in these common
having ratings from different users on the same set of items new coordinates. If a user cluster in X(1) and another one
as shown in Fig. 4. By factorizing X(2) , we can group users in X(2) are geometrically close in the new coordinates,
with blue circle and those with yellow triangle preferences it indicates these user clusters have similar interests on
(2)
in columns of user factor U(2) or U0 . items of the same characteristics. Hence, we match the S(1)
Although the behaviors of users in X(1) and those in and S(2) weights that are close geometrically in the new
(2)
X are highly correlated, finding a match between user coordinates defined by V(1) and V(2) .
clusters in X(1) and X(2) is not an obvious task because of We use this principle to find matches of user clusters
(1) (2)
two reasons. Firstly, users in X(1) and X(2) are different. It across domains. A user cluster U∗,i matches with U∗,j if
means that we do not know about the users nor their order (1) (2)
Si,∗ is the closest to Sj,∗ . Thus, we compare the Euclidean
in the rating matrices. Secondly, as alternating least square
distances between rows of S(1) and those of S(2) to match
method [28] is often used to find the factors, factors are
user clusters whose distances are the shortest.
initiated randomly. Therefore, a particular user cluster may
never be fixed in a particular column of the user factor. In 4.2.4 Aligning correlated clusters
other words, we have no information of whether users with
Once matches between user clusters across domains are
brown circles are to be captured in the first or the second
found, we align them to be as close as possible in optimiza-
column of the user factor of X(1) . In the same way, users
tion processes. This alignment is to enforce similar users in
with blue circles can be captured in the first or the second
latent feature space to rate new items similarly. Here, we
column of the user factor of X(2) . Therefore, there are four
propose to regularize the difference between centroids of
possible cases for matching user clusters of X(1) and X(2) as
the closely related user clusters across U(1) and U(2) .
demonstrated in Fig. 5. In the event the input datasets have
N clusters, the situation is more complex as there will be N2 (0) T 2 (0) T 2
possible cases. (1) (1) V + X(2) − U(2) S(2) V
(1)
L= X − U S (1)

V V(2)
r 2
4.2.3 Matching correlated clusters across domains

X (1) T (2) T
+ m
l − m∗ + λθ
To find matches of users clusters across domains, we re- l=1
visit our intuition about similar user groups. We assume (2)
that users have similar behavior on the common item set. (1) T
When we extract user clusters by matrix tri-factorization as where ml denotes a row vector of the l-th user cluster’s
(2) T
discussed in Sect. 4.2.1, item factors V(1) and V(2) have centroid in U(1) and m∗ denotes a row vector of the
6
xy
u1,x+u2,x+u3,x u1,y+u2,y+u3,y
m T1 = ( , )
3 3
U u4,x+u5,x u4,y+u5,y
m T2 = ( , )
2 2
Fig. 6: An illustration of how centroid of a cluster is computed. U

captures two user clusters: the first three users are in the first cluster
and the last two are in the second one. x and y denote the first and the
second column of U, respectively.
matched user cluster’s centroid in U(2) ; an example of how

(1) T (2) T
ml and m∗ are computed is shown in Fig. 6.
When the user clusters across domains are matched as
Fig. 5a and Fig. 5d, then:
r 2 q
X (1) T (2) T (1) (2) 2 (1) (2) 2
m
l − m∗ = (m1,x − m1,x ) + (m1,y − m1,y )
l=1
q
(1) (2) (1) (2)
+ (m2,x − m2,x )2 + (m2,y − m2,y )2
However, if the user clusters across domains are matched

(2) T
as Fig. 5b and Fig. 5c, the order of columns in m∗ has to
(1) T
be adjusted to match with the order of columns ml such
that
r 2 q
(1) T (2) T
= (m(1) − m(2) )2 + (m(1,y) − m(2) )2
X
m
l − m ∗ 1,x 2,y 1 2,x
l=1
q
(1) (2) (1) (2)
+ (m2,x − m1,y )2 + (m2,y − m1,x )2
Fig. 8: (Best viewed in color) An illustration of how well our cluster
alignment method works. Sci-fi fans and comedy ones from X(1) and
sci-fi sci-fi sci-fi comedy comedy comedy sci-fi sci-fi sci-fi comedy comedy comedy X(2) in Fig. 7 are captured in user factors U(1) and U(2) . Circles are
#1 #2 #3 #1 #2 #3 #1 #2 #3 #1 #2 #3 users captured in the first column of U, and triangles are for those
captured in the second column of U. As U(1) and U(2) are randomly
initiated, users are scattered. User clusters are formed and aligned
over iterations. From iteration 2, two user clusters in X(1) are clearly
separated and gradually aligned with those in X(2) . There is a change
in cluster matching in iteration 3: the first column of U(1) is matched
(a) X(1) (b) X(2) with the second column of U(2) and vice versa. Since then, centroids of
clusters are aligned to be as close as possible from iteration five till the
Fig. 7: Generated ratings of two domains X(1) and X(2) . Blank entries last iteration.
are unobserved ratings. Although users in X(1) and X(2) are different,
they share implicit similarities in their preferences: those who like sci-fi
genre rate sci-fi movies highly and those who are comedy fans do so
for comedy movies. first column and y is the second column of U(2) . Otherwise,
we change the order of columns in U(2) before plotting its
To elaborate how this alignment works, we run Algo- users so that matching columns between U(1) and U(2) are
rithm 1 with a synthetic data in Fig. 7 where X(1) and aligned, i.e., y is the first and x is the second column of U(2) .
X(2) each has six unique users who rate the same list of Initially, the algorithm randomizes the user factors. So
six movies (three sci-fi and three comedy movies). Some the users are scattered randomly in the x-y space. Iteration
users in X(1) and X(2) like sci-fi genre and provide high 1 starts grouping users to different clusters and aligning
ratings for sci-fi movies; some other prefer comedy category them. After iteration 2, users with similar behaviors in U(1)
and rate comedy movies five stars. Thus, users in each and those in U(2) are formed, and we can see the first and
domain can be grouped by their preferences: one cluster the second user clusters of X(1) are aligned with the first
of sci-fi and another one of comedy fans. Furthermore, each and the second user clusters of X(2) , respectively. In other
group of users in X(1) shares its implicit preference with words, sci-fi fans in the first column of U(1) (blue circles) are
a corresponding group of users in X(2) . As a result, when matched with those in the first column of U(2) (red circles)
algorithm 1 factorizes X(1) and X(2) with rank r = 2, each whereas comedy fans in the second column of U(1) (blue
user cluster is captured in a column of the user factors. In triangles) are matched with those in the second column of
addition, their user clusters are aligned so that those with U(2) (red triangles). Iteration 3 makes a correction on cluster
similar preferences are close to each other. alignment: sci-fi fans, now in the second column of U(2) ,
We then plot the resulted user clusters of several iter- are matched with those in the first column of U(1) (blue
ations into an x-y coordinates in Fig 8 where the x-axis is circles) whereas comedy fans (in the first column of U(2)
the first column of U(1) and the y-axis is its second column. (red circles) are matched with ones in the second column
If the first and the second column of U(2) match with the of U(1) (blue triangles). In this case, the order of columns
first and the second column of U(1) respectively, x is the U(2) is corrected so that users in U(1) and U(2) are in
7
Algorithm 1: HISF: Utilizing both explicit and implicit partial derivative of L with respect to it is set to zero.
similarities from two matrices m
δL X (1) (1) T (01) (01) T
Input : X(1) , X(2) , E =−2 Xi,j − ui vj vj
(1) T
δui j
Output: U(1) , S(1) , V(0) , V(1) , U(2) , S(2) , V(2)
(1) T (1) T
− bT + 2λui

+ 2 ui
1 Randomly initialize all factors
(1) T (1) T T
2 Initialize L by a small number = − 2xi,∗ V(01) + 2ui V(01) V(01)
(1) T (1) T
3 repeat + 2ui − 2bT + 2λui
4 PreL = L
(1) T (2) T (1) T
where bT = −ml + m∗ + ui , l is the cluster
5 Find matches between clusters in U(1) and U(2)
(1) T
user i belongs to and xi,∗
is a row vector of all observed
6 Solve U(1) by (5) (1)
7 Solve U(2) by (6) Xi,j , ∀j ∈ [1, m].
8 Solve common V(0) by (7) By setting δL (1) T
= 0, we achieve updating rule for
δui
9 Solve domain-specific V(1) by (8) (1) T
ui :
10 Solve domain-specific V(2) by (9)
+
Solve S(1) by (10)

11 (1) T T (1) T
ui = V(01) V(01) + (λ + 1)I xi,∗ V(01) + bT (5)
12 Solve S(2) by (11)
13 Compute L following (2) (2) T
In the same way, optimal uk can be derived from:
14 until ( PreL−L
PreL < E ) +
(2) T T (2) T
uk = V(02) V(02) + (λ + 1)I xk,∗ V(02) + bT (6)
(1) T (2) T (2) T

the same coordinates. In the next iterations, centroids of where bT = −ml + m∗ + uk , l is the cluster
corresponding clusters in X(1) and X(2) are adjusted to be (2) T
matched with the one user k belongs to and xk,∗ is a row
close as observed in iteration 5’s result or the last iteration’s. (2)
vector of all observed Xk,j , ∀j ∈ [1, m]. I is the identity
matrix.
4.3 Optimization
4.3.2 Solving common V(0)
Equation (2) is not a convex function. However, it is convex (1) T (10) (11) T (1) (2) T
with respect to a factor when the others are fixed. Therefore, Let ui = ui |ui S and uk =
we use the alternating least square framework to take turns
(20) (22) T (2) (10) T (20) T 1∗c
uk |uk S where ui , uk ∈ R and
optimizing one factor in function (2) while fixing the others, (11) T (22) T
as shown in Algorithm 1. Moreover, as data for updating ui , uk ∈ R1∗r−c , then Equation (2) can be rewritten
rows of the factors are independent, our algorithm com- as:
n,m 2
putes factors in a row-wise manner instead of full matrix X (1) (10) T (0) (11) T (1)
L= Xi,j − ui vj − ui vj
operations so that the computation can be done in parallel i,j
with multiple CPU cores or a distributed system. To this p,m
X 2
(2) (20) T (0) (22) T (2)
end, we rewrite equation (2) as the following: + Xk,j − uk vj − uk vj
k,j
n,m (0)
" #
X 2 r 2
(1) (1) T (1) vj X (1) T (2) T
L= Xi,j − ui S (1) + m
l − m ∗
+ λθ
i,j
vj l=1

p,m
" (0) # r
2 2
X (2) (2) T vj X (1) T (2) T Similar to the case of U(1) and U(2) , we now optimize
+ Xk,j − uk S(2) (2) + m
l − m∗ + λθ
vj (0)
k,j l=1 vj while fixing all other parameters. Again, by setting the
(3) (0)
partial derivative of L with respect to vj to zero, we can
(0)
achieve optimal value of vj .
4.3.1 Solving U(1) and U(2)
" (0) # " (0) # n
δL (10) T (0) (10)
(1) vj (2) vj
(1)
X
(01) (02) =−2 Yi,j − ui vj ui
Let vj = S (1) and vj = S (2) , then δvj
(0)
vj vj i
p
Equation (3) becomes X (2) (20) T (0) (20) (0)
−2 Yk,j − uk vj uk + 2λvj
n,m 2 p,m 2 k
X (1) (1) T (01)
X (2) (2) T (02) T (1) T (0)
L= Xi,j − ui vj + Xk,j − uk vj = − 2U(1) y∗,j + 2U(1) U(1) vj
i,j k,j T (2)
(4) (2) T (0) (0)
r 2 − 2U(2) y∗,j + 2U U(2) vj + 2λvj
X (1) T (2) T
+ m
l − m∗ + λθ
(1) (1) (11) T (1) (2) (2)
l=1 where Yi,j = Xi,j − ui vj and Yk,j = Xk,j −
(22) T (2) (1) (2)
By fixing all vj , Equation (4) is a convex function with uk vj ; y∗,j and y∗,j are column vectors of all observed
(1) T (1) T (1) (2)
respect to ui . As a result, ui is optimal when the Yi,j , ∀i ∈ [1, n] and Yk,j , ∀k ∈ [1, p], respectively.
8
(0) (1) T
Our updating rule for vj is: xi,∗ V(01) , and performs Cholesky decompositions for the
ΩX(1)
(0)

T T
+
T (1) T (2)

pseudo inverse. Preparing V(01) requires O( n r) opera-
vj = U(1) U(1) + U(2) U(2) + λI U(1) y∗,j + U(2) y∗,j (1) T
(01) T (01) Ω
(7)
tions while V V and takes O( Xn r2 )
xi,∗ V(01) (1)
each. Cholesky decomposition of r × r matrices is O(r3 ).

4.3.3 Solving domain-specific V(1) and V(2) Thus, solving U(1) takes O(nr3 + ΩX(1) r2 ).
(1) In a similar way, solving U(2) takes O(pr3 + ΩX(2) r2 );
We now perform the similar operations with respect to vj
(2)
V(0) takes O(mc3 + (ΩX(1) + ΩX(2) )c2 ), V(1) and V(2)
and vj . take O(m(r − c)3 + (ΩX(1) + ΩX(2) )(r − c)2 ), S(1) takes
δL T (1) T (1) (1) O(rr3 + ΩX(1) r2 ) and S(2) takes O(rr3 + ΩX(2) r2 ). As
= −2U(1) z∗,j + 2U(1) U(1) vj + 2λvj =0
δvj
(1) a result, the complexity of the Algorithm 1 is
O((n + p + m + r)r3 + (ΩX(1) + ΩX(2) )r2 ).
(1) (1) (10) T (0) (2) (2)
where Zi,j = Xi,j − ui vj and Zk,j = Xk,j −
(20) T (0) (1) (2)
uk vj ; z∗,j and z∗,j are column vectors of all observed 5 E XTENSION TO THREE OR MORE MATRICES
(1) (2)
Zi,j , ∀i ∈ [1, n] and Zk,j , ∀k ∈ [1, p], respectively. We can find more than two correlated matrices in many
(1)
Then we derive the updating rule for vj : cases. For example, Amazon has ratings from different do-
+ mains, e.g., Books, Movies and TV, Electronics, Digital Mu-
T T (1)
vj
(1)
= U(1) U(1) + λI U(1) z∗,j (8) sics, etc. These rating matrices may have a close relationship
that can help to collaboratively improve recommendation
(1) (2) accuracy. In this case, suppose we have three correlated
Analogy to optimizing vj , optimal vj is obtained by:
matrices X(1) , X(2) and X(3) . They are coupled in their
second dimension, i.e. X(1) is a rating matrix from n users
+
(2) T T (2)
vj = U(2) U(2) + λI U(2) z∗,j (9) for m items, X(2) is another rating matrix from k users
for the same m items and X(3) is another rating matrix
4.3.4 Solving weighting factor S(1) and S(2) from l users for the same m items. Our idea above can
Let be extended to utilize common parts of coupled factors of
the three matrices (explicit similarities) and align clusters
T (1) (1) (1) (1) (1) (1)
s(1) =(S1,1 ,S2,1 ,...,Sr,1 ,S1,2 ,S2,2 ,...,Sr,2 ,...,S1,r ,S2,r ,...,S(1)
r,r )
(1) (1)
among non-coupled factors (implicit similarities). Thus, we
propose the following extension for utilizing explicit and
and implicit similarities among three or more matrices:
(1) (1) (1) (1) (1) (1) (1) (1) (1) (1) (1) (1)
aT =(Ui,1 Vj,1 ,Ui,2 Vj,1 ,...,Ui,r Vj,1 ,Ui,1 Vj,2 ,Ui,2 Vj,2 ,...,Ui,r Vj,2 ,
3 (0) T 2
X − U(i) S(i) V
(1) (1) (1) (1) (1) (1)
X (i)
...,Ui,1 Vj,r ,Ui,2 Vj,r ,...,Ui,r Vj,r ) L= + simimplicit + λθ (12)
V(i)
i=1
Equation (3) then becomes:
n,m
where simimplicit is the regularization term of implicit simi-
X − s(1) T a 2 + λ(ks(1) T k2 ) + const larities across domains and is defined by:
X (1)
L= i,j
i,j
r 2 Xr 2
X (1) T (2) T (2) T (3) T
where const is the remaining regularization terms. simimplicit = m
l − m ∗
+

m
l − m ∗

T l=1 l=1
We can achieve optimal s(1) when ∂s∂L (1) T
=0 r
(3) T
2
(1) T
X
+ m
l − m ∗

∂L (1) T T T
T
= − 2x∗,∗ A + 2s(1) AT A + 2λs(1) =0 l=1
∂s(1) (10)
⇔ s(1)
T
= AT A + λI
−1 (1) T
x∗,∗ A
Updating rules for optimizing U(1) , U(2) , U(3) , S(1) ,
S , S(3) , V(1) , V(2) and V(3) can be achieved similarly to
(2)
(1) T (1) (1) T
where x∗,∗ contains observed Xi,j ; x∗,∗ ∈ R1xΩX(1) ; A ∈ the case of two matrices in Section 4.3. Moreover, following
RΩX(1) x(rxr) . the same derivations for common parts of two matrices in
T (0)
In an analogy way, the update rule for s(2) is: Section 4.3.2, we reach an updating rule for vj , which is
−1 (2) T the common parts of three matrices:
T
⇔ s(2) = AT A + λI x∗,∗ A (11)
+
(0) T T T
vj = U(1) U(1) + U(2) U(2) + U(3) U(3) + λI ×
4.4 Complexity analysis (13)
T (1) T (2) T (3)
U(1) y∗,j + U(2) y∗,j + U(3) y∗,j
The computational complexity of Algorithm 1 is
O((n + p + m + r)r3 + (ΩX(1) + ΩX(2) )r2 ) By optimizing Equation (12), we leverage both explicit
similarities (in a form of common V(0) parts) and implicit
Proof. Algorithm 1 computes U(1) in line 6, U(2) in line 7,
similarities (in a form of aligned user groups in U(1) , U(2)
common V(0) in line 8, domain specific V(1) and V(2) in line
and U(3) ) among X(1) , X(2) and X(3) . Utilization of corre-
9 and 10, S(1) in line 11 and S(2) in line 12. Finding U(1)
(1) T
lations from four or more matrices can be easily extended
involves solving n ui in Equation 5. Equation 5 prepares in the same way. Therefore, our proposed method does not
(01) ΩX(1) T
V of size n r (average), computes V(01) V(01) and limit itself to a certain number of correlated matrices.
9
TABLE 2: Dimension and number of known entries for training, valida-
tion and testing of census data on New South Wales (NSW)(X(1) ) and
• Total Household Income (weekly) by Rent (weekly)
Victoria (VIC) (X(2) ) states as well as crime statistics of NSW (X(3) ). • Family Composition by Mortgage Repayment
• Family Composition by Income Comparison for Par-
Characteristics X(1) X(2) X(3) ents/Partners
Dimension 154 × 7,889 81 × 7,889 154 × 62
• Family Composition and Social Marital Status by
Training 91,069 47,900 661 Number of Dependent Children
Validation 4,793 2,521 34 • Selected Labour Force, Education and Migration
Testing 23,965 12,605 173 Characteristics
• Family Composition and Labour Force Status of Par-
TABLE 3: Dimension and number of known entries for training, vali-
dation and testing of Amazon dataset on books (X(4) ), Movies (X(5) )
ent(s)/Partners by Total Family Income (weekly)
and Electronics (X(6) ). • Non-School Qualification: Level of Education by Age
by Sex
Characteristics X(4) X(5) X(6) • Non-School Qualification: Field of Study by Age by
Dimension 5,000 × 5,000 5,000 × 5,000 5,000 × 5,000 Sex
Training 158,907 94,665 41,126 • Labour Force Status by Age by Sex
Validation 8,363 4,982 2,164 • Industry of Employment by Sex
Testing 18,585 11,071 4,809 • Occupation by Sex
We form a matrix X(1) (LGA by population and family
6 E XPERIMENTS AND A NALYSIS profile) for NSW and another matrix X(2) for VIC. Values of
We evaluate our proposed HISF in comparison with existing the matrix are normalized by row. We then randomly select
algorithms, including CMF [19], CST [4], CBT [23] and 10% of the data and use its 80% for training and 20% for
CLFM [24]. This evaluation thoroughly studies two test testing.
cases: one with two matrices and another one with three
matrices. Our goal is to evaluate how well these algorithms 6.1.2 Dataset #2
suggest unknown information based on the observed cross- Bureau of Crime Statistics and Research (BOCSAR)2 pro-
domain ratings. For this purpose, we compare them based vides a statistical record of criminal incidents within 154
on the commonly used root mean squared error (RMSE) LGAs of New South Wales. There are 62 specific offence cat-
metric. egories. We form a matrix X(3) (LGA by offence categories)
v whose entries represents how many cases were reported for
u 2
an offence category in an LGA. Values of the matrix are
uP
u n,m U × S × VT − X
u i,j i j i,j normalized by row. We randomly select 10% of the data for
RMSE =
t
ΩX the experiments. Among them, 80% of X(3) are for training
and the rest are for testing.
where ΩX is the number of observations of X.
6.1.3 Dataset #3
6.1 Data for the experiments Three matrices of ratings for books, movies and electronics
We use three publicly available datasets for our experi- are extracted from Amazon website [29]. The book data
ments. Their characteristics are summarized in Table 2 and contains ratings from 305,475 users on 888,057 books; the
Table 3. movie data is ratings from the same 305,475 users on 128,097
movies and TV programs; the electronics data is of the same
6.1.1 Dataset #1 users on 196,894 items. All ratings are from 1 to 5. For this
Australian Bureau of Statistics (ABS)1 publishes comprehen- experiment, the data is constructed as the following:
sive census data of New South Wales (NSW) and Victoria - We first adopt the same sub-sampling approach as in
(VIC) states. The dataset for NSW comprises of populations [4] by randomly extracting 104 × 104 dense rating matrices
and family profiles of 154 areas, so-called “local government from these three matrices. Then, we take three sub-matrices
areas” (LGA), and of 81 LGAs in VIC. 7,889 aspects of of 5,000×5,000 each, as summarized in Table 3. All sub-
population and family profile are extracted from the below matrices share the same users, but no common items.
census categories: - All ratings are normalized by 5 so that their values are
from 0.2 to 1.
• Rent (weekly) by Landlord Type
• Rent (weekly) by Family Composition for Couple
6.2 Experimental settings
Families
• Rent (weekly) by Family Composition for One Parent We performed factorization with different ranks. Also, each
Families algorithm was run five times. The mean and standard
• Total Family Income (weekly) by Number of Chil- deviation of the results are reported in the next Section.
dren for Couple Families Furthermore, we assume small changes across consecutive
• Total Family Income (weekly) by Number of Chil- iterations indicate an algorithm’s convergence. Thus, we
dren for One Parent Families stopped the algorithms when changes were less than 10−5 .
1. ABS: http://www.abs.gov.au/websitedbs/censushome.nsf/ 2. BOCSAR: http://www.bocsar.nsw.gov.au/Pages/bocsar crime

home/datapacks stats/bocsar crime stats.aspx
10
TABLE 4: Mean and standard deviation of tested RMSE on ABS NSW and ABS VIC data with different algorithms. For CST, when X(1) is the
target, X(2) is used as an auxiliary data and vice versa. HISF1 denotes the algorithm with only explicit similarities whereas HISF2 uses both
implicit and explicit similarities. Best results for each rank are in bold. The Hotelling’s T squared tests row presents p-value of Hotelling’s T
squared tests between each algorithm and our proposed HISF.
Rank Dataset CMF CBT CLFM CST HISF1 HISF2

ABS NSW X(1) 0.0226 ±0.0026 0.0839 ±0.0002 0.0838 ±0.0002 0.0248 ±0.0005 0.0183 ±0.0028 0.0132 ±0.0002
5
ABS VIC X(2) 0.0364 ±0.0031 0.0844 ±0.0003 0.0845 ±0.0004 0.0271 ±0.0015 0.0335 ±0.0091 0.0266 ±0.0030
ABS NSW X(1) 0.0222 ±0.0009 0.0836 ±0.0004 0.0842 ±0.0006 0.0244 ±0.0005 0.0192 ±0.0032 0.0131 ±0.0003
7
ABS VIC X(2) 0.0428 ±0.0020 0.0845 ±0.0004 0.0849 ±0.0004 0.0288 ±0.0025 0.0318 ±0.0038 0.0239 ±0.0025
ABS NSW X(1) 0.0241 ±0.0011 0.0841 ±0.0002 0.0848 ±0.0009 0.0231 ±0.0009 0.0192 ±0.0020 0.0143 ±0.0004
9
ABS VIC X(2) 0.0476 ±0.0040 0.0852 ±0.0003 0.0848 ±0.0003 0.0283 ±0.0024 0.0304 ±0.0036 0.0221 ±0.0019
ABS NSW X(1) 0.0265 ±0.0026 0.0846 ±0.0007 0.0841 ±0.0007 0.0223 ±0.0006 0.0186 ±0.0012 0.0143 ±0.0003
11
ABS VIC X(2) 0.0501 ±0.0029 0.0858 ±0.0005 0.0851 ±0.0007 0.0267 ±0.0020 0.0302 ±0.0025 0.0242 ±0.0015
ABS NSW X(1) 0.0237 ±0.0024 0.0851 ±0.0002 0.0850 ±0.0005 0.0222 ±0.0010 0.0177 ±0.0023 0.0151 ±0.0004
13
ABS VIC X(2) 0.0489 ±0.0032 0.0860 ±0.0003 0.0852 ±0.0003 0.0290 ±0.0026 0.0286 ±0.0038 0.0227 ±0.0015
ABS NSW X(1) 0.0229 ±0.0029 0.0853 ±0.0005 0.0847 ±0.0005 0.0219 ±0.0007 0.0180 ±0.0031 0.0150 ±0.0000
15
ABS VIC X(2) 0.0514 ±0.0041 0.0862 ±0.0002 0.0854 ±0.0006 0.0285 ±0.0023 0.0274 ±0.0054 0.0215 ±0.0005
Hotelling’s T squared tests 1.15×10−6 3.97×10−17 1.13×10−18 4.91×10−8 -
6.3 Empirical results columns are 3, 5, 5, 7, 7 and 7 for the decomposition rank
Three scenarios are tested as the following: of 5, 7, 9, 11, 13 and 15, respectively. It means that common
parts together with domain-specific parts better capture the
true correlation nature of datasets. These explicit similarities
6.3.1 Case #1. Latent demographic profile similarities and
together with implicit similar group alignments allow better
latent LGA groups similarities can help to collaboratively
knowledge leveraging between datasets, thus, improving
suggest unknown information in these states
recommendation accuracy.
We use X(1) and X(2) as described in Sect. 6.1 which are To confirm the statistical significance of our method,
from different LGAs of two states. Nevertheless, both of we perform Hotelling’s T squared tests [30] which is the
them are ratings for the same demographic categories. They multivariate version of the t-tests in univariate statistics.
share some common explicit demography similarities as Our objective is to validate if our proposed algorithm dif-
well as implicit LGAs’ latent ones. We would like to assess fers significantly from baselines. We use this multivariate
how well both explicit similarities in demography dimen- Hotelling’s T squared tests here because each population
sion and implicit ones in LGA dimension collaboratively involves observations from two variables: NSW (ABS NSW
suggest unknown information. X(1) ) and VIC (ABS VIC X(2) ). For testing the null hypoth-
Table 4 shows mean and standard deviation of RMSE esis that each pair of algorithms (CMF vs. HISF, CBT vs.
of all algorithms on tested ABS data for New South Wales HISF, CLFM vs. HISF and CST vs. HISF) has identical mean
(X(1) ) and Victoria (X(2) ) states. Both CBT and CLFM that RMSE vectors, let
assume two states have similar demography patterns in H0 : population mean RMSEs are identical for all of the variables
latent sense clearly perform the worst. The results demon- (µN SW 1 = µN SW 2 and µV IC 1 = µV IC 2 )
strate that these explicit similarities (in the form of latent H1 : at least one pair of these means is different
demography patterns) do not fully capture the correlation
(µN SW 1 6= µN SW 2 or µV IC 1 6= µV IC 2 )
nature of two datasets in this case. Thus, they do not help
both CBT and CLFM to improve their performance signif- Because all p-values are smaller than α (0.05), we re-
icantly. CMF applies another approach to take advantages ject the null hypothesis. Therefore, the observed difference
of explicit correlations between NSW state’s population and between the baselines and our proposed algorithm is statis-
family profile and those of VIC state. Specifically, CMF’s as- tically different.
sumption on the same population and family profile factor
between NSW and VIC helps improve its performance over 6.3.2 Case #2. Latent LGA similarities and latent similari-
that of CBT and CLFM. CST allows a more flexible utiliza- ties between demography and crime can help to collabora-
tion of explicit correlations between NSW’s population and tively suggest unknown crime and state information
family profile and those of VIC state than CMF does. As a The advantages of both explicit and implicit similarities
result, CST achieves a little bit higher accuracy than CMF in are further confirmed in Table 5. In this case, they are
recommending NSW’s and VIC’s missing information. applied to other cross domains: ABS NSW demography
Nevertheless, the prediction accuracy can be improved (X(1) ) and NSW Crime (X(3) ). These datasets have explicit
even more as illustrated with our proposed idea of explicit similarities in their LGA latent factors. At the same time,
and implicit similarities discovery. Utilizing them helps our implicit similarities in demography profile and criminal
proposed HISF to achieve about two times higher accuracy behaviors are also utilized. Our proposed HISF leveraging
compared with CMF, about up to 47% improvement for both similarities outperforms existing algorithms. It is worth
NSW and up to 25% for VIC compared to CST. These im- to note here that performance of CST, in this case, is worse
pressive results are achieved when the numbers of common than that of CMF (about two times worse for NSW and a
11
TABLE 5: Mean and standard deviation of tested RMSE on ABS NSW demography and BOCSAR NSW crime data with different algorithms. For
CST, when X(1) is the target, X(3) is used as an auxiliary data and vice versa. HISF1 denotes the algorithm with only explicit similarities whereas
HISF2 uses both implicit and explicit similarities. Best results for each rank are in bold. The Hotelling’s T squared tests row presents p-value of
Hotelling’s T squared tests between each algorithm and our proposed HISF.
Dataset Rank CMF CBT CLFM CST HISF1 HISF2

Demography X(1) 0.0209 ±0.0016 0.0840 ±0.0001 0.0840 ±0.0001 0.0304 ±0.0080 0.0166 ±0.0015 0.0174 ±0.0015
5
Crime X(3) 0.2796 ±0.0204 0.3411 ±0.0035 0.3422 ±0.0071 0.3216 ±0.0052 0.3316 ±0.0175 0.2697 ±0.0073
Demography X(1) 0.0223 ±0.0024 0.0840 ±0.0002 0.0855 ±0.0006 0.0324 ±0.0040 0.0151 ±0.0021 0.0143 ±0.0004
7
Crime X(3) 0.2907 ±0.0265 0.3432 ±0.0021 0.3912 ±0.0188 0.3440 ±0.0055 0.3155 ±0.0100 0.2716 ±0.0029
Demography X(1) 0.0199 ±0.0027 0.0838 ±0.0002 0.0850 ±0.0008 0.0337 ±0.0050 0.0152 ±0.0013 0.0143 ±0.0003
9
Crime X(3) 0.2813 ±0.0261 0.3562 ±0.0134 0.3722 ±0.0249 0.3434 ±0.0087 0.3192 ±0.0094 0.2648 ±0.0058
Demography X(1) 0.0212 ±0.0049 0.0839 ±0.0001 0.0843 ±0.0004 0.0402 ±0.0073 0.0168 ±0.00147 0.0146 ±0.0003
11
Crime X(3) 0.2689 ±0.0143 0.3539 ±0.0061 0.3712 ±0.0199 0.3495 ±0.0051 0.3089 ±0.0118 0.2618 ±0.0012
Demography X(1) 0.0194 ±0.0022 0.0837 ±0.0001 0.0837 ±0.0003 0.0383 ±0.0113 0.0169 ±0.0008 0.0149 ±0.0003
13
Crime X(3) 0.2700 ±0.0150 0.3481 ±0.0070 0.3500 ±0.0135 0.3599 ±0.0024 0.2957 ±0.0070 0.2623 ±0.0024
Demography X(1) 0.0173 ±0.0014 0.0835 ±0.0001 0.0834 ±0.0002 0.0503 ±0.0084 0.0172 ±0.0015 0.0149 ±0.0002
15
Crime X(3) 0.2647 ±0.0031 0.3485 ±0.0038 0.3580 ±0.0099 0.3610 ±0.0011 0.3114 ±0.0097 0.2625 ±0.0015
Hotelling’s T squared tests 8.22×10−4 1.37×10−15 2.51×10−15 1.38×10−6 -
c parameter in Figure 9, 10. These figures illustrate the case

where decomposition rank equals 11. All the curves of HISF
have an identical trend. As the number of common row
HISF HISF of explicit similarities (c parameter) increases, recommen-
dation accuracy improves to reach its minimum rate. Then,
the performance reduces as c increases. This pattern con-
firms the significance of preserving domain-specific parts of
C C
each dataset. In other words, common and domain-specific
(a) (b) parts better capture the correlation nature of datasets across
Fig. 9: Tested mean RMSEs under different number of common row domains, thus, achieve a higher recommendation accuracy.
(c) in coupled factor of HISF with Rank=11. a) Results on ABS NSW Moreover, when c equals rank 11, HISF and CMF utilize
dataset; b) Results on ABS VIC dataset. CMF and CBT do not have c the same explicit similarities (identical coupled factors). Yet,
parameter, thus, their results are used as references. Lower is better.
HISF uses an extra correlation from implicit similarities.
RMSE of HISF outperforms its competitors when c equals 7.
This additional knowledge enables HISF to achieve a lower
RMSE compared to CMF (Figure 9b). All of these suggest
the advantages of both explicit and implicit similarities in
joint analysis of two matrices.
HISF
6.3.3 Case #3. Latent users’ taste similarities when buying
HISF
books, movies and electronics devices and latent item group
similarities can help to collaboratively improve recommenda-
C C tion accuracy
(a) (b) Joint analysis of all three rating matrices from Amazon
website, X(4) for books, X(5) for movies and X(6) for
Fig. 10: Tested mean RMSEs under different common row (c) in coupled electronics, is studied. All of them are explicitly from the
factor of HISF with Rank=11. a) Results on ABS NSW demography;
b) Results on BOCSAR NSW Crime. CMF and CBT do not have c
same users. Nevertheless, their hidden similarities in per-
parameter, thus, their results are used as references. Lower is better. sonal preferences and items’ characteristics can also help to
RMSE of HISF outperforms its competitors when c equals 9. better suggest missing information. We want to assess how
these explicit together with implicit similarities can help to
collaboratively improve recommendation accuracy of each.
bit lower for VIC). The results suggest that flexibility of Table 6 summarizes the performance of CMF, CBT,
utilizing explicit similarities does not work here. On the CLFM and our proposed algorithm on Books, Movies and
contrary, our proposed idea of utilizing both explicit and Electronics datasets from Amazon. CST is not applied due
implicit similarities eventually achieves more stable and to its limit on one auxiliary dataset for one factor; in this
better results (demonstrated in Table 4 and 5). case, we have two side data for user dimension. The re-
We perform similar Hotelling’s T squared tests as of the sults in this three datasets are consistent with those for
case #1. Again, as all p-values are smaller than α (0.05) here, two matrices above. In particular, CBT and CLFM which
we reject the null hypothesis. It is therefore convincing to assume shared rating patterns among the three perform the
conclude that the mean RMSEs between the baselines and worst. CMF improves accuracy a lot in comparison with
our proposed idea differ significantly. CBT and CLFM. Nevertheless, its performance is once more
We also show how HISF works with different values of time outperformed by our proposed ideas. The results, in
12
TABLE 6: Mean and standard deviation of tested RMSE on book, movie and electronics data with different algorithms. CST is not applied here
as it does not support two or more principal coordinates on one factor. Best results for each rank are in bold. The Hotelling’s T squared tests row
presents p-value of Hotelling’s T squared tests between each algorithm and our proposed HISF.
Rank Dataset CMF CBT CLFM HISF-N

Books X(4) 0.2169 ±0.0005 0.2187 ±0.0026 0.2161 ±0.0023 0.1951 ±0.0005
5 Movies X(5) 0.2273 ±0.0022 0.3342 ±0.0039 0.3474 ±0.0063 0.2212 ±0.0010
Electronics X(6) 0.2375 ±0.0017 0.4206 ±0.0036 0.4596 ±0.0066 0.2642 ±0.0011
Books X(4) 0.2170 ±0.0009 0.3068 ±0.0031 0.3059 ±0.0024 0.1977 ±0.0010
7 Movies X(5) 0.2279 ±0.0009 0.4368 ±0.0031 0.4337 ±0.0055 0.2213 ±0.0007
Electronics X(6) 0.2348 ±0.0009 0.6607 ±0.0135 0.6389 ±0.0116 0.2497 ±0.0016
Books X(4) 0.2185 ±0.0010 0.3163 ±0.0004 0.3150 ±0.0026 0.2001 ±0.0005
9 Movies X(5) 0.2302 ±0.0010 0.4661 ±0.0046 0.4594 ±0.0065 0.2218 ±0.0008
Electronics X(6) 0.2354 ±0.0019 0.7098 ±0.0232 0.6909 ±0.0175 0.2411 ±0.0019
Books X(4) 0.2219 ±0.0014 0.3207 ±0.0044 0.3204 ±0.0021 0.2012 ±0.0004
11 Movies X(5) 0.2319 ±0.0021 0.4865 ±0.0019 0.4795 ±0.0043 0.2247 ±0.0005
Electronics X(6) 0.2417 ±0.0010 0.7390 ±0.0090 0.7223 ±0.0123 0.2420 ±0.0018
Books X(4) 0.2244 ±0.0018 0.3267 ±0.0020 0.3291 ±0.0027 0.2021 ±0.0010
13 Movies X(5) 0.2349 ±0.0025 0.5014 ±0.0057 0.4908 ±0.0028 0.2263 ±0.0009
Electronics X(6) 0.2422 ±0.0012 0.7695 ±0.0118 0.7684 ±0.0155 0.2410 ±0.0022
Books X(4) 0.2260 ±0.0014 0.3303 ±0.0024 0.3353 ±0.0025 0.2022 ±0.0008
15 Movies X(5) 0.2365 ±0.0018 0.5132 ±0.0044 0.5070 ±0.0065 0.2270 ±0.0009
Electronics X(6) 0.2447 ±0.0013 0.8067 ±0.0126 0.7688 ±0.0122 0.2417 ±0.0024
Hotelling’s T squared tests 2.80×10−7 1.15×10−5 1.11×10−6 -
this case, demonstrate our idea can be generalized to the

situation of more than two datasets.
We perform Hotelling’s T squared tests for three vari-
ables, i.e. Books, Movies and Electronics, between each HISF HISF
baseline and HISF-N to confirm our method’s statistical
significance. Our Hotelling’s T squared tests this time also
result in very small p-values. This once again confirms C C
the performance between them is significantly different.
We, thus, can conclude that their observed difference is (a) (b)
convincingly significant.
We again show how HISF-N works with different values
of c parameter in Figure 11. There are two conclusions we
can draw by observing the figures. Firstly, domain-specific HISF
parts are quite substantial in addition to the common parts.

Secondly, our model reduces X(4) ’s RMSE the most, then
X(5) ’s and X(6) ’s one in comparison with other algorithms. C
This can be explained by their observed data size. Our (c)
model, which optimizes the loss function (12), gives more
preference on optimizing the domain with more data while Fig. 11: Tested mean RMSEs under different common row (c) in coupled
factor of HISF-N with Rank=11. a) Results on Amazon Books; b) Results
preserving comparable performance on the other domains. on Amazon Movies; c) Results on Amazon Electronics. CMF and CBT
do not have c parameter, thus, their results are used as references. RMSE
of HISF outperforms its competitors when c equals 7.
7 R ELATED WORK
Joint analysis of two datasets explores the possibility of
using their correlations to better understand the underly- their coupled dimension. The concept of explicit similarities
ing nature of the datasets. This thorough understanding as the common factor has been widely used. Bhargava et al.
provides more meaningful insights to improve recommen- [31] proposed location, activity and time would provide a
dation accuracy. Many methods have been proposed to complete picture of users. Thus, different data sources were
achieve this goal. In this section, we summarize popular modeled to have common coupled factors to fuse their ex-
fundamental ideas and some of their extensions. plicit similarities. Transfer by Collective Factorization (TCF)
[32] was proposed to use ratings and binary like/ dislike
from the same users and items. TCF assumed that both user
7.1 Collective Matrix Factorization (CMF) and item factors would be the same between the inputs.
CMF was proposed to jointly analyze two datasets having There are also several ideas to extend this CMF’s coupled
one common dimension [19]. CMF was based on the as- loss function. Acar et al. [20] introduced the concept for
sumption that they have a common low-rank subspace in the joint analysis of a tensor [26] and another matrix. The
13
authors proposed Coupled Matrix and Tensor Factorization 8 C ONCLUSION

(CMTF) to capture explicitly similarities (common coupled In this research, we propose to discover both explicit and im-
factor) between the tensor and the matrix. Weighted Non- plicit similarities between datasets across domains. We pro-
negative Matrix Co-Tri-Factorization (WNMCTF) [34] ex- pose to share common and preserve specific latent variables
tended CMF’s idea to non-negative matrix tri-factorization in coupled factors to better capture the true explicit char-
[35]. acteristics of datasets across domains. Besides, we present
an idea to discover non-coupled factors’ latent similarities
7.2 CodeBook Transer (CBT) and make use of them to further improve the recommenda-
tion. Moreover, on the non-shared dimension, we propose
Targeting on improving recommendation on one domain to use the middle matrix of the tri-factorization to match
by utilizing similar latent rating patterns from another do- the unique factors, and align the matched unique factors
main, Li et al. [23] suggested one as a source domain and to further transfer cross-domain implicit similarities and
another one as a target domain. In this research, the authors thus further improve the recommendation. The advantages
employed matrix tri-factorization to decompose the source of our ideas, validated with real-world datasets, suggest
into user, item and weighting factors. This weighting factor combining both explicit and implicit similarities has a signif-
defined rating patterns of the users for the items, being icant impact on improving the performance of cross-domain
used as the codebook to be transferred to the target. In recommendation. The empirical results also encourage us
this case, the explicit similarities were defined as the rating to extend our idea to a generalized model which enables
patterns transferred from the source to the target to improve joint analysis of both similarities from multiple correlated
recommendation of the target. An extension, named Rating- datasets. In the future, we would like to extend our model to
Matrix Generative Model (RMGM) [36], combined both the very high dimensional tensors using distributed processing.
steps of extracting the codebook and transferring it. In
case there are more than two datasets, Moreno et al. [37]
introduced multiple explicit codebooks shared among them. ACKNOWLEDGMENT
Each codebook was also assigned a different weight to uti- This research is partially funded by the Natural Science
lize the correlations better. Even though codebook transfer Foundation of China (Grant No. 61602141) and the Sci-
worked on any datasets, they are only effective when all ence and Technology program of Zhejiang Province (No.
datasets have the same rating patterns. This condition limits 2018c04001).
their applications to a broader available data.
R EFERENCES
7.3 Cluster-Level Latent Factor Model (CLFM)
[1] Y. Koren, R. Bell, and C. Volinsky, “Matrix factorization techniques
The assumption that two datasets from different domains for recommender systems,” Computer, 2009.
have the same rating patterns would be unrealistic in [2] P. Lops, M. de Gemmis, and G. Semeraro, “Content-based rec-
ommender systems: State of the art and trends,” in Recommender
practice. This intuition motivated Gao et al. [24] to pro- Systems Handbook, 2011.
pose a generalized codebook sharing model by introducing [3] Y. Koren and R. Bell, “Advances in collaborative filtering,” in
CLFM. In specific, CLFM only shared common parts of Recommender Systems Handbook, 2011.
[4] W. Pan, E. W. Xiang, N. N. Liu, and Q. Yang, “Transfer learning in
the rating patterns between datasets. By doing so, CLFM collaborative filtering for sparsity reduction,” in Proceedings of the
model learned the only explicitly shared latent parts be- AAAI Conference on Artificial Intelligence, 2010.
tween datasets. When the number of commonly shared [5] W. Chen, W. Hsu, and M. L. Lee, “Making recommendations from
parts equals decomposition rank, CLFM is CBT. In all other multiple domains,” in Proceedings of the ACM SIGKDD International
Conference on Knowledge Discovery and Data Mining, 2013.
cases, CLFM model learns the only the shared latent space [6] C.-Y. Li and S.-D. Lin, “Matching users and items across domains
while preserving the unique characteristics of each data. to improve the recommendation quality,” in Proceedings of the ACM
This model, thus, can be used to joint analyze multiple SIGKDD International Conference on Knowledge Discovery and Data
Mining, 2014.
datasets to overcome the sparsity of each of them. [7] M. Jiang, P. Cui, X. Chen, F. Wang, W. Zhu, and S. Yang, “So-
cial recommendation with cross-domain transferable knowledge,”
IEEE Transactions on Knowledge and Data Engineering, 2015.
7.4 Coordinate System Transfer (CST) [8] Y. Wei, Y. Zheng, and Q. Yang, “Transfer knowledge between
cities,” in Proceedings of the ACM SIGKDD International Conference
All of the above methods were based on an equally shared on Knowledge Discovery and Data Mining, 2016.
correlation between datasets. Basing on observations that [9] D. Yang, J. He, H. Qin, Y. Xiao, and W. Wang, “A graph-based
the rating matrix for recommendation is often sparse while recommendation across heterogeneous domains,” in Proceedings of
the ACM International on Conference on Information and Knowledge
auxiliary data of the users and items can sometimes be Management, 2015.
found, Weike et al. [4] suggested that common coordinates, [10] Y.-F. Liu, C.-Y. Hsu, and S.-H. Wu, “Non-linear cross-domain
e.g., users tastes, would exist between the main matrix collaborative filtering via hyper-structure transfer,” in Proceedings
and its side data. As a result, these coordinates could be of the International Conference on Machine Learning, 2015.
[11] T. Iwata and T. Koh, “Cross-domain recommendation without
transferred from the auxiliary data to the main rating matrix. shared users or items by sharing latent vector distributions,” in
The author proposed Coordinate System Transfer (CST) Proceedings of the International Conference on Artificial Intelligence and
model to construct principal coordinates in low-dimensional Statistics, 2015.
space for the users and items. It then utilized these explicit [12] H. Jing, A. C. Liang, S. D. Lin, and Y. Tsao, “A transfer probabilistic
collective factorization model to handle sparse data in collabora-
similarities on principal coordinates as regularization terms tive filtering,” in Proceedings of the IEEE International Conference on
on the target data. Data Mining, 2014.
14
[13] L. Zhao, S. J. Pan, E. W. Xiang, E. Zhong, Z. Lu, and Q. Yang, [37] O. Moreno, B. Shapira, L. Rokach, and G. Shani, “Talmud: Transfer
“Active transfer learning for cross-system recommendation,” in learning for multiple domains,” in Proceedings of the ACM Interna-
Proceedings of the AAAI Conference on Artificial Intelligence, 2013. tional Conference on Information and Knowledge Management, 2012.
[14] L. Hu, J. Cao, G. Xu, L. Cao, Z. Gu, and C. Zhu, “Personalized
recommendation via cross-domain triadic factorization,” in Pro-
ceedings of the International Conference on World Wide Web, 2013.
[15] B. Wang, M. Ester, Y. Liao, J. Bu, Y. Zhu, Z. Guan, and D. Cai, “The
million domain challenge: Broadcast email prioritization by cross-
domain recommendation,” in Proceedings of the ACM SIGKDD Quan Do is a PhD student in Advanced An-
International Conference on Knowledge Discovery and Data Mining, alytics Institute, University of Technology Syd-
2016. ney, Australia. He conducts research on data
[16] C. C. Hsu, M. Y. Yeh, and S. d. Lin, “A general framework for mining and publishes papers on tensor fac-
implicit and explicit social recommendation,” IEEE Transactions on torization, recommendation systems and cross-
Knowledge and Data Engineering, 2018. domain learning.
[17] F. Wu, Z. Yuan, and Y. Huang, “Collaboratively training sentiment
classifiers for multiple domains,” IEEE Transactions on Knowledge
and Data Engineering, 2017.
[18] J. D. Zhang, C. Y. Chow, and J. Xu, “Enabling kernel-based
attribute-aware matrix factorization for rating prediction,” IEEE
Transactions on Knowledge and Data Engineering, 2017.
[19] A. P. Singh and G. J. Gordon, “Relational learning via collective
matrix factorization,” in Proceedings of the ACM SIGKDD Interna-
tional Conference on Knowledge Discovery and Data Mining, 2008.
Wei Liu is a Senior Lecturer at the Advanced
[20] E. Acar, T. G. Kolda, and D. M. Dunlavy, “All-at-once optimiza-
Analytics Institute, School of Software, Faculty
tion for coupled matrix and tensor factorizations,” arXiv preprint
of Engineering and Information Technology, the
arXiv:1105.3422, 2011.
University of Technology, Sydney. Before joining
[21] K. Shin, L. Sael, and U. Kang, “Fully scalable methods for dis- UTS, he was a Research Fellow at the Univer-
tributed tensor factorization,” IEEE Transactions on Knowledge and sity of Melbourne and then a Machine Learning
Data Engineering, 2017. Researcher at NICTA. He obtained his PhD from
[22] W. Pan, “A survey of transfer learning for collaborative recom- the University of Sydney. He works in the areas
mendation with auxiliary data,” Neurocomput., 2016. of machine learning and data mining and has
[23] B. Li, Q. Yang, and X. Xue, “Can movies and books collaborate?: been publishing papers in research topics of
Cross-domain collaborative filtering for sparsity reduction,” in tensor factorization, adversarial learning, graph
Proceedings of the International Joint Conference on Artifical Intelli- mining, causal inference, and anomaly detection.
gence, 2009.
[24] S. Gao, H. Luo, D. Chen, S. Li, P. Gallinari, and J. Guo, “Cross-
domain recommendation via cluster-level latent factor model,”
in Proceedings of the European Conference on Machine Learning and
Knowledge Discovery in Databases, 2013.
[25] C. Ding, T. Li, W. Peng, and H. Park, “Orthogonal nonnegative Jin Fan received the B.Sc. degree from Xian
matrix t-factorizations for clustering,” in Proceedings of the ACM Jiaotong University, China, the M.Sc. and Ph.D.
SIGKDD International Conference on Knowledge Discovery and Data degrees from Loughborough University, U.K.
Mining, 2006. She is currently an Associate Professor with the
[26] T. G. Kolda and B. W. Bader, “Tensor decompositions and applica- Department of Computer Science and Technol-
tions,” SIAM Rev., 2009. ogy, Hangzhou Dianzi University, China. Her re-
[27] C. Bauckhage, “k-means clustering is matrix factorization,” in search interests lie in the general areas of mo-
arxiv: 1512.07548, 2015. bile sensing and related data analytics aspects.
[28] Y. Hu, Y. Koren, and C. Volinsky, “Collaborative filtering for
implicit feedback datasets,” in Proceedings of the IEEE International
Conference on Data Mining, 2008.
[29] R. He and J. McAuley, “Ups and downs: Modeling the visual
evolution of fashion trends with one-class collaborative filtering,”
in Proceedings of the International Conference on World Wide Web,
2016.
[30] H. Hotelling, “The generalization of student’s ratio,” The Annals of Dacheng Tao (F’15) is Professor of Computer
Mathematical Statistics, 1931. Science and ARC Laureate Fellow in the School
of Information Technologies and the Faculty of
[31] P. Bhargava, T. Phan, J. Zhou, and J. Lee, “Who, what, when, and
Engineering and Information Technologies, and
where: Multi-dimensional collaborative recommendations using
the Inaugural Director of the UBTECH Sydney
tensor factorization on sparse user-generated data,” in Proceedings
Artificial Intelligence Centre, at the University of
of the International Conference on World Wide Web, 2015.
Sydney. He mainly applies statistics and mathe-
[32] W. Pan, N. N. Liu, E. W. Xiang, and Q. Yang, “Transfer learning to
matics to Artificial Intelligence and Data Science.
predict missing ratings via heterogeneous user feedbacks,” in Pro-
His research interests spread across computer
ceedings of the International Joint Conference on Artificial Intelligence,
vision, data science, image processing, machine
2011.
learning, and video surveillance. His research
[33] Y. Shi, M. Larson, and A. Hanjalic, “Mining contextual movie
results have expounded in one monograph and over 500+ publications
similarity with matrix factorization for context-aware recommen-
at prestigious journals and prominent conferences, such as IEEE T-
dation,” ACM Trans. Intell. Syst. Technol., 2013.
PAMI, T-NNLS, T-IP, JMLR, IJCV, NIPS, ICML, CVPR, ICCV, ECCV,
[34] J. Yoo and S. Choi, “Weighted nonnegative matrix co-tri- ICDM; and ACM SIGKDD, with several best paper awards, such as the
factorization for collaborative prediction,” in Proceedings of the best theory/algorithm paper runner up award in IEEE ICDM’07, the best
Asian Conference on Machine Learning: Advances in Machine Learning, student paper award in IEEE ICDM’13, the distinguished student paper
2009. award in the 2017 IJCAI, the 2014 ICDM 10-year highest-impact paper
[35] D. D. Lee and H. S. Seung, “Algorithms for non-negative matrix award, and the 2017 IEEE Signal Processing Society Best Paper Award.
factorization,” in Advances in Neural Information Processing Systems He received the 2015 Australian Scopus-Eureka Prize, the 2015 ACS
13, 2001. Gold Disruptor Award, and the 2015 UTS Vice Chancellor’s Medal for
[36] B. Li, Q. Yang, and X. Xue, “Transfer learning for collaborative Exceptional Research. He is a Fellow of the IEEE, AAAS, OSA, IAPR,
filtering via a rating-matrix generative model,” in Proceedings of and SPIE.
the International Conference on Machine Learning, 2009.

Unveiling Hidden Implicit Similarities For Cross-Domain Recommendation

Uploaded by

Document Informationclick to expand document information

Copyright:

Available Formats

Unveiling Hidden Implicit Similarities For Cross-Domain Recommendation

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Unveiling Hidden Implicit Similarities For Cross-Domain Recommendation

Uploaded by

Copyright:

Available Formats

This article has been accepted for publication in a future issue of this journal, but has not been

Unveiling Hidden Implicit Similarities for

Index Terms—Cross-Domain Learning; Recommendation System; Matrix Factorization.

1 I NTRODUCTION require adequate ratings to analyze the latent similarities of

High income Low crime

Low income High crime

idea can be extended to multi-sources (Sect. 6). 2.2 Our notation

3 P ROBLEM DEFINITION 4.1 Explicit similarities - Sharing common and preserv-

Fig. 4: (Best viewed in color) Suppose X(2) is obtained from another

Fig. 6: An illustration of how centroid of a cluster is computed. U

matched user cluster’s centroid in U(2) ; an example of how

However, if the user clusters across domains are matched

(1) T (2) T (2) T

each. Cholesky decomposition of r × r matrices is O(r3 ).

1. ABS: http://www.abs.gov.au/websitedbs/censushome.nsf/ 2. BOCSAR: http://www.bocsar.nsw.gov.au/Pages/bocsar crime

Rank Dataset CMF CBT CLFM CST HISF1 HISF2

Dataset Rank CMF CBT CLFM CST HISF1 HISF2

c parameter in Figure 9, 10. These figures illustrate the case

Rank Dataset CMF CBT CLFM HISF-N

this case, demonstrate our idea can be generalized to the

parts are quite substantial in addition to the common parts.

authors proposed Coupled Matrix and Tensor Factorization 8 C ONCLUSION

You might also like