# Communication activity in a social network: relation between long-term correlations and inter-event clustering

## Abstract

The timing patterns of human communication in social networks is not random. On the contrary, communication is dominated by emergent statistical laws such as non-trivial correlations and clustering. Recently, we found long-term correlations in the user’s activity in social communities. Here, we extend this work to study collective behavior of the whole community. The goal is to understand the origin of clustering and long-term persistence. At the individual level, we find that the correlations in activity are a byproduct of the clustering expressed in the power-law distribution of inter-event times of single users. On the contrary, the activity of the whole community presents long-term correlations that are a true emergent property of the system, i.e. they are not related to the distribution of inter-event times. This result suggests the existence of collective behavior, possible arising from nontrivial communication patterns through the embedding social network.

Various constituents of social systems have been found to follow remarkable statistical regularities. Only the recent availability of relevant data made it possible to unravel such features. Tracking bank notes or cell phones it has been shown that humans follow simple and reproducible mobility patterns [1]. The communication via e-mails occurs in bursts, exhibiting a broad distribution of times between successive messages of individuals (inter-event times) [3]. Recently, we have found that the act of sending messages of individual users in two online communities present long-term correlations [5] characterized by power-law correlation functions obtained via standard Detrended Fluctuation Analysis.

In the present work we examine the relation between the two empirical findings of broad inter-event time distributions [3] and the long-term persistence identified in the communication activity [5]. Therefore, we investigate the communication activity of actants in a social online community with special consideration of the timing and study long-term correlations in the communication as well as clustering of successive messages.

Long-term correlations have been found in the dynamics of many physical, technological, and natural systems. They are characterized by a divergent correlation time, i.e. a power-law decaying auto-correlation function (for a review see [6]). Such correlations lead to a pronounced mountain-valley-structure on all time scales – comprising indeterministic epochs of small and large values [7]. This type of persistence represents a surprising regularity since it is present in many different data such as DNA-sequences, human heartbeat, climatological temperature, etc. [8]. Long-term persistence in human related data has been reported for highway traffic [11], Wikipedia access [13], Ethernet traffic [14], finance and economy [15], written language [18], as well as physiological records [9], Human brain activity [22] and human motor activity [25] also comprise long-term correlations as well as city growth [26], biological networks [30] and the spreading of disease [31].

The distributions of inter-event times (times between successive messages) have been found to be rather broad, described by power-laws [3]. If many short intervals are separated by few long ones, the activity as messages per unit time comprises persistence, i.e. epochs of large and small activity. Since such distributions have been described with power-laws, we wish to investigate the relation between the long-term correlations in activity [5] and the broad (power-law) distribution of inter-event times [3]. We will test two possible scenarios: (i) In the first scenario, the long-term correlations found in the communication activity [5] result from Levy type distributions, i.e. correlations are only due to the power-law inter-event time distribution (with exponents in the specific range) [32]. In the second scenario, (ii) the activity comprises ’real’ correlations, i.e. the inter-event time distributions do not follow a power-law, but the communication activity is temporally not independent, namely long-term correlated.

We study the activity of sending messages based on detailed temporal data from a social online community and obtain the long-term correlation exponent via DFA. The exponent depends on the overall activity of the members; the more active the members the larger the fluctuation exponents. This exponents reaches a value for the most active users from an uncorrelated value for the less active ones. Then, we compare the value of with the corresponding exponents of randomized data and a theoretical prediction relating correlations with clustering in the inter-event times. From the consistency of the comparison of this three measures, we conclude that the long-term correlations found in the activity of sending messages for single users is a direct consequence of the power-law distributed inter-event time of the individuals. Thus, the burstiness in the user activity explains the long-range correlations.

More interesting results are found when we consider the activity of the whole community as a sum of the activity of its members. Again we find non-trivial long-range correlations with exponents in the same range as the individual users. However, the origin of this correlations is not related to the inter-event activity. This is probed by shuffling the activity data but preserving the distribution of inter-event times. In this case, this shuffling destroys the long-range correlations, implying that the correlations are not a byproduct of the broad distribution of inter-event times. We conclude that the whole system acts as a true long-term correlated system where correlations are not directly related to the Levy distributions of events.

We analyze the data of an online community (www.pussokram.com, POK, [33]) covering the complete lifetime of the community over days from February 2001 until June 2002. We record the activity among almost members with more than 500,000 messages sent. This internet-site has been used for general social interactions and dating. The data consist of the time when the messages are sent and anonymous identification numbers of the senders and receivers. The data has been analyzed by us in [5]. In contrast to similar network data sets consisting only of snapshots, i.e. temporally aggregated social networks expressing who sent messages to whom, the advantage of this data set is that it provides the exact time when the messages were sent. For a discussion see [37].

