in the Agent-Artifact Space: Emergence of Hierarchical Order in the 2000–2014 Photonics Techno-Economic Complex System

. In this work we represent a techno-economic complex system based on the agent-artifact space theoretical model. The objective is to structure a methodology to statistically investigate the presence of hierarchical order, as an emerging property of this system. To analyse the agent-artifact space, two statistical methodologies are initially employed. The ﬁrst is a community detection method, employed with the objective to detect groups of agents that are likely to intensively exchange information within the considered complex system. The second is a natural language processing method, the LDA topic model, employed with the objective of identifying types of artifacts as technological subdomains through textual information that describes the activities of agents. After this initial part, we address the investigation of the structure of the agent-artifact space by estimating the involvement of each community in the detected topics. This is eﬀectuated by means of a statistic that considers the information ﬂow percentage of agents, the fractional count of activities, and the probability of agents’ activities to belong to topics. We then estimate the hierarchical order of the topics’ distribution in communities, by computing its nestedness temperature, which is adopted by studies on ecological systems. This statistic’s signiﬁcance is ﬁnally evaluated with z -scores based on homogeneous systems. The case study is a system consisted of economic agents (e.g. ﬁrms, universities, governmental institutions) patenting in the technological domain of photonics. The analysis is eﬀectuated over ﬁve time


Introduction
The present work aims to investigate the emergent presence of physical order in an observed techno-economic complex system. Initially, two parallel methodologies are implemented with the respective objectives of (i) detecting group of agents that are likely to interact among themselves, and (ii) identifying the technological thematic subdomains that emerge under the initially considered technology, i.e. 'photonics'. As described in the relevant literature regarding the conceptualization of the agent-artifact space [19][20][21], two elements have to be considered when analyzing evolutionary dynamics regarding innovation and technological developments. The first element is the presence of interactive mesostructures of agents, which in this work are investigated as communities of agents (i.e. groups) intensively interacting among themselves. The second element is the presence of different kinds of artifacts belonging to the same technological domain. These are investigated in this study as thematic topics (i.e. technological subdomains). From a methodological point of view, the first part of the analysis is developed based on a community detection over a multilayer network (MLN), while the second analysis is developed based on an unsupervised natural language processing method, namely topic modeling.
After the detection of communities and topics, the distribution of topics in communities is analyzed with an approach based on ecology [3,27]. In order to look for the existence of any hierarchical structure, we investigate nestedness. In ecological systems, one of the interpretations of the nestedness's presence is that: (i) species that are located in the largest number of different habitats, are the only ones located in the habitat presenting the minimum diversity, and (ii) species that are located in few habitats, are located in habitats that are populated by a large diversity of species. When either of these cases occurs, the system is hierarchically ordered. In order to measure this order, the nestedness is measured as the temperature of the matrix describing the presence of species in habitats [3], where a lower temperature indicates a more nested system. Finally, to assess the presence of order in the considered agent-artifact space, i.e. a complex system of economic institutions patenting in photonics from 2000 to 2014, we refer to the previously introduced concepts. In the matrix describing the involvement of communities of agents in the technological subdomains, we compute the nestedness temperature. Then its statistical significance is evaluated with computed homogeneous systems. The obtained results prove that the detected levels of nestedness are statistically significant in the considered case study.
The originality of this work is the development of a method to analyze technoeconomic complex systems towards two directions. The first is the use of the conceptualization of the agent-artifact space in order to structure the analysis of a techno-economic complex system. The second is the investigation of its emerging properties and, more specifically for this work, the presence of physical order, which is here measured in terms of nestedness. Therefore, we direct the investigation of agents' interactive groups (the communities) and of the existing types of artifacts (the topics), towards the measurement of their intertwining.
The economic implications of the implemented analysis and a prediction model based on the obtained results, are not developed in this work. Both these will be the subject of subsequent studies.
In Sect. 2, we introduce the methodology developed to detect communities of agents and technological topics. In Sect. 3 we refer to nestedness temperature to assess the involvement of the communities in the topics, and we measure its significance. Finally, in Sect. 4 we implement the methodology in a real observed system and we discuss the obtained results.

