Pattern theory


Pattern theory, formulated by Ulf Grenander, is a mathematical formalism to describe knowledge of the world as patterns. It differs from other approaches to artificial intelligence in that it does not begin by prescribing algorithms and machinery to recognize and classify patterns; rather, it prescribes a vocabulary to articulate and recast the pattern concepts in precise language. Broad in its mathematical coverage, Pattern Theory spans algebra and statistics, as well as local topological and global entropic properties.
In addition to the new algebraic vocabulary, its statistical approach is novel in its aim to:
The Brown University Pattern Theory Group was formed in 1972 by Ulf Grenander. Many mathematicians are currently working in this group, noteworthy among them being the Fields Medalist David Mumford. Mumford regards Grenander as his "guru" in Pattern Theory.

Example: Natural Language Grammar

We begin with an example to motivate the algebraic definitions that follow. If we want to represent language patterns, the most immediate candidate for primitives might be words. However, set phrases, such as “in order to" immediately indicate the inappropriateness of words as atoms. In searching for other primitives, we might try the rules of grammar. We can represent grammars as finite state automata or context-free grammars. Below is a sample finite state grammar automaton.
The following phrases are generated from a few simple rules of the automaton and programming code in pattern theory:
To create such sentences, rewriting rules in finite state automata act as generators to create the sentences as follows: if a machine starts in state 1, it goes to state 2 and writes the word “the”. From state 2, it writes one of 4 words: prince, boy, princess, girl, chosen at random. The probability of choosing any given word is given by the Markov chain corresponding to the automaton. Such a simplistic automaton occasionally generates more awkward sentences:
From the finite state diagram we can infer the following generators that creates the signal. A generator is a 4-tuple: current state, next state, word written, probability of written word when there are multiple choices. That is, each generator is a state transition arrow of state diagram for a Markov chain.
Imagine that a configuration of generators are strung together linearly so its output forms a sentence, so each generator "bonds" to the generators before and after it. Denote these bonds as 1x,1y,2x,2y,...12x,12y. Each numerical label corresponds to the automaton's state and each letter "x" and "y" corresponds to the inbound and outbound bonds. Then the following bond table is equivalent to the automaton diagram. For the sake of simplicity, only half of the bond table is shown—the table is actually symmetric.
1x1y2x2y3x3y4x4y5x5y6x6y7x7y8x8y9x9y10x10y11x11y12x12y
1x1
1y1
2x1
2y1
3x1
3y11
4x
4y11
5x
5y1
6x
6y1
7x1
7y
8x
8y1
9x
9y1
10x
10y1
11x1
11y1
12x
12y

As one can tell from this example, and typical of signals that are studied, identifying the primitives and bond tables requires some thought. The example highlights another important fact not readily apparent in other signals problems: that a configuration is not the signal that is observed; rather, its image as a sentence is observed. Herein lies a significant justification for distinguishing an observable from a non-observable construct. Additionally, it provides an algebraic structure to associate with hidden Markov models. In sensory examples such as the vision example below, the hidden configurations and observed images are much more similar, and such a distinction may not seem justified. Fortunately, the grammar example reminds us of this distinction.
A more sophisticated example can be found in the link grammar theory of natural language.

Algebraic foundations

Motivated by the example, we have the following definitions:

  1. A generator, drawn as
    is the primitive of Pattern Theory that generates the observed signal. Structurally, it is a value with interfaces, called bonds, which connects the 's to form a signal generator. 2 neighboring generators are connected when their bond values are the same. Similarity self-maps s: G -> G express the invariances of the world we are looking at, such as rigid body transformations, or scaling.
  2. Bonds glue generators into a configuration, c, which creates the signal against a backdrop Σ, with global features described locally by a bond coupling table called. The boolean function is the principal component of the regularity 4-tuple , which is defined as
    seems to capture the notion of allowable generator neighbors. That is, Regularity is the law of the stimulus domain defining, via a bond table, what neighbors are acceptable for a generator. It is the laws of the stimulus domain. Later, we will relax regularity from a boolean function to a probability value, it would capture what stimulus neighbors are probable.
    Σ is the physical arrangement of the generators. In vision, it could be a 2-dimensional lattice. In language, it is a linear arrangement.
  3. An image captures the notion of an observed Configuration, as opposed to one which exists independently from any perceptual apparatus. Images are configurations distinguished only by their external bonds, inheriting the configuration's composition and similarities transformations. Formally, images are equivalence classes partitioned by an Identification Rule "~" with 3 properties:
    1. ext = ext whenever c~c'
    2. sc~sc' whenever c~c'
    3. sigma ~ sigma whenever c1~c1', c2~c2' are all regular.
    A configuration corresponding to a physical stimulus may have many images, corresponding to many observer perception's identification rule.
  4. A pattern is the repeatable components of an image, defined as the S-invariant subset of an image. Similarities are reference transformations we use to define patterns, e.g. rigid body transformations. At first glance, this definition seems suited for only texture patterns where the minimal sub-image is repeated over and over again. If we were to view an image of an object such as a dog, it is not repeated, yet seem like it seems familiar and should be a pattern.
  5. A deformation is a transformation of the original image that accounts for the noise in the environment and error in the perceptual apparatus. Grenander identifies 4 types of deformations: noise and blur, multi-scale superposition, domain warping, and interruptions.
    With the benefit of being given generators and complete bond tables, a difficult part of pattern analysis is done. In tackling a new class of signals and features, the task of devising the generators and bond table is much more difficult.
    Again, just as in grammars, identifying the generators and bond tables require some thought. Just as subtle is the fact that a configuration is not the signal we observe. Rather, we observe its image as silhouette projections of the identification rule.
