Individual position diversity in dependence socioeconomic networks increases economic output

TL;DR

This study reveals that individuals' diversity in their positions within dependence networks from MMORPGs correlates positively with their economic output, highlighting the importance of network structure in socioeconomic systems.

Contribution

It introduces a novel analysis of dependence network motifs and classifies individuals based on structural similarity, linking position diversity to economic performance.

Findings

Node position diversity is positively correlated with economic output.

Leaf nodes have significantly lower economic outputs than other nodes.

Structural similarity of motif positions helps classify individuals into meaningful groups.

Abstract

The availability of big data recorded from massively multiplayer online role-playing games (MMORPGs) allows us to gain a deeper understanding of the potential connection between individuals' network positions and their economic outputs. We use a statistical filtering method to construct dependence networks from weighted friendship networks of individuals. We investigate the 30 distinct motif positions in the 13 directed triadic motifs which represent microscopic dependences among individuals. Based on the structural similarity of motif positions, we further classify individuals into different groups. The node position diversity of individuals is found to be positively correlated with their economic outputs. We also find that the economic outputs of leaf nodes are significantly lower than that of the other nodes in the same motif. Our findings shed light on understanding the influence of network structure on economic activities and outputs in socioeconomic system.

Equations

$${\alpha }_{i,m}=1-(k_i-1){\int }_0^{w_{i,m}}(1-x{)}^{k_i-2}\mathrm{d}x<\alpha ,$$

$$Z_{i,j}=\frac{p_{i,j}-⟨p_{i,j}⟩}{\sigma (p_{i,j})},$$

$$s_{i,m}=\frac{E[(Z_{i,j}-⟨Z_i⟩)(Z_{m,j}-⟨Z_m⟩)]}{E[(Z_{i,j}-⟨Z_i⟩{)}^2{]}^{1/2}E[(Z_{m,j}-⟨Z_m⟩{)}^2{]}^{1/2}},$$

$$d_{i,m}=1-s_{i,m}.$$

$$d_i=-\sum _{j=1}^{30}{p}_{i,j}\ln p_{i,j}.$$

Full text

††preprint: PRE

I Introduction

Considerable studies have offered us a deep understanding of the influence of network structures on the dynamics of complex systems, such as the spreading of diseases and information Wang et al. (2009) and emerging of collaborations Newman (2001). However, the connection between network position and economic output is less studied. It is reported that specific network structures may enhance economic outputs Bhattacharya and Dugar (2014); Xie et al. (2014); Eagle et al. (2010); Bettencourt et al. (2014); Ortman et al. (2015). Eagle et al argue that the diversity of individual relationships within a community strongly correlates with economic development of communities Eagle et al. (2010). Furthermore, Bettencourt et al and Ortman et al find that the diversity of relationships is positively correlated with the productivity of individuals and communities Bettencourt et al. (2014); Ortman et al. (2015).

Network motifs are building blocks of complex networks Milo et al. (2004a, 2002, b). It is found that network motifs of social networks may reflect the driving forces for forming social structures Kovanen et al. (2013); Klimek and Thurner (2013). Similar to the friendship networks of US students Ball and Newman (2013), Xie et al. studied triadic motifs in dependence networks of virtual societies and found that low level individuals have preference of forming links to high level individuals Xie et al. (2014). Their findings in virtual world are consistent with empirical findings in real society that “collaboration is easier when both partners share the same social status, and the probability of partnership formation decreases significantly as the status gap between the partners increases” Bhattacharya and Dugar (2014).

