Object embedding using an information geometrical perspective

Acquiring vector representations of objects is essential for applying machine learning, statistical inference, and visualization. Although various vector acquisition methods have been proposed considering the relationship between objects in target data, most of them are supposed to use only a specific relevance level. In real-world data, however, there are cases where multiple relationships are contained between objects, such as time-varying similarity in time-series data or various weighted edges on graph-structured data. In this paper, a vector acquisition method which assigns vectors in a single coordinate system to objects preserving the information given by multiple relations between objects is proposed. In the proposed method, a logarithmic bilinear model parameterized by representation vectors is utilized for approximating relations between objects based on a stochastic embedding idea. The inference algorithm proposed in this study is interpreted in terms of information geometry: the m-projection from the probability distribution constructed from observed relations on the model manifold and the e-mixture in the model manifold are alternately repeated to estimate the parameters. Finally, the performance of the proposed method is evaluated using artificial data, and a case study is conducted using real data.


Introduction
The acquisition of vector representations for objects is useful for information processing.Vector representations allow the visualization of relationships between objects and can be used as a preprocessing tool for statistical inference and machine learning algorithms.However, there are many cases, such as images and sentences, that are difficult to represent in fixed-length vectors or cases where only relationships between objects, such as relevance level, are observed.In particular, we propose a method for obtaining an appropriate vector representation in situations where multiple relevance degrees are given for each object.For example, consider the case where a movie is rated in each country at year.In this case, multiple relevance degrees are given for the movie.Each country has different preferences for movies, and these preferences change over time.We propose a method to obtain a vector representation that takes into account multiple degrees of relevance for data with multiple observations over time.In another example, we can discuss cultural differences and trends in scientific papers.It is assumed that each journal or international conference has a different culture, as different words are used in different venues, and that they also change over time.The proposed method is expected to be applied to the analysis of such data.This example case study will be discussed in detail in Sect. 4. Such as graph data and image data, there are cases where only distance and similarity between objects are observed.In this paper, such distances and similarities are treated as special cases of relevance and discussed comprehensively.The proposed method can be applied to obtain vector representations for various data, including time-series data for which similarities are given at different times and graph-structured data for edges of which multiple pattern weights are assumed.
In the following, a set of objects is denoted by O, and individual objects are denoted by o, with the subscript i where necessary to distinguish them.The problem of acquiring a vector x to an object o to approximate the relationships between the objects is generally called object embedding or vectorization.This problem can be divided into two types: the first is "constellation", in which a vector representation is assigned to an individual object, And the second is "embedding", in which a function is constructed to give the vectors to the entire set of target objects, For the latter embedding problem, the problem of learning the mapping φ from the data using a function approximator, such as a neural network, is considered.In this paper, we deal in particular with the former constellation problem.
For most of the methods that construct the vectors independently for each given relevance level is proposed and the comparison or manipulation of estimated vectors from multiple relevance levels is not assumed.In the example of movies and journals mentioned earlier, they could be used to analyze a single year or each year separately, but they are not intended for analysis of data that is acquired across years, i.e., multiple times.The individual constellations for multiple relevance levels will result in vectors that do not preserve the data characteristics of the given relevance levels.
Therefore, we propose a constellation method for multiple relevance levels.In this paper, we assign the vectors to objects using the stochastic embedding approach represented by Stochastic Neighborhood Embedding (SNE) [1], introduced in Sect. 2. If we reformulate the problem setting in terms of stochastic embedding, we can view the current problem as a problem of estimating the appropriate probability distribution to represent them in a situation where multiple empirical distributions are available for each object.The approach of simply estimating the appropriate distribution from each empirical distribution cannot capture the relationship between the empirical distributions.It is necessary to express the relationship between empirical distributions by imposing constraints on the target space.A proposed logarithmic bilinear model (logbilinear model) approximates the relevance levels using two vectors, one for the objects and the other for the features as parameters.One of the properties of the log-bilinear model is that by fixing one parameter, the space stretched from the remaining parameters becomes the space of exponential families.In this paper, we use this property to represent the relationship between multiple empirical distributions by sharing some of the parameters of the probability distributions to be estimated.We also show that the weighted sum of the vectors obtained can represent a mixture of probability distributions in the space of probability distributions and can be easily visualized after estimation.
The contributions of this paper are described as follows.
• Our proposed framework can provide vectors for even multiple relevance levels.
• Our algorithm can be interpreted as a simple operation in terms of information geometry.• The proposed method enables vectorization even when additional data are generated.
The structure of this paper is as follows.Section 2 provides an overview of existing mapping methods based on similarity and discrepancy, comparing classical and probabilistic methods.In Sect.3, we first introduce auxiliary features for objects and then propose approximate representations using a log-bilinear model with two vectors of the features and objects, and its optimization, on the basis of an information geometry interpretation.We evaluate the effectiveness of the proposed method and its visualization using artificial and real data in Sects.4, and 5, we summarize the proposed method and discuss the remaining issues on its possible applications.

