Data mining of audiology patient records: factors influencing the choice of hearing aid type

Background This paper describes the analysis of a database of over 180,000 patient records, collected from over 23,000 patients, by the hearing aid clinic at James Cook University Hospital in Middlesbrough, UK. These records consist of audiograms (graphs of the faintest sounds audible to the patient at six different pitches), categorical data (such as age, gender, diagnosis and hearing aid type) and brief free text notes made by the technicians. This data is mined to determine which factors contribute to the decision to fit a BTE (worn behind the ear) hearing aid as opposed to an ITE (worn in the ear) hearing aid. Methods From PCA (principal component analysis) four main audiogram types are determined, and are related to the type of hearing aid chosen. The effects of age, gender, diagnosis, masker, mould and individual audiogram frequencies are combined into a single model by means of logistic regression. Some significant keywords are also discovered in the free text fields by using the chi-squared (χ2) test, which can also be used in the model. The final model can act a decision support tool to help decide whether an individual patient should be offered a BTE or an ITE hearing aid. Results The final model was tested using 5-fold cross validation, and was able to replicate the decisions of audiologists whether to fit an ITE or a BTE hearing aid with precision in the range 0.79 to 0.87. Conclusions A decision support system was produced to predict the type of hearing aid which should be prescribed, with an explanation facility explaining how that decision was arrived at. This system should prove useful in providing a "second opinion" for audiologists.


Background
This research looks for factors influencing the choice between two common hearing aid types: BTE (worn behind the ear) or ITE (worn in the ear). This choice is typically made by audiology technicians working in outpatient clinics, on the basis of audiogram results and consultation with the patient. In many cases, the choice is clear cut, but at other times the technicians might benefit from a second opinion given by an automatic system with an explanation of how that second opinion was arrived at. The production of such a decision support system is the main goal of this paper. Our data set is unusual in that ITE hearing aids are not generally available on the British National Health Service in England, as they are more expensive than BTE hearing aids. However, both types of aid are prescribed at James Cook University Hospital in Middlesbrough, UK. The data, collected between 1992 and 2001, consists of the following types of records: • Audiograms (graphs of the auditory thresholds, or faintest sounds audible to the patient at six different pitches or frequencies, where 0 shows perfect hearing and higher thresholds show impaired hearing), e.g., 40,35,35,35,85,70,15,20,20,30,55, where the first six values are AC (air conduction) and the last five are for BC (bone conduction). AC is measured by placing headphones over the ears, and determines the overall level of hearing. BC is measured by placing the sound source tightly on the mastoid bone behind the ear, and measures the level of hearing of the inner part of the ear. A constraint on the data is that BC must always be the same or better than AC. The difference between the AC and the BC is called the airbone gap, and measures the hearing ability of the middle and outer parts of the ear.
• Brief free text notes made by the technicians, e.g., IMPS. TAKEN FOR BINAURAL AIDS., where IMPS is an abbreviation for "impressions", and BINAURAL means "worn in both ears".

Principal component analysis on audiograms
This section describes how a Principal Component Analysis (PCA) was performed on the set of 11,462 audiograms where all AC and BC thresholds for the right ear were recorded, to determine the main audiogram types found among hearing aid users. The rows of our input matrix were the individual audiograms for the right ear, while the 11 columns were for six air conduction and five bone conduction thresholds. Although the patients were originally tested at 11 frequencies, the principle of PCA is that certain frequencies tend to vary together, and thus can be grouped into a smaller number of underlying variables called principal components (PC). Each PC has a set of coefficients in the range -1 to +1, corresponding to the degree of influence of each of the original thresholds on that PC, given in Table 1.
A method was devised for converting PCA scores to typical audiogram types, which is an approximation only. It assumes that the PCA scores are directly related to audiogram thresholds, whereas in reality they measure different things: PCA measures the importance of a threshold in distinguishing audiograms, while audiograms measure the degree of hearing loss at the same frequency. The range of PCA scores is -1 to +1, while the range of audiogram thresholds is 0 to 120 dB. If we assume that the relationship is linear, then for each frequency a PCA value of -1 corresponds to an audiogram threshold of 0, and a PCA value of 0 = an audiogram threshold of 60, and a PCA value 1 = an audiogram threshold of 120. The formula relating the two was: Audiogram threshold = 60 + (60 × PCA score). (1) Each individual patient audiogram was classified into one of the main audiogram types identified, according to least Euclidean distance. A chi-squared test was then performed to determine whether there was any association between the audiogram class of each patient and the type of hearing aid worn. This test was done on the set of 7,437 records where all AC and BC thresholds were available for the right ear, and either a BTE or an ITE aid was specified. In the final logistic regression model, rather than simply using the identified broad audiogram types, each individual hearing threshold was used.
Use of the chi-squared test to discover other factors related to hearing aid type In the previous section, it was shown that the choice of hearing aid type was related to the shape of the audiogram. This section describes how the simple chi-squared test was used to discover which of the category data fields were significantly associated with the choice of hearing aid type, and also to discover free-text keywords which were significantly associated with either BTE or ITE hearing aids.
To discover those free-text keywords which were significantly associated with either BTE or ITE hearing aids, first a large contingency table was created where the rows stood for hearing aid type, while each column stood for a candidate keyword (one of 664 distinct words found to have occurred at least once in the free-text fields of the entire record set). The observed value in each cell was the number of times that word had been found in the free-text fields of patients with that type of hearing aid. This table (Table 2) had an overall chi-squared value of 5421.84 for 663 degrees of freedom, giving p < 0.001. This data showed, with 99.9% confidence that these free text words were not randomly distributed, but some text  words are associated with hearing aid type. To find the association between individual free-text words and hearing aid type, the quantity (O -E) 2 / E was examined to rank the keywords according to importance. Some of the keywords in Table 3, were stemmed forms such as 'reshel' for 'reshell' and 'tinnitu' for 'tinnitus', since all the text was passed through Porter's stemmer [1] for the removal of grammatical endings. Some significant keywords were abbreviations, such as IMP for 'impression'. Dictionaries were not used to group terms with the same meaning (synonyms), instead using only the surface form of words because the procedure of keyword selection to was made as automatic as possible. It was assumed that the text notes used in the database were homogeneous, where the various technicians were consistent in their use of terminology and abbreviations. Although it was not done, consistency in writing conventions could have been verified manually, since the identity of the technician is given in one of the record fields. Similarly, it was assumed that two technicians treating the same patient would produce identical audiograms and choose the same hearing aid type. To do otherwise would have resulted in excessive subdivision of the data set.