Understanding the structure and function of social networks are of great importance to investigate economic activities and outputs in social systems Kovanen et al. (2013); Palla et al. (2007); Onnela et al. (2007); Kumpula et al. (2007); Eagle et al. (2009); Jo et al. (2011); Jiang et al. (2013). In real world, social networks are not large and samples are biased Ball and Newman (2013); Currarini et al. (2009, 2010), which hinders the empirical investigation of relationship between network structures and economic activities in social system. Some other empirical studies have shown that individuals’ or firms’ network positions are closely related to success measured in terms of economic production Uzzi (1996); Guimera et al. (2014); Cantner and Joel (2011); Garas et al. (2014). In the era of big data, information technology provides us alternative methods to collect data of social relationships and economic activities, for example, massively multiplayer online role-playing games (MMORPGs) Jiang et al. (2009, 2010); Thurner et al. (2012); Szell et al. (2012); Szell and Thurner (2012), which enable us to study complex social and economic behaviours of human in online social systems Bainbridge (2007); Papagiannidis et al. (2008); Williams (2010). Bainbridge et al also emphasized the scientific potential of virtual world for future research Bainbridge (2007). Empirical studies show that social behaviours in MMORPGs are representative of human behaviours in many aspects in real society Chesney et al. (2009); Szell et al. (2010); Szell and Thurner (2010); Klimek and Thurner (2013); Szell and Thurner (2013); Grabowski and Kosiński (2008); Grabowski and Kruszewska (2007).

In this paper, we construct dependence networks based on weighted friendship networks of individuals and identify 30 distinct motif positions in 13 directed triadic motifs which represent local dependence among individuals. Using the $k$-means algorithm, we further classify individuals into $k$ clusters based on the motif position profiles. Our results indicate that the motif position of individuals do have great influence on their economic outputs.

II Materials and Methods

II.1 Data description

We use a huge database recorded from $124$ servers to investigate the potential connection between network structure and economic output for individuals within the virtual world of a popular Massively Multiplayer Online Role-Playing Game (MMORPG) in China. There are two opposing camps or societies in a virtual world residing in a server, thus giving us $248$ virtual societies. There exist great differences about the numbers of avatars among different virtual societies. The populations of virtual worlds vary from thousands up to fifty thousands. The distribution of the number of avatars in each virtual world is drawn in Fig. 1G. In each society, three professions have different skills. The advantage of warrior’s skill is the power of attack, the advantage of mage’s skill is the power of defense, and the advantage of priest’s skill is the ability to cure illness. An avatar can be a warrior, a priest, or a mage. To improve their skills, the avatars cooperate with friends to accomplish tasks. The more friends, the higher efficiency. Two individuals $i$ and $m$ are allowed to establish social ties to satisfy their desire of making friends and enhance their utility of collaborations. The strength of social ties is measured by the intimacy $I_{i,m}$, which increases according to the collaborative activities of $i$ and $m$ if they belong to the same society; otherwise, $I_{i,m}$ remains zero if $i$ and $m$ belong to two different societies. Hence the friendship networks of the two camps are essentially separated. As a measure of closeness to each friendship, the values of intimacy are recorded every day. When two individuals in the same society form a team and collaborate to accomplish a task, their intimacy increases. The evolving intimacy allows us to track the evolution of the cooperation behavior in the socioeconomic networks. Each individual can maintain a friendship list, denote as $\mathcal{F}_{i}$ for individual $i$. The social tie is symmetric: if $i\in\mathcal{F}_{m}$, then $m\in\mathcal{F}_{i}$. In addition to the friendship network, the game data contains other socioeconomic networks, such as the face to face trading networks between initiators and receivers, the vendor trading networks between vendors and costumers, the mail networks between senders and receivers, the mentor networks between students and mentors, the kill networks between killers and victims. In this paper, our focus is the friendship network. More details of the database can be found in our earlier works about the triadic motifs in dependence networks Xie et al. (2014) and skill complementarity in collaboration networks Xie et al. (2016).

II.2 Economic outputs of individuals

We measure the economic output by converting the virtual money and items into a standardized currency for each individual. There are two virtual currencies, Xingbi and Jinbi. The Xingbi and Jinbi can be exchanged in the built-in platform in each virtual world. The virtual system has an approximately stable exchange rate between Xingbi and the Chinese currency Renminbi. Jinbi is produced by the economic activities of the individuals when they form a team and collaborate to accomplish the tasks. The currency Jinbi and virtual items, such as weapons, clothes, and medicines, are awarded to individuals when monsters are killed and tasks are accomplished.