A Statistical Approach for the Representation of the Agent-Artifact Space
The informational basis needed for the implementation of the proposed methodology is a bipartite network made of economic agents (e.g. companies, universities, research centers, local government institutions) participating in economic activities related to innovation processes (in this work, for instance, we consider patents). Apart from the aforementioned network structure, the requested information is dual: the geographic location of the agents, and the textual information describing the economic activities in which they are involved. On this basis, two methodologies are implemented in parallel [25]. These are a (i) community detection based on a MLN structure, and (ii) a topic modelling based on the textual technological information.
The first, i.e. the community detection, aims to detect groups of agents that are likely to exchange information, due to proximity determined by interactions in multiple dimensions. The second, i.e. the topic modelling, aims to identify technological subdomains, in order to categorize the innovative economic activities based on their content. These two analyses are selected according to the relevant literature regarding the conceptualization of the agent-artifact space [19][20][21], in which innovation processes are explained to be the results of a series of interactions that (i) occur among the agents involved in the system, and (ii) whose core objective is the discussion about the content of a specific technological artifact. In order, to pertinently develop the methodology with respect to the considered theoretical approach, we investigate how agents interacted, and about what they interacted. In the rest of the section, we present the two parallel analyses that constitute the first part of our work.

A MLN Representing a Multi-dimensional Space of Interactions.
The initial bipartite network of agents and activities, conceptually agents and artifacts, is transformed in its one-mode projection. The result is a one-mode network in which agents are connected among themselves based on the common activities they performed together. Then, two additional networks are generated. In the first, the same agents of the aforementioned network are connected only when belonging to the same geographical area, which in this work is described as sub-regions. In the second, the same agents are connected on the basis of the use of similar keywords in the patents they developed. These three networks, which are formed by the same agents, are subsequently conceived as layers, and then a single MLN is generated. The MLN can be formally described as a graph G = (V, E, J), where V = {1, 2, ..., x, ..., N } is the set of N agents (i.e. the nodes of the graph), E = {1, 2, ..., m, ..., M } is the set of M edges, and J = {j 1 , j 2 , j 3 } is the set of the three considered layers. As any edge is located in one and only one layer, mutually exclusive subsets of E can be created so as to separate the edges depending on the layer to which they belong. Formally, we have that where Γ is the function that univocally assigns any edge m to its corresponding layer j.
The scope of this MLN is to represent three dimensions of interaction in which the agents are involved. These three dimensions acknowledge socio-economic and biological complex system theories regarding the formation of bottom-up mesostructures [5,20,21,23,25]. In particular, the three aspects that these theories address as fundamental are: processes, structures and functions. With respect to our work, the processes are determined by observed agents' co-participations in patents, which identify occurred exchanges of information. The structures are determined by the agents' location, which can determine concentration of specific technologies and related know-how, according to theories regarding economic districts [4,9,26]. The functions are determined by the agents' technological orientations, which can reveal agents' semantic proximities and convergence towards similar economic processes [11,13,16]. The objective of this MLN is to represent agents in a multi-dimensional space and, more specifically, to infer their potential interaction intensity based on their proximity.

Community Detection of Interactive Economic Agents.
In order to identify groups of multi-dimensional interactive agents, a community detection information theory based algorithm, namely Infomap [18,24,28], is implemented on the MLN. The objective is the detection of groups of nodes, i.e. communities of agents, that are likely to exchange information. The community detection algorithm is selected not only for its implementability in MLNs [12], but also for its core and basic functioning. The Infomap algorithm performs the detection of subsets of agents by simulating a spread of a flow throughout the MLN, and minimizing the information needed to describe the circulation of the flow. This allows the identification of meso-structures within which intense information exchanges occur. As information management is an essential organizing principle in the initial formation and in the life-cycle of economic and biological systems [2,6,10,14,15,17,29], the groups detected by Infomap are likely to activate and/or join adaptive bottom-up dynamics.
In order to balance the weight of the three layers, the weights of the connections of each layer are scaled based on the following a ratio. This ratio is computed between the sum of the weights of the connections that belong to the considered layer, and the sum of the connections that belong to a selected layer that is used as reference. For any edge m ∈ E, a weight ϕ(m) that defines the intensity of the connection, is considered. The sum of the weight of the edges of any layer can be computed as ϕ(E j ) = m∈Ej ϕ(m). In addition, with ϕ * = max j∈J ϕ(E j ) we refer to the weight of the 'heaviest' layer. Then, in order to normalize the layers' weight, we compute a new edges' weight, namelyφ(m), as followsφ where m ∈ E j . As a result, we have thatφ(E j1 ) =φ(E j2 ) =φ(E j3 ), wereφ(E j ) is defined specularly to ϕ(E j ). Qualitatively speaking, the sum of the edges' new weight of any layer is equal to the sum of the new edges' weight of any other layer. ϕ(m) is the set of intra-layer edges' weight that is considered when implementing the community detection algorithm 1 . Thanks to this normalisation, the flow that is simulated to circulate throughout the entire MLN structure has the same probability to move to any layer. The result of this part of analysis is that each agent x is assigned univocally to a community c, as the algorithm is set to search for a hard-partition of G (i.e. nonoverlapping subsets of V ). In addition, Infomap provides also the information regarding the percentage of simulated flow crossing each agent, namely q.

