Closeness of human behavior to model.

From Phil Runkel on 3 March 97.

Some time ago some people on the CSGnet were proposing various
numerical indicators of the closeness of a trace of human
behavior to a trace of model behavior. The motivation was that a
correlation tells the extent to which one trace parallels or
makes the same bends as another, but it does not tell the
distance between the traces. Instead of the word "trace," you
could use "profile" or "graph" or even "curve."

It seems to me a good measure could be made by some integrating
calculation of the area between the profiles. But my calculus is
not good enough to propose a formula. Anyway, maybe Cronbach's
d-statistic comes close to that. Did anybody propose Cronbach's
d?

The distance between two points p1 and p2 in 2-space is [(x1 -
x2)^2 + (y1 - y2)^2]^1/2. You can think of the trace or profile
produced by the human as a single point in a multidimensional
space, and of the profile produced by the model as another point
in that same space. The number of dimensions is the number of
pairs of points that are compared in the run of the experiment.
If, in every pair, one point falls exactly upon the other, the
distance in the pair is zero, and all the distances sum to zero.
The core of the measure, if I am remembering this correctly, is
[(x1 - x2)^2 + (y1 - y2)^2 + ... + (n1 - n2)^2]^1/2. I forget
the precise formula that Cronbach recommended. If anybody has
any further interest in this, I'll be glad to go find a
statistics book and look it up. (I gave away all my statistics
books a few years ago.)

Cronbach developed his measure for much the same reason the
question arose on our net. Researchers were wanting to assess
person 1 on several dimensions and person 2 on the same
dimensions and then compare the two profiles. Most measures told
about the degree of parallelism or the elevation of the one
compared to the other, but it was hard to find one that partook
of both.

I have seen the "stability factor" mentioned in several letters
on the CSGnet. Maybe that measure does as well.

And maybe my memory of Cronbach's d-statistic is faulty. It has
been a long time since I have had to think about such matters.

[From Bill Powers (970303.1618 MST)]

Phil Runkel, 3 March 97.

Some time ago some people on the CSGnet were proposing various
numerical indicators of the closeness of a trace of human
behavior to a trace of model behavior.

...
Did anybody propose Cronbach's d?

The distance between two points p1 and p2 in 2-space is [(x1 -
x2)^2 + (y1 - y2)^2]^1/2.

...
If anybody has

any further interest in this, I'll be glad to go find a
statistics book and look it up. (I gave away all my statistics
books a few years ago.)

What we're doing is comparing two records of handle position, point by
point. This means that x is always the same for both points being compared,
so the only difference we care about is in y. If you compute the squared
distance between the points, add all of the squares up, take the average,
and then take the square root, what have you got? Right. Cronbach's
whatchamacallit in one dimension, otherwise known as the RMS difference
between the records.

And maybe my memory of Cronbach's d-statistic is faulty. It has
been a long time since I have had to think about such matters.

Count your blessings.

Best,

Bill P.

[Martin Taylor 970304 11:00]

Phil Runkel on 3 March 97.

Some time ago some people on the CSGnet were proposing various
numerical indicators of the closeness of a trace of human
behavior to a trace of model behavior. The motivation was that a
correlation tells the extent to which one trace parallels or
makes the same bends as another, but it does not tell the
distance between the traces.

As a main contributor to that discussion, my motivation was a bit
different. If there is control of a certain precision, as assessed
by the stability factor, the correlation, or whatever parameter one
uses, there is a multitude of possible control system structures that
could achieve this precision. The stability factor (etc., etc.) shows
_that_ control is happening, but not _which control structure_ is a
better model than which other that has the same degree of fit.

You can think of the trace or profile
produced by the human as a single point in a multidimensional
space, and of the profile produced by the model as another point
in that same space. The number of dimensions is the number of
pairs of points that are compared in the run of the experiment.