We convert the produced items and Jinbi to Xingbi to obtain the real economic output for each individual on each day. On average, the normal life span of virtual societies is close to 5 months Xie et al. (2016). Therefore, we calculate the output of each individual for a fixed period of 145 days for all virtual societies, denoted as $y_{i}=\ln\sum_{t=1}^{t=145}y_{i,t}$.

II.3 Construction of dependence network

For a given friendship network, we construct a dependence network by removing the insignificant edges based on statistical validation Serrano et al. (2009). First, we define the relative intimacy of individual $m$ in reference to all friends of individual $i$ as $w_{i,m}=I_{i,m}/\sum_{m=1}^{N}I_{i,m}$. Obviously, $w_{i,m}\neq w_{m,i}$. Following the statistical validation Serrano et al. (2009), a directed tie $i\rightarrow m$ is significant at the level of $\alpha$ if

[TABLE]

where $k_{i}$ is the degree of individual $i$. If $\alpha_{m,i}<\alpha$, the directed tie $m\rightarrow i$ is significant. For each society, we set the significant level $\alpha$ and remove the insignificant links, resulting a directed dependence network. If link $i\to m$ is significant, it means that individual $m$ is relatively important to individual $i$ in $i$’s friends. In other words, individual $i$ depends individual $m$. By the disparity filter Serrano et al. (2009), the dependence networks are the backbones of the original friendship networks.

Fig. 1 (A, B, C) show the topological structure of the dependence network constructed from a virtual society with significant level $\alpha=0.01,0.05,0.1$ respectively. We analysis the out-degree distribution $P(k_{\rm{out}})$ and in-degree distribution $P(k_{\rm{in}})$ of dependence networks for different significant level $\alpha=0.01,0.05,0.1$ in Fig. 1 (D, E). The degree distribution of the dependence network is obviously different from the original friendship network ($\alpha=1$). Using the smaller significance level $\alpha$, the disparity filter reduces more edges. For different significant level $\alpha$, the fraction of nodes, edges, weight kept in the dependence network increase with the significant level $\alpha$ as shown in Fig. 1 (F). And the average degree of the 248 dependence networks is monotone increasing with the significant level $\alpha$ in Fig. 1 (I). For different significant level $\alpha=0.01,0.05,0.1$, the distribution $P(N)$ of the 248 dependence networks’ sizes $N$ have nearly the same shape in Fig. 1 (G). The numbers of the avatars range from thousand to fifty thousands. Comparing with the global threshold filter, the disparity filter considers the relevant edges and ensures that the edges with small intimacy are not neglected. So we can find that the average intimacy kept in the dependence network is not monotone decreasing with the significant lever $\alpha$ in the Fig. 1 (H).

II.4 Quantifying position ratio profile

Following the method Milo et al. (2002) to identify motifs, we are able to identify 13 different directed triadic motifs, as shown in Fig. 2, in dependence networks. These motifs uncover the dependence structures among individuals at the microscopic level. For example, motif \scriptsize{2}⃝$\leftarrow$\scriptsize{1}⃝$\rightarrow$\scriptsize{2}⃝ stands for the situation that one individual depends on the other two individuals, motif \scriptsize{3}⃝$\rightarrow$\scriptsize{4}⃝$\rightarrow$\scriptsize{5}⃝ means that one individual depends on another individual which in turn depends on the third individual, whereas motif \scriptsize{a}⃝$\rightarrow$\scriptsize{b}⃝$\rightarrow$\scriptsize{c}⃝$\rightarrow$\scriptsize{a}⃝ represents the situation that individual $a$ depends on individual $b$, individual $b$ depends on individual $c$, and individual $c$ also depends on individual $a$. We can further locate 30 distinct motif positions Stouffer et al. (2012) within the 13 different directed triadic motifs, as shown in Fig. 2. We directly enumerate the relative frequency $p_{i,j}$ that individual $i$ appears in position $j$ across the 13 motifs, which gives the motif position ratio profile $p_{i}=(p_{i,1},p_{i,2},...,p_{i,30})$ for individual $i$. Note that we have $\sum_{j=1}^{30}p_{i,j}=1$.

