Abstract:
A confusability tool generates a confusability cost associated with two phonemic transcriptions. The confusability cost measures the likelihood that a human or machine hearing the first word will mistakenly hear the second word. The cost calculation is based on a weighting of the Levinstein distance between the transcription pair.

Description:
BACKGROUND OF THE INVENTION 
     The present invention relates generally to speech processing and recognition, and specifically, to methods and systems of predicting the likelihood of a human or machine confusing two different words or phrases. 
     Words that sound alike tend to be confused with one another. For example, a listener hearing the word &#34;brush&#34; may mistakingly believe she is hearing &#34;rush,&#34; while the same listener hearing the word &#34;falcon&#34; is not likely to believe that she is hearing the word &#34;peckham.&#34; 
     Knowing that two words are likely to be confused can be useful in a number of applications. For example, two applications that may benefit from knowing the likelihood of confusion of two words are: speech recognition algorithms, which attempt to associate a spoken word with its intended orthography (i.e., the English written version of the word), and vocabulary rejection algorithms, which determine whether a spoken word is present in a predefined dictionary. 
     One method to determine the likelihood of confusion between spoken words is to simply have a human expert rate the potential confusability of the words based on intuition and experience. However, this method is laborious and results in subjective and inconsistent confusability measurements. 
     Thus, there is a need to be able to automatically generate, from a written version of a word, an objective metric of the likelihood of confusing a spoken word with another spoken word. 
     SUMMARY OF THE INVENTION 
     It is an object of the present invention to automatically generate an objective metric of the likelihood of confusing two words. 
     A method consistent with the present invention determines the likelihood of confusing a first word for a second word. The method comprises the steps of (1) phonemically transcribing the first and second words into first and second transcriptions; (2) calculating a Levinstein distance between the first and second transcriptions as the number of edit operations required to transform the first transcription into the second transcription; (3) obtaining a phonemic transformation weight for each edit operation of the Levinstein distance; and (4) summing the weights to generate a value indicating the likelihood of confusion between the first and second words. 
     Additionally, a computer system consistent with the present invention comprises a memory, a processor, and a display. The memory stores a confusability tool program and a transcripted phonemic lexicon. The confusability tool program includes means for calculating a Levinstein distance between first and second phonemic transcriptions from the lexicon as the number of edit operations required to transform the first transcription into the second transcription, obtaining a phonemic transformation weight for each edit operation of the Levinstein distance, and summing the weights to generate a value indicating the likelihood of confusion between the first and second transcriptions. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate several embodiments of the invention and, together with the description, serve to explain the principles of the invention. In the drawings, 
     FIG. 1 is a block diagram of an exemplary computer system used to execute methods consistent with the present invention; 
     FIG. 2 is high level block diagram illustrating the functional relationships of software components consistent with the present invention; 
     FIG. 3 is a flowchart illustrating a method of confusability cost calculation for phonemic transcription pairs; and 
     FIG. 4 is a flow chart illustrating the use of the present invention in a vocabulary diagnostic tool. 
    
    
     DETAILED DESCRIPTION 
     Reference will now be made in detail to the embodiments of the invention, examples of which are illustrated in the accompanying drawings. 
     Generally, methods and systems consistent with the present invention generate, for each of multiple pairs of phonemic transcriptions, a numeric indication of the potential confusability of the spoken version of the phonemic transcriptions. The numeric indication is based on the Levinstein distance, which determines the number of edit operations required to transform a first phonemic transcription into a second. The number of edit operations can be roughly equated to the degree of similarity, or the confusability, of the two phonemic transcriptions. The final numeric indication of confusability is derived from the Levinstein distance by multiplying each required edit operation by a weight corresponding to the significance of the edit operation, and summing the resultant values. 
     A &#34;phonemic transcription&#34; of a word is a representation of the word comprising a series of phonemes, which are the basic units of sound in a given language that distinguish one word from another. For example, the initial sound in &#34;cat&#34; and &#34;kick&#34; may be represented by the phonemic symbol /k/ while the one in &#34;circus&#34; may be represented by the symbol /s/. Throughout this disclosure, a phonemic alphabet of 40 symbols is used, although other phonemic alphabets could equivalently be used. Further, backslashes will be used to distinguish a symbol as a phonemic one. In contrast to a phonemic transcription of a word, an orthographic transcription of the word refers to the typical spelling of the word. 
     Phonemic transcriptions are important because English spelling is not phonetic; that is, there is often no direct correspondence between the spelling and the sounds we actually produce. One sound can be represented by a number of different letters or combinations of letters. For example, the first sound in the words &#34;cat,&#34; &#34;kick,&#34; &#34;quick,&#34; and &#34;chemistry&#34; is the same. 
     FIG. 1 is a block diagram of an exemplary computer system used to execute methods consistent with the present invention. Computer system 102 includes a chassis 110, which holds the computer&#39;s main processor and main memory; an input device such as keyboard 112; a storage device such as floppy or hard disk drive 114; and a display such as monitor 116. Preferably, the methods described herein are executed by the processing section of computer system 102. Computer system 102 is shown connected to network 118, and may be operated directly by a user or through the network. 
     Many variations of computer system 102 are possible. For example: storage device 114 may additionally include storage media such as optical disks, and user input device 112, instead of or in addition to keyboard 112, may include any type of user input device such as: an electronic mouse, a trackball, a lightpen, a touch-sensitive pad, a digitizing tablet, or a joystick. 
     FIG. 2 is high level block diagram illustrating a functional relationship of the data operated on by computer system 102. As shown, a confusability tool 202 receives, as inputs, a phonemically transcribed lexicon 204 and a phonemic transformation weight matrix 206. Lexicon 204 may be, for example, a pre-selected corpus of orthographies that have been phonemically transcribed either manually or through a known automated transcription method. The phonemic transformation weight matrix 206, W(/x/, /y/), specifies a weight value for each pair of phoneme symbols (/x/, /y/), as well as the phoneme null pairs (/x/,` `) and (` `,/y/). 
     Based on phonemic lexicon 204 and weight matrix 206, confusability tool 202 generates an ordered set of predicted auditory confusions 208. The generated output set is preferably in the form {orthography 1, orthography 2, cost}, where the cost is a number between zero and one. A small value of cost indicates a high level of potential confusability and a large value of cost indicates a high level of distinctiveness. 
     In addition to the ordered set of predicted auditory confusions, the tool may output a global confusability parameter. This parameter measures the confusability of the entire input lexicon and is preferably generated by inverting the confusability cost for the most confusable pairs of orthographies and summing the inverted values. 
     FIG. 3 is a flowchart illustrating the operation of the confusability tool 202 for phonemic transcription pairs from lexicon 204. For each pair of phonemic transcriptions T 1 . and T 2 , the tool calculates a Levinstein distance and the edits needed to transform T 1  into T 2  (described in more detail below) and uses this information to calculate a weighted Levinstein distance (steps 302, 304, and 306). A weighted Levinstein distance is calculated for each possible Levinstein distance, and the final weighted Levinstein distance, D(T 1 , T 2 ), is defined as the minimum of the possible weighted Levinstein distances (step 308). The final cost value is a normalized version of D(T 1 , T 2 ) (step 310). The normalization is preferably performed by dividing the final weighted Levinstein distance by the length of T 1  plus T 2 . 
     The final weighted Levinstein distance, determined in step 308, can be represented in a single equation as: ##EQU1## Here, T 1  is the source phonemic transcription, T 2  is the target phonemic transcription, N is the number of transformations between T 1  and T 2 , i is the index representing the edit operations of the Levinstein distance set (i.e., the particular substitutions, deletions, or insertions), and /x/ and /y/ are the phonemes involved in the edit, and may be the null phoneme ` ` for an insertion or deletion. The function W is the phonemic transformation weight 206 associated with the phonemic pair (/x/, /y/). As described above with reference to step 308, the minimization indicated by &#34;min&#34; is over the possible phonemic transformations between T 1  and T 2 . 
     Some orthographies have more than one phonemic transcription. To calculate the minimum weighted Levinstein distance between a pair of orthographies having different transcriptions, steps 302, 304, 306, and 308 are repeated for each possible pairing of the phonemic transcriptions of the two orthographies. If the weighted Levinstein distance between a pair of orthographies is denoted D(O 1 , O 2 ), where O 1  and O 2  are the two distinct orthographies, then, in terms of equation (1), D(O 1 , O 2 ) can be written as: 
     
         min D(T.sub.1 (i), T.sub.2 (j))i=[1 . . . n]; j=[1 . . . m](2) 
    
     Here, the weighted Levinstein distance between the pair of orthographies is defined as the minimum weighted Levinstein distance between the n possible transcriptions of T 1  and m possible transcriptions of T 2 . 
     Calculation of Levinstein Distance 
     As discussed above, calculation of the weighted Levinstein distance requires calculation of the Levinstein distance between two phonemic transcriptions (step 304). Suitable methods are known for computing the Levinstein distance between two strings, and include, for example, a method using maximum substring matching and a method using the Viterbi search algorithm across a transcription graph. For purposes of explanation, an exemplary Levinstein distance calculation will now be described. 
     Calculating the Levinstein distance involves &#34;editing&#34; a first string to transform it into a second string. Three different types of edits are preferably used to perform the transformation: (1) substitution, (2) deletion, and (3) insertion. With these three edit operations, any word can be transformed into another. For example, the substitution transformation {`W`, `f`} may indicate that &#34;wish&#34; should be transformed into &#34;fish&#34;; the deletion transformation {`t`, ` `} may indicate that &#34;tent&#34; is to become &#34;ten&#34;; and the insertion transformation {` `, `h`} may indicate that &#34;tent&#34; is to become &#34;tenth.&#34; Each substitution, deletion, or insertion is defined to involve only a single character symbol. 
     The Levinstein distance is the number of substitutions, deletions, or insertions required to transform the first phonemic transcription, T 1 , into the second, T 2 . Multiple transformation sets may be possible between T 1  and T 2 . For example, &#34;fishing&#34; to &#34;fashion&#34; can be transformed as: (1) {`i`, `a`}, {` `, `o`}, and {`g`, ` `}; or as (2) {`i`, `a`}, {`n`, `o`}, and {`g`, `n`}. In this example, each transformation has a Levinstein distance of three. Note, however, that because W(` `, `o`) may have a value different than W(`n`, `o`), the weighted Levinstein distance of (1) and (2) may vary. Accordingly, between (1) and (2), the chosen weighted Levinstein distance will be the smaller of the two. 
     Computation of Phonemic Transformation Weights W i  (x, y) 
     Consistent with the present invention, either one of two methods is used to calculate the phonemic transformation weights. Preferably these weights are precalculated and stored in a phonemic transformation weight matrix. In the case of a forty symbol phonemic alphabet (plus a null phoneme), a 41×41 phonemic transformation weight matrix would result. 
     Computation of Phonemic Transformation Weights Based on Phonological Features 
     A first method usable with the present invention for calculating the phonemic transformation weight for a phoneme pair (/x/, /y/) is based on descriptive values assigned by an expert to each phoneme. The phonemic transformation weights are given by: ##EQU2## where N is the number of phonological feature types and α i  is the relative significance of the feature type i. The functions f(/x/) and f(/y/) are the feature values of feature type i for phonemes /x/ and /y/, and t(f(/x/), f(/y/)) is the transformation penalty between feature values f i  (/x/) and f i  (/y/). The meaning of feature values, feature types, and transformation cost will be further explained. 
     As mentioned previously, forty phonemes may be used to describe English language words. Each phoneme sound is described by seven feature types, each of which may take on a number of values. The feature types, along with their permissible values, are shown below in Table I. 
     
                       TABLE I______________________________________Feature Types and ValuesFeature Type       Feature Values______________________________________voicing     voiced, unvoicedplace       bilabial, labio-dental, dental, alveolar,       palatal, palato-alveolar, velar, labio-velar,       labio-palatal, uvular, pharyngeal, glottalmanner      stop, fricative, nasal, approximant, vowel,       lateral, trill, flap, click, ejective,       implosive, affricaterealization aspirated, unexploded, syllabic, murmured,       long, velarized, labialized, palatalized,       rhoticized, nasalized, pharyngealized, tensefrontness   front, center, backheight      high, semi-high, upper-middle, middle, lowroundness   rounded, unrounded______________________________________ 
    
     The phoneme /a/, for example, as in the word &#34;cat,&#34; is assigned the feature values: &#34;low,&#34; &#34;front,&#34; and &#34;unrounded&#34;. These feature values correspond to the feature types &#34;height,&#34; &#34;frontness,&#34; and &#34;roundness,&#34; respectively. The feature types &#34;voicing,&#34; &#34;place,&#34; &#34;realization,&#34; and &#34;manner&#34; are not specifically defined for /a/, and accordingly, /a/ is said to have a null feature value for these feature types. Table II lists the preferred 40 phonemes for English, and their corresponding feature values. 
     
                       TABLE II______________________________________Phoneme   Feature values______________________________________A         low back unroundeda         low front unroundedaj        low front unrounded/semi-high front unroundedaw        low front unrounded/semi-high back roundedb         voiced bilabial stopD         voiced dental fricatived         voiced alveolar stopdZ        voiced alveolar stop/voiced palato-alveolar fricativeE         middle front unroundede         upper-middle front unroundedf         voiceless labio-dental fricativeg         voiced velar stopg-        voiced velar nasal stoph         voiceless glottal fricativeI         semi-high front unroundedi         high front unroundedj         voiced palatal approximantk         voiceless velar stopl         voiced alveolar lateral approximantm         voiced bilabial nasal stopn         voiced alveolar nasal stopO         back middle roundedOj        back middle rounded/semi-high front unroundedo         back upper-middle roundedp         voiceless bilabial stopr         voiced alveolar approximants         voiceless palato-alveolar fricatives         voiceless alveolar fricativeT         voiceless dental fricativet         voiceless alveolar stoptS        voiceless alveolar stop/     voiceless palato-alveolar fricativeU         back semi-high roundedu         back high roundedv         voiced labio-dental fricativeW         voiceless labio-velar approximantw         voiced labio-velar approximantZ         voiced palato-alveolar fricativez         voiced alveolar fricative     back middle unrounded*         center unrounded______________________________________ 
    
     Each feature type relates loosely to the physiological action used to make the phoneme sound. &#34;Place&#34; refers to the position of the primary articulator (e.g., tongue blade, tongue dorsum, lips), &#34;manner&#34; refers to the articulatory method used to generate the sound, &#34;realization&#34; refers to allophonic attributes of the production of the sound, &#34;frontness&#34; and &#34;height&#34; refer to tongue tip position in the production of the vowels, and &#34;roundness&#34; refers to the lip shape during production of the vowels. 
     The feature types are not all equally significant in the production of a sound. For example, the significance of &#34;manner&#34; is stronger than that of &#34;rounding.&#34; The relative significance value, α i , is used to quantify the significance of the various feature types. The values of α i  were assigned by a speech expert and are shown in Table III. 
     
                       TABLE III______________________________________Feature Type Numeric Significance, α.sub.i______________________________________voicing      0.3place        0.7manner       0.8realization  0.25frontness    0.4height       0.6roundness    0.2______________________________________ 
    
     In addition to the relative significance value, α i , transformation penalties are used to calculate the phonemic transformation weights. The transformation penalty is a function of two feature values of a feature type, and loosely measures the acoustic consequences of changing from an articulation of one feature type to another. In a manner similar to the significance values, the transformation penalties were determined by an expert. Tables IV-IX illustrate the preferred transformation penalties for the feature values except the feature values for the &#34;realization&#34; feature type. 
     
                       TABLE IV______________________________________(Feature Type = Height)from/to high     semi-high upper-mid                              mid   low______________________________________high    0.0      0.2       0.3     0.4   0.6semi-high   0.2      0.0       0.2     0.3   0.5upper-mid   0.3      0.2       0.0     0.2   0.4mid     0.4      0.3       0.2     0.0   0.3low     0.6      0.5       0.4     0.3   0.0______________________________________ 
    
     
                                           TABLE V__________________________________________________________________________(Feature Type = Place)from/to    bila  l-d     dent        alve           p-a              pala                 vela                    l-y                       l-p                          uvul                             phar                                glot__________________________________________________________________________bila    0.0  0.05     0.1        0.3           0.35              0.25                 0.2                    0.15                       0.2                          0.55                             0.7                                0.8l-d 0.05  0.0     0.1        0.25           0.3              0.35                 0.2                    0.3                       0.25                          0.5                             0.7                                0.8dent    0.1  0.1     0.0        0.05           0.1              0.15                 0.25                    0.3                       0.35                          0.45                             0.6                                0.7alvc    0.3  0.2     0.05        0.0           0.05              0.15                 0.25                    0.3                       0.35                          0.4                             0.5                                0.6p-a 0.35  0.3     0.1        0.05           0.0              0.05                 0.15                    0.2                       0.1                          0.35                             0.4                                0.45pala    0.25  0.35     0.15        0.15           0.05              0.0                 0.05                    0.15                       0.1                          0.3                             0.4                                0.45vela    0.2  0.2     0.25        0.25           0.15              0.05                 0.0                    0.05                       0.1                          0.2                             0.4                                0.45l-v 0.15  0.25     0.3        0.2           0.2              0.15                 0.05                    0.0                       0.05                          0.3                             0.45                                0.5l-p 0.2  0.25     0.35        0.1           0.0              0.05                 0.1                    0.05                       0.0                          0.35                             0.45                                0.5uvul    0.55  0.5     0.45        0.4           0.35              0.3                 0.2                    0.25                       0.3                          0.0                             0.3                                0.4pbar    0.7  0.7     0.6        0.5           0.4              0.4                 0.4                    0.45                       0.45                          0.3                             0.0                                0.15glot    0.8  0.8     0.7        0.6           0.45              0.45                 0.45                    0.5                       0.5                          0.4                             0.15                                0.0__________________________________________________________________________ (Key: bila = bilabial; ld = labiodental; dent = dental; alve = alveolar; pa = palatoalveolar; pala = palatal; vela = velar; lv = labiovelar; lp = labiopalatal; uvul = uvular; phar = pharyngeal; glot = glottal). 
    
     
                       TABLE VI______________________________________(Feature Type = Roundness)from/to         rounded unrounded______________________________________rounded         0.00    1.0unrounded       1.0     0.0______________________________________ 
    
     
                       TABLE VII______________________________________(Feature Type = Frontness)from/to   front        center  back______________________________________front     0.0          0.2     0.4center    0.2          0.0     0.2back      0.4          0.2     0.0______________________________________ 
    
     
                                           TABLE VIII__________________________________________________________________________(Feature Type = Manner)from/to    stop  fric     nas        appr           vow              lat                 tril                    flap                       clic                          ejec                             impl                                affr__________________________________________________________________________stop    0.0  0.3     0.4        0.8           0.9              0.7                 0.6                    0.5                       0.7                          0.5                             0.5                                0.2fric    0.3  0.0     0.3        0.6           0.8              0.6                 0.5                    0.35                       0.8                          0.6                             0.6                                0.15nas 0.4  0.4     0.0        0.3           0.6              0.2                 0.4                    0.3                       0.8                          0.6                             0.6                                0.35appr    0.8  0.4     0.3        0.0           0.2              0.2                 0.4                    0.25                       0.8                          0.8                             0.8                                0.7vow 0.9  0.8     0.6        0.2           0.0              0.6                 0.7                    0.7                       0.9                          0.9                             0.9                                0.85lat 0.7  0.6     0.2        0.2           0.3              0.0                 0.35                    0.3                       0.8                          0.8                             0.8                                0.7tril    0.6  0.5     0.4        0.4           0.7              0.35                 0.0                    0.2                       0.8                          0.8                             0.8                                0.4flap    0.5  0.35     0.3        0.25           0.7              0.3                 0.2                    0.0                       0.8                          0.8                             0.8                                0.4clic    0.7  0.8     0.8        0.8           0.9              0.8                 0.8                    0.8                       0.0                          0.4                             0.5                                0.7ejec    0.5  0.6     0.6        0.9           0.9              0.8                 0.8                    0.8                       0.4                          0.0                             0.5                                0.3impl    0.5  0.6     0.6        0.9           0.9              0.8                 0.8                    0.8                       0.5                          0.5                             0.0                                0.4affr    0.2  0.15     0.35        0.7           0.85              0.7                 0.4                    0.4                       0.7                          0.3                             0.4                                0.0__________________________________________________________________________ (Key: stop = stop consonant; fric = fricative; nas = nasal consonant; app = approximant; vow = vowel; lat = lateral; tril = trill; flap = flap; cli = click; ejec = ejective; impl = implosive; affr = affricate) 
    
     
                       TABLE IX______________________________________(Feature Type = Voicing)from/to         voiced  voiceless______________________________________voiced          0.0     1.0voiceless       1.0     0.0______________________________________ 
    
     The &#34;realization&#34; feature type is preferably a sequence of binary features. The particular binary values for the sequence are zero when a &#34;realization&#34; feature value is transformed into the identical feature and one when the &#34;realization&#34; feature value is transformed into a second feature value. In more precise terms, let x be one of the realization attributes X, where X is the set {aspirated, unexploded, syllabic, murmured, long, velarized, labialized, palatalized, rhoticized, pharyngealized, tense}. If x changes from x to not-x (e.g., from aspirated to unexploded), the transformation penalty is one, while if x does not change (e.g., the feature value of both phonemes is aspirated), the transformation penalty is zero. 
     Using equation (3) and tables II-IX, the phonemic transformation weights can be calculated for any phoneme symbol pair (/x/, /y/). Additionally, insertion and deletion costs, W(/x/, &#34;&#34;) and W(&#34;&#34;, /y/), respectively, are determined by using an expert to estimate the phonemic transformation weights of each phoneme and the null value. As an example transformation weight calculation, consider the phoneme symbol pair (/b/) {voiced bilabial stop}, and /d/ {palato-alveolar}. Applying equation (3), and using the appropriate values from the tables, the feature weight for place is 0.7, and the transformation penalty from &#34;bilabial&#34; to &#34;pal-alv&#34; is 0.3, which leads to a phoneme transformation penalty of D(/b/,/d/)=0.21. Table X, below, lists additional exemplary phonemic transformation weight pairs calculated using phonological features. 
     
                       TABLE X______________________________________  Phoneme Pair           Weight______________________________________  /b/, /d/ 0.21  /b/, /g/ 0.18  /b/, /p/ 0.3  /b/, /n/ 0.61  /b/, /w/ 0.84______________________________________ 
    
     Although the values assigned to the phonemes, as shown above in Tables IV-IX, are preferred, other values may also be used. Indeed, the number of phonemes used to describe a language, the feature types, and the feature values were all picked based on the best judgement of language experts, but are, of course, not the only possible representation of a language. 
     Computation of Confusability Weights Based on Acoustic Models 
     A second method usable with the present invention for calculating the phonemic transformation weights uses the relative entropy between an acoustic model, such as the Hidden Markov Model (HMM). Relative entropy is first found by aligning the states in the HMMs, and then considering the dominant pair of components of the Gaussian mixtures. 
     Using HMMs to acoustically model phonemes is known. Generally, in a HMM model, phonemes are assigned multiple states. Each state is represented by two probability distributions: the state probability distribution function and the transition probability distribution function. The state probability represents the probability that the state is occupying a certain value, and the transition probability represents the probability that the stochastic process--a derivation of the speech signal--will leave the state. 
     Consistent with the present invention, each phoneme in the HMM is assigned an equal number of states. The states of the phonemes are aligned in a one-to-one relationship. 
     The state probability is a mixture of multi-variate Gaussian probability distributions, and is given by ##EQU3## where G(μ i , Σ i ) is a multi-variate Gaussian distribution with mean, μ i , and covariance matrix, Σ i , N is the number of distributions in the mixture, and η i  is the mixture weight. 
     The distance between HMM state x i  and state y i  of phoneme /x/ and /y/, respectively, is defined to be the relative entropy between mixtures P x (i) and P y (i), denoted H(x i , y i ). Methods for calculating H(x i , y i ) are known. One such description is given in the conference report &#34;State Clustering of HMMs with Continuous Mixture Densities,&#34; by P. Stubley and J. F. Crespo. IEEE Workshop on Automatic Speech Recognition, Snowbird, Colorado, 1995. 
     The phonemic transformation weight between phonemes x and y, in terms of H(x i , y i ), is: ##EQU4## Here, n is the number of HMM states used to model the phoneme (e.g., 3), and v i  is a state weight, set heuristically, with the following guiding principle: because the steady state behavior of a phoneme is usually represented by the central states in the acoustic model, these states are usually assigned a higher weight. Conversely, the initial and final states are usually assigned a smaller weight. Table XI, below, lists exemplary confusability weight pairs calculated using equation (5). In the table, three HMM states were used and the weights, v i , were set at {0.2, 0.6, 0.2}. Note that the specific value of the state weights V 1  are not restricted to the example weights given in this disclosure. 
     
                       TABLE XI______________________________________  Phoneme Pair           Weight______________________________________  /b/, /d/ 0.18  /b/, /g/ 0.21  /b/, /p/ 0.25  /b/, /n/ 0.54  /b/, /w/ 0.56______________________________________ 
    
     Application of the Confusability Tool: Vocabulary Diagnostic Tool 
     Computer recognition of human speech is valuable in many applications requiring human/computer interaction. In designing such a speech recognition system, there is some flexibility in the selection of the lexicon the computer is trained to recognize. The confusability tool consistent with the present invention may be used to automatically find word pairings that tend to cause high confusability. The designer can then replace the relevant orthographies with alternate synonyms. In general, simplification of a lexicon by replacing confusable words with non-confusable synonyms can be useful in many different communication areas, such as simplification of a lexicon to facilitate understanding across a communication medium. Alternatively, the creation of a highly confusable lexicon may also be valuable, in, for example, creating a confusable test lexicon to rigorously test a speech recognizer. 
     FIG. 4 is a flow chart illustrating the use of confusability tool 202 as a vocabulary diagnostic aide. The proposed lexicon of the speech recognizer is phonemically transcribed by a human expert or by an automated transcription tool (step 402) For each pair of orthographies, the weighted Levinstein distance is computed pursuant to equation (1) (step 404). Orthographies may have more than one transcription. For these, the minimum normalized weighted Levinstein distance is computed, pursuant to equation (2), to obtain D(O 1 , O 2 ) (step 406). The vocabulary diagnostic tool then apples a predefined threshold value to the normalized minimum values of each orthography pair, D(O 1 , O 2 ), of each triplet {O 1 , O 2 , D(O 1 , O 2 )}(step 408). Finally, the remaining triplets are sorted in ascending order by weight and presented to the designer (step 410). 
     Accordingly, the designer is automatically presented with a list of the orthographies in the proposed lexicon, sorted so that those orthographies most likely to be confused are presented first. Based on the list, the designer may modify the orthographies in the lexicon, and, if he desires, re-runs the vocabulary diagnostic tool (step 412). 
     A confusability tool for objectively and automatically generating confusability ratings between pairs of orthographies or phonemic transcriptions has been described. The confusability tool has many potential applications, including use in a vocabulary diagnostic tool and in a speech recognition algorithm. Other uses are possible. For example, the confusability tool can be used to help speakers acquire a second language by indicating word pairings that the speaker should spend extra time learning to pronounce because of their potential confusion, such as the words &#34;tin&#34; and &#34;teen&#34; for speakers whose first language does not differentiate between the sounds /i/ and /l/. As a second example, knowledge of potential confusable words could be used in the speech pathology field to help speakers differentiate or pronounce difficult word pairings. 
     It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the scope or spirit of the invention. For example, many-to-many phoneme transformations can be used instead of one-to-one transformations (e.g., T 1  and T 2  in D(T 1 , T 2 ) may include pairs such as /xy/ and /ab/, where /x/, /y/, /a/, and /b/ are all individual phonemes. Other possible modification include: (1) combining the weight and feature transformations using a non-linear combination; (2) the transformation weights may take into account bundles of features rather than a single feature at a time (e.g., mapping from {voiced stop} to {voiceless fricative} pairs rather than cumulating for {voiced} to {voiceless} and {stop} to {fricative}; (3) the table of 41×41 phoneme transformation penalties can be assigned heuristically by an expert; (4) and the table of 41×41 phoneme penalties can be assigned using results from a speech recognition experiment. Other embodiments of the invention will be apparent to those skilled in the art from consideration of the specification and practice of the invention disclosed herein. It is intended that the specification and examples be considered as exemplary only, with the true scope and spirit of the invention being indicated by the following claims.