Logistic regression (LR) model for ITE/BTE right ear hearing aids
Having determined that audiogram frequencies, gender, presence of tinnitus, use of tinnitus masker, age, mould and certain keywords were all associated with the decision between fitting a BTE or an ITE aid, all these factors were combined into a single logistic regression model. The input to the model was the matrix of patient data, where columns corresponded to attributes and rows corresponded to individual patient records. The output was a formula [2] in the form: In our case p is the probability that the patient should be fitted with an ITE aid, while (1-p) is the probability that the patient should be given a BTE aid. b 0 is a constant, and b 1 to b k are called the coefficients of the model. The values x 1 to x k are all either 1 or 0, depending on whether a given attribute in the patient's record is present or absent. The overall value L is greater than 0 if it is more likely that the patient should be given a BTE aid, while it is less than 0 if it is more likely that the patient should be given an ITE aid.
Before performing the logistic regression, the actual set of records was divided into two parts, one containing 80 percent (5,736) of the records and the other containing the remaining 20 percent (1,433) of these records. The 80 percent subset was used as the training set for model construction, and the remaining records for testing the model. The sampling method was to extract every fifth record for testing. The logistic regression was performed on the records which had all fields filled for the right ear: AC (air conduction) and BC (bone conduction) thresholds, gender, age and text keywords (5,736 records), of which 128 also had non-null entries for diagnosis, 98 had non-null entries for masker, and 3983 had non-null entries for mould. This data was converted into discrete numeric values as inputs to the model, as follows: For AC thresholds below the first quartile (40 db) a value of 0; for thresholds between the first and second quartile (55 dB) a value of 1; for thresholds between the second and third quartile (75 dB) a value of 2; and for thresholds above the third quartile a value of 3. The same method was used to assign values for BC thresholds and age, except in that the quartile thresholds were 25, 40 and 55 dB and 60, 70 and 78 years respectively. The values for diagnoses were 0 for no tinnitus diagnosis and 1 for tinnitus. Finally, for gender, 0 was assigned for male and 1 for female. Regression coefficients and associated p values were found for all the model variables and those variables with p values more than 0.05 for the constant were discarded. Thus, BC4000 (bone conduction at 4000 Hz), age, diagnosis and masker were not considered in the final model as the p values of their constants was more than 0.05, as shown in Tables 4 to 7. Due to data sparseness, it was not possible to incorporate all the keywords discovered by an analysis of the free text into a single model. However, a few significant keywords (meaningful words producing the highest chisquared values, shown in Table 8), were included into the final model, along with all the categorical and numeric fields. These keywords were all acronyms where APPT stands for appointment, FTA for first time appointment, GP for general practitioner, MAP for processor amplification map (associated with cochlear

Principal component analysis (PCA)
The coefficients of the first PC (PC1) were all negative and approximately equal. This suggests that the main source of variation between the patients was simply the overall degree of hearing loss. The coefficients of the second PC (PC2) were negative for frequencies at or below 1000 Hz, but positive for higher frequencies, for both air and bone conduction, and thus differentiated patients according to whether they have a predominantly high frequency or low frequency hearing loss. The coefficients of the third PC (PC3) were positive for air conduction at all frequencies, but negative for bone conduction, showing a contrast between patients with and without an air-bone gap. The fourth component (PC4) was similar to the third, but corresponded to an air-bone gap at low frequencies. No clear patterns were seen for the fifth or subsequent principal components. The first four PCs corresponded to audiogram types frequently encountered in audiology clinics. The percentage of the overall variability in the data explained by the first four principal components respectively was 59.5, 13.4, 9.7, and 5.2, giving a total of 87.8%.
The thresholds corresponding to the first four PCs by using formula (1) are given in Table 9. The top row of each PC refers to the air conduction frequencies at 250, 500, 1000, 2000, 4000 and 8000 Hz respectively, while the second row refers to the corresponding bone conduction frequencies in the range 250 to 4000 Hz. The results are shown in Tables 10 to 12. The overall χ 2 (chi-squared value), calculated as the sum of the cells in Table 12, is 548.07, which for one degree of freedom gives p < 0.001, so audiogram type is clearly related to hearing aid type. Also, in Table 12, the (O -E) 2 / E values which make the greatest contribution to the overall χ 2 value are those in the PCA1 and PCA3 columns. Thus, flat hearing loss (PCA1) audiograms without air-bone were associated with ITE aids and flat audiograms with additional airbone gaps (PCA3) were associated with BTE aids. This result is in accordance with Stephens [3], who found that the fitting of ITE aids was limited in cases of severe hearing loss. Audiograms were also clustered using K-means clustering [4] which produced similar results to PCA, in that the mild to moderate hearing loss cluster was associated with ITE aids and the severe hearing loss cluster was associated with BTE aids. Thus, it was demonstrated that the audiogram is a factor influencing the choice of hearing aid type.

Chi-squared test
The contingency table showing the relationship between gender and hearing aid type is shown in Table 12. The raw counts are given at the top of each cell, where for example there were 3196 male patients who wore BTE hearing aids. In each cell the Observed frequencies (O) are not enclosed in brackets, Expected frequencies (E) are in () and the quantity (O -E) 2 / E is in []. The overall chisquared value (the sum of the values in [] for all four cells) was 33.68, which for one degree of freedom is significant at p < 0.001. Males tended more to use ITE hearing aids and females tended more to use BTE hearing aids. For the relationship between hearing aid type and a diagnosis of     tinnitus (ringing in the ear), the overall chi-squared value was 31.75, again significant at p < 0.001 for one degree of freedom. Patients with tinnitus tended more to wear ITE hearing aids. The relationship between the wearing of a tinnitus masker (a soothing sound source designed to drown out tinnitus) and hearing aid type, among patients diagnosed with tinnitus, had the overall chi-squared value of 17.16, which for one degree of freedom, was also significant at p < 0.001. The data for the cross-tabulation of hearing aid type and age produced the overall chi-squared value of 10.53, which for one degree of freedom, showed significance at p < 0.001. Mould type was also cross-tabulated with hearing aid type and the overall chi-squared value was 9844.18, which for 30 degrees of freedom was significant at p < 0.001. Thus all the category data types were significantly associated with hearing aid type. All the data in the patient records was used without considering confounding effects, where for example it might have been the choice of hearing aid type affecting the choice of mould, rather than vice versa. It is believed that this may have been the case, since many mould types never occurred in conjunction with one or the other hearing aid type. The set of free-text keywords which tended to occur significantly more and less often (called positive and negative keywords respectively) in records where the patient wore either BTE or ITE aids are shown in Table  3. The association between these keywords and one or other type of hearing aid suggests the following: BTE aids were associated with high gain (amplification), e.g., be34, be36 and be52, and cases where changes had been made to the ear mould. ITE hearing aid types tended to use lacquer, had vents, required reshelling of ear impressions, had changes made to the hearing aid itself, were reviewed and the wearers were making progress.

Logistic regression (LR) model
In Table 8, the part of the model which takes into account the occurrence or otherwise of the selected keywords in deciding which type of hearing aid to suggest is given. Using keywords alone, the relative likelihoods of the patient needing an ITE or a BTE aid are given by equation (3): (3) As shown in Table 13, the logistic regression coefficient for gender was calculated. In Table 14, logistic regression values for air conduction (AC) at 250 dB are given for each quartile and in Table 15, predicted log odds are calculated using the regression coefficient values (b) from Table 14. Similarly, predicted log odds were calculated for AC500, AC1000, AC2000, AC4000, AC8000, BC250, BC500, BC1000, BC2000, gender and mould.
To show how the model works, a sample record from the test set is taken, as shown in Table 16. The attributes of this record are shown in the first column, and their values are shown in the second. Starting with a predicted log odds of 0 (meaning a BTE and ITE aid are assumed equally likely), the values of the attributes are examined in the record one by one, and add on the regression coefficient corresponding to that value of that attribute to the running total. The value of the first attribute, age is disregarded, but the next attribute gender has value "male". For "male", the regression coefficient is calculated to be -0.23, so the running total becomes 0 -0.23 = -0.23. The relevant regression coefficients for each attribute are added in turn, ending with adding on -0.04 for the presence of the keyword "REV"   in the free-text field. The final total of the regression coefficients is 10.1, which is the final log odds value, suggesting that it is much more likely that this patient would benefit most from a BTE hearing aid as opposed to an ITE aid. Testing of these logistic regression models showed that overall there was 81.64% agreement between the predictions of our model and the actual hearing aid chosen by the audiologist (as given in the "type" field) as shown in Table 17. The agreement rate was higher for patients fitted with ITE aids (86%) than for those fitted with BTE aids (76%). The results were analyzed according to precision, recall and F-measures [5] using equations (5), (7) and (8) respectively, as shown in Table 18. For comparison, a similar analysis using a Naïve Bayesian approach was performed, and 0.67, 0.76, and 0.71 were obtained as the precision, recall and F-score respectively for ITE and 0.66, 0.56 and 0.60 as the precision, recall and F-score respectively for BTE.
In Table 19, 782 and 651 are the counts of ITE and BTE aids respectively in the human-annotated test data, while 833 and 600 are the counts of ITE and BTE aids respectively in the machine predicted results. The overall agreement is much better than random (50%), but the performance of a classifier should also be compared against the "simplest possible algorithm" [6]. In our case, this would be to assume that all the patients should be assigned the more commonly prescribed type of hearing aid. In our test set 782 out of 1433 patients in the test set were given ITE aids, so simply assigning all the patients this type of aid would provide 54.6% agreement, which is referred to as the ZeroR baseline.
The theoretical upper bound of classifier performance is the inter-annotator agreement [2], in our case the rate at which two expert audiologists would assign the same hearing aid to the same patient. Unfortunately, we do not have data on this. Five-fold cross validation (repeated subsampling of the data to produce five nonoverlapping test sets for an unbiased estimation of model accuracy) was performed. The overall similarity was in the range 82 to 85%, precision was in the ranges 0.79 to 0.87 for ITE and 0.82 to 0.85 for BTE, recall was 0.84 to 0.88 for ITE and 0.74 to 0.85 for BTE, and the F measure was 0.83 to 0.86 for ITE and 0.79 to 0.83 for BTE. For most of the cross-validation runs BC2000, BC4000, Diagnosis and Masker were discarded from the model, since these variables have p values of more than 0.05 for their constants. These results show that for each run, both the final model and the success rates were similar.

Discussion
Although this LR model did not find age as a significant factor, Meredith and Stephens [7] have found that the ITE hearing aid presents handling problems only in subjects over 75 years of age. Dillon [8] also found that BTE aids are easier to operate as they are larger in size, and thus would be more popular with older people. The literature shows men and women preferring the two types for different reasons. Martin, et al. [9] found that more males choose ITE aids than females, because they perceive them to be a more advanced technology -though in reality the same makes and specifications are available in both styles, and neither model is more advanced than the other. They also found more females reporting that ITE aids are easier to handle than BTE. Mueller, et al. [10] found no difference in how embarrassed males and females feel about using a BTE aid. This LR model did not include diagnosis (as mentioned above for Table 6), although the authors previously found [11] that there was a significant association between the choice of BTE hearing aid and a diagnosis other than tinnitus (ringing in the ear), by using the chi-squared test. We also found, by using the chi-squared test that BTE hearing aids were atypical of tinnitus-withmasker. Other factors mentioned in the literature which could not be tested with this data are the greater cosmetic acceptability of the smaller ITE aids, comfort in wear, ease of use with spectacles, and sound quality [12].

Conclusions
The associations between hearing aid type and audiogram type were confirmed by both the PCA/chi-squared and LR experiments described in this paper, and also by the authors' previous work on associations between words found in the database and hearing aid type, and the previous findings by audiologists [3]. These approaches will form the basis for an audiology decision support system, where unseen patient records would be presented to the system, and the relative probability that the patient should be fitted with an ITE aid as opposed to a BTE hearing aid would be returned. The advantage of these techniques for the combination of evidence is that it is easy to see which variables contributed to the final decision. It is planned to validate these results by obtaining feedback from a professional audiologist, and by using an approach (Bayesian networks) which constructs model with interaction between variables. A major advantage of both Naïve Bayes and Logistic Regression is that they enable an explanation facility to be incorporated into any decision support tool, since it is easy to read back and