Hence, we first define the $z$-sore of occurrence frequency for position $j$:

[TABLE]

where $\langle p_{i,j}\rangle=\sum_{i=1}^{N}p_{i,j}/N$ and $\sigma(p_{i,j})$ are the mean and standard deviation of $p_{i,j}$ over all individuals in the dependence network of a given virtual society Milo et al. (2004a, b). The structural similarity between individual $i$ and $m$ are thus defined as the correlation coefficients between $Z_{i}$ and $Z_{m}$, such that,

[TABLE]

where $E[x]$ is the mathematical expectation of $x$.

II.5 Classifying individuals based on their position ratio profile

We employ the $k$-means algorithm to classify $N$ individuals based on their position ratio profiles $\{Z_{i}\}$ in the dependence social networks. The position ratio profile of individual $i$ is wrote as $Z_{i}=(Z_{i,1},Z_{i,2},...,Z_{i,30})$. For each virtual society, we can get a matrix of position ratio profile, denoted by $\cal{Z}$, in which rows correspond to individuals and columns correspond to positions in triadic motifs. By adopting the $k$-means algorithm, we partition the individuals into $k$ clusters in terms of minimizing the sum, over all clusters, of the within-cluster sums of point-to-cluster-centroid distances. The partition algorithm is implemented in an iterative way.

To evaluate the optimal number of clusters, we calculate the clustering evaluation object containing Davies-Bouldin index values and find the minimum criterion value. Then we can get the optimal number of clusters $k$. The correlation distance defined as one minus the sample correlation between points (treated as $Z_{i}$) is used to capture the closeness between individuals in $k$-means algorithm, such as,

[TABLE]

The $k$-means algorithm returns the cluster indices of each individual and the $k$ cluster centroid locations in a $k$-by-$30$ matrix. Each centroid is the component-wise mean of the points in that cluster. We denote the $k$-th cluster centroid locations as $P_{k}=(P_{k,1},P_{k,2},...,P_{k,30})$, where $P_{k,j}=\sum_{i\in{\cal{C}}_{k}}p_{i,j}/c_{k}$ and $j\in\{1,2,...,30\}$. $c_{k}$ is the number of individuals in cluster ${\cal{C}}_{k}$. We have $\sum_{j=1}^{30}P_{k,j}=1$, and $P_{k}$ is the position ratio profile of cluster $k$.

III Results

III.1 Quantifying position ratio profile

Individuals inside the virtual societies are involved in all kinds of social and economic activities. The theory of structural holes already tells us that the agent occupies special position in social network will have great influence on his economic performance Burt (2009). Here, we empirically investigate the underlying connections between network structure and economic output for individuals in 248 virtual societies. Our data are game logs from a popular MMORPG in China (see data description in section II). For each society, we construct a dependence network of individuals, which only keeps the multiscale backbone of original friendship network (see construction of dependence network in section II). Fig. 1 illustrates the dependence network from a virtual society and some basic statistics about the topological structure. One can see that there are lots of in-stars in dependence network, which is in accordance with the findings Xie et al. (2014) that the in-degree distribution is much fatter than out-degree distribution. Such phenomena can be explained by that in virtual societies the individuals with high levels may play a relatively important role in friend lists of low level individuals, which is supported by the strong preference of connecting high level individuals for low level individuals Xie et al. (2014).