Topic Modeling of Technological Artifacts.
In parallel, a second methodology is developed to investigate the agent-artifact space from the artifacts' perspective. In particular, in order to explore their semantic content, an unsupervised learning algorithm is implemented. The unsupervised generative model that is used is the Latent Dirichlet Allocation (LDA) [8]. LDA classifies the collected discrete textual information to a finite number of thematic topics, with the conjecture that they represent combinations of technological subdomains. In this scope, topic per document and words per topic models are established to obtain the most probable thematic groups, or else topics, with Dirichlet multinomial distributions. The assumptions that are made are the following: -textual information, stored as documents, is a mixture of one or multiple topics simultaneously (as in natural language words may belong to multiple topics), which generate relevant words based on probability distributions -a topic is a mixture of words from several documents [7,22,30] -words' order is not considered (exchangeability and bag-of-word assumptions) [1,8,30].
This analysis provides a set of θ k (y), each of which represent the probability that the activity y belongs to topic k, using the textual information provided by the activity. For each identified topic, the corresponding θ k is computed.

Analysis of Hierarchical Order in the Distribution of Communities and Topics
The second part of the outlined methodology is made of three stages, which are described in the paragraphs contained in this section.

The Involvement of Communities in Topics:
w c,k . In order to assess the involvement of interactive communities of agents in activities' technological subdomains, we compute the following statistic: where x is an agent belonging to community c, Y x is the set of activities in which agent x is involved, y is an activity for which y ∈ Y x , ξ is the fractional counting of the activities, q is the Infomap flow associated with the agents by the analysis of the MLN, and θ k is the probability that an activity belongs to topic k, as computed by the topic modelling. The fraction counting, i.e. ξ, is a function that computes the reciprocal of the number of agents involved in that activity 2 . This function is used (i) to equally distribute the weight of the activity among all the agents involved and, (ii) to ensure that all the activities have the same weight in the system. In fact, if this was not implemented, the entire weight of an activity would have been repeated for all the agents involved in it. Finally, the values of w c,k are linearly scaled in the interval [0, 1]. In this way, the values w c,k are obtained.

The Binary Matrix B h and the Nestedness Temperature T .
The W matrix is determined, with c indicating the rows, k the columns, andw c,k the value of the cells. Then, the presence of topics in communities, as described by matrix W , is computed in a binary way using a threshold. If the presence of a topic in a community, i.e.w c,k , is above a threshold h, where h ∈ R and 0 ≤ h ≤ 1, then the topic is considered to belong to this community. Communities with all topics below the threshold are discarded, hence not considered to be part of the matrix. This allows to focus on the communities with a minimum strength in terms of involvement in a topic. Based on this step, the binary matrix B h is generated from the matrix W , depending on the selected value of h. Formally, the value of the elements of the matrix B h is 1 if the corresponding w c,k ≥ h, otherwise is 0. Subsequently, the nestedness temperature T of the matrix B h is computed, as defined by Atmar and Patterson [3]. Rows and columns of the matrices are sorted in decreasing order by the row-sums and the column-sums, respectively. T which is within the range of [0, 100] degrees, measures the unexpectedness of the nonempty cells that are detected below the anti-diagonal, and the unexpectedness of the empty cells that are detected above the anti-diagonal. These cells are considered as unexpected, as they contribute to increase the disorder in the matrix 3 . The larger the disorder, the higher the nestedness temperature T .
The Statistical Significance of T . In order to assess the statistical significance of the computed nestedness temperature T , homogeneous systems are used. Starting from the binary matrix B h , 1,000 homogeneous matrices B h are computed. In order to be homogeneous with respect to B h , each B h presents the following characteristic: the number of communities to which each specific topic belongs is the same as observed in B h . This means that B h are matrices that are randomly generated, with the only constraint that the column-sums are equal to the column-sums of B h . The matrices B h represent homogeneous systems because (i) each topic belongs to the same number of communities as in the original system (B h ), and (ii) the communities to which the topics belong are randomly determined. This stage allows the creation of a sample of matrices to be used to investigate the significance of the nestedness temperature of the original system. The distribution of the values of the nestedness temperature T of the matrices B h is used to compute the z-scores of the nestedness temperature T of the matrix B h . This statistic, namely z T (B h ), is calculated as follows: where T (B h ) is the average temperature of B h matrices, and σ(T (B h )) is the standard deviation of the temperature of B h matrices.