Before shutdown, the members could log in and meet virtually. In such communities, there are different ways of interacting. Usually, it is possible to choose favorites, i.e. certain members, that a person somehow feels committed to. Such platforms also offer the possibility to discuss in groups with other members about specific topics. We focus on messages sent among the members – they are similar to e-mails but have the advantage that they are sent within a closed community where there are no messages coming from or going outside. Figure ? illustrates patterns of sending messages for typical single users [a-d] and for the whole community [e]. The data is publically available^{1}

## Results

### Study of correlations in individual activity

Applying DFA [38] we have found in [5] that the individual activity records, , i.e. messages per unit time (records of messages per day or per week), exhibit long-term correlations. The fluctuation function provided by DFA scales as

where the exponent is also known as the Hurst exponent. In the case of long-term correlations – which are characterized by a power-law decaying auto-correlation function:

where denotes the average, is standard deviation of , and is the correlation exponent () – one finds , whereas larger exponents correspond to more pronounced long-term correlations. For uncorrelated or short-term correlated records (, or in general , is the substrate dimension) the asymptotic fluctuation exponent is . In the range both exponents are related via

For an overview, we refer to [39]. DFA removes polynomial trends of the order from the original record , i.e. DFA2 copes with linear trends.

It is important to note that the DFA fluctuation function Eq. (Equation 1) is not applied to the activity , but to the integrated signal . Thus would be the analogous to the steps in a random walk and the displacement. DFA incorporates an additional detrending of the data. The integration leads to the appearance of long range correlation when the interval between each step is power-law distributed. We will come back to this result when explaining the long-term correlations in terms of the burstiness.

We have measured the fluctuation exponents by applying least squares fits to vs. on the scales weeks conditional to the member’s activity level, e.g. their total number of messages, [5]. Figure ? depicts the DFA results. We find that the less active members, sending only one or two messages in the period of data acquisition, exhibit uncorrelated behavior. The more messages the members send, the more correlated is their activity. The fluctuation exponent increases with and reaches values up to (value obtained for sending messages, we disregard the last points, , which have too large errors bars). The uncorrelated behavior– – for small activity can be understood since when there is not enough time in the data acquisition window to capture long-term correlations. In [36] we propose a model which reproduces the dependence of the fluctuation exponents on the activity level of the members. For receiving messages we find almost identical results [36]. We use weekly resolution in order to cope with possible weekly oscillations [40].

Similar long-term correlations have been found in [43] in traded values of stocks and e-mail communication. The fluctuation exponent increases with the mean trading activity of the corresponding stock or with the average number of e-mails similarly as in our results.

### Study of clustering in individual activity

The timing of human communication activity has been found to comprise bursts where many events occur in relatively short periods which are separated by long periods with few or no events at all. Such patterns can be characterized with the inter-event times, i.e. the times, , between successive messages. For e-mail communication it has been argued that their probability density follows a power-law,

with exponent [3]. As an origin for such heavy tails in human dynamics a queuing model has been suggested [3] according to which each individual performs tasks from a priority list. It has been confirmed that such a process can reproduce bursts of activity or clustering, see e.g. [47]. In contrast, analyzing the same e-mail data, a log-normal distribution has been found to be more appropriate to describe the inter-event time distribution [49]. We would like to remark that fitting fat tailed distributions is disputed [51]. There is neither a consensus on a typical functional form nor on a proper fitting technique. Recently, a cascading Poisson process based on daily and weekly cycles has been proposed as origin of slower-than-exponential decays of [4]. We studied the cascading Poisson process in [36].

In [55], memory in the sequences of has been studied for different data sets, characterizing the inter-event times in terms of a burstiness parameter, which is based on the distribution, and in terms of a memory coefficient, which is the auto-correlation function at lag . In addition, the authors locate the corresponding data sets in a phase diagram defined by these two quantities. Nevertheless, we would like to note that the quantification of long-term correlations in the can be hindered by noise [56].

Next, we study the POK data, i.e. the inter-event times between successive messages of individual members, and relate their statistics to the long-term correlations. The finding of long-term correlations opens the question of the origin of such a persistence pattern in the social communication. From a statistical physics point of view, we consider two possible scenarios:

In the first scenario, the intervals between the messages follow a power-law [3]. Accordingly, the activity pattern comprises many short intervals and few long ones, implying persistent epochs of small and large activity. This fractal-like clustering in the activity can – depending on the exponent – lead to long-term correlations with (see the analogous problem of the origin of long-term correlations in DNA sequences as discussed in [59]). This scenario implies a direct link between the correlations in the activity and the distribution of inter-event times which can be obtained analytically [60]. We call this scenario “Levy correlations” since the actual activity may not be correlated per-se, but correlations arise as a byproduct of integrating a signal with a power-law distribution of inter-events in the DFA formalism.

In the second scenario, the intervals between the messages may or may not follow a power-law distribution, but the values of the inter-event times are not independent of each other and comprise ’real’ long-term persistence. For example, the distribution of inter-event times could be stretched exponential (see recent work on the study of extreme events of climatological records exhibiting long-term correlations [61]) and then the only way to explain long-term correlations in the activity are correlations in the inter-event times. We call this scenario “true correlations” since the correlations are not related to the distribution of inter-events but they reflect ’real’ correlations in the dynamics of the communication activity.

A possible way to discern between these two scenarios is to shuffle the temporal activity, keeping the inter-event distribution intact. While in the case of Levy type correlations shuffling the inter-event times should not influence the long-term correlation properties of , in the case of ’real’ long-term correlations shuffling the inter-event times should destroy the (asymptotic) long-term correlations since the memory is due to the arrangement of the inter-event times. In what follows, we investigate the activity of individual members and the activity of the whole POK community.

### Study of inter-event distribution of individual members

Figure ? exhibits the fluctuation exponents for individual members when we shuffle the data but preserve the distribution of inter-event times. The corresponding exponents also reach high values, almost as high as for the original data, and do not drop for very active members. This agreement is a first indication of Levy correlations in single user activity.

Further evidence is found by studying the distribution of inter-event times in the activity of each individual. Figure ? shows the probability density, , of times between messages of the same users sent in the online community. A power-law regime of approximately two decades can be seen with an exponent , which differs from the exponent reported for e-mail communication [3], i.e. . A reason for these different findings might be that in the case of [3] only one user is considered and that depends on the activity level of the users, as we show below. In addition, here we study all messages from a closed community. The exponent we find is closer to the one reported for reply times (waiting times), i.e. the time individuals spend between receiving and sending to the same communication partner. For reply times of e-mails and land mail has been reported [3].

Since we found a dependence of the fluctuation exponent on the activity level , i.e. the total number of messages each member sends, we suspect that also might depend on . Thus, in Fig. ? we plot for sending messages in POK (daily resolution) the for groups of different activities, i.e. different total number of messages . We find that for the most active members decays rather steeply, while for the least active members decays much slowly. Due to the finite size of the data it is not quite clear which functional form the curves follow. If one assumes a power-law decay then the exponents are roughly in the range .

As discussed above, the power-law distribution of inter-event times, Eq. (Equation 3), can lead to long-term correlations in activity, without requiring temporal dependencies between the intervals themselves. It can be shown that the long-term persistence properties of this point process are characterized by the fluctuation exponent which theoretically depends on according to [32]:

see Fig. ?. Apart from detrending, DFA provides an integration of the original record. So if there are long periods of no activity due to power-law inter-event times, then, this is reflected in long persistence in the signal calculated by DFA. Thus, the existence of long-term correlations is due to the long periods distributed via Levy distributions as expressed by the direct relation between correlations and Levy inter-event activity, Eq. (Equation 4).

Applying least squares fits (in the straight range) to the for sending in POK (Fig. ?) we obtain values for as a function of the activity level and determine the corresponding fluctuation exponents, , as expected from Eq. (Equation 4). We would like to note that the curves in Fig. ? are not always straight lines leading to large uncertainty regarding the estimated values of .

Figure ? depicts the fluctuation exponents from Eq. (Equation 4) in comparison with the values obtained from DFA. We find for a big part of the range. The exponents are also close to of the shuffled records where the inter-event times are preserved. The fact that when we shuffle the signal, respecting the corresponding distribution of inter-event, gives rise to the same correlation function, indicates that the origin of the long-term correlation obtained in DFA are due to the Levy correlations. This is further corroborated by the agreement between from DFA and the prediction . From Fig. ? we see that the three curves are in a reasonable agreement. This supports that the correlations in single user activity can be due to the power-law distribution of the inter-event times, which is in favor of Levy type correlations.

### Study of whole community activity

Next, we investigate the activity of the community as a whole. While we have studied the activity of single users, it is of interest to investigate the activity of the whole community by considering the number of messages sent by all members in a specified period of time. Figure ?(e) shows such activity temporally aggregated to one day. The interest arises since we would like to test the existence of correlations arising from collective behavior in the communication patterns at the level of the whole community.

For this study, we disregard who sends the messages to whom and only consider the instants when any message was sent. In order to have a sufficiently long record to apply DFA, we aggregate the data to messages per hour (instead of daily or weekly resolution). As can be seen in Fig. ?(e), the record contains oscillations [4]. Since such periodicities lead to erratic fluctuation functions [39], we subtract the hourly averages over all days: .

The DFA fluctuation functions are shown in Fig. ?. The hump on scales around hours in the results of DFA1 and DFA2 are residual oscillations, i.e. they were not completely removed. On larger scales this effect vanishes and we find a fluctuation exponent . The straight line in the case of DFA0 is due to the fact that the maximum exponent is 1 [39]. More importantly, when the record of the whole community is shuffled but preserving the inter-event distribution, the asymptotic scaling is . That is, in contrast to the result for individual activity, when the shuffle the signal of the whole community, we obtain the uncorrelated exponent: (dashed lines in Fig. ?). The fact that the correlations vanish ( ) when the data is shuffled indicates that the long-term correlations found in the activity of the community as a whole are not due to Levy correlations. Instead, correlations in the whole community are “true correlations” appearing as a manifestation of collective behavior of the scale of the entire community.

Another surprise appears when we calculate the distribution of inter-event times for the whole community. Here we define inter-event the time between the sending of two consecutive messages of any member in the community. This contrasts with the same study done at the single user level (Fig. ?) when inter-event is defined as the time between two events of the same user. In a sense, for the entire community captures the collective behavior emerging from the entire community as the information travels through the whole network.

In Fig. ? the resulting probability density is displayed. We find a plateau up to seconds followed by a power-law decay according to Eq. (Equation 3) with . Thus, the distribution of inter-event activity of the community as a whole is also a Levy type like the single user activity, albeit with a larger exponent. Such a larger exponent reflects the fact that is narrower for the community than for the individuals, as expected.

When we convert the exponent to the through the Levy distribution model, Eq. (Equation 4), we find . Thus surprisingly, Eq. (Equation 4) may also explain the persistence as in the individual activity. However, the main evidence of Fig. ?, that is, the fact that the correlations vanish when we shuffle the data, probe that, even if Eq. (Equation 4) provides a good estimation of , the long-term correlations are due to ’real’ correlations and are not an artifact of the integration of a Levy type activity with DFA.

The long-term correlations found in the behavior of the entire community is more understandable than in the activity of single members, since the activity of the community is based on the communication patterns of the messages and information flowing through the whole system. The existence of at the whole level and the indications that the correlations are real ones is an interesting instance of the emergence of critical behavior in the collective dynamics of the system as a whole.

We conclude that while at the individual level we find Levy correlations, the activity of the whole community comprises ’real’ correlations, which is due to the (possibly correlated) superposition of the individuals activity into a collective self-organized information flow in the system. Such a behavior is reminiscent of critical systems in phase transitions.

## Discussion

We have studied the timing of communication in a social online community and find long-term persistence in the activity of sending messages at the single user level and the whole community level. Furthermore, we have addressed the question of the origin of these long-term correlations and whether these are Levy type or ’real’ correlations. While in the case of Levy type correlations the inter-event times need to be power-law distributed, ’real’ long-term correlations are independent of the distributions, since they are due to interdependencies in the activities.

Our work, then, still leaves unanswered the question of the cause of the long-term persistence in the communication patterns at the whole level. One possibility is that the temporal correlations are related to correlations in the network structure [64]. The persistence could also be due to social effects, i.e. the dynamics in the social network [66] induces persistent fluctuations, such as cascades. An example could be that a group of friends tries to make an appointment and therefore sends many subsequent messages in a relatively short time [67]. After agreeing, the communication activity among the group drops. The activity patterns of individuals could be understood as a superposition of many such cascades. On the other hand, it could be purely due to a state of mind [23], solipsistic, emerging from moods. More research is needed to thoroughly understand the interesting properties of human activity and its motives.

In conclusion, we have determined exponents to characterize communication activity: (i) , the fluctuation exponent of the original data, (ii) , the fluctuation exponent when the data is shuffled preserving the inter-event times, (iii) , the fluctuation exponent which is expected from power-law distributed inter-event times. We find that which supports the hypothesis of Levy correlations in the single user activity, while we find for the collective behavior of the whole community revealing non-trivial long-term correlations and self-organization at the level of the whole system.

We should mention a third scenario which we leave for future work. It is possible that the correlations comprise more complex features. It has been shown that nonlinear correlations in multifractal data sets lead to power-law distributed inter-event times (of peaks over threshold) [68]. In fact, the authors of [68] find in their Fig. 1(c) a similar dependence of on the total number of events as we do for in our Fig. ?. Additional analysis is needed to fully characterize the multifractal properties [69] of communication activity via e-mails or messages in online communities.

## Acknowledgments

We thank C. Briscoe, J.F. Eichner, L.K. Gallos, and H.D. Rozenfeld for useful discussions. This work was supported by National Science Foundation Grants NSF-SES-0624116 and NSF-EF-0827508 and ARL. F.L. acknowledges financial support from The Swedish Bank Tercentenary Foundation. S.H. thanks the European EPIWORK project, the Israel Science Foundation, ONR and DTRA for financial support.

## Author contributions

All authors contributed equally to the work presented in this paper.

## Additional information

The authors declare no competing financial interests.

FIG. ?. Examples of activity of sending messages and overall activity in POK. The vertical lines in (a) and (c) represent the instants when the messages have been sent by two arbitrary members. The panels (b) and (d) show the records of number of messages per day, , of the same two members. The record of the total number of messages sent by all members per day within POK is depicted in (e). (a) and (b): member 326 (); (c) and (d): member 9414 ().

FIG. ?. Fluctuation exponents of the communication activity sending messages. The exponents are plotted as a function of the activity level (final number of messages) for the original data (green circles), shuffled data preserving the individual inter-event times (blue squares), as well as the exponents expected from Eq. (Equation 4) (brown triangles down) and from the distribution of inter-event times as characterized by power-law fits to the curves of Fig. ? providing the exponent .

FIG. ?. Probability density of inter-event times between successive messages sent by a single member of POK, in daily resolution. The values are extracted by considering every single individual sending messages in the period of data acquisition and then joined from all members. The dotted line in the top corresponds to the exponent .

FIG. ?. Probability density of inter-event times between successive messages sent by all individual members of POK in daily resolution. The values are extracted for the individuals and unified among members according to their activity level . The curve for the most active members is in the bottom, while the one for the least active is in the top. The dotted lines correspond to the exponents (top) and (bottom).

FIG. ?. Levy correlations and persistence. Theoretical relation between the inter-event time distribution exponent and the fluctuation exponent according to Eq. (Equation 4).

FIG. ?. Fluctuation function of the record of messages sent by any member of POK. The record is the same as in Fig. ?(e) but in hourly resolution. Prior to applying DFA, the record has been deseasoned according to . The different curves differ in the DFA-order (DFA0-DFA2, from top to bottom), which determines the capability of detrending. DFA2 eliminates linear trends in [39]. The dotted line in the bottom corresponds to a power-law with exponent and serves as guide to the eye while the continuous line at the top represents .

FIG. ?. Probability density of inter-event times between successive messages sent by any member of POK in seconds. The dotted straight line has corresponds to a power-law with exponent and serves as guide to the eye.

### Footnotes

- http://lev.ccny.cuny.edu/hmakse/soft_data.html

### References

**The scaling laws of human travel.**

D. Brockmann, L. Hufnagel, and T. Geisel.*Nature*, 439(7075):462–465, 2006.**Understanding individual human mobility patterns.**

M. C. Gonzalez, C. A. Hidalgo, and A.-L. Barabási.*Nature*, 453(7196):779–782, 2008.**The origin of bursts and heavy tails in human dynamics.**

A.-L. Barabási.*Nature*, 435(7039):207–211, 2005.**A poissonian explanation for heavy tails in e-mail communication.**

R. D. Malmgren, D. B. Stouffer, A. E. Motter, and L. A. N. Amaral.*Proc. Nat. Acad. Sci. U.S.A.*, 105(47):18153–18158, 2008.**Scaling laws of human interaction activity.**

D. Rybski, S. V. Buldyrev, S. Havlin, F. Liljeros, and H. A. Makse.*Proc. Nat. Acad. Sci. U.S.A.*, 106(31):12640–12645, 2009.*Encyclopedia of Complexity and System Science*, chapter entry 00620: Fractal and Multifractal Time Series.

J. W. Kantelhardt. Springer, 2009.**Method for generating long-range correlations for large systems.**

H. A. Makse, S. Havlin, M. Schwartz, and H. E. Stanley.*Phys. Rev. E*, 53:5445–5449, 1996.**Long-range correlations in nucleotide sequences.**

C.-K. Peng, S. V. Buldyrev, A. L. Goldberger, S. Havlin, F. Sciortino, M. Simons, and H. E. Stanley.*Nature*, 356(6365):168–170, 1992.**Long-range anticorrelations and non-gaussian behavior of the heartbeat.**

C.-K. Peng, J. Mietus, J. M. Hausdorff, S. Havlin, H. E. Stanley, and A. L. Goldberger.*Phys. Rev. Lett.*, 70(9):1343–1346, 1993.**Indication of a universal persistence law governing atmospheric variability.**

E. Koscielny-Bunde, A. Bunde, S. Havlin, H. E. Roman, Y. Goldreich, and H.-J. Schellnhuber.*Phys. Rev. Lett.*, 81(3):729–732, 1998.**Power-law fluctuation in expressway traffic flow: Detrended fluctuation analysis.**

S. Tadaki, M. Kikuchi, A. Nakayama, K. Nishinari, A. Shibata, Y. Sugiyama, and S. Yukawa.*J. Phys. Soc. Jpn.*, 75(3):034002, 2006.**Detrended fluctuation analysis of traffic data.**

Z. Xiao-Yan, L. Zong-Hua, and T. Ming.*Chin. Phys. Lett.*, 24(7):2142–2145, 2007.**Burst event and return interval statistics in wikipedia access and edit data.**

M. Kämpf, S. Tismer, J. W. Kantelhardt, and L. Muchnik.*submitted*, 2011.**On the self-similar nature of ethernet traffic (extended version).**

W. E. Leland, M. S. Taqqu, W. Willinger, and D. V. Wilson.*IEEE/ACM Trans. Networking*, 2(1):1–15, 1994.**Statistical properties of the volatility of price fluctuations.**

Y. Liu, P. Gopikrishnan, P. Cizeau, M. Meyer, C.-K. Peng, and H. E. Stanley.*Phys. Rev. E*, 60(2):1390–1400, 1999.*An Introduction to Econophysics: Correlations and Complexity in Finance*.

R. N. Mantegna and H. E. Stanley. Cambridge University Press, Cambridge, 1999.*The Science of Disasters*, chapter 13. Market Fluctuations I: Scaling, Multiscaling, and Their Possible Origins, pages 373–409.

F. Lux and M. Ausloos. Springer-Verlag, Berlin, 2002.**Long range correlations in human writings.**

A. Schenkel, J. Zhang, and Y.-C. Zhang.*Fractals*, 1(1):47–57, 1993.**Language time series analysis.**

K. Kosmidis, A. Kalampokis, and P. Argyrakis.*Physica A*, 370(2):808–816, 2006.**Sleep-wake differences in scaling behavior of the human heartbeat: Analysis of terrestrial and long-term space flight data.**

P. Ch. Ivanov, A. Bunde, L. A. N. Amaral, S. Havlin, J. Fritsch-Yelle, R. M. Baevsky, H. E. Stanley, and A. L. Goldberger.*EPL*, 48(5):594–600, 1999.**Correlated and uncorrelated regions in heart-rate fluctuations during sleep.**

A. Bunde, S. Havlin, J. W. Kantelhardt, T. Penzel, J.-H. Peter, and K. Voigt.*Phys. Rev. Lett.*, 85(17):3736–3739, 2000.**Long-range temporal correlations and scaling behavior in human brain oscillations.**

K. Linkenkaer-Hansen, V. V. Nikouline, J. M. Palva, and R. J. Ilmoniemi.*J. Neurosci.*, 21(4):1370–1377, 2001.**Spontaneous brain activity as a source of ideal noise.**

P. Allegrini, D. Menicucci, R. Bedini, L. Fronzoni, A. Gemignani, P. Grigolini, B. J. West, and P. Paradisi.*Phys. Rev. E*, 80(6):061914, 2009.**A small-world of weak ties provides optimal global integration of self-similar modules in functional brain networks.**

L. K. Gallos, H. A. Makse, and M. Sigman.*Proc. Nat. Acad. Sci. USA*, 109:2825–2830, 2012.**Endogenous circadian rhythm in human motor activity uncoupled from circadian influences on cardiac dynamics.**

P. Ch. Ivanov, K. Hu, M. F. Hilton, S. A. Shea, and H. E. Stanley.*Proc. Nat. Acad. Sci. U.S.A.*, 104(52):20702–20707, 2007.**Modelling urban growth patterns.**

H. A. Makse, S. Havlin, and H. E. Stanley.*Nature*, 377:608–612, 1995.**Modeling urban growth patterns with correlated percolation.**

H. A. Makse, J. S. Andrade, M. Batty, S. Havlin, and H. E. Stanley.*Phys. Rev. E*, 58:7054–7062, 1998.**Laws of population growth.**

H. D. Rozenfeld, D. Rybski, J. S. Andrade Jr., M. Batty, H. E. Stanley, and H. A. Makse.*Proc. Nat. Acad. Sci. USA*, 105:18702–18707, 2008.**The area and population of cities: New insights from a different perspective on cities.**

H. D. Rozenfeld, D. Rybski, and H. A. Gabaix, X. Makse.*American Economic Review*, 101:560–580, 2011.**Modularity map of the network of humn cell differentiation.**

G. Galvao, J. G. V. Miranda, R. F. S. Andrade, J. S. Andrade Jr., L. K. Gallos, and H. A. Makse.*Proc. Nat. Acad. Sci. USA*, 107:5750–5755, 2010.**Collective behavior in the spatial spreading of obesity.**

L. K. Gallos, P. Barttfeld, S. Havlin, M. Sigman, and H. A. Makse.*Sci. Rep.*, 2012.**Lévy dynamics of enhanced diffusion: Application to turbulence.**

M. F. Shlesinger, B. J. West, and J. Klafter.*Phys. Rev. Lett.*, 58(11):1100–1103, 1987.**Network dynamics of ongoing social relationships.**

P. Holme.*EPL*, 64(3):427–433, 2003.**Network bipartivity.**

P. Holme, F. Liljeros, C. R. Edling, and B. J. Kim.*Phys. Rev. E*, 68(5):056107, 2003.**Structure and time evolution of an internet dating community.**

P. Holme, C. R. Edling, and F. Liljeros.*Soc. Networks*, 26(2):155–174, 2004.**Communication activity in social networks: growth and correlations.**

D. Rybski, S. V. Buldyrev, S. Havlin, F. Liljeros, and H. A. Makse.*Eur. Phys. J. B*, 84(1):147–159, 2011.**How people interact in evolving online affiliation networks.**

L. K. Gallos, D. Rybski, F. Liljeros, S. Havlin, and H. A. Makse.*submitted*, 2011.**Mosaic organization of dna nucleotides.**

C.-K. Peng, S. V. Buldyrev, S. Havlin, M. Simons, H. E. Stanley, and A. L. Goldberger.*Phys. Rev. E*, 49(2):1685–1689, 1994.**Detecting long-range correlations with detrended fluctuation analysis.**

J. W. Kantelhardt, E. Koscielny-Bunde, H. H. A. Rego, S. Havlin, and A. Bunde.*Physica A*, 295(3–4):441–454, 2001.**Rhythms of social interaction: messaging within a massive online network.**

S. Golder, D. M. Wilkinson, and B. A. Huberman.*online-arXiv*, arXiv:cs/0611137v1 [cs.CY], 2006.**Planetary-scale views on an instant-messaging network.**

J. Leskovec and E. Horvitz.*online-arXiv*, arXiv:0803.0939v1 [physics.soc-ph], 2008.**On universality in human correspondence activity.**

R. D. Malmgren, D. B. Stouffer, A. S. L. O. Campanharo, and L. A .N Amaral.*Science*, 325(5948):1696–1700, 2009.**Scaling theory of temporal correlations and size-dependent fluctuations in the traded value of stocks.**

Z. Eisler and J. Kertész.*Phys. Rev. E*, 73(4):046109, 2006.**Fluctuation scaling in complex systems: Taylor’s law and beyond.**

Z. Eisler, I. Bartos, and J. Kertész.*Adv. Phys.*, 57(1):89–142, 2008.**Probing human response times.**

A. Johansen.*Physica A*, 338(1–2):286–291, 2004.**Comment on A.-L. Barabasi, Nature****435**207-211 (2005).

A. Johansen.*online-arXiv*, arXiv:physics/0602029v1 [physics.soc-ph], 2006.**Exact results for the barabasi model of human dynamics.**

A. Vázquez.*Phys. Rev. Lett.*, 95(24):248701, 2005.**Modeling bursts and heavy tails in human dynamics.**

A. Vázquez, J. G. Oliveira, Z. Dezsö, K. I. Goh, I. Kondor, and A.-L. Barabási.*Phys. Rev. E*, 73(3):036127, 2006.**Comment on “The origin of bursts and heavy tails in human dynamics” by Barabasi, Nature****435**, 207 (2005).

D. B. Stouffer, R. D. Malmgren, and L. A. N. Amaral.*online-arXiv*, arXiv:physics/0510216v1 [physics.data-an], 2005.**Reply to comment on “the origin of bursts and heavy tails in human dynamics”.**

A.-L. Barabási, K.-I. Goh, and A. Vazquez.*online-arXiv*, arXiv:physics/0511186v1 [physics.data-an], 2005.**Power laws, Pareto distributions and Zipf’s law.**

M. E. J. Newman.*Contemp. Phys.*, 46(5):323–351, 2005.**Log(rank-1/2): a simple way to improve the ols estimation of tail exponents.**

X. Gabaix and R. Ibragimov. Discussion Paper 2106 (26), Harvard Institute of Economic Research, Cambridge, Massachusetts, February 2006.**Power-law distributions in empirical data.**

A. Clauset, C. R. Shalizi, and M. E. J. Newman.*SIAM Rev.*, 51(4):661–703, 2009.**Gibrat’s law for cities: uniformly most powerful unbiased test of the Pareto against the lognormal.**

Y. Malevergne, V. Pisarenko, and D. Sornette.*online-arXiv*, arXiv:0909.1281v1, 2009.**Burstiness and memory in complex systems.**

K.-I. Goh and A.-L. Barabási.*EPL*, 81(4):48002, 2008.**Statistics of return intervals in long-term correlated records.**

J. F. Eichner, J. W. Kantelhardt, A. Bunde, and S. Havlin.*Phys. Rev. E*, 75(1):011128, 2007.**Eliminating finite-size effects and detecting the amount of white noise in short records with long-term memory.**

S. Lennartz and A. Bunde.*Phys. Rev. E*, 79(6):066101, 2009.**Random walk models for spike activity of single neuron.**

G. L. Gerstein and B. Mandelbrot.*Biophys. J.*, 4(1P1):41–68, 1964.**Generalized Lévy-walk model for DNA nucleotide sequences.**

S. V. Buldyrev, A. L. Goldberger, S. Havlin, C.-K. Peng, M. Simons, and H. E. Stanley.*Phys. Rev. E*, 47(6):4514–4523, 1993.*Encyclopedia of Complexity and System Science*, volume Fractals and multifractals, chapter Fractals in Biology.

S. V. Buldyrev. Springer, 2010.**Long-term memory: A natural mechanism for the clustering of extreme events and anomalous residual times in climate records.**

A. Bunde, J. F. Eichner, J. W. Kantelhardt, and S. Havlin.*Phys. Rev. Lett.*, 94(4):048701, 2005.**Darwin and Einstein correspondence patterns.**

J. G. Oliveira and A.-L. Barabási.*Nature*, 437(7063):1251–1251, 2005.**Analysis, synthesis, and estimation of fractal-rate stochastic point processes.**

S. Thurner, S. B. Lowen, M. C. Feurstein, C. Heneghan, H. G. Feichtinger, and M. C. Teich.*Fractals*, 5(4):565–595, 1997.**Mechanisms and models of human dynamics.**

A. Kentsis.*Nature*, 441(7092):E5–E5, 2006.**Quantifying long-range correlations in complex networks beyond nearest neighbors.**

D. Rybski, H. D. Rozenfeld, and J. P. Kropp.*EPL*, 90(2):28002, 2010.**Dynamical and bursty interactions in social networks.**

J. Stehle, A. Barrat, and G. Bianconi.*Phys. Rev. E*, 81(3):035101, 2010.**Quantifying social group evolution.**

G. Palla, A.-L. Barabási, and T. Vicsek.*Nature*, 446(7136):664–667, 2007.**Effect of nonlinear correlations on the statistics of return intervals in multifractal data sets.**

M. I. Bogachev, J. F. Eichner, and A. Bunde.*Phys. Rev. Lett.*, 99(24):240601, 2007.**Multifractal detrended fluctuation analysis of nonstationary time series.**

J. W. Kantelhardt, S. A. Zschiegner, E. Koscielny-Bunde, S. Havlin, A. Bunde, and H. E. Stanley.*Physica A*, 316(1–4):87–114, 2002.**Multifractality of river runoff and precipitation: comparison of fluctuation analysis and wavelet methods.**

J. W. Kantelhardt, D. Rybski, S. A. Zschiegner, P. Braun, E. Koscielny-Bunde, V. Livina, S. Havlin, and A. Bunde.*Physica A*, 330(1–2):240–245, 2003.**Long-term persistence and multifractality of precipitation and river runoff records.**

J. W. Kantelhardt, E. Koscielny-Bunde, D. Rybski, P. Braun, A. Bunde, and S. Havlin.*J. Geophys. Res.-Atmos.*, 111(D1):D01106, 2006.