For a given dependence network, we can estimate a position ratio profile for each individual (see quantifying position ratio profile in section II). Xie et al. revealed that the open motifs have higher occurrence frequency than close motifs Xie et al. (2014), suggesting that the positions in the open motifs may have higher occurrence frequency than positions in close motifs. Fig. 3 (A) shows the position ratio profile for all individuals in dependence networks in Fig. 1. There are more zeros in the position ratio profiles of individuals with low degrees than those with high degrees. More extremely, there are lots of individuals occupying only one specific position ($j$ = 3, 9, 10, 16, 17, and 18) of triadic motifs in dependence networks, which leads to $p_{i,j}=1$. The individuals with lager value of $p_{i,3},p_{i,9},p_{i,16}$ in their position ratio profile may correspond to leaf nodes in dependence networks. Note that the occurrence frequency of motif \scriptsize{a}⃝$\rightarrow$\scriptsize{b}⃝$\rightarrow$\scriptsize{c}⃝$\rightarrow$\scriptsize{a}⃝ is 0, resulting in that the position 21 is impossible to be observed.

III.2 Similarity of position ratio profiles between individuals

The triadic motifs are regarded as building blocks of complex networks Milo et al. (2002), suggesting that the 30 positions in triadic motifs may contain important structural information of nodes in complex networks. Here we utilize the occurrence frequency of the 30 different positions in dependence networks to represent the network structure profile for each individual. Based on these profiles, we can assess the structural similarity between individuals, which allows us further to group individuals. In dependence networks, it is observed that some motifs, for example, \scriptsize{3}⃝$\rightarrow$\scriptsize{4}⃝$\rightarrow$\scriptsize{5}⃝, \scriptsize{9}⃝$\rightarrow$\scriptsize{10}⃝$\leftarrow$\scriptsize{9}⃝, \scriptsize{16}⃝$\rightarrow$\scriptsize{17}⃝$\rightleftharpoons$\scriptsize{18}⃝, appears more frequency than other motifs Xie et al. (2014).

By ordering the individuals according to the rule that nodes with large similarity are close to each other, we plot the structural similarity in Fig. 3 (D). The color bar stands for the value of similarities $s_{i,m}$. One intriguing observation is that there is a block-diagonal structure, strongly indicating the function of grouping individuals for position ratio profiles. This inspires us to further classify the individuals into clusters by $k$-means algorithm (see classifying individuals based on their position ratio profile in section II). In Fig. 3(A) and Fig. 3(D), there are some classes which are made up of some special individuals, such as individuals with $p_{i,3}=1$ in motif \scriptsize{3}⃝$\rightarrow$\scriptsize{4}⃝$\rightarrow$\scriptsize{5}⃝, or individuals with $p_{i,9}=1$ or $p_{i,10}=1$ in motif \scriptsize{9}⃝$\rightarrow$\scriptsize{10}⃝$\leftarrow$\scriptsize{9}⃝, or individuals with $p_{i,16}=1$, $p_{i,17}=1$, or $p_{i,18}=1$ in motif \scriptsize{16}⃝$\rightarrow$\scriptsize{17}⃝$\rightleftharpoons$\scriptsize{18}⃝.

The $k$-means algorithm is an iterative partitioning algorithm, which maximizes the similarity of the within-cluster. To evaluate the optimal number of clusters, we create a clustering evaluation object containing Davies-Bouldin index values in Fig. 3 (E). From the minimum criterion value, we can get the optimal number of clusters $k=41$. Based on the definition of position ratio profiles, there is a correlation among the 30 positions. From the dendrogram of the 30 position in Fig. 3 (B), we can find that a avatar at position 22 has a higher chance to be also in position 23 or 24 in the same motif. But there is no obvious classification among the 30 position. The block-diagonal structure can not be explained by this relation among the 30 position. An alternative method to elucidate the position significant profiles is to calculate the difference of position ratio profiles between the dependence networks and the reference randomized dependence networks Milo et al. (2002). But this method is too computationally intensive for our 248 dependence networks and some of the networks’ sizes are larger than fifty thousands.

III.3 Relationship between position ratio profiles and economic outputs

