Now that we all agree, or should agree, that Martin’s objection is without merit, and that the analyses of Marken & Shaffer (M&S) and Maoz et al (Maoz) are identical, I will explain why M&S and Maoz came to diametrically different conclusions about the power law.
First, let me summarize the results of the M&S and Maoz analysis. First, both note that the test for fit of a power function to a movement is done by regressing measures of curvature, C, on measures of velocity, V, of the movement using the following regression equation:
log(V) = log(k) + Beta * log (C) … (1)
If the regression results in a value of Beta close to -1/3 with an R^2 greater than .7 then the movement is consided to be power law compliant. M&S and Maoz point of that equation (1) leaves out a covariate caalled affine velocity, Alpha, which is part of the equation that relates C to V. So the regression equation for power ower law analysis should be:
log(V) = log(k) + Beta * log (C) + BetaA*log (Alpha) … (2)
If you run this regression for ANY movement trajectory, regardless of how it is produced, you will always find that log(k) =0, Beta = -1/3 and BetaA = 1/3. That’s because the true mathematical relationship between V and C includes Alpha:
log(V) = -1/3 * log (C) + 1/3 * log (Alpha) … (3)
Both M&S and Maoz found that when Alpha is left out of the regression used to find the fit of a power law to a movement, the value of the exponent of C that is found by that regression will deviate from its true mathermatical value, -1/3, by an amount proportional to:
Delta = (Covariance [log(C), log(Alpha)] /Variance(log(C)))/3
The size of Delta, for any randomly (or systematically) produced smooth movement is generally very small. Maoz ran several different simulations, generating 100 randomly generated movement trajectories made up of 1,000,000 points in three different ways. The typical value of Beta, averaging over the different types of simulations, was -.31 ± .03 and the average size of R2 was 0.61 ± 0.07. As noted by Maoz “Both Beta and its R2 magnitudes are within what is considered by experimentalists to be the range of applicable values for the power-law”. In other words, you see the power law whether the movement trajectory was generated randomly by a computer algorithm or by a living control system.
I have found essentially the same thing in my own simulations, though I didn’t report them in either M&S paper. What I did show in the M&S paper is that you will see a power law whether the movement was generated by a living organism (human) or a non-living one (a helcopter).
Maoz et al understood the unpleasant ramifications of their simulation studies, as can be seen in the fact that they begin the Dscussion of their results by saying “We do not suggest that the power-law, which stems from analysis of human data, is a bogus phenomenon…”. They say this because they showed that noise added to a non-power law conforming trajectory will produce a bogus power law conforming result. Thus, they conclude that since power law researchers typically low pass filter trajectories produced by living organisms it is unlikely that the power law conforming results they find were produced by contaminating noise.
I find this conclusion very puzzling. Maoz et al have just shown that non-power law conforming trajectories are very rare, as indicated by the very small value (.03) of the standard deviation around the average power coefficient value of .31. Assuming that these Beta values are normally distributed, 95% of the time, for ANY movement trajectory, Beta will be found to lie between .25 and .36. And this range is probably even smaller. As Maoz et al point out “standard outlier detection and removal techniques as well as robust linear regression make Beta approach closer to −.33 and increase the R2 value”. Since observing a Beta close to -.33 is highly probable for any trajectory, filtered or not, it would be impossible to tell when an observed power coefficient is “real” (whatever that means) or a consequence of what Maoz demonstrated earlier in the paper: that the observed Beta is a consequence of omitting affine velocity from the regression analysis.
I also find their conclusion puzzling because they had just shown that the true (mathematical) value for Beta is .33 but that the observed value will differ depending on the covariance between curvature and affine velocity. The filtering procedures described by Maoz can’t protect the researcher from falsely concluding that a non-power law conforming trajectory is power law conforming becuase filtering simply changes the covariance between curvature and affine velocity, and does so in (until some mathemetician figures it out) ann unpredicable way, sometime increasing and sometimes decreasing this covariance (as Maoz et al found out).
Rather than concluding, as Maoz et al did, that the power law is not a bogus phenomenon, M&S concluded that the power law is simply a statistical artifact; a result of omitting afffine velocity from the regression used to determine whether or not a trajectory is power law conforming. Per M&S, the power law is a consistent, but irrelevant, side effect of controlling for the changing position of the moved entity, the change resulting from variations in the agent’s reference for the position of that entity.
The power law is exactly analogous to the word “hello” that is written as a side effect of a person controlling a cursor relative to a target in a demo by Powers (does anyone know where the video of that demo is?). The explanation of that consistent, irrelevant side effect is that the disturbance has been created to vary as the inverse of “hello”.
The explanation of the irrelevant side effect that is the power law is the fact that the method of determining conformity to that law omits one of the variables that determines that value of velocity; this omission leads to a slight deviation from -1/3, the true value of the power relationship between curvature and velocity, but this deviation is typicaly fairly small since in almost all smooth curves the covariance between curvature and the omitted variable, affine velocity – which is the cause of the deviation – is small (as demonstrated by Maoz et al)
Best, Rick