Yes, that's the approach I used. Then my idea was to compare three
points in this space, not two. One is the "ideal" of perfect control,
one the actual human trace, and one the simulation trace. We have
chosen the simulation model as the best fit to the human data that
we can get by varying the parameters for control systems of that
structural type. This means that we are concerned with possible
control systems of other structural types that would fit equally well.

But the simulated control system is just that, a _control_ system. So
some part of the fit will be due to the fact that the human controls,
and the simulation controls. The question is how much of the fit is
due to that fact, and how much is due to the simulation having a
structure that matches the human's control structure reasonably
well.

Now we think about the layout of points in this high-dimensional space.
Two observations apply: (1) the higher the dimensionality, the more
likely it is that two vectors chosen at random will be very close to
orthogonal, and (2) we are dealing with only three points, so we can
always fit a plane--a two-dimensional flat surface--through them. We
are therefore dealing with a triangle and the angles of the triangle.
Call the corners of the triangle C (centre--the ideal), H (human) and
S (simulation). High dimensionality means having many _independent_
sample points. If the control bandwidth is, say, 5 Hz, there are 10
independent sample points per second, no matter how many measurements
are made.

Let's consider extreme cases:

(1) The structure of the simulation is unrelated to that of the human,
and the fit is entirely due to the fact that both the human and the
simulation control equally well. The possible structures for which the
best fit is exactly this good provide data that fall on a circle around
point H. If the fit is largely due to the fact that both human and
simulation control, then the point S will be closer to C than to H, and
in the extreme will be directly between H and C. The simulation will
control too well. The angle CHS approaches zero.

(2) The structure of the simulation is perfect so far as it goes (by which
I mean stuff like being a one-level model, or not having adaptive gain
when the human does--if the human does), and mismatch is due to
things like mouse-slips, momentary failures of attention, and other
things not modelled. Now we can say that the simulation controls because
it has the structure of the human being modelled (so far as it goes)
and the human controls. If the human didn't control, neither would the
simulation, in this case. Of all the possible structures we might have
chosen equally like the human structure, some will control better than
the human, some worse, but we can't say anything about those, since we
don't know what they are. In this case, the angle CHS approaches 90 degrees.

(3) There's a rather absurd possibility, but it has to be included, and that
is that the simulation fits the human data _in spite of_ the fact that the
human controls. This case corresponds to values of the angle CHS appreciably
greater than 90 degrees. If this case ever showed up, I'd worry about the
programming, the modelling, or something. Occasional high values of CHS
are to be expected simply from statistical considerations, but for long
tracks they should be _very_ rare.

All of these cases depend on there being a distance measure of some kind,
whether it be the RMS measure you propose, a correlation measure, or some
other measure of fit. The key is to measure the angle CHS. The noisier
the data, the nearer that angle will be to 45 degrees, where the data are
equally due to the fact that the simulation control and to the fact that
it controls the same way that the human controls (that's an intuitive
statement, not an analytic one, and it could well be wrong).

To recap: values of angle CHS near zero imply inappropriate structures
for the simulation, values near 90 degrees imply good structures, values
much greater than 90 degrees should lead to soul-searching on the part of
the programmer or the analyst, and values near 45 degrees should suggest
that the data may be poor (or that the structure lacks something but also
has something right about it).

Martin

[Hans Blom, 970305d]

(Philip Runkel <runk@OREGON.UOREGON.EDU>, 3 March 97)

Some time ago some people on the CSGnet were proposing various
numerical indicators of the closeness of a trace of human
behavior to a trace of model behavior. The motivation was that a
correlation tells the extent to which one trace parallels or
makes the same bends as another, but it does not tell the
distance between the traces. Instead of the word "trace," you
could use "profile" or "graph" or even "curve."

It seems to me a good measure could be made by some integrating
calculation of the area between the profiles.

Mathematics has invented a huge number of different distance measures. All
have the same purpose: to be able to compare things, i.e. to provide strict
definitions of concepts such as "more" or "less". Which distance measure is
most suited under which circumstances remains an unsolved problem which has
only heuristic solutions. The measure that you discuss is the Euclidean
distance. The Manhattan distance is a different measure; the Manhattan
distance between two points (x1,y1) and (x2,y2) is defined as

        abs (x1-x2) + abs (y1-y2)

which can also be extended to higher dimensions. One group of distance
measures has the form

        abs (x1-x2)^N + abs (y1-y2)^N

where you recognize the Manhattan distance if N=1 and the Euclidean distance
if N=2. N=infinity is also frequently used.

In Artificial Intelligence and game theory, another frequently used distance
measure is the number of moves from the current position of a game to its
solution. This distance is usually unknown but assumed to be knowable in
principle. Solving a game is then the attempt to perform moves such that at
each move this distance is decremented by one.

So my answer to your question what the most desirable distance measure would
be to compare curves or responses would be: that depends on _what_ you
want to measure the "more" or "less" of. Is it important, for instance, that
the curves were traversed at the same speed (as is often the case when one
compares controllers), or is it sufficient that they have the same shape?

Well, I don't think that this answer helps much. Except, maybe, to
demonstrate that there are no easy answers. One general feature of all
distance measures, however, is that when a distance is zero, we have found
the "optimum" or, in case of a game, the "solution".

Greetings,

Hans

[From Rick Marken (970305.1010 PST)]

Martin Taylor (970304 11:00) --

The stability factor (etc., etc.) shows _that_ control is
happening, but not _which control structure_ is a better model
than which other that has the same degree of fit.

If models have the same degree of fit to data then the only way to tell
which is the better model is to look for predictions that distinguish
the models and then _test_ (collect _more_ data) to see which
predictions match what is actually _observed_. This is called "doing
science" and it used to be considered a very cool way to try
to understand the nature of reality.

Best

Rick

[Martin Taylor 970305 15:00]

Rick Marken (970305.1010 PST)]

Martin Taylor (970304 11:00) --

> The stability factor (etc., etc.) shows _that_ control is
> happening, but not _which control structure_ is a better model
> than which other that has the same degree of fit.

If models have the same degree of fit to data then the only way to tell
which is the better model is to look for predictions that distinguish
the models and then _test_ (collect _more_ data) to see which
predictions match what is actually _observed_. This is called "doing
science" and it used to be considered a very cool way to try
to understand the nature of reality.

I can only conclude that you never read the message you quote in part.

There are two ways, at least, of interpreting your comment. I'll deal
with them in sequence.

My message showed that your statement "The only way ..." is false in
this case (and in many others), if you don't believe that examining
the data is a way of testing predictions.

The message had to do with the degree to which the "goodness" of a model
of human behaviour was generic or fortuitous, not with the degree of
"goodness" itself. A model that fits because it has the right structure
is a more useful model than one that happens to fit because both it and
the human happen to control equally well.

It's a bit like the difference between saying (of a physical causal system,
not a control system) that A correlates with B, and that A influences B.
You can have an amazing fit between A and B even though they have no physical
connection, because C causes both. Or you can have the same amazing fit
because B depends on A. Your statement is the equivalent of saying

"If A and B have a high correlation, the only way to tell whether A influences
B is to look for changes in A and see if they are reflected in B." Of
course they will be.

I could, however, be charitable, and interpret your statement the second way.
If I did that, I would point out that what I was doing was to "look for
predictions that distinguish the models and then _test_ to see which
predictions match what is actually observe." That actually is what I
discussed in the message that you criticize--ways of finding differences
between proposals (I don't use "models" because I don't mean "simulations").

Once again, although I have no evidence that it does any good, I repeat
my suggestion that you read the messages you criticize, at least once,
and possibly even read them all the way through.

Martin

From Phil Runkel on 13 March 97.

        Thanks to Bill Powers and Martin Taylor (and perhaps others I have
not yet read) for replying to my poor algebra. I do not know enough to
comment further.

        Best wishes.