In social network, the individuals in key positions may have great impacts of network dynamics if some controlling strategies are applied on these key nodes. This leads to the conjecture that the position ratio profiles may have potential influence on the economic outputs for individuals. We estimate the economic output $y_{i}$ for each individual in a given world (see economic outputs of individuals in section II). To make these economic outputs comparable between different virtual worlds, we standardize the economic output within each society. The standardization of economic performance of individual $i$ is defined as output per capita, $Y_{i}=[y_{i}-\langle y_{i}\rangle]/\sigma(y_{i})$, where $\langle y_{i}\rangle$ and $\sigma(y_{i})$ are the mean and standard deviation of $y_{i}$ over all individuals in a given society. Therefore, the mean of economic performance of the class ${\cal{C}}_{k}$ can be estimated via $Y_{k}=\sum_{i\in{\cal{C}}_{k}}Y_{i}/c_{k,}$, where $c_{k}$ is the number of individuals in class ${\cal{C}}_{k}$.

Motivated by social diversity Eagle et al. (2010), we define the individual position diversity $d_{i}$ as a function of Shannon entropy to quantify how individual appears in the 30 positions:

[TABLE]

Fig. 3 shows the results from the analysis of position ratio profiles for individuals in Fig. 1 with significant level $\alpha=0.05$. The position ratio profiles of 4331 individuals are presented in Fig. 3 (A). To improve the visualization, the color bar represents the log value of position ratio profiles $p_{i,j}$. The position diversity $d_{i}=0$ means that the individual only occupy one triadic motif position. As shown in Fig. 3 (A), there are individuals with $p_{i,3}=1$ in motif \scriptsize{3}⃝$\rightarrow$\scriptsize{4}⃝$\rightarrow$\scriptsize{5}⃝, $p_{i,9}=1$, $p_{i,10}=1$ in motif \scriptsize{9}⃝$\rightarrow$\scriptsize{10}⃝$\leftarrow$\scriptsize{9}⃝, $p_{i,16}=1$, $p_{i,17}=1$, and $p_{i,18}=1$ in motif \scriptsize{16}⃝$\rightarrow$\scriptsize{17}⃝$\rightleftharpoons$\scriptsize{18}⃝. Fig. 3 (C, F) show the economic output $y_{i}$ and the individual position diversity $d_{i}$ respectively. The individuals with position diversity $d_{i}=0$ account for a large part of avatars in the dependence network of a given society. Fig. 3 (L) shows the scatter plot of economic output $y_{i}$ and position diversity $d_{i}$. The correlation coefficient $\rho$ is 0.63 and $p$-value is less than $10^{-6}$, implying that the correlation between $y_{i}$ and $d_{i}$ is positive and highly significant. Our results indicate that the individuals who appear in more triadic motif positions have higher economic outputs. It is easy to find that more active players have more friends (more in-degree and out-degree) and appear in more triadic motif positions. Here we show the plots of the relation between the economic output $y_{i}$ and the individuals’ in-degree $k_{\rm{in}}$, out-degree $k_{\rm{out}}$ in Fig. 3 (G, H) respectively. This agrees with Fuchs’ result Fucha and Thurner (2014) that the output is correlated to both in- and out-degree. To avoid the impact of the in- and out-degree, we investigate the relation between the economic output $y_{i}$ and the position diversity $d_{i}$ of individuals with fixed $k_{\rm{in}},k_{\rm{out}}=1,3,5$ in Fig. 3 (J, K, L). And we can get the same conclusion that the individual position diversity increases economic output of individuals with fixed $k_{\rm{in}}$ or $k_{\rm{out}}$. What is more, we also find that the economic output is susceptible to asymmetries between individuals’ in- and out-degree in our dependence socioeconomic network. The individuals with high in-degree have higher economic outputs than with high out-degree. Because it is difficult to get the records and measure the activity of individuals accurately, we assume that the individuals with the same in- and out-degree have almost the same activity.