The Analysis over the Complex System of Photonics Patents in 2000-2014
Based on a tech-mining approach, a set of 4,926 patents in the field of 'photonics' of the period 2000-2014, are identified. This allows the initial definition of our system, in which the agents are the 1,313 economic institutions (e.g. firms, research institutes and governmental institutions) that were filing at least one the detected patents. As the considered time period is sufficiently large to be cut in smaller time spans, five distinct systems are defined, each of them referring to a different three-year period. According to Eq. 2, the computation of the valuesw c,k is performed for each considered system, so as to obtain the corresponding W matrices. Then, in order to generate the binary matrices B h , five possible values of the threshold h are considered, namely 0.01, 0.02, 0.05, 0.1 and 0.2. For each of them, nestedness temperature T is calculated. In Fig. 1, the five matrices B h , with h equal to 0.05, are represented with empty cells in white (so as to represent a community non intensively involved in the corresponding topic), and non-empty cells in red (so as to represent a community intensively involved in the corresponding topic). In the same figure, in each matrix B h also the corresponding value of T is reported.
Finally, in order to assess the statistical significance of the obtained T , zscores are calculated as described in Eq. 3, based on 1,000 B h for each B h . The z T (B h ) that are obtained for the five systems and the five currently considered thresholds, which are presented in Table 1, are all lower than −2, except from two cases. These are for the systems referred to period '2000-2002' and '2012-2014', when h = 0.2 (the corresponding z T (B h ) equals 1.097 and −1.367, respectively). The T (B h ) are proven to be significantly far from the mean values of T (B h ) (at least two standard deviations), which means that the nestedness temperature is statistically significantly low. Therefore, the distribution of topics and communities expresses a significantly high level of hierarchical order, independently from the adopted threshold's value (with only two non significant results both associated to h = 0.2, i.e. the largest threshold considered).

Conclusions
For the developed work, an emerging property of the considered complex system is detected and its statistical significance is confirmed. The analyses reveal that (i) the topics with the minimum diffusion are associated to the communities with the largest set of topics, and (ii) the topics with the maximum diffusion are the only ones to be included in mono-topic communities. Given the obtained results, the methodological approach here outlined to model the agent-artifact space is proven to be pertinent to detect an emergent property of it, at least for the considered case study. More specifically, this work provides statistical proofs to support the presence of hierarchical order in the distribution of technological subdomains (i.e. the topics) representing distinct types of artifacts, over interactive groups of agents (i.e. the communities) representing areas of the system hosting intense exchanges of information.
The developed analysis has to be discussed in the context of economic theories, regarding the distribution of different technological topics over different communities belonging to the same agent-artifact space. Further analyses will consider the dynamic of the system, as the communities detected in each period do not necessarily have continuity with those detected in the following period. In fact, the presented work addresses the static analysis of five instances of the same system. In addition, the 'unexpected' combinations observed in matrices B h is inviting to investigate the innovative dynamics. In this perspective, the topological proximities of communities, as observed in the MLN, and the thematic proximities of topics will be considered in a subsequent analysis. Finally, other case studies will be evaluated following the same methodological approach outlined in this work.