Transition analysis of boundary-based active configurations in temporal simplicial complexes for ingredient co-occurrences in recipe streams

Aiming at knowledge discovery for temporal sequences of cooking recipes published in social media platforms from the viewpoint of network science, we consider an analysis of temporal higher-order networks of ingredients derived from such recipe streams by focusing on the framework of simplicial complex. Previous work found interesting properties of temporal simplicial complexes for the human proximity interactions in five different social settings by examining the configuration transitions before and after triplet interaction events corresponding to 2-simplices. In this paper, as an effective extension of the previous work to the case of higher dimensional n-simplices corresponding to newly published recipes, we propose a novel method of configuration transition analysis by incorporating the following two features. First, to focus on changes in the topological structure of temporal simplicial complex, we incorporate analyzing the transitions of boundary-based configurations. Next, to focus on the temporal heterogeneity in usage activities of ingredients, we incorporate analyzing the transitions of active configurations by introducing the activity degree of configuration. Using real data of a Japanese recipe sharing site, we empirically evaluate the effectiveness of the proposed method, and reveal some characteristics of the temporal evolution of Japanese homemade recipes published in social media from the perspective of ingredient co-occurrences.

The increasing popularity of social media platforms for sharing cooking recipes is enabling us to investigate creative homemade recipes of ordinary people. Recently, the interest in food science and computing has been growing (Min et al. 2019), and the data analysis methods from network science have also been used to explore the co-occurrence properties of ingredients in recipes (Ahn et al. 2011; Teng et al. 2012; Jain et al. 2015). Here, networks provide a fundamental framework for modeling complex systems, and network-based analysis has been successfully applied to various areas including social media analytics and behavioral and social sciences (Barabási 2016). In traditional network-based models, graphs have been widely used, where vertices represent the elementary units of the underlying system, and edges encode their pairwise interactions and relationships. However, in many real-world systems, it can be important to analyze interactions among more than two units (Benson et al. 2019. For instance, people mostly communicate in groups, more than two species can interact in an ecosystem, social topics can often be explained by multiple keywords, and a cooking recipe usually consists of more than two ingredients. For modeling such higher-order interactions, attention has recently been devoted to using hypergraphs and simplicial complexes since traditional graphs cannot give faithful representations (Preti et al. 2021). Hypergraphs are one of the most natural generalizations of graphs, where the concept of an edge is generalized to a hyperedge described as a subset of vertices. Simplicial complexes are a special case of hypergraphs, satisfying the downward closure property (Preti et al. 2021), where the generalization of an edge is called a simplex. Theory of simplicial complexes has a close relation with algebraic topology (Croom 2007), and modeling frameworks based on simplicial complexes have been successful in developing new insights into several research problems including brain organization (Saggar et al. 2018), protein interaction (Estrada and Ross 2018) and social influence spreading (Iacopini et al. 2019). In this paper, by focusing on the framework of a temporal simplicial complex, we consider the dynamical properties of higher-order relationships among ingredients appearing within a recipe stream, consisting of recipes with time-stamps, that is generated through diverse user interactions on a social media platform dedicated to sharing homemade recipes (see the “Datasets and experimental settings” section in the “Experiments” section for more details on such recipe stream data).

For an analysis of growing simplicial complexes, Benson et al. (2019) explored a higher-order link prediction problem for 19 datasets from various domains in terms of simplicial closure, and revealed fundamental differences between higher-order networks (i.e., simplicial complexes) and traditional dyadic networks (i.e., traditional graphs). Unlike link prediction such as simplicial closure, Cencetti et al. (2021) found interesting properties of the temporal simplicial complexes derived from human proximity interactions for five different social settings in terms of configuration transitions. More specifically, for each event corresponding to a 2-simplex (i.e., an interaction among exactly three people), they analyzed the configurations among the three people involved one step before the interaction and one step after it, and revealed the characteristics for transitions of the configurations. However, in the case of n-simplices with \(n \ge 3\), such configuration transition analysis was not performed. In fact, this analysis method can intrinsically suffer from combinatorial explosion, i.e., the number of possible configurations can become very large as the dimension n of simplices to be considered gets increasingly larger. Thus, it is desirable to develop an effective and efficient approach for configuration transition analysis.

In this paper, we effectively extend the previous work to higher dimensional n-simplices, and consider the problem of analyzing configuration transitions for a temporal simplicial complex of ingredient co-occurrences derived from a recipe stream. Here, first assuming that nice homemade foods can be made through fascinating combinations of ingredients, we focus on each combination of the n ingredients for a recipe consisting of \((n+1)\) ingredients as a first-step. Next, assuming that ingredients used in recipes can strongly depend on season and trend, we also take temporal usage patterns of ingredients into account. We thus present a novel method of configuration transition analysis from the following perspectives: First, aiming at effectively capturing changes in the basic qualitative features of ingredient combinations in recipes, we focus on changes in the topological structures of simplicial complexes before and after occurrences of n-simplices corresponding to newly published recipes, and propose analyzing the transitions of boundary-based configurations according to homology theory (Croom 2007) from algebraic topology, instead of examining all possible configurations. Next, for the n-simplices corresponding to new recipes, we propose analyzing the transitions of active configurations around them by introducing the activity degree of configuration to focus on the temporal heterogeneity in usage activities of ingredients in the recipe stream. Using real data from a social media site for sharing Japanese recipes, we empirically evaluate the proposed method for configuration transition analysis, and extract characteristics of the temporal evolution of Japanese homemade recipes in terms of ingredient combinations.

This paper extends our conference paper (Fujisawa et al. 2023) in which we presented the basic idea for analyzing the transitions of boundary-based active configurations between the preceding and forthcoming periods for occurrences of n-simplices, and showed the analysis results for a part of the dataset of Japanese recipe sharing site “Cookpad”^{Footnote 1} only in the cases of 3-simplices (i.e., \(n = 3\)) and 4-simplices (i.e., \(n = 4\)) by setting one month as the preceding and forthcoming periods. In this paper, we first provide the explanation of our configuration transition analysis method in more detail. By improving the analysis method, we further propose such a framework that is able to effectively analyze the whole picture of the configuration transitions around n-simplices over all dimensions n. We also evaluate the proposed analysis framework in other parts of the Cookpad dataset instead of repeating the experimental evaluation for the same dataset as in the conference paper (Fujisawa et al. 2023). To show the significance of incorporating the activity degree of configuration, we newly compare it with a baseline method derived as a straightforward extension of the previous work (Cencetti et al. 2021). Moreover, we examine the effect of the length of the preceding and forthcoming periods, which is an important parameter for analysis.

The paper is organized as follows. The “Related work” section briefly mentions related work. As preliminaries, the “Preliminaries” section briefly describes the previous work (Cencetti et al. 2021) for configuration transition analysis of human social interactions, explains the temporal simplicial complexes derived from recipe streams, and introduces several notations used in the later sections. The “Analysis framework” section explains the proposed analysis framework in detail, and the “Experiments” section reports the experimental results for real data. The “Conclusion” section summarizes the main achievement and future plans.

Food science and computing have been attracting attention in recent years (Min et al. 2019). As an analysis from the perspective of complex network science, pairwise relationships between ingredients or cuisines were explored. Several researchers investigated dyadic networks of ingredients for various cuisines in terms of flavor compounds (Ahn et al. 2011; Jain et al. 2015; Park et al. 2021; Makinei and Hazarika 2022). For example, Jain et al. (2015) analyzed ingredient networks in terms of spices for Indian cuisines. To make a comparative analysis of cuisines in the world, Sajadmanesh et al. (2017) examined cuisine networks in terms of ingredients and flavors. On the other hand, West et al. (2013) investigated population-wide dietary preferences from recipe queries on the Web, Jiang et al. (2017) compared food cultures in the world by jointly visualizing recipe density and ingredient categories, and Min et al. (2018) made a cross-region analysis of culinary culture by automatically extracting cuisine-course topics of recipes on the basis of ingredient combinations. Also, a variety of food-oriented applications are being explored. For example, complement and substitute networks of ingredients were analyzed toward recipe recommendation (Teng et al. 2012), and it was investigated how algorithmic solutions relate to the healthiness of recipes (Trattner and Elsweiler 2017). However, as far as we know, there have been few attempts to analyze higher-order relationships among ingredients for home recipes. In this paper, we explore the characteristics for the temporal evolution of higher-order relationships among ingredients for Japanese homemade recipes published in social media.

This paper is also related to link prediction tasks. Much effort has been devoted to studying link prediction for traditional dyadic-networks (Liben-Nowell and Kleinberg 2007; Leskovec et al. 2010; Lu and Zhou 2011). However, only a few attempts have been made at higher-order link prediction for large real-world data since it can be computationally challenging. For example, Xu et al. (2013) provided HPLSF which is a method of predicting higher-order links of an arbitrary order for social networks using latent features. By restricting a set of candidate higher-order links, Zhang et al. (2018) presented Coordinated Matrix Minimization (CMM) which is a higher-order link prediction method based on adjacency space, and empirically showed that CMM outperforms several baselines including HPLSF. As mentioned before, to reduce the computational load of analyzing growing simplicial complexes, Benson et al. (2019) focused on simplicial closure phenomena and examined higher-order link prediction in terms of simplicial closure. On the other hand, there is an increasing interest in analyzing simplicial complexes from the perspective of data science. In fact, such an analysis framework has been successfully applied to various real problems (Saggar et al. 2018; Estrada and Ross 2018; Iacopini et al. 2019), and novel generalizations of methods for traditional graphs are being developed for simplicial complexes. For instance, Shaub et al. (2020) studied diffusion processes, random walks and Laplacians on simplicial complexes, and Preti et al. (2021) developed an efficient algorithm for truss decomposition of simplicial complexes. Also, random simplicial complex models have been mathematically investigated (Bobrowski and Krioukov 2022). In this paper, as with Cencetti et al.’s work (2021), we analyze the configuration transitions in a temporal simplicial complex that are associated with new higher-order events.

We consider a recipe stream \({\mathcal R}\) from a social media platform of sharing recipes, where \({{\mathcal {R}}}\) is a set of recipes with time-stamps. For each recipe \(r \in {\mathcal R}\), let t(r) denote the time-stamp of r. Here, we measure t(r) by using one day as unit of time. We explore the evolution of higher-order relationships among ingredients obtained from the discrete-time event history \({\mathcal R}\) through a configuration transition analysis before and after the publication of a new recipe in \({\mathcal R}\).

Previous work

We begin with revisiting the previous work (Cencetti et al. 2021) for configuration transition analysis of temporal evolution of higher-order social interactions.

Unlike link prediction tasks such as simplicial closure (Benson et al. 2019), Cencetti et al. (2021) explored the temporal evolution of the higher-order structure of human proximity interactions in terms of configuration transition. They said that a group of n individuals is formed if there are all pairwise interactions among them during some time period. Here, higher-order interactions among people happen in a group of them. For each event corresponding to the formation of a group of n individuals, they addressed the relational structures around the n people at the preceding and forthcoming periods. In particular, they focused on group formations involving exactly three individuals (i.e., triplet interaction events), and investigated the transition rates of the relational structures around three people before and after a triplet interaction event. Here, as for possible configurations around three individuals, there are the following four classes. In the first class, there are no pairwise interactions among them, i.e., there are no groups including two of them. In the second class, there is only a single pairwise interaction among them. In the third class, there are only two pairwise interactions among them. In the fourth class, there is a group including all of them. For five different social settings, they empirically demonstrated that the configuration transitions towards and from social triplet interaction events are characterized by the above second and third classes (Cencetti et al. 2021).

Example for topological structures of simplicial complexs. Consider 2-simplex \(S = \{ v_0, v_1, v_2\}\), and its boundary-faces \(S_0 = \{ v_1, v_2 \}\), \(S_1 = \{ v_0, v_2 \}\), \(S_2 = \{ v_0, v_1 \}\). Suppose that simplicial complexes at time-step 0 and time-step 1 are \({\mathcal K}_0 = \{ S_0, S_1, S_2, \{ v_0 \}, \{ v_1 \}, \{ v_2 \} \}\) and \({\mathcal K}_1 = \{ S, S_0, S_1, S_2, \{ v_0 \}, \{ v_1 \}, \{ v_2 \} \}\), respectively. Then, the first homology groups are given by \(H_1 ({\mathcal K}_0) = {\mathbb Z}\) and \(H_1 ({\mathcal K}_1) = \{ 0 \}\). Namely, \([{\mathcal K}_0]\) has one 1-dimensional hole, while \([{\mathcal K}_1]\) has no 1-dimensional holes. This implies that the topological structures of \([{\mathcal K}_0]\) and \([{\mathcal K}_1]\) are different, and \({\mathcal K}_0\) and \({\mathcal K}_1\) are qualitatively distinct

Full size image

Temporal simplicial complex from recipe stream

Given a finite set of vertices V, a k-simplex S is a subset of V with \(|S| = k + 1\) for a non-negative integer \(k < |V|\), and a simplicial complex \({\mathcal K}\) is a set of simplices which satisfies the following downward closure property (see Schaub et al. 2020; Preti et al. 2021): For every simplex \(S \in {\mathcal K}\), all subsets of S belong to \({\mathcal K}\). Here, for each subset \(S'\) of simplex S with \(|S'| = k' + 1\), \(S'\) is called a \(k'\)-face of S, and S is called a coface of \(S'\). Also, for a k-simplex S, k is referred to as the dimension of S. Simplicial complex \({\mathcal K}\) corresponds to a polytope \([{\mathcal K}]\) called its geometric realization, and is studied in the field of algebraic topology (Croom 2007). For an integer k with \(0< k < |V|\) and a k-simplex \(S = \{ v_0, v_1, \dots , v_k \}\), the boundary \(\partial S\) of S consists of its \((k-1)\)-faces

which are \((k-1)\)-simplices in \({\mathcal K}\) (see Fig. 1). We refer to \(S_j\) as boundary-face of S for each \(j = 0, 1, \dots , k\). By algebraically examining the simplices in \({\mathcal K}\) and their boundaries, the \(k'\)th homology group of \({\mathcal K}\), \(H_{k'}({\mathcal K})\), is defined for any integer \(k'\) with \(0 \le k' \le |V|\) as a global geometric structure of the polytope \([{\mathcal K}]\) corresponding to \({\mathcal K}\) (Croom 2007). It is known that \(H_0 ({\mathcal K})\) algebraically expresses information of the connected components in \([{\mathcal K}]\), and \(H_{k'}({\mathcal K})\) algebraically represents information of the \(k'\)-dimensional holes (i.e., the holes with a \(k'\)-dimensional boundary) in \([{\mathcal K}]\) for each \(k' \ge 1\) (see Fig. 1). In this paper, for configurations around new simplices appearing in a temporal simplicial complex, we consider focusing on their boundary-faces from the viewpoint of the topological structures of simplicial complexes through the homology theory.

As for recipe stream \({\mathcal R}\), we choose a set of ingredients as V, and consider the set of recipes at each time-step t,

We focus on a simplicial complex \({\mathcal K}_t\), which is defined from \({\mathcal R}_t\) through ingredient co-occurrences in the following way: A subset \(S = \{ v_0, \dots , v_k \}\) of V with \(|S| = k+1\) is called a k-simplex at time-step t if there is an \(r \in {\mathcal R}_t\) such that ingredients \(v_0, \dots , v_k\) are included in recipe r. \({\mathcal K}_t\) is defined as the set of all simplices at time-step t. We refer to \(\{ {\mathcal K}_t \}\) as the temporal simplicial complex derived from recipe stream \({\mathcal R}\).

For an arbitrary recipe \(r \in {\mathcal R}\), let \(\sigma (r)\) denote the set of all the ingredients included in recipe r. We refer to \(\sigma (r)\) as the simplex corresponding to recipe r, or a recipe simplex. Note that \(\sigma (r) \in {\mathcal R}_{t(r)}\). We investigate the configuration transitions associated with occurrences of recipe simplices in temporal simplicial complex \(\{ {\mathcal K}_t \}\).

Illustration of temporal change for recipe simplex \(\sigma (r)\)

Full size image

Boundary-based configuration for recipe simplex

For a recipe simplex

of dimension n in temporal simplicial complex \(\{ {\mathcal K}_t \}\), we focus on a period of time just before the occurrence of \(\sigma (r)\),

and a period of time just after it,

where a positive integer \(\tau\) is a parameter of our analysis, and indicates the length of investigation periods \(I^B(r)\) and \(I^A(r)\) (see Fig. 2). Let \({\mathcal K}^B(r)\) and \({\mathcal K}^A(r)\) be the simplicial complexes during periods \(I^B(r)\) and \(I^A(r)\), respectively,

For temporal simplicial complex \(\{ {\mathcal K}_t \}\), we consider the configurations around \(\sigma (r)\) within investigation periods \(I^B(r)\) and \(I^A(r)\), respectively. Note that if we straightforwardly apply the idea of the previous work (Cencetti et al. 2021) to this case, we must examine the k-faces of \(\sigma (r)\) for any integer k with \(0 \le k < n\) and all the cofaces of \(\sigma (r)\). To avoid this combinatorial issue, we choose not to examine all the possible configurations. Alternatively, aiming to make a connection with the topological structures of \({\mathcal K}^B(r)\) and \({\mathcal K}^A(r)\) through homology theory, we focus on its boundary-faces \(\sigma _0(r)\), \(\sigma _1(r)\), \(\dots\), \(\sigma _{n}(r)\) for the configurations around \(\sigma (r)\) (see Fig. 3), where

To this end, we introduce a configuration feature vector within \(I^B(r)\) (see Fig. 2), which is defined as

and a configuration feature vector within \(I^A(r)\) (see Fig. 2), which is defined as

Now, we define \({\varvec{p}}^B(r)\) and \({\varvec{p}}^A(r)\) in explicit detail. First, \(p^B_h (r)\) and \(p^A_h(r)\) are the probabilities that a recipe simplex including \(\sigma (r)\) occurs within \(I^B(r)\) and \(I^A(r)\), respectively. Next, for each \(j = 0, 1, \dots , n\), \(p^B_j (r)\) and \(p^A_j(r)\) are the probabilities that a recipe simplex including only one \(\sigma _j(r)\) among \(\sigma _0(r), \sigma _1(r), \dots , \sigma _n(r)\) occurs within \(I^B(r)\) and \(I^A(r)\), respectively (see Fig. 3). Finally, \(p^B_\ell (r)\) and \(p^A_\ell (r)\) are the probabilities that a recipe simplex including none of \(\sigma _0(r), \sigma _1(r), \dots , \sigma _n(r)\) occurs within \(I^B(r)\) and \(I^A(r)\), respectively. As for the probabilities described above, we empirically estimated them in our experiments.

Example for the boundary-faces of a recipe n-simplex \(\sigma (r)\) (\(n=3\))

Full size image

To take into account the temporal heterogeneity in usage activities of ingredients in recipe stream \({\mathcal R}\), we also introduce a concept of activity degree. To this end, we consider a sufficiently long period of time before the occurrence of \(\sigma (r)\),

and introduce a configuration feature vector within \({{\tilde{I}}}^B(r)\), which is defined as

(see Fig. 2). Here, the length \({\tilde{\tau }}\) of investigation period \({{\tilde{I}}}^B(r)\) is a positive integer with \({\tilde{\tau }} \gg \tau\), and \({\tilde{\varvec{p}}}^B(r)\) is defined within \({{\tilde{I}}}^B(r)\) in the same way as \({\varvec{p}}^B(r)\) within \(I^B(r)\). Note that unlike \({\varvec{p}}^B(r)\), \({\tilde{\varvec{p}}}^B(r)\) represents the ordinary feature for the co-occurrence of ingredients \(v_0(r)\), \(v_1(r)\), \(\dots\), \(v_n(r)\) in recipes before time-step t(r) if \({{\tilde{I}}}^B(r)\) is a sufficiently long period. In the experiments, we set \({\tilde{\tau }}\) to 365 days (i.e, one year). We are also interested in how configuration feature vectors \({\varvec{p}}^B(r)\) within \(I^B(r)\) and \({\tilde{\varvec{p}}}^B(r)\) within \({{\tilde{I}}}^B(r)\) affect configuration feature vector \({\varvec{p}}^A(r)\) within \(I^A(r)\) (see Fig. 2).

For temporal simplicial complex of ingredient co-occurrences \(\{ {\mathcal K}_t \}\) derived from recipe stream \({\mathcal R}\), we propose an analysis framework of configuration transition.

Transition analysis of active configurations

For each positive integer n with \(0< n < |V|\), let \({\mathcal R}(n)\) denote the subset of \({\mathcal R}\) such that \(\sigma (r)\) is an n-simplex for any \(r \in {\mathcal R}(n)\). We first provide an analysis method of configuration transition associated with occurrence of \(\sigma (r)\) in \(\{ {\mathcal K}_t \}\) for any \(r \in {\mathcal R}(n)\) under fixed n, and by appropriately improving it, we propose an analysis framework of configuration transition for the whole \({\mathcal R}\).

For any recipe n-simplex \(\sigma (r) = \{ v_0(r), v_1(r), \dots , v_n(r) \}\) (\(r \in {\mathcal R}(n)\)), we explore boundary-based configurations around \(\sigma (r)\) within investigation periods \(I^B(r)\) and \(I^A(r)\) by measuring the activity degree of configuration on the basis of the configuration feature vectors \({\varvec{p}}^B(r)\), \({\varvec{p}}^A(r)\) and \({\tilde{\varvec{p}}}^B(r)\). First, for \(I^B(r)\) which is the time period just before the occurrence of \(\sigma (r)\), we consider the decomposition of \({\mathcal R}(n)\),

and define the boundary-based active configuration \({\mathcal C}(\sigma (r); I^B(r))\) around \(\sigma (r)\) within \(I^B(r)\) as follows: If \(p^B_h(r) > {{\tilde{p}}}^B_h(r)\), then \(r \in {\mathcal R}^B_h (n)\) and we say that the boundary-based active configuration around \(\sigma (r)\) within \(I^B(r)\) is \(\sigma (r)\) itself, i.e., \({\mathcal C}(\sigma (r); I^B(r))\) \(=\) \(\{ \sigma (r) \}\). Note that \(r \in {\mathcal R}^B_h (n)\) means that simplices including \(\sigma (r)\) actively occur within \(I^B(r)\), and thus \(\sigma (r)\) itself is active within \(I^B(r)\). In the case of \(r \notin {\mathcal R}^B_h (n)\), we examine the inequality \(p^B_j(r) > {{\tilde{p}}}^B_j(r)\) for each boundary-face \(\sigma _j(r)\) (\(j = 0, 1, \dots , n\)), and set

Let \(\lambda\) denote the number of the boundary-faces belonging to \({\mathcal C}(\sigma (r); I^B(r))\), i.e., \(\lambda\) \(=\) \(| {\mathcal C}(\sigma (r); I^B(r)) |\). If \(\lambda > 0\), then \(r \in {\mathcal R}^B_\lambda (n)\) and we say that the boundary-based active configuration around \(\sigma (r)\) within \(I^B(r)\) is the set of \(\lambda\) active boundary-faces \({\mathcal C}(\sigma (r); I^B(r))\) obtained by Eq. (2). Note that \(r \in {\mathcal R}^B_\lambda (n)\) implies that \(\lambda\) boundary-faces of \(\sigma (r)\) actively occur within \(I^B(r)\). If \(\lambda = 0\), then \(r \in {\mathcal R}^B_\ell (n)\) and we say that the boundary-based active configuration around \(\sigma (r)\) within \(I^B(r)\) is the empty set, i.e., \({\mathcal C}(\sigma (r); I^B(r))\) \(=\) \(\emptyset\). Note that \(r \in {\mathcal R}^B_\ell (n)\) means that none of the boundary-faces of \(\sigma (r)\) actively occur within \(I^B(r)\).

Next, for \(I^A(r)\) which is the time period just after the occurrence of \(\sigma (r)\), we define the decomposition of \({\mathcal R}(n)\),

and the boundary-based active configuration \({\mathcal C}(\sigma (r); I^A(r))\) around \(\sigma (r)\)within \(I^A(r)\) in the same way as the case for \(I^B(r)\). Here, if \(p^A_h(r) > {{\tilde{p}}}^B_h(r)\), we define \(r \in {\mathcal R}^A_h (n)\) and \({\mathcal C}(\sigma (r); I^A(r))\) \(=\) \(\{ \sigma (r) \}\). In the case of \(r \notin {\mathcal R}^A_h (n)\), we set

and \(\lambda\) \(=\) \(| {\mathcal C}(\sigma (r); I^A(r)) |\). If \(\lambda > 0\), we define \(r \in {\mathcal R}^A_\lambda (n)\), and say that the boundary-based active configuration around \(\sigma (r)\) within \(I^A(r)\) is the set of \(\lambda\) active boundary-faces \({\mathcal C}(\sigma (r); I^A(r))\) obtained by Eq. (3). If \(\lambda = 0\), we define \(r \in {\mathcal R}^A_\ell (n)\), and \({\mathcal C}(\sigma (r); I^A(r))\) \(=\) \(\emptyset\).

Five classes of the boundary-based active configurations for a recipe 2-simplex

Full size image

Therefore, for any recipe n-simplex \(\sigma (r) = \{ v_0(r), v_1(r), \dots , v_n(r) \}\), we divide each of boundary-based active configurations \({\mathcal C}(\sigma (r); I^B(r))\) and \({\mathcal C}(\sigma (r); I^A(r))\) into the following \((n+2)\) classes in the proposed analysis framework (see Fig. 4): First, we say that \({\mathcal C}(\sigma (r); I^B(r))\) is in class h if \(r \in {\mathcal R}^B_h (n)\), and \({\mathcal C}(\sigma (r); I^A(r))\) is in class h if \(r \in {\mathcal R}^A_h (n)\). Next, for each \(k = 1, \dots , n+1\), we say that \({\mathcal C}(\sigma (r); I^B(r))\) is in class \(\lambda = k\) if \(r \in {\mathcal R}^B_k (n)\), and \({\mathcal C}(\sigma (r); I^A(r))\) is in class \(\lambda = k\) if \(r \in {\mathcal R}^A_k (n)\). Finally, we say that \({\mathcal C}(\sigma (r); I^B(r))\) is in class \(\ell\) if \(r \in {\mathcal R}^B_\ell (r)\), and \({\mathcal C}(\sigma (r); I^A(r))\) is in class \(\ell\) if \(r \in {\mathcal R}^A_\ell (r)\). Note first that \(\sigma (r)\) itself is active within \(I^B(r)\) and \(I^A(r)\) if \({\mathcal C}(\sigma (r); I^B(r))\) and \({\mathcal C}(\sigma (r); I^A(r))\) are in class h, respectively. For each \(k = 1, \dots , n+1\), \(\sigma (r)\) has k active boundary-faces within \(I^B(r)\) and \(I^A(r)\) if \({\mathcal C}(\sigma (r); I^B(r))\)and \({\mathcal C}(\sigma (r); I^A(r))\) are in class \(\lambda = k\), respectively. Moreover, \(\sigma (r)\) has no active boundary-faces within \(I^B(r)\) and \(I^A(r)\) if \({\mathcal C}(\sigma (r); I^B(r))\) and \({\mathcal C}(\sigma (r); I^A(r))\) are in class \(\ell\), respectively.

For a fixed integer n with \(0< n < |V|\), we make a detailed analysis for the transition of boundary-based active configurations associated with occurrences of recipe n-simplices by using the transition matrix T(n),

To effectively analyze the whole picture of the configuration transitions around recipe n-simplices across all dimensions n, we propose an analysis framework of using the aggregated transition vector

where

Note that \(AT_b(n)\) indicates the occurrence probability of recipe n-simplices around which the configuration transitions are comprehensively characterized by their active boundary-faces. We examine each component of AT(n) as a function n.

Influence analysis

We arbitrarily fix the dimension n of recipe simplices to be considered, and explore the relationship of configuration feature vector \({\varvec{p}}^A(r)\) within \(I^A(r)\) with configuration feature vectors \({\varvec{p}}^B(r)\) within \(I^B(r)\) and \({\tilde{\varvec{p}}}^B(r)\) within \({{\tilde{I}}}^B(r)\) for each recipe n-simplex \(\sigma (r) = \{ v_0(r), v_1(r), \dots , v_n(r) \}\).

In view of situations before and after occurrence of \(\sigma (r)\), we naturally speculate that we might have \({\varvec{p}}^A(r) \sim {\varvec{p}}^B(r)\), \({\varvec{p}}^A(r) \sim {\tilde{\varvec{p}}}^B(r)\), or \({\varvec{p}}^A(r) \propto {\varvec{p}}^B(r) / {\tilde{\varvec{p}}}^B(r)\) (i.e., \(p^A_h(r) = \nu \, p^B_h(r)/ {{\tilde{p}}}^B_h(r)\), \(p^A_j(r) = \nu \, p^B_j(r)/ {{\tilde{p}}}^B_j(r)\) (\(j = 0, 1, \dots , n\)), and \(p^A_\ell (r) = \nu \, p^B_\ell (r)/ {{\tilde{p}}}^B_\ell (r)\) for some \(\nu >0\)). In this paper, we examine the influence of \({\varvec{p}}^B(r)\) and \({\tilde{\varvec{p}}}^B(r)\) on \({\varvec{p}}^A(r)\) by modeling

where \(c > 0\), and

is the model parameter called the weight vector

By assuming a multinomial model within \(I^A(r)\) based on the probability vector \({\varvec{p}}^A(r)\) defined by Eq. (7), we estimate the weight vector \({\varvec{w}}\) from the observed data vector within \(I^A(r)\),

for any \(\sigma (r) = \{ v_0(r), v_1(r), \dots , v_n(r) \} \in {\mathcal R}(n)\) according to the MAP estimation framework. Here, let \({\mathcal S}^A(r;n)\) be the set of such recipes in \({\mathcal R} (n)\) that are published within \(I^A(r)\), i.e.,

Then, \({\varvec{m}}(r)\) is defined as follows:

First, \({\mathcal S}^A_h(r;n)\) is the set of recipes \(r' \in {\mathcal S}^A(r;n)\) such that recipe simplex \(\sigma (r')\) includes \(\sigma (r)\), i.e.,

Next, for each \(j = 0, 1, \dots , n\), \({\mathcal S}^A_j(r;n)\) is the set of recipes \(r' \in {\mathcal S}^A(r;n)\) such that \(r' \notin {\mathcal S}^A_h(r;n)\) and recipe simplex \(\sigma (r')\) includes boundary-face \(\sigma _j(r)\), i.e.,

Finally, \({\mathcal S}^A_\ell (r;n)\) is the set of recipes \(r' \in {\mathcal S}^A(r;n)\) such that \(r' \notin {\mathcal S}^A_h(r;n)\) and \(r' \notin {\mathcal S}^A_j(r;n)\) for any \(j = 0, 1, \dots , n\), i.e.,

In the assumed generative model, the probability \(P({\varvec{m}}(r) \, | \, {\varvec{w}})\) of observing \({\varvec{m}}(r)\) \(=\) \((m_h(r)\), \(m_0(r)\), \(m_1(r)\), \(\dots\), \(m_n(r)\), \(m_\ell (r))\) within \(I_A(r)\) is given by

We assume a Gaussian prior for \({\varvec{w}}\) and estimate weight vector \({\varvec{w}}\) by maximizing the function

where \(\mu > 0\) is a hyper-parameter, \(\Vert {\varvec{w}}\Vert\) stands for the Euclidean norm of \({\varvec{w}}\), and \({\mathcal L}({\mathcal R}(n); {\varvec{w}})\) is the log-likelihood of \({\mathcal R}(n)\) for weight vector \({\varvec{w}}\),

(see Eq. 8). Here, by Eq. (7), we have

where \(c' > 0\). Note that maximizing \(F ({\varvec{w}})\) reduces to a sort of softmax optimization problem in neural network [see Eqs. (9), (10) and (11)]. In the experiments, we employed a gradient-based method (Bishop 1995).

Trends in the number of recipes posted and the number of Cooksnaps received for the Dessert, Meat-dish and Vegetable dish datasets. The number of such posts within each week is displayed for the period from Jan 1, 2010 to Feb 28, 2013

Full size image

We conducted an empirical evaluation of the proposed analysis framework using actual temporal simplicial complexes of ingredients obtained from Cookpad, a Japanese recipe-sharing service. Our analysis focuses on the dynamical properties of higher-order relationships involving more than two ingredients.

Number of recipe n-simplices for the Dessert (2012), Meat-dish (2012) and Vegetable-dish (2012) datasets

Full size image

Datasets and experimental settings

We utilized Cookpad data from Jan 1, 2010 to Feb 28, 2013, and constructed datasets for recipe stream \({\mathcal R}\), which correspond to the recipe categories of Cookpad, a social meidia platform dedicated to sharing Japanese homemade recipes. In particular, we focused on the three main categories of Cookpad, “Dessert”, “Meat-dish” and “Vegetable-dish”, and explored the temporal evolution of higher-order relationships among ingredients for recipes appearing within \({\mathcal R}\). In Cookpad, users can post their original recipes. Also, for each posted recipe, other users can send a “Thank You” message with a dish’s photo when they actually cooked and loved it, and this type of message is called a Cooksnap. Since the number of Cooksnaps for a recipe is seen as a measure of its popularity, it is natural for Cookpad users to aspire to create their own original recipes that receive a large number of Cooksnaps. Thus, they should try to create unique and fascinating recipes by taking into account seasonal variations and current trends. Moreover, when they come across a recipe that has received several Cooksnaps, they may be inspired to create similar ones. Consequently, it is considered that Cookpad users are likely to interact in this manner. To confirm the presence of such interactions among users for each of the Dessert, Meat-dish and Vegetable-dish datasets, we examined the temporal nature of recipe stream \({\mathcal R}\). Figure 5 shows the number of recipes posted and the number of Cooksnaps received within each week as a function of week for the three datasets. We see that the number of Cooksnaps received is much larger than the number of recipes posted for the three datasets. These results suggest that there are active interactions among Cookpad users through Cooksnaps for all three datasets. Therefore, we can assume that the three datasets for \({\mathcal R}\) are generated through such interactions among Cookpad users, and it can be important to investigate the characteristics of recipe streams for the three datasets in terms of revealing the evolution of Japanese homemade recipes.

We separately investigated the configuration transitions for the recipes r published in 2011 and the recipes r published in 2012. In our conference paper (Fujisawa et al. 2023), we presented the analysis results for the 2011 Cookpad datasets in the cases of \(n = 3, 4\) for the dimension n of recipe simplex \(\sigma (r)\) and \(\tau =30\) days for the length \(\tau\) of investigation periods \(I^B(r)\) and \(I^A(r)\), while in this paper, we have conducted more detailed investigations. To avoid duplication, we only report the results for the 2012 Cookpad datasets, “Dessert (2012)”, “Meat-dish (2012)” and “Vegetable-dish (2012)”, in this paper. For each recipe stream \({\mathcal R}\), we adopted the set of its major ingredients^{Footnote 2} as a set of vertices V. Then, we had \(|V| = 976\), \(|V| = 753\), and \(|V| = 1,023\) for the Dessert (2012), Meat-dish (2012) and Vegetable-dish (2012) datasets, respectively. Fig. 6 displays the number of recipe n-simplices as a function of n for the three datasets.

On the other hand, the length \(\tau\) of investigation periods is an important parameter for our configuration transition analysis since varying the value of \(\tau\) can generally affect the emergence of active boundary-faces within \(I^B(r)\) and \(I^A(r)\) for any recipe simplex \(\sigma (r)\). However, we only investigated the case of \(\tau = 30\) days in our conference paper (Fujisawa et al. 2023). In this paper, we examine the effect of \(\tau\) for the cases of \(\tau = 14\) days, \(\tau = 30\) days and \(\tau = 60\) days. Also, in the previous work of Cencetti et al. (2021), social interactions among three people (i.e., for the case of \(n = 2\)) were only explored for configuration transition analysis. In this paper, we thus focus on the case of \(n \ge 3\).

Results of aggregated transition vector AT(n) for the Dessert (2012), Meat-dish (2012) and Vegetable-dish (2012) datasets. The results for \(\tau = 14\) days, \(\tau = 30\) days and \(\tau = 60\) days are displayed, where \(\tau\) is the length of each of the investigation periods \(I^B(r)\) and \(I^A(r)\) for any recipe \(r \in {\mathcal R}\)

Full size image

Results for transition analysis of active configurations

We first analyze the whole picture of the configuration transitions around recipe simplices by using the aggregated transition vector AT(n) over all dimensions n (see Eq. 6). Next, compared with a baseline method straightforwardly derived from the previous work (Cencetti et al. 2021), we describe the detailed analysis results for a particular dimension n in terms of the transition matrix T(n) (see Eq. 4).

Analysis results for aggregated transition vectors

For the real recipe stream \({\mathcal R}\), we examined the aggregated transition vectors \(\{ AT(n) \}\). Fig. 7 shows the results for the three datasets, where each component of AT(n) is displayed as a function of the dimension n of a recipe simplex. As for the length \(\tau\) of each of the investigation periods \(I^B(r)\) and \(I^A(r)\) for any recipe \(r \in {\mathcal R}(n)\) (see Fig. 2), we examined the cases of \(\tau = 14\) days, \(\tau = 30\) days and \(\tau = 60\) days.

We first see from Fig. 7 that as n increases, \(AT_b(n)\) tends to decrease and \(T_{\ell , \ell }(n)\) tends to increase. In particular, \(AT_b(n)\) and \(T_{\ell , \ell }(n)\) are dominant when n is relatively small, and \(T_{\ell , \ell }(n)\) is only dominant when n becomes larger. Note that the exact value of n that marks this shift in dominance can vary depending on the specific datasets and the value of \(\tau\). We also observe that \(T_{\ell , h}(n)\), \(T_{h, \ell }(n)\) and \(T_{h, h}(n)\) are usually small. These imply that when n is relatively small, there are a significant number of recipe n-simplices around which the configuration transitions are comprehensively characterized by their active boundary-faces. However, when n is relatively large, the configuration transitions around recipe n-simplices are no longer able to be comprehensively characterized by their active boundary-faces (i.e., classes \(\lambda = 1, \dots , n+1\)) and themselves active (i.e., class h). These results suggest that for such a Japanese recipe that consists of relatively small number of ingredients and is published on a social media platform, there is a high probability that some of the ingredient combinations obtained by excluding only one ingredient from the recipe can actively appear in recipes within the period before or after its occurrence. This demonstrates the effectiveness of the proposed framework for analyzing boundary-based active configurations.

Next, we consider the effect of \(\tau\). In view of the activity degree, \(T_{\ell , \ell }(n) = 1\) and the other components of AT(n) are zero if \(\tau\) becomes equal to \({\tilde{\tau }}\) (i.e., one year). On the other hand, as \(\tau\) increases, the simplicial complexes \({\mathcal K}^B(r)\) and \({\mathcal K}^A(r)\) for each \(r \in {\mathcal R}(n)\) (see Eq. 1) should become including more diverse simplices since there are an increasing number of recipes contained within \(I^B(r)\) and \(I^A(r)\). From Fig. 7, we see that as \(\tau\) increases up to 60 days, \(AT_b(n)\) tends to increase and \(T_{\ell , \ell }(n)\) tends to decrease. Note that the magnitude of the changes observed can depend on the datasets. These results suggest that in a reasonable range of \(\tau\) such as 14 days \(\le\) \(\tau\) \(\le\) 60 days, it is highly likely that the configuration transitions around recipe simplices of small dimensions can be comprehensively characterized by their active boundary-faces.

Analysis results for transition matrices

For the real recipe stream \({\mathcal R}\), we made a detailed analysis of the configuration transitions for each dimension n in terms of the transition matrix T(n). We also evaluated the significance of incorporating the concept of activity degree through an ablation study, where we consider a baseline method of only examining the presence or absence of its boundary-faces and cofaces within investigation periods \(I^B(r)\) and \(I^A(r)\) for each recipe n-simplex \(\sigma (r)\). Note that the baseline method is obtained by ignoring the activity degree of configuration through only setting \({\tilde{\varvec{p}}}^B(r)\) to the zero vector \({\textbf { 0}}\) in the proposed analysis framework, and is regarded as a straightforward extension of the previous work (Cencetti et al. 2021). In this paper, we only report the results for the Dessert (2012) dataset in the case of \(n = 3\) and \(\tau = 30\) days. Fig. 8 shows the visualization results of transition matrix T(3), where the results of the proposed and baseline methods are displayed in Figs. 8a and b, respectively.

As for the proposed method, we see from Fig. 8a that various entries of T(3) contributed to \(AT_b(3)\). This implies that the configuration transitions characterized by active boundary-faces can have various types. For this dataset, those configuration transitions were mainly related to one or two active boundary-faces (i.e., classes \(\lambda = 1\) and \(\lambda = 2\).). Thus, it is possible to extract interesting phenomena of configuration transitions in the recipe stream as illustrated later in this section. This also indicates the effectiveness of the proposed analysis framework.

From Fig. 8b, we observe that the entry \(T_{h, h}(3)\) obtained by the baseline method was large. This means that for the Dessert (2012) dataset, there were many recipe 3-simplices satisfying the condition that there exist recipe simplices including them before and after they occur. Note that this result of configuration transitions for ingredient co-occurrences is different from the properties of human interactions found by the previous work (Cencetti et al. 2021).

Visualization results of transition matrix T(3) for the Dessert (2012) dataset (\(\tau = 30\) days)

Full size image

We also see from Figs. 8b and b that the transition matrix T(3) obtained by the proposed method was completely different from that obtained by the baseline method. For example, we focus on the entry \(T_{h,4}(3)\) \(=\) \(|{\mathcal R}^B_h(3) \cap {\mathcal R}^A_4(3)| / |{\mathcal R}(3)|\) (see Eq. 4). By the baseline method, we have a recipe \(r' \in {\mathcal R}^B_h(3) \cap {\mathcal R}^A_4(3)\) such that \(\sigma (r')\) \(=\) \(\{v_0(r')=``\textrm{cocoa}''\), \(v_1(r')=``\textrm{butter}''\), \(v_2(r')=``\textrm{chocolate}''\), \(v_3(r')=``\textrm{honey}''\}\). However, the boundary-face \(\sigma _3(r')\) of \(\sigma (r')\) is not active within \(I^A(r')\) since \(p^A_3(r') = 0.33\%\) and \({{\tilde{p}}}^B_3(r') = 1.97\%\). This means that it is quite usual to employ the combination of cocoa, butter and chocolate in a desert recipe. Thus, we have \(r' \notin {\mathcal R}^A_4(3)\) by the proposed method. This example clarifies the difference between the proposed and baseline methods for \(T_{h,4}(3)\). On the other hand, by the proposed method, we can extract a recipe \(r'' \in {\mathcal R}^B_h(3) \cap {\mathcal R}^A_4(3)\) such that \(\sigma (r'')\) \(=\) \(\{v_0(r'')=``\textrm{egg}''\), \(v_1(r'')=``\textrm{cake flour}''\), \(v_2(r'')=``\mathrm{cream~cheese}''\), \(v_3(r'')=``\mathrm{fresh~cream}''\}\). This example reveals an interesting phenomenon as mentioned below. In fact, the configuration transition around \(\sigma (r'')\) is described as follows: Just before \(\sigma (r'')\) occurs, recipe simplices including \(\sigma (r'')\) actively occurred in simplicial complex \({\mathcal K}^B (r'')\). Furthermore, just after \(\sigma (r'')\) occurred, the 3-simplex \(\sigma (r'')\) decays into its four active boundary-faces \(\sigma _0 (r'')\), \(\sigma _1 (r'')\), \(\sigma _2 (r'')\) and \(\sigma _3 (r'')\) in simplicial complex \({\mathcal K}^A (r'')\) (see Fig. 3). Namely, all the four boundary-faces actively occurred in \({\mathcal K}^B (r'')\), but recipe simplices including \(\sigma (r'')\) did not actively occur in \({\mathcal K}^A (r'')\). These results demonstrate that incorporating the concept of activity degree can be significant for exploring the ingredient co-occurrences in the recipe stream, and imply that we can find several interesting properties of temporal simplicial complex \(\{ {\mathcal K}_t \}\) by examining boundary-based active configurations.

Evaluation results of the proposed influence analysis model for the Dessert (2012), Meat-dish (2012) and Vegetable-dish (2012) datasets (\(n = 3, 4\) and \(\tau = 30\) days)

Full size image

Evaluation of influence analysis model

For each recipe n-simplex \(\sigma (r) = \{ v_0(r), v_1(r), \dots , v_n(r) \}\) (\(r \in {\mathcal R}(n)\)), we assess the influence of the configuration feature vectors \({\varvec{p}}^B(r)\) within \(I^B(r)\) and \({\tilde{\varvec{p}}}^B(r)\) within \({{\tilde{I}}}^B(r)\) on the configuration feature vector \({\varvec{p}}^A(r)\) within \(I^A(r)\). We first conducted an empirical evaluation of the proposed influence analysis model [see Eqs. (7) and (8)] in terms of predictive performance. We divided \({\mathcal R}(n)\) into a training set \({{\mathcal {R}}}^{\text{train}}(n)\) and a test set \({{\mathcal {R}}}^{\text{test}}(n)\) along the time-axis, with a 7 : 3 ratio, and evaluated the predictive performance of the learned model on \({{\mathcal {R}}}^{\text{train}}(n)\) using a prediction log-likelihood ratio PLR defined by

(see Eq. 10). Here, PLR indicates the difference in log-likelihood between the learned model and a uniformly random model, and quantifies the relative performance of the learned model versus the random guessing on \({{\mathcal {R}}}^{\text{test}}(n)\). Note that the uniformly random model is obtained by setting \({\varvec{w}}= {\textbf { 0}}\), i.e., \(w_b = w_h = w_\ell = {{\tilde{w}}}_b = {{\tilde{w}}}_h= {{\tilde{w}}}_\ell = 0\) (see Eq. 7).

We compared the proposed generative model against four baseline models referred to as baselines 1, 2, 3, and 4, which are defined as follows: First, the baselines 1 and 2 are obtained by “\(w_b = w_h = w_\ell = 1\), \({{\tilde{w}}}_b = {{\tilde{w}}}_h= {{\tilde{w}}}_\ell = 0\)”, and “\({{\tilde{w}}}_b = {{\tilde{w}}}_h = {{\tilde{w}}}_\ell = 1\), \(w_b = w_h = w_\ell = 0\)”, respectively (see Eq. 7). Namely, the baseline 1 model supposes \({\varvec{p}}^A(r) = {\varvec{p}}^B(r)\), and the baseline 2 model supposes \({\varvec{p}}^A(r) = {\tilde{\varvec{p}}}^B(r)\). Next, the baselines 3 and 4 are obtained by “\(w_b = w_h = w_\ell = 1\), \({{\tilde{w}}}_b = {{\tilde{w}}}_h = {{\tilde{w}}}_\ell = -1\)”, and “\(w_b = w_h = w_\ell = -1\), \({{\tilde{w}}} = {{\tilde{w}}}_h = {{\tilde{w}}}_\ell = 1\)” (see Eq. 7). Namely, by considering active and inactive situations, the baseline 3 model assumes \({\varvec{p}}^A(r) \propto {\varvec{p}}^B(r) / {\tilde{\varvec{p}}}^B(r)\), and the baseline 4 model assumes \({\varvec{p}}^A(r) \propto {\tilde{\varvec{p}}}^B(r) /{\varvec{p}}^B(r)\), where the elementwise division of two vectors is used. In this paper, we only present the results of influence analysis for small n and \(\tau = 30\) days according to the analysis results of aggregated transition vectors (see Fig. 7).

Figure 9 indicates the evaluation results of the proposed, baseline 1 and baseline 2 models for the Dessert (2012), Meat-dish (2012) and Vegetable-dish (2012) datasets in terms of the PLR metric. Here, the baseline 3 and 4 models are excluded since they were on par with the uniformly random model. First, it is evident that the proposed, baseline 1 and baseline 2 models significantly outperformed the random guessing, while the value of PLR can vary depending on datasets. It should be noted that the baseline 2 model exhibited marginally better prediction performance than the baseline 1 model. However, the baseline 3 and 4 models derived from the simple ratios of \({\varvec{p}}^B(r)\) and \({\tilde{\varvec{p}}}^B(r)\) in view of active and inactive situations were quite ineffective. In contrast, the proposed model consistently outperformed the other models. These findings show the effectiveness of the proposed model, prompting us to utilize it for our subsequent influence analysis.

Results of influence analysis for the Dessert (2012), Meat-dish (2012) and Vegetable-dish (2012) datasets (\(n=4\) and \(\tau = 30\) days)

Full size image

Results for influence analysis

By applying the proposed influence analysis model [see Eqs. (7) and (8)], we examined the influence of \({\varvec{p}}^B(r)\) and \({\tilde{\varvec{p}}}^B(r)\) on \({\varvec{p}}^A(r)\) for any recipe n-simplex \(\sigma (r) = \{ v_0(r), v_1(r), \dots , v_n(r) \}\) (\(r \in {\mathcal R}(n)\)) in terms of the weight vector \({\varvec{w}}\). In this paper, we in particular describe the analysis results for \(n = 4\) and \(\tau = 30\) days.

Figure 10 indicates the results of \({\varvec{w}}\) for the Dessert (2012), Meat-dish (2012) and Vegetable-dish (2012) datasets. Regarding the five boundary-faces \(\sigma _j(r)\) (\(j = 0, 1, \dots , 4\)), the features \({{\tilde{p}}}^B_j(r)\) (\(j = 0, 1, \dots , 4\)) within long-term period \({{\tilde{I}}}^B(r)\) had a stronger influence on \({\varvec{p}}^A(r)\) than the features \(p^B_j(r)\) (\(j = 0, 1, \dots , 4\)) within short-term period \(I^B(r)\) (see the results of \(w_b\) and \({{\tilde{w}}}_b\)). For the ingredient relationships including \(\sigma (r)\) (called higher dimensional relationships), the influence of the feature \({{\tilde{p}}}^B_h(r)\) within \({{\tilde{I}}}^B(r)\) was stronger than or comparable to that of the feature \(p^B_h(r)\) within \(I^B(r)\) (see the results of \(w_h\) and \({{\tilde{w}}}_h\)). Moreover, we see that which of the boundary-faces and the higher dimensional relationships had a stronger influence on \({\varvec{p}}^A(r)\) depended on the datasets. On the other hand, as for the ingredient relationships lower than boundary-faces (called lower dimensional relationships), both the feature \({{\tilde{p}}}^B_\ell (r)\) within \({{\tilde{I}}}^B(r)\) and the feature \(p^B_\ell (r)\) within \(I^B(r)\) had a relatively weak influence on \({\varvec{p}}^A(r)\), compared to the other features (see the results of \(w_\ell\) and \({{\tilde{w}}}_\ell\)). In general, this can be attributed to the smallness of dimension n. It is worth noting that the observed properties can necessarily vary depending on the datasets to be considered. These findings suggest that the proposed influence analysis method has a potential for uncovering intriguing properties of ingredient co-occurrences in the recipe stream of social media

We addressed the problem of exploring configuration transitions associated with occurrences of recipe n-simplices in a temporal simplicial complex derived from a recipe stream in social media. By extending the previous work of Cencetti et al. (2021), we have proposed a novel framework of analyzing the transitions of boundary-based active configurations for the preceding and forthcoming investigation periods. First, we gave an analysis method of using the transition matrices, and by improving it, we further devised a method of employing the aggregated transition vectors so as to effectively analyze the whole picture across all dimensions n. Next, by introducing a probabilistic generative model for influence analysis, we provided a method of examining how the configuration feature vectors within the preceding periods affect the configuration feature vector within the forthcoming investigation period.

Using real data from a Japanese recipe sharing site, we empirically evaluated the effectiveness of the proposed analysis framework. We first showed that for the investigation periods of a reasonable length, the configuration transitions around recipe simplices of relatively small dimensions n are comprehensively characterized by their active boundary-faces. Next, compared with a baseline method straightforwardly derived from the previous work (Cencetti et al. 2021), we demonstrated the significance of incorporating the concept of active degree. Moreover, we revealed some interesting category-specific properties for the temporal evolution of Japanese homemade recipes published in social media from the perspective of ingredient co-occurrences.

In this paper, we focused on temporal simplicial complexes obtained from recipe streams. Clearly, the proposed framework for analyzing configuration transitions has the potential to be applied to temporal simplicial complexes for datasets from other domains. Our immediate future work is to evaluate it for other kinds of datasets. Inspired by homology theory in algebraic topology, we focused on the boundary-based active configurations before and after occurrences of simplices. Our future work includes examining changes in the topological structures of simplicial complexes in terms of homology groups for temporal simplicial complexes in various domains.

The Cookpad dataset we used in this paper was provided by Cookpad Inc. via IDR Dataset Service of National Institute of Informatics: https://www.nii.ac.jp/dsc/idr/cookpad/

1. https://cookpad.com/https://cookpad.com/

2. First, general-purpose ingredients typically used in Japanese cuisines such as soy sauce, salt, sugar, water and edible oil were excluded. Next, the ingredients that appeared in at least five recipes were identified as a set of major ingredients.

This work was supported in part by JSPS KAKENHI Grant Number JP21K12152.

Authors and Affiliations

1. Faculty of Advanced Science and Technology, Ryukoku University, Otsu, Japan

   Koudai Fujisawa, Masahito Kumano & Masahiro Kimura

Contributions

All authors designed research, contributed new analysis methods, analyzed data, and wrote the paper.

Corresponding author

Correspondence to Masahiro Kimura.

Competing interests

The authors declare that they have no competing interests.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions