[5] Effect of binning and flipping

The following relates to [2]-[4]:

Martin has pointed out that the natural binning of the raw data and the different flipping operation of Gonzalez et al. could be partly responsible for the marked differences between their and our results. In particular, the double-peak structure might be reduced if our raw data are artificially binned and the same flipping rule (most occupied cell at right hand side) is used.

When applied to the experimental microbead data, the results are considerably closer to the single-peaked, triangular distributions of Gonzalez:

[3] Universal PDF of 2D Brownian and Persistent Trajectories

I am grateful to Martin Reichelsdorfer, who has programmed and applied the CRFSA-method (center, rotate, flip, scale and average), as described in post [2], to various synthetic and experimental trajectories in 2D.

In contrast to the original method of Gonzalez et al., his raw data trajectories consisted of real points in 2D. Cells were only introduced in the final averaging step. Therefore, the flipping operation had to be redefined: It was done in such a way that the last point of the trajectory was always located right of the first point.

He first tested the method for a standard Brownian Random Walk. For a given time period dt, the PDF of a Random Walk is known to be a Gaussian, centered at the origin, with a variance proportional to dt. This is also what one naively expects to see as the ensemble averaged spatial distribution. However, after CRFSA one obtaines a double-peak structure:


Second, Martin applied CRFSA to measured trajectories of micro-beads bound to living MEVO cells. On long time scales, those trajectories show superdiffusive behaviour with fractional powerlaw exponents in the mean squared displeacement (MSD), corresponding to directional persistence in the bead motion. This persistence is a qualitative difference to the Brownian Random Walk and one might expect new features for the ensemble averaged spatial distribution function (SDF). However, the results for the Persistent Walks are not much different from before.


No significant changes are observed when the ensemble of trajectories is divided into groups of similar powerlaw exponent.

These results lead to a couple of questions:

• Why is the CRFSA-SDF of a Random Walk not a Gaussian ?
  
• What exactly is the origin of the double-peak ?
• What makes the trajectories of mobile phone users so qualitatively different from Brownian or peristent walks ?
  

[2] Universal Spatial Distribution of Mobile Agents

In a recent paper "Understanding individual human mobility patterns" published by Marta C. Gonzalez et al. in Nature, the authors tracked the coordinates of mobile phone users over an extended period.

In this study, space was divided into cells (the range of a single mobile phone antenna) and the actual raw data consisted of the sequence of cells in which each user was found. In most of the cells, the same user was found several times. In particular, there was a cell of maximum occupancy for each user.

For each individual, the spatial distribution pattern (formally equivalant to a mass distribution) was artificially centered around the origin. Next, it was rotated around its center of mass, so that the main axis of the inertia tensor became alligned along the x-axis. Next, the pattern was flipped, such that the cell of maximum occupancy was located at the right side. It was then scaled to a fixed standard deviation in both x- and y-directions. The centered, rotated, flipped and scaled distributions were finally averaged over all individual trajectories (CRFSA method) . The result was a "universal" and surprisingly non-trivial distribution pattern:



It seems worthwhile to apply a similar procedure to the trajectories of microbeads attached to the living cytoskeleton. I shall name this project "Universal PDF".