Related works
In this section, we provide an overview of existing methods of mapping vectors based on relevance level.First, let us organize the symbols.The set of objects is denoted by O, and the number of objects is represented by n.Individual objects are distinguished by subscript i as Each object is given a d-dimensional vector x ∈ R d , and a matrix X is defined as the matrix of vectors x corresponding to object o, arranged in rows as column vectors, Here, we introduce a matrix R ∈ R n×m to estimate the matrix X that represents the relevance between objects and features.We assume that object o and feature u have an relevance rel(o, u) ∈ R between them, and estimate the vector x.Let U be a set of features.
where (R) i j is a (i, j)th entry of matrix R. As mentioned above, similarity and distance can be regarded as special cases of relevance R, in which case U = O and (R) i j can be regarded as indicating their values.
In the following, we describe methods for estimating the matrix X from the relevance R. First, Torgerson's classical metric multidimensional scaling (MDS) is a technique to obtain a low-rank approximation X using the similarity matrix as a following loss function, where • F denotes the Frobenius norm [2].Note that R is a symmetric matrix where each element represents the degree of similarity between objects.On the other hand, Kruskal's non metric MDS [3] uses the Minkowski distance, a generalization of the Euclidean distance, to find the matrix X that approximates the relevance matrix R, where • p denotes the Minkowski distance of order p and φ(•) is a suitable monotonic transformation.In this case, the relevance R is a matrix such that each element represents the distance between objects.There are also other methods such as robust Euclidean embedding [4], which improves robustness to outliers by evaluating differences between relevance levels and estimated vectors with an L 1 -loss function, and non metric MDS, which optimizes ordinal statistics [3].Since essential information is often located in a low-dimensional space in practical problems using the high-dimensional data, various techniques have been actively studied to reconstruct it appropriately in the low-dimensional space.The difference among various methods is how they incorporate information from a neighborhood that is assumed to contain important information.For example, Isomap [5], which recalculates the geodetic distances from a neighborhood, semidefinite embedding (SDE) [6] using kernel matrix, locally linear embedding (LLE) [7], and laplacian eigenmaps [8].In the practical problem of acquiring a vector representation of a word, a low-rank approximation of the word-document matrix is proposed to acquire the vectors [9,10].
The approximation error was evaluated using the norm of the matrix in the classical method described above.On the other hand, a method called stochastic embedding has been proposed, which evaluates the reconstruction error of vectors based on the statistical distance between probability distributions.
Stochastic neighborhood embedding (SNE) is the pioneer of stochastic embedding.In this method, the vectors are mapped by expressing the probability of an object o being in the neighborhood of o as a conditional probability mass function p o (o ), Note that f is a function of the similarity between objects.Let x be the vector of an object o, and we denote the conditional probability q that is, an object o is a neighborhood of an object o , where g is a function that transforms such that the smaller the distance, the higher probability.
The statistical distance is employed as the evaluation function between probability distributions in a stochastic embedding.The Kullback-Leibler divergence (KL divergence) is generally defined as one of the statistical distances and defined as ( The statistical distance between the two probability distributions is used as one of loss functions with respect to the matrix X in a stochastic embedding,

123
Various methods have been proposed for the stochastic embedding, depending on a choice of statistical distances and similarity calculation methods to capture the relevance levels between objects.t-SNE [11], a variant of SNE, uses a function whose distribution after transformation is the t-distribution to improve estimation robustness and speed up computation.Other proposed methods include UMAP [12], which applies the theory of algebraic topology, and word2vec [13] which estimates the dual vectors using the asymmetry of the conditional probability defined from the similarity.GloVe [14] is also a highly relevant and key technique that is also related to the proposed method, which proposes vector acquisition by focusing on the ratio of conditional probabilities of objects i, j.The problem setting of this paper is related to embedding for time-series data.Dynamic t-SNE is an extension of t-SNE for mapping vectors to time-series data [15].In addition to the conventional t-SNE loss function, this method enables the mapping of objects by adding penalties to differences between times.Several methods have been proposed to enable the application of Word Embedding to time-series data, called "Dynamic Word Embedding" [16][17][18].Bunte et al. describes the these methods from a unified perspective [19].
Various techniques have been proposed to construct vectors that approximate relevance, although this problem can be viewed as a low-rank approximation problem of the relevance matrix R by a matrix X , R X X T .
The classical method, metric MDS uses the evaluation of approximation errors by the Frobenius norm and results in a singular value decomposition of a symmetric matrix, In this case, approximation errors are evaluated irrespective of the similarity owing to the nature of the Frobenius norm.On the other hand, the stochastic embedding approach considers the probabilistic models defined by relevances and vectors as and evaluates the approximation error on the basis of the statistical distance as follows, In this case, the nature of the statistical distance gives more weight to relevance levels between objects with high similarity, i.e., high probability values.Various frameworks have been proposed for geometric approaches, such as Poincaré embedding [20], which is a method to improve performance by changing the calculation method of the distance between two objects and embedding them in the hyperbolic space.Related research in information geometry is described below, along with the proposed method.

Proposed methods
Since the vectors are obtained in existing methods for a single relevance matrix, as described in Sect.2, they are not designed for use of multiple matrices.Therefore, it is difficult to compare the obtained vectors when the existing method applies to multiple matrices.In this section, we propose a method of mapping vectors that can be used when multiple relevance levels are obtained by assuming two types of vectors: object vectors and feature vectors.We can compare the object vectors using shared feature vectors among multiple relevance levels.
The symbols are described here.Let O = {o i } n i=1 be the set of objects and U = {u j } m j=1 be the set of features.For example, in the example of movie ratings by country described in the introduction, O is the set of movies and U is the set of countries.The function rel(o, u) gives the relevance of an object o ∈ O to a feature u ∈ U.In the case of the movie data example, the rating of the movie o in country u is defined as rel(o, u).
In this section, we propose a method to estimate object vectors when k relevance matrices R (1) , . . ., R (K ) are given.First, consider the case of a single relevance matrix.The proposed method approximates the relevance of an object and a feature by the inner product of vectors, given vector x i of the object o i and the vector y j of the feature u j , Here, considering the two matrices X and Y with each vector x i and y j aligned, Equation ( 2) can be treated as a problem, which approximates the relevance matrix R by the matrices X and Y .In the following, R assumes that the value of (R) i j takes on larger values when o i and u j are strongly related.
Next, we define the probability distributions P and Q using the relevance matrix R and the matrices X and Y to evaluate the approximation accuracy of the decomposition in the framework of stochastic embedding.They indicate the probability that an object o and a feature u are observed simultaneously.The proposed method approximates the probability mass function p(o, u) of the empirical distribution P using a log-bilinear model on the basis of the inner product of the vectors x and y to correspond to the object o and the feature u, respectively, where φ(R) and ζ(X , Y ) are the normalization terms for P(R) and Q(X , Y ) to be the probability distributions.The KL divergences of the joint probability distribution is used as the evaluation function for this problem, This concept of relying on a log-bilinear model is also the idea used in word2Vec [13] and GloVe [14].However, minimizing the KL divergence for multiple relevance matrices does not allow us to estimate the appropriate vectors.If A is a regular matrix, then Since R of Eq. ( 3) has the same value as both a pair of X and Y , and a pair of X and Y , the value of Eq. ( 6) give the same solution, This means that the solution to this problem is arbitrary in terms of regular matrices.It is difficult for given multiple relevance levels R (k) K k=1 to discuss the relationship between their solutions X (k) for this reason.We can derive an operation such that the convex combination of matrices X (k) can be regarded a mixture of probability distributions Q (k) by using the shared matrices Y for any k, We can extend the method to estimate the matrices {X (k) } for multiple relevance matrices by this shared matrix Y .Let us consider a dividing point x at vectors x (1)  and x (2) estimated using the relevance R (1) and R (2) , x = αx (1) + (1 − α)x (2) , 0 ≤ α ≤ 1.
The probability mass function q(o, u; X , Y ) of the joint probability distribution Q( X , Y ) by x can be written as (1) (o, u; X (1) , Y ) α q (2) (o, u; X (2) , Y ) where q (k) is a probability mass function of Q (k) , and exp(−ζ ( X , Y )) is the normalization term.The distribution Q is a convex combination of two distributions Q (1) and Q (2) under log transformation, Therefore, the distribution Q( X , Y ) is a mixture distribution called the e-mixture.Since the e-mixture emphasizes regions where both distributions have a high probability, it is intuitively an AND-like operation that takes the common part of the distributions.In this way, the mixture of the vectors corresponds to an operation of the distributions.The details of "e-mixture" are described in the next section.In the proposed method, we can estimate vectors for multiple relevance matrices by fixing Y .In this case, the matrices {X (k) }, Y are not uniquely determined.In this method, we reduce the degree of freedom of the solution by imposing the following constraints, where I d ∈ R d×d is the identity matrix.Note that even with this condition, the estimation of X and Y remains arbitrary in the orthogonal matrices.Given multiple relevance levels, The optimization problem is written as follows for the matrices is written as with respect to X (k) , ∀k and Y In this paper, we imposed the constraint Eq. ( 8), but it could be any other constraint to reduce the degree of freedom of the solution.

Preliminary on information geometry
Information geometry is a framework for discussing the behavior of statistical inference and machine learning using a geometric perspective.A wide variety of methods and algorithms have been understood and extended by information geometry.For example, Akaho proposed a method for obtaining vectors by dimension reduction in [21].As another an information geometry approach to object embedding, Riccard et al. propose alpha embedding, which computes between two stochastic distributions using the computed alpha-divergence [22].In both studies, the information geometry approach is used to extend the conventional methods for embedding and dimension reduction of a single data.As previously discussed, when multiple datasets are obtained, it is difficult to directly apply the object embedding approach to consider the relationship between multiple data, so it is necessary to extend the method to consider the relationship.As previously mentioned, directly applying the object embedding approach to deal relationships between multiple datasets is challenging.In this section, we describe the basic ideas necessary to explain the proposed method, especially the handling of statistical models and their estimation in information geometry [23].
Let Z be the sample space and Z be the random variable.And let S be the space containing an arbitrary probability distribution, and P be the probability distribution whose probability mass function is p(z), The set of probability distributions with parameter θ is defined as which is called the model manifold M(θ ).Let z 1 , . . ., z N be the observed data sampled from a probability mass function p with appropriate parameters θ * .The observed data can be regarded as a probability distribution using a delta function called the empirical distribution, Estimating the parameter θ refers to finding the closest probability distribution Q on the subspace M from the empirical distribution P using the statistical distance D in information geometry [23].The statistical distance D refers to the pseudo-distance between two probability distributions, P and Q, and here we employ the KL divergence D KL (P Q) in Eq. ( 1).
We define a set G m and G e of interior points for the two probability distributions P (1) and P (2) on the space S of probability distributions, where α ∈ [0, 1].The subspace G m expressed in Eq. ( 10) is called the m-geodesics, and G e expressed in Eq. ( 10) is called as the e-geodesics, given K probability distributions, the subspace M m is defined as the m-flat subspace spanned by their convex combination, and the subspace M e is defined as the e-flat subspace spanned by a convex combination of logarithmic representations.
where α = α (1) , . . ., α (K ) T ∈ R K , α (k) ≥ 0 and K k=1 α (k) = 1.And p(x) is called the m-representations, log ( p(z)) is called the e-representations, and the operations to construct P m and P e by the convex combination for each representation are called the m-mixture and the e-mixture, respectively.
We now consider the KL divergence for the three probability distributions P, Q, and R, where r is the probability mass function of R. When the right-hand side of Eq. ( 11) is zero, that is, when the m-geodesics from probability distribution P to R and the e-geodesics from R to Q are orthogonal, the Pythagorean theorem holds This is useful for describing parameter estimation in the stochastic models.We now introduce a set of probability distributions called an exponential family to describe the geometric properties of parameter estimation in statistical models.An exponential family is a set of probability distribution whose probability mass function can be written in the following form: where θ are parameters of P, and ζ(θ) is the normalization term as Probability distributions in the exponential family of distributions are e-flat [23].
For an e-flat subspace M e , if there exists probability distributions P / ∈ M e , Q and R ∈ M e .Let these three points satisfy Eq. ( 12).Then the point Q that is the shortest point from P to M e is equivalent to R. because the m-geodesics from P to R and e-geodesics from R to Q are orthogonal by the above Pythagorean theorem.And this process of estimating R from P is called the m-projection, R = arg min Q∈M e D KL (P Q).
In the next section, we explain how the proposed method can be interpreted in terms of information geometry for the case where multiple empirical distributions P are observed.

Algorithms for information geometric interpretation
This section describes an information geometrical interpretation of the proposed model.We use a vec: R N ×M → R N M and a Kronecker product operator ⊗ : where B is a U × V matrix.The probability that an object o and a feature u are observed simultaneously, can be generalized as Equation ( 13) is equivalent to Eq. ( 5) given o i and u j as a concrete event, Here, δ o i , δ u j , δ o , and δ u are Using the vec operator and the Kronecker product, we can transform Eq. ( 13) as where is an nd × md matrix and is an n × m matrix defined as We see that the probability mass function q is expressed in Eq. ( 14) is an exponential family with respect to θ with ξ fixed, and we see that it is an exponential family with respect to ξ with θ T fixed as well.Here, we write v = (o, u)ξ and v = θ T (o, u), respectively, where θ and ξ are We propose an optimization process of Eq. ( 9) on the basis of these results.The proposed algorithm estimates vectors from multiple relevance levels by a combination of the m-projection to the model manifold and the e-mixture on the model manifold.
Probability distributions P (1) , P (2) , . . ., P (K ) computed from multiple relevance levels are represented by a single point in the space of probability distributions S as empirical distributions.Similarly, it is represented as a single point on S, given the parameters X and Y .In other words, this problem is equivalent to finding an appropriate mapping from the empirical distribution P (k) to the model manifold M(X , Y ), where ξ † = vec(Y † ).The subspace M Y † (X ) is e-flat with respect to the parameter X because it is a member of the exponential family.This means that the optimal distribution in Q from the empirical distribution P is given by the m-projecteion.
Similarly let M X † (Y ) be a subspace of M(X , Y ) where X is fixed to X † , and Y is freely chosen, where θ † = vec(X † ).
The convex combination of the matrices on the model manifold is an e-mixture of probability distributions as shown in Eq. (7).Therefore, the model manifold M Ȳ computed from Ȳ is expected to preserve the properties of Y (k) , The feature matrix Y is constrained by Eq. ( 8) to limit the degree of freedom.The matrix satisfying this constraint is called a Stiefel manifold and which it defined as The mixture matrix Ȳ given in Eq. ( 15) is not necessarily a matrix on this manifold.Therefore, we introduce a retraction function π : Ȳ → Y ∈ St.
There are various types of retraction for matrix A, and we use the method based on Q R decomposition [24], where D is a orthogonal matrix, and U is an upper triangular matrix.The proposed algorithm estimates parameters separately for each data set and integrates them, which is computationally less expensive than estimating all parameters simultaneously.However, efficient estimation when the amount of objects in a dataset is large is a future challenge.In particular, depending on the retraction method, it may become a bottleneck in computational complexity and should be chosen carefully.For the QR decomposition selected this time, O(n 3 ), it is necessary to consider other options when the data size is large.
We revisit the problem formulation from an information geometric perspective.The task involves projecting multiple empirical distributions onto a large model manifold while preserving the relationships between these empirical distributions.
To provide an intuitive explanation of the proposed algorithm, our approach involves constructing a large model manifold parameterized by X and Y .By applying constraints that fix a subset of parameters, we project the empirical distributions onto a subspace within the large model manifold.By iteratively performing this process, we expect to accurately estimate the parameters while retaining the characteristics of each empirical distribution.
The concept of an algorithm of the proposed method is shown in Fig. 1.Algorithm 1 is the proposed algorithm for estimating {X (k) } and Y .These figures describe the algorithm as follows.First, it is assumed that {P (k) } 3 k=1 are obtained, which are computed from multiple relevance levels {R (k) } 3 k=1 , respectively.In addition, the parameters of the iteration index t are given as {X t+1 } gives the subspace M Ȳt+1 (X ).In Fig. 1d, Y t+1 is derived by the retraction of Ȳt+1 onto a Stiefel manifold to satisfy the constraints of Eq. ( 8).
The matrices {X (k) } and Y estimated by this process maintains the neighborhood relations of {P (k) }.The shared Y also allows the matrices of objects to be compared across matrices.

Algorithm 1 an algorithm of the proposed mapping method
Require:  In addition to introducing two coordinates as in the proposed method when multiple relevance levels are given, another possible approach is to estimate them collectively, for instance, by vertically concatenating them and treating them as a single dataset.In this case, X (k) can be estimated directly, but the computational cost becomes substantial when the number of steps and data points increases.Therefore, the proposed method, which allows estimating the problem with a smaller size for each relevance level, offers an advantage.The weighted sum of the coordinates estimated by the proposed method also represents a mixture of probability distributions, enabling adjustments weights to accommodate subsequent tasks and visualizations after the vectors are estimated.

Appearance and disappearance of data
The proposed method can naturally estimate the vectors of time-series data, even when the data appear with the transition of time.Let P (1) , P (2) , and P (3) be the empirical distributions at time k = 1, 2, and 3, respectively.And the set of objects at each time be denoted by O , O , and O .In addition, the number of objects at each time is The proposed method works equally well for the time-series data regardless of the number of objects by assuming that the number of the features The estimated model are the subspaces M Y fix (X ) ⊂ M (X , Y ), and M Y fix (X ) ⊂ M (X , Y ) respectively, with the additional constraint that Y = Y fix at each time.Therefore, the proposed method can be applied similarly to Algorithm 1 in estimating the object matrices X (k) .The case of k = 3, the mapping is processed in the same way.

Experiments
Generally, when considering embedding real data, only the relevance is given, not the actual vector that is the correct answer.We evaluate the effectiveness of the proposed method using an artificial data with the true matrices X (k) * and Y * , and also show as a case study that the proposed method works and behaves intuitively on real data.We compare the proposed method with three other methods: t-SNE [11] applied independently at each time(ind t-SNE), t-SNE applied to all data at once(concat t-SNE), and Dynamic t-SNE [15], which is a probabilistic embedding method that considers the time series.The evaluation criteria will be employed Spearman's ranking correlation and MAP@H .We introduce some symbols used for MAP@H .Let N (o i , X , X , h) be a set of h neighbors of object o using the object matrix X , Here, o(X (k) , h) is the h th nearest object from o using the object matrix X (k) .MAP@H evaluates the match between the neighborhoods of the true vectors and the estimated vectors as We define the criteria to evaluate whether the information is preserved between time series.Normally, MAP@H is calculated inside the kth relevance matrix (k = k), but in this case, we also calculate MAP@H using the kth relevance matrix and the k − 1 th relevance matrix to verify whether the information between relevance is retained, which we call MAP@H in this paper.

Mapping using artificial time-series relevance levels
We axevaluate the mapping of time-series data from multiple relevance levels.Each object belongs to one of the clusters and moves toward the center of each cluster over time.We map the vectors to this dataset and evaluate how well the characteristics hold.
Let K be the time length and c(i) be the cluster to which each object o i belongs.The speed at which the object moves toward the center is defined for each cluster and is written as v(c(i)).And artificial time-series data are generated using Algorithm 2.
The quantitative evaluation are shown in Table 1.Experimental settings not listed in the Table 1 include K = 8, m = 10, and also H = 50 for evaluation.The results of the proposed method show that for d ≥ d * , the ranking correlation and the MAP@H are high, indicating that the original relationship is reconstructed.When the true data dimension exceeds the estimated dimension, i.e., in d * ≥ d, the values of each evaluation index are low.The evaluation metric, MAP@H , represents the relationships among the estimated vectors in each relevance matrix.When t-SNE is applied, we observe a significant decline in this value, suggesting a compromised representation of these relationships.In contrast in the case of our proposed method, this value does not decrease, indicating that our method successfully maintains the relationships among the estimated vectors within each relevance matrix.Thus, it is suggested that our proposed method effectively considers the relationships among the estimated vectors in each relevance matrix.However, it has been observed that the proposed method is less accurate than other methods when the dimension of the true data is large.Other methods show that MAP@H is higher, and ranking correlation is lower than those obtained by the proposed method, because they expresses the neighborhood relationship of each object rather than the whole relationship.The results show that even time-series data can be used to estimate vectors that reproduce each relevance levels.The mapping results for time-series data generated as true dimension d * = 2 and time length K = 8 are shown in Fig. 3: We have shown that the proposed method preserves the relationship between times for the true vectors, while reproducing each relevance level.

Case study: DBLP dataset
In this case study, we will use the DBLP dataset [25] for real data.This dataset consists of meta and citation information on 4,107,340 papers.
This section presents the results of a case study on providing vectors using this dataset.The data for providing the vectors is created by aggregating this dataset using the fields in Table 2. Let O be the set of terms in the title and U be the venue on Fig. 3 The color-coded matrices for each cluster at each time k.In the independent t-SNE, which does not consider the relationships among different relevance levels, the vectors are mapped to distinct positions at each time step, even when they belong to the same cluster.By comparing the results with two other methods that account for the relationships among relevance levels, it becomes evident that the proposed method produces outcomes most similar to the ground truth data.While other methods converge at a uniform speed, the proposed method in particular captures differences in the speed of movement to the center and appropriately represents the cohesiveness of the clusters at varying times which the paper was published.The preprocess to create a dataset whose parameters are n = 300, m = 30, K = 12 and {R (k) } K k=1 , the start year is 2006, and the last year is 2017 shown in Fig. 4.
The mapping results for d = 30 are shown in Fig. 5. Here, we selected the most commonly used words for each venue at each time point in 2006 and 2017, applied PCA on them, and then visualized the results in two dimensions.
This figure shows that "deep" has significantly changed its position near "neural" and "network" from 2010 to 2017, while the positions of other words have not changed much.The reason for the movement of these words is assumed to be that "deep learning" has been discussed in various papers in recent years.This method can capture the trend changes given multiple relevance levels.We showed that it is possible to simultaneously acquire high-dimensional vectors that retain relevance levels and enable the visualization of specific objects of interest in low dimensions by using PCA after the vectors are acquired in this experiment.Furthermore, it also enables the comparison of the mapping as time-series data by fixing the feature matrix Y for the features at all times.

Conclusion
In this paper, we proposed a method that provides vectors to an object on the basis of multiple relevance levels defined between the objects and features using a stochastic embedding.It is difficult to interpret the results of existing SNE and t-SNE, which are used for the stochastic embedding when applied independently, because they are not designed to acquire vectors to multiple relevance levels (Fig. 3).
We solved these problems using the log-bilinear model combined vectors of the objects and features.In this paper, we showed that the estimation of the object vectors using the common feature vectors is equivalent to mapping to the e-flat subspaces on the model manifold.In this method, the manipulation of the vectors corresponds to the manipulation of the probability distributions, since the convex combination of vectors is the e-mixture of the probability distributions.
We evaluated the effectiveness of the proposed method using artificially generated data and the DBLP dataset as a real-world time series data.The results show that the proposed method can provide the appropriate vectors to the objects while preserving the trend changes for the multiple relevance levels to the artificial data and the real-world data.This method enables the visualization of relationships that change over time, such as texts and videos.The establishment of a baseline for evaluating embedding methods across multiple datasets will be facilitated by defining versatile classification tasks utilizing multiple datasets in the future.When the goal is visualization itself, it is also expected that quantitative evaluations of the estimated vectors will be conducted through means such as questionnaires and interviews.From here on, we will describe the gradients for X and Y .

A.1 Gradient in the X direction
We first calculate for x h and then show the gradient for the whole X , Thus, the gradient of X is induced as

A.2 Gradient in the Y direction
We similarly derive the gradients of y h and Y , Thus, the gradient of X is similarly induced as Y t .Figure 1a firstly shows the matrices {X (k) t+1 } are estimated by the m-projection from {P (k) } to a subspace M Y t (X ) on the model manifold with Y t fixed.{Y (k) t+1 } is estimated by the m-projection from {P (k) } to the model manifold with {X (k) t+1 } fixed in Fig. 1b. Figure 1c describes that a convex combination of {Y (k)

Fig. 1
Fig. 1 An optimization process.The first step a shows that is the m-projection from P (k) to M Y t under with it fixed Y t .In the b step, each X (k) t+1 is fixed, and Y (k) t+1 is estimated with the second m-projection.c Next, we mix {Y t+1 } to calculate Ȳt+1 by using the e-mixture.d Finally, Y t+1 is derived of Ȳt+1 by the retraction onto the Stiefel manifold

Fig. 2
Fig. 2 Geometric conceptual figure of data appearance.At time k = 1, M is used as the model manifold, but at times k = 2, 3, it is extended to M and M owing to the increase in the amount of objects.In addition, the vectors are estimated with the model manifold constrained by fixing Y = Y fix

Fig. 5
Fig. 5Two-dimensional representation of the relationship between the terms in the titles of papers published between 2010 and 2017

∂
∂ y h L KL (P Q(X , Y )) i y h x i .

Table 1
Mapping results using multiple relevance level

Table 2
Data schema of DBLP dataset(except)