The $k$-means algorithm also gives cluster centroid locations denoted as $P_{k}=(P_{k,1},P_{k,2},...,P_{k,30})$ for each class ${\cal{C}}_{k}$. Similar to Eq. 5, we define the cluster position diversity as $D_{k}=-\sum_{j=1}^{30}P_{k,j}\ln P_{k,j}$. In Fig. 4(A-F), the clusters are sorted by the cluster position diversity $D_{k}$ in ascending order. Considering the impact of different significant level $\alpha$ on the relation between the economic output and position diversity, we analysis the dependence network from all the virtual society with different significant level $\alpha=0.01,0.05,0.1$. In Fig. 4, the three columns from left to right correspond to the dependence networks with significant level $\alpha=0.01,0.05,0.1$ respectively. Fig. 4 (A, B, C) show the plots of the average economic output $Y$ of classes and Fig. 4 (D, E, F) illustrate the plots of the value of $P_{k}$ for different significant level $\alpha=0.01,0.05,0.1$ respectively. The clusters are sorted by the cluster position diversity $D_{k}$ in ascending order, so one can see that the economic output $Y$ increases with the cluster position diversity $D_{k}$. We can find that lots of clusters are comprised of individuals with $P_{k,3},P_{k,9},P_{k,16}$ equalling to one. Although these individuals have the same individual position diversity $d_{i}=0$, their economic outputs could be different. For each dependence network with a given $\alpha$, we calculate the Davies-Bouldin index values to evaluate the optimal number of clusters and get the optimal number of clusters $k$. The distributions of the optimal number $k$ are shown in Fig. 4 (G, H, I). The average optimal number $k$ is between 50 and 60.

To compare the relation coefficients between the economic output and the individuals’ in-degree $k_{\rm{in}}$, out-degree $k_{\rm{out}}$, position diversity $d_{i}$, we denote the coefficients as $\rho_{k_{\rm{in}}}$, $\rho_{k_{\rm{out}}}$ and $\rho_{d_{i}}$ respectively and draw the distribution of the three kinds of coefficients for different significant level $\alpha$. When the significant levels $\alpha$ is smaller, there are greater difference between the three kinds of coefficients. At the same time, the asymmetries between individuals’ in- and out-degree in our dependence socioeconomic network is obvious with the significant level $\alpha=0.01$.

III.4 Economic output of individuals with the position diversity $d_{i}=0$

Here, we conduct a comparison of the economic output of individuals with individual position diversity $d_{i}=0$ in the dependence networks from all the virtual society with significant level $\alpha=0.01,0.05,0.1$. The individual position diversity $d_{i}=0$ means that the individual appears only in one triadic motif position.

We perform $t$-tests on pairs of economic output of the individuals at two position in the same motif. Fig. 5 shows the box plots of economic output of individuals with position diversity $d_{i}=0$ for 30 unique positions $j\in\{1,2,...,30\}$ in 13 directed triadic motifs.

In Fig. 5, the superscript *, **, *** stand that the relationship is acceptable at the significant level $0.05$, $0.001$ and $0.0001$ respectively. For all the dependence networks with significant level $\alpha=0.01,0.05,0.1$, the results of the $t$-tests are similar.

In Fig. 5, we perform $t$-tests on pairs of individual economic output of the individuals at two position in the dependence networks. It is worthy stressing that the mean economic output of individual at position $1$ is greater than that of individuals within position $2$ in motif \scriptsize{2}⃝$\leftarrow$\scriptsize{1}⃝$\rightarrow$\scriptsize{2}⃝. This relationship can be expressed as $1>2$ in Fig. 5. The mean economic output of individuals at position $4$ and $5$ are greater than that of individuals at position $3$ in motif \scriptsize{3}⃝$\rightarrow$\scriptsize{4}⃝$\rightarrow$\scriptsize{5}⃝. These relationship can be expressed as $4>3$ and $5>3$ in Fig. 5. The mean economic output of individuals at position $10$ is greater than that of individuals at position $9$ in motif \scriptsize{9}⃝$\rightarrow$\scriptsize{10}⃝$\leftarrow$\scriptsize{9}⃝. This relationship can be expressed as $10>9$ in Fig. 5. The mean economic output of individuals at position $17$ and $18$ are greater than that of individuals at position $16$ in motif \scriptsize{16}⃝$\rightarrow$\scriptsize{17}⃝$\rightleftharpoons$\scriptsize{18}⃝. These relationship can be expressed as $17>16$ and $18>16$ in Fig. 5.