Bond
Values
012345
011
111
21
3
4
5

Entropy

Pattern Theory defines order in terms of the feature of interest given by p.

Statistics

Grenander's Pattern Theory treatment of Bayesian inference in seems to be skewed towards on image reconstruction. That is given image I-deformed, find I. However, Mumford's interpretation of Pattern Theory is broader and he defines PT to include many more well-known statistical methods. Mumford's criteria for inclusion of a topic as Pattern Theory are those methods "characterized by common techniques and motivations", such as the HMM, EM algorithm, dynamic programming circle of ideas. Topics in this section will reflect Mumford's treatment of Pattern Theory. His principle of statistical Pattern Theory are the following:
Statistical PT makes ubiquitous use of conditional probability in the form of Bayes theorem and Markov Models. Both these concepts are used to express the relation between hidden states and observed states. Markov Models also captures the local properties of the stimulus, reminiscent of the purpose of bond table for regularity.
The generic set up is the following:
Let s = the hidden state of the data that we wish to know. i = observed image. Bayes theorem gives:
The following conditional probability examples illustrates these methods in action:

Conditional probability for local properties

N-gram Text Strings: See Mumford's Pattern Theory by Examples, Chapter 1.
MAP ~ MDL

Conditional probability for hidden-observed states

Bayes Theorem for Machine translation

Supposing we want to translate French sentences to English. Here, the hidden configurations are English sentences and the observed signal they generate are French sentences. Bayes theorem gives p'p = p = p'p and reduces to the fundamental equation of machine translation: maximize p = pp over the appropriate e. This reduces the problem to three main calculations for:
  1. p for any given e, using the N-gram method and dynamic programming
  2. p for any given e and f, using alignments and an expectation-maximization algorithm
  3. e that maximizes the product of 1 and 2, again, using dynamic programming
The analysis seems to be symmetric with respect to the two languages, and if we think can calculate p, why not turn the analysis around and calculate p directly? The reason is that during the calculation of p the asymmetric assumption is made that source sentence be well formed and we cannot make any such assumption about the target translation because we do not know what it will translate into.
We now focus on p in the three-part decomposition above. The other two parts, p and maximizing e, uses similar techniques as the N-gram model. Given a French-English translation from a large training data set :
NULL And the program has been implemented
Le programme a ete mis en application
the sentence pair can be encoded as an alignment that reads as follows: the first word in French comes from the second English word, the second word in French comes from the 3rd English word, and so forth. Although an alignment is a straightforward encoding of the translation, a more computationally convenient approach to an alignment is to break it down into four parameters:
  1. Fertility: the number of words in the French string that will be connected to it. E.g. n = probability that “implemented” translates into three words – the word's fertility
  2. Spuriousness: we introduce the artifact NULL as a word to represent the probability of tossing in a spurious French word. E.g. and its complement will be.
  3. Translation: the translated version of each word. E.g. t = translation probability that the English "has" translates into the French "a".
  4. Distortion: the actual positions in the French string that these words will occupy. E.g. d = distortion of second position French word moving into the fifth position English word for a four-word English sentence and a six-word French sentence. We encode the alignments this way to easily represent and extract priors from our training data and the new formula becomes
For the sake of simplicity in demonstrating an EM algorithm, we shall go through a simple calculation involving only translation probabilities t, but needless to say that it the method applies to all parameters in their full glory. Consider the simplified case without the NULL word where every word has fertility 1 and there are no distortion probabilities. Our training data corpus will contain two-sentence pairs: bcxy and by. The translation of a two-word English sentence “b c” into the French sentence “x y” has two possible alignments, and including the one-sentence words, the alignments are:
b c b c b
| | x |
x y x y y
called Parallel, Crossed, and Singleton respectively.
To illustrate an EM algorithm, first set the desired parameter uniformly, that is
Then EM iterates as follows
The alignment probability for the “crossing alignment” got a boost from the second sentence pair b/y. That further solidified t, but as a side effect also boosted t, because x connects to c in that same “crossing alignment.” The effect of boosting t necessarily means downgrading t because they sum to one. So, even though y and c co-occur, analysis reveals that they are not translations of each other. With real data, EM also is subject to the usual local extremum traps.

HMMs for speech recognition

For decades, speech recognition seemed to hit an impasse as scientists sought descriptive and analytic solution. The sound wave p below is produced by speaking the word “ski”.
Its four distinct segments has very different characteristics. One can choose from many levels of generators : the intention of the speaker's brain, the state of the mouth and vocal cords, or the ‘phones’ themselves. Phones are the generator of choice to be inferred and it encodes the word in a noisy, highly variable way. Early work on speech recognition attempted to make this inference deterministically using logical rules based on binary features extracted from p. For instance, the table below shows some of the features used to distinguish English consonants.
In real situations, the signal is further complicated by background noises such as cars driving by or artifacts such as a cough in mid sentence. The deterministic rule-based approach failed and the state of the art is to use a family of precisely tuned HMMs and Bayesian MAP estimators to do better. Similar stories played out in vision, and other stimulus categories.
ptkbdgmnfsvz
Continuant
Voiced
Nasal
Labial
Coronal
Anterior
Strident