Comparing the output for individuals at different positions in the dependence networks with significant level $\alpha=0.05$, the 30 positions can be classified into two groups. The first group contains the positions: \scriptsize{3}⃝, \scriptsize{16}⃝, \scriptsize{9}⃝, \scriptsize{10}⃝, \scriptsize{18}⃝, \scriptsize{4}⃝, \scriptsize{17}⃝, \scriptsize{5}⃝, \scriptsize{7}⃝, \scriptsize{20}⃝, \scriptsize{8}⃝, \scriptsize{13}⃝, \scriptsize{2}⃝, \scriptsize{26}⃝ . The second group contains positions: \scriptsize{25}⃝, \scriptsize{23}⃝, \scriptsize{12}⃝, \scriptsize{24}⃝, \scriptsize{11}⃝, \scriptsize{19}⃝, \scriptsize{6}⃝, \scriptsize{15}⃝, \scriptsize{29}⃝, \scriptsize{22}⃝, \scriptsize{1}⃝, \scriptsize{27}⃝, \scriptsize{30}⃝, \scriptsize{14}⃝, \scriptsize{28}⃝. All the positions are sorted by means of individual economic output in ascending order. The mean economic output of individuals in the first group is less than 0, while in the second group the mean output is greater than 0. Such differences could result from that most of positions of the first group are in open motifs while most of positions in second group belong to close motifs.

IV Discussion

In this paper we analyzed the relationship between friendship structure and economic output for individuals in a massively multiplayer online role-playing game. We found that the individual position diversity of individuals is positively correlated with their economic output and social status. We have developed a new approach to study structural similarity of individuals in networks and classify individuals into different clusters based on their position ratio profiles. It is found that the economic output of leaf nodes is significantly lower than the nodes at the other triadic motif positions. The individual position diversity positively correlates with the economic output. Our results are in consistent with the results from the physical society Bhattacharya and Dugar (2014); Ball and Newman (2013). Furthermore, Eagle et al. found that the diversity of individual relationships within a community strongly correlates with the economic development of the community Eagle et al. (2010). In a recent analysis, Xie et al. found the degree of skill complementarity is positively correlated with their output Xie et al. (2016).

Much evidence shows that the behaviour of individuals in virtual society is representative in many aspects of human behaviour in physical society. This is rational because individuals’ decision process is determined by players that are real human being. The invisible hand operates not only in modern societies but also in ancient societies, not only in real societies but also in virtual societies. Yee et al. argued that online environments such as MMORPGs could potentially be unique research platforms for the social sciences and clinical therapy, but it is crucial to firstly establish that social behavior and norms in virtual environments are comparable to those in the physical world Yee et al. (2007). To investigate the relation between friendship (or socioeconomic) networks in virtual and real worlds, Grabowski and Kruszewska conducted a survey among the players of an online game and construct the off-line network Grabowski and Kruszewska (2007). They showed that the structure of the friendship network in virtual world was very similar to the structure of different social networks in real world.

We believe that our interdisciplinary work represents significant scientific evidence for understanding the behaviour of social systems. It involves topics that range from social science, network science, and economics to human dynamics. It also enriches our understanding on the formation of socio-economic networks and proposes a new method to classify nodes of complex networks for understanding of people’s economic behaviour from the big data of massive players. Our work sheds new light on the scientific research utility of virtual worlds for studying human behaviour in complex socio-economic systems.

Acknowledgements.

This work was supported by the National Natural Science Foundation of China [11505063, 11375064 and 11605062]; the Shanghai Chenguang Program [15CG29]; the Ph.D. Programs Foundation of Ministry of Education of China [20120074120028]; and the Fundamental Research Funds for the Central Universities.