AN EXPLORATION OF DIABETIC FOOT SCREENING PROCEDURES DATA BY A MULTIPLE CORRESPONDENCE ANALYSIS Vilma URBANČIČ ROVAN1*, Jože ROVAN2 1University Medical Centre Ljubljana, Department of Endocrinology, Zaloška 7, 1000 Ljubljana, Slovenia 2University of Ljubljana, Faculty of Economics, Deptartment of Statistics, Kardeljeva ploščad 17, 1000 Ljubljana, Slovenia Received: Feb 16, 2016 Accepted: Aug 16, 2016 Original scientific article Aims. Gangrene and amputation are among most feared complications of diabetes mellitus. Early detection of patients at high risk for foot ulceration can prevent foot complications. Regular foot screening (medical history, foot examination and classification into risk groups) was introduced at the out-patient diabetes clinic in Ljubljana in November 1996. We aimed to explore the relationships between the observed variables, check the appropriateness of the risk status classification and of the post-screening decisions. Methods. The data of 11.594 patients, obtained in 18 years, were analysed by multiple correspondence analysis (MCA). Most of the observed variables were categorical. Results. The majority of the screened population was free of foot complications. We demonstrated an increasing frequency and severity of foot problems with an increasing age, as well as the association between the loss of protective sensation and the history of foot ulceration, foot deformity and callus formation, the history of foot ulcer or amputation and acute foot ulceration. A new finding was that the location of foot deformity points was closer to female than male gender, indicating the possible role of fashionable high-heel footwear. The appropriateness of therapeutic decisions was confirmed: the points representing absent foot pulses and referral to vascular specialist were close together, as well as points representing foot deformity and special footwear prescription or callus formation and referral to pedicurist. Conclusions. MCA was applied to the data on foot pathology in the population attending the out-patient diabetes clinic. The method proved to be a useful statistical tool for analysing the data of screening procedures. Namen. Gangrena in amputacija sodita med najhujše zaplete sladkorne bolezni. Zaplete na nogah lahko preprečimo z zgodnjim odkrivanjem ogroženih bolnikov. V Diabetološki ambulanti Ljubljana izvajamo presejalni test za diabetično stopalo od novembra 1996. Test obsega anamnezo, klinični pregled in klasifikacijo glede na ogroženost. Želeli smo raziskati povezanost med opazovanimi spremenljivkami, preveriti pravilnost klasifikacije bolnikov in ustreznost odločitev po pregledu. Metode. Podatke 11.594 bolnikov, dobljene v 18 letih, smo analizirali z multiplo korespondenčno analizo (MCA). Večina opazovanih spremenljivk je bila kategoričnih. Rezultati. Večina opazovane populacije ni imela zapletov na nogah. Dokazali smo, da pogostnost in izraženost težav z nogami narašča s starostjo ter da obstajajo povezave med izgubo zaščitne občutljivosti in razjedo v anamnezi, med deformacijo nog in tvorbo kalusa, med ulkusom ali amputacijo v anamnezi in akutno razjedo. Nova ugotovitev je lokacija točke, ki označuje deformacijo stopala, bliže ženskemu kot moškemu spolu, kar kaže na možen vpliv modnih čevljev z visokimi petami. Potrdili smo ustreznost terapevtskih odločitev: točke, ki označujejo odsotne stopalne pulze, in tiste, ki označujejo napotitev k angiologu, so bile blizu skupaj, prav tako deformacija nog in predpis posebne obutve ali tvorba kalusa in napotitev k pedikerju. Zaključek. MCA smo uporabili za analizo podatkov o patologiji stopal pri populaciji bolnikov iz naše diabetološke ambulante. Metoda se je izkazala kot uporabno statistično orodje za analizo podatkov o presejalnih testih. ABSTRACT Keywords: diabetic foot, screening, statistics IZVLEČEK Ključne besede: diabetično stopalo, presejanje, statistika *Corresponding author: Tel: ++ 386 41 792 771; E-mail: vilma.urbancic@kclj.si 10.1515/sjph-2017-0009 Zdr Varst 2017; 56(1): 65-73 65 ANALIZA PODATKOV PRESEJALNEGA TESTA ZA DIABETIČNO NOGO Z MULTIPLO KORESPONDENČNO ANALIZO Urbančič Rovan V, Rovan J. An exploration of diabetic foot screening procedures data by a multiple correspondence analysis. Zdr Varst 2017; 56(1): 65-73 © Nacionalni inštitut za javno zdravje, Slovenija. 1 INTRODUCTION Diabetes mellitus is a chronic metabolic disorder affecting 8-11% of the general population (1). It is one of the major public health problems worldwide, including Slovenia (2). The disease is characterized by elevated blood sugar levels, which damage blood vessel walls, ultimately leading to end-organ injury. The financial burden of diabetes treatment and its complications is enormous (3). Gangrene and amputation are among the most feared complications of diabetes. More than 50% of all non-traumatic amputations are performed on diabetic patients (4, 5). The initial event is usually a foot ulcer, which affects about 15% of all patients with diabetes mellitus. The major risk factors for diabetic foot ulceration and gangrene are neuropathy, impaired blood supply, trauma and infection. Other contributing factors include a history of previous foot ulcer, improper footwear, low-quality podiatry service, poor metabolic control, psychological factors, tobacco smoking, old age and low social status. Foot complications usually develop due to an interplay of several component causes (4, 6, 7). Diabetic neuropathy can involve sensory, motor and autonomic nerves. The most prevalent form is distal symmetric sensory neuropathy (8, 9, 10), which affects vibration, light touch, pain, temperature and proprioceptive sensation. Motor neuropathy leads to muscle atrophy, foot deformity (claw-hammer toes), altered biomechanics of walking, and redistribution of foot pressures during standing and walking. Abundant callus formation on pressure points, together with thinning of the submetatarsal head fat pads, additionally increases the magnitude of plantar pressure, and ultimately results in foot ulceration (11, 12). Neuropathy of sympathetic sudomotoric nerves results in diminished or even absent sweating, which causes dry foot skin, prone to cracks and fissures. The distribution of neuropathic changes is usually symmetrical (9, 10). Due to impaired or absent pain sensation, patients with sensory neuropathy often neglect warning signs of foot disease. Early detection of high-risk patients by systematic screening enables timely implementation of preventative measures to avoid gangrene and amputation (13). Screening means examination of asymptomatic people to detect those with a high probability of having a given disease (14). In the past, there used to be mass screening for lung tuberculosis; nowadays, all blood products are screened for the presence of HIV. Foot screening protocol is a set of simple and inexpensive procedures, which reveal the patients at risk for the development of foot ulceration. It consists of one’s medical history (previous foot ulceration, symptoms of neuropathy), foot examination and classification into risk groups. In Slovenia, there is no national diabetes registry. In 2014, there were 104.550 patients on antidiabetic medications, with the estimated prevalence of diabetes mellitus in the population over 15 years of age being 6.9% (confidence interval 6.3 - 7.6 %) (15). The out-patient diabetes unit at the University Medical Centre Ljubljana is the biggest in Slovenia (with over 11.000 registered patients). Here, we started to perform foot screening in November 1996. The procedure is repeated in every patient at least once a year. In this retrospective survey, we decided to analyse the results of the first screening test. The main analytical tool was multiple correspondence analysis (MCA). The aims of the study were: 1. To reveal the possible new causal relationships among the observed variables and to confirm the causal pathways to foot complications which have already been described before. 2. To evaluate the appropriateness of the classification into risk groups. 3. To evaluate the appropriateness of the post-screening decisions (referral, footwear prescription). 2 METHODS 2.1 Patients and Data Collection Between November 1996 and December 2014, 11.594 patients with diabetes mellitus were screened (55.2% men, average age 60.81 years): all patients registered at the clinic in November 1996, and afterwards, all newly-diagnosed patients with type 2 diabetes, as well as the patients with type 1 diabetes five years after the diagnosis. For every patient, only the data of the first foot examination were included in the analysis. The protocol was adopted by the Slovenian Working Group on the Diabetic Foot in 1996, and it is coherent with international recommendations (13). The examination is done by a specialist nurse and evaluated by a diabetologist. The protocol (Appendix 1) includes demographic data (ID, age, sex), medical history (previous foot ulcer, amputation, various symptoms) and foot examination (various deformities, hard skin, ulcer, dry skin, redness, arterial pulses). Risk status classification (groups 1 – 4, Table 1) is done according to the data from medical history and the findings upon foot examination. Therapeutic measures taken after the examination are also recorded: education, footwear prescription, referrals (foot clinic, angiologist, surgeon and chiropodist). The duration of diabetes and the level of metabolic regulation (HbA1c level) are not included in the protocol. 10.1515/sjph-2017-0009 Zdr Varst 2017; 56(1): 65-73 66 10.1515/sjph-2017-0009 Zdr Varst 2017; 56(1): 65-73 67 2.2 Multiple Correspondence Analysis The relationship between variables was summarized by cross-tabulations and analysed by multiple correspondence analysis (MCA). MCA is a powerful descriptive statistical technique for handling larger, more complex datasets, including high-dimensional categorical data often encountered in social sciences, marketing, health economics and biomedical research (17). The primary goal of this exploratory statistical method is to transform numerical information into graphical displays (‘maps’) and related numerical statistics. The position of the category-points in MCA maps is the basis for revealing the relationship among the investigated variables (17, 18). The original data set consists of 56 variables. We included the data of both feet about previous and acute ulceration, amputation and foot pulses. As we have previously demonstrated that the data on neuropathy and foot deformity did not differ between the two feet (16), we decided to include only the data of one (right) foot in the analysis. Among the remaining 24 variables under consideration, all except for age and loss of protective sensation were nominal, mostly dichotomous. The list of the analysed variables is shown in Appendix 2. Age has been recoded to 3 age groups (0: under 51 years, 1: 51 – 70, 2: more than 70 years). The loss of protective sensation has been recoded to 3 categories (0: no loss of protective sensation, 1: unable to feel the monofi lament on 1 – 3 points, 2: unable to feel the monofi lament on 4 or more points). In this way, there were 52 categories to be analysed altogether. Statistical analysis was performed with the SAS 9.4, procedures Corresp and IML. The correspondence maps were designed with XLSTAT 2015, procedure Multiple Correspondence Analysis (MCA). Table 1. Risk status classifi cation. 1 2 3 4 Normal sensation, palpable pedal pulses, no foot deformity Impaired sensation, no foot deformity Absent pedal pulses, normal sensation, no foot deformity • Combination: foot deformity and impaired sensation and / or absent foot pulses • Previous ulcer or amputation • Charcot foot Category Description 2.3 Burt Matrix MCA can be defi ned as the correspondence analysis of the so-called ‘Burt matrix’ B. It is a partitioned symmetric matrix containing all pairs of crosstabulations among a set of categorical variables. Each crosstabulation Fqq (q=1,2,...,Q) on the main diagonal of the Burt matrix is a diagonal matrix of the marginal frequencies (i.e., a crosstabulation of a variable with itself). Each off-diagonal crosstabulation is an ordinary two-way contingency table Fqq’ (q,q’=1,2,...,Q, q≠q’).Each contingency table above the diagonal has a transposed counterpart below the diagonal. Table 2. The Burt matrix, the presentation of crosstabulations reduced to the fi rst two variables and the last variable (the complete matrix at (19)). SEXF SEXM NSY0 NSY1 AMP0 AMP1 Total 5193 0 2519 2674 5167 26 124632 2519 3876 6395 0 6376 19 153480 5167 6357 6376 5148 11524 0 276576 0 6401 3876 2525 6357 44 153624 2674 2525 0 5199 5148 51 124776 26 44 19 51 0 70 1680 124632 153624 153480 124776 276576 1680 6678144 SEXF NSY0 AMP1SEXM AMP0NSY1 Total As mentioned before, we have analysed a data set of 24 variables (Q=24), containing altogether 52 categories (J=52). Therefore,the Burt matrix is a symmetric matrix consisting of 52 rows and columns (Table 2). 10.1515/sjph-2017-0009 Zdr Varst 2017; 56(1): 65-73 68 On the basis of the Burt matrix, one can calculate rows of relative frequencies (the division of row frequencies by the marginal total), called profi les, representing each variable category (Table 3). The marginal row of relative frequencies at the bottom of the table presents the average profi le (Total). Categories with similar profi les have similar characteristics for the variables considered. Unfortunately, due to numerous categories, direct inspection of all profi les and formation of homogenous profi le groups is impossible. For that reason, the numerical information contained in the Burt matrix has been analysed with MCA. Table 3. The matrix of row profi les (source: Table 2). 2.4 Dimensionality of MCA Solution Based on the Burt matrix one can form the matrix where DB is a diagonal matrix with the elements of the average profi le on the main diagonal, PB is a matrix of relative frequencies of the Burt matrix (division of Burt matrix by the grand total) and pB is an average profi le (Table 3). The expression in brackets represents the deviations of observed relative frequencies of the Burt matrix from the model of independence, while DB-1/2 is a weighting factor matrix. Basic calculations are made in three steps: 1. Find the eigenvalue-eigenvector decomposition of the matrix (2) where is an orthogonal matrix of eigenvectors and is a diagonal matrix of eigenvalues The number of nontrivial singular values is equal to L = J - Q (L=28). In case of the analysis of Burt table, the principal inertias are defi ned as the squared values of singular values . 2. Calculate the matrix of standard coordinates 3. The corresponding principal coordinates are given by SEXF SEXM NSY0 NSY1 AMP0 AMP1 Total 4.1667 0.0000 1.6413 2.1430 1.8682 1.5476 1.8663 2.0212 2.5230 4.1667 0.0000 2.3053 1.1310 2.2982 4.1458 4.1380 4.1543 4.1258 4.1667 0.0000 4.1415 0.0000 4.1667 2.5254 2.0236 2.2985 2.6190 2.3004 2.1455 1.6436 0.0000 4.1667 1.8613 3.0357 1.8684 0.0209 0.0286 0.0124 0.0409 0.0000 4.1667 0.0252 100 100 100 100 100 100 100 SEXF NSY0 AMP1SEXM AMP0NSY1 Total 10.1515/sjph-2017-0009 Zdr Varst 2017; 56(1): 65-73 69 Table 5. The adjustment of inertias. 1 2 3 4 5 6 7 8 9 0.20583 0.11983 0.07324 0.06369 0.04969 0.04743 0.04398 0.04196 0.04174 0.04237 0.01436 0.00536 0.00406 0.00247 0.00225 0.00193 0.00176 0.00174 66.72182 15.12592 2.46806 1.20083 0.15938 0.08224 0.01325 0.00021 0.00001 0.02934 0.00665 0.00109 0.00053 0.00007 0.00004 0.00001 0.00000 0.00000 66.72182 81.84773 84.31580 85.51663 85.67600 85.75824 85.77149 85.77170 85.77171 ********************** *** * Singular values Burt principal inertia Percentage of adjusted inertia Adjusted inertia Cumulative percentage of adjusted inertia 16 32 48 64 80 ----+----+----+----+----+--- with the condition that these be calculated for only. Table 4. The fi rst two principal coordinates of the matrix Y. SEXF SEXM NSY0 NSY1 HVA0 HVA1 AMP0 AMP1 0.045545 -0.036950 -0.337288 0.414879 -0.071098 0.340644 -0.011657 1.919002 0.264007 -0.214184 -0.219527 0.270028 -0.165217 0.791587 -0.010940 1.801034 2. principal coordinate 1. principal coordinate MCA includes the fi tting of the diagonal submatrices Fqq (q=1,2,...,Q) of the Burt matrix. As a result, the total inertia is infl ated and thus the proportions of the fi rst few principal inertias as parts of the total inertia are reduced (21, p.155). One way to address this problem, proposed by Benzécri (22), is to consider only those principal axes whose eigenvalues are higher than 1/Q (i.e. 1/24=0.041667 in this application, see Table 5). These adjusted inertias can be calculated according to Benzécri’s formula The number of inertias has been reduced from 28 to only 9, with strongly dominating values of the fi rst two principal inertias. The values of the adjusted inertias are shown in Table 5. 0.04398 10.1515/sjph-2017-0009 Zdr Varst 2017; 56(1): 65-73 70 The next question is the quality of the presentation of the position of the profi les based on the fi rst few principal coordinates. M. Greenacre (21) calculates the percentage of inertia as follows where is an average of the off diagonal inertias , i.e. where is the off-diagonal inertia; e.g., chi-square measure of association (based on particular crosstabulation in the Burt matrix) divided by the number of units observed. The values of percentages of adjusted principal inertias and the cumulative percentages of adjusted principal inertia are also shown in Table 5. According to the aforementioned Greenacre’s approach to the calculation of percentages of inertia, 81.85% of the total inertia is explained by the fi rst two principal axes. Therefore, without losing too much information, the representation of the positions of the profi les (category points) based on the fi rst two principal coordinates (on a two-dimensional map) can serve as a good basis for the thorough analysis of the relationship between considered variables. 3 RESULTS: MAPS AND ANALYSIS The spread of the profi les is evaluated as total inertia. When the cumulative percentage of inertia of the fi rst two dimensions is relatively high, then most of the profi les are well represented in a two-dimensional map (by their projections onto a plane). Thus the two-dimensional solution of our example (Table 4) explaining 81.85% of the total inertia (Table 5) can serve as a good basis for the display of the profi les (Figure 1). In correspondence analysis, the average profi le is represented as a point in the origin of the coordinate system. The deviation of a particular variable category point from the origin is determined by its share in the sample: the lower the share, the farther the category point from the origin. Figure 1. A two-dimensional map – a display of all categories (extended labels in Appendix 2). In our case (Figure 1), the categories with high shares (younger age group - AGE0, no history of foot ulcer – HUL0, no loss of protective sensation - LPS0, absence of foot deformities, palpable pulses of pedal arteries) are concentrated in the left lower quadrant, close to the origin. On the other side, the points representing the presence of particular foot deformities, loss of protective sensation, and acute foot ulceration are distributed in the right upper quadrant. There is a separate group of category points in the right lower quadrant, which represents the patients with absent foot pulses. Close to them is the point representing referral to vascular specialist. The point representing group 3 (pure ischaemia and no neuropathy) is a little away. 10.1515/sjph-2017-0009 Zdr Varst 2017; 56(1): 65-73 71 Figure 2. A zoom-in display of the upper part of Figure 2. The zoomed-in upper part of Figure 1, where the majority of the points are located, is shown in Figure 2. The points representing the ordinal variables are connected with broken lines. The solid line connecting the age groups (AGE1, AGE2, AGE3) is stretched in the direction from the left lower to the right upper quadrant. The youngest age group (under 50, AGE0) is situated at the left lower end. The dotted line, representing the ordinal character of protective sensation loss, from LPS0 to LPS2, follows almost the same direction as age groups, with no loss of sensation (LPS0) in the left lower quadrant and loss of sensation (LPS2) in the right upper quadrant. 4 DISCUSSION AND CONCLUSIONS The decision to conduct the study was based on the belief that it could contribute to improved diabetic foot care in Slovenia. This would be important from the public health aspect, since Slovenia had the second highest age-sex standardised amputation rate per 100 000 population among the observed countries, and the highest age-sex standardised amputation rate per 100 000 people with diabetes, according to the data from the OECD health statistics report published in 2015 (23). Regular preventative foot examinations are the cornerstone of effi cient amputation prevention strategies. With the introduction of a new organisation at the primary health care level, called model practices, the frequency of foot examinations has increased signifi cantly (24). Our retrospective analysis of the large foot screening data set has confi rmed most of the already known relationships in the development of foot pathology, and revealed a stronger association of foot deformity with female than with male gender. The survey has confi rmed that the risk status classifi cation and the post-screening decisions were concordant with the data obtained by screening. We have decided to use MCA because it has proved to be a strong tool for examining the association between signs and symptoms described by mostly nominal data. In 1989, Crichton described its use on two big data sets of patients with chest pain and acute abdominal pain (25). MCA has been used in the analysis of data in psychiatry (26), rheumatology (27), infectious diseases (malaria) (28) and oral health (29). Recently, Sacco et al.(30) used MCA to reveal early indicators of diabetic polyneuropathy. Their study included 193 patients without partial or total foot amputation, major vascular complications or severe nephropathy. On the contrary, our study was performed on an unselected population of patients with diabetes. To the best of our knowledge, no analysis of such large data set has been published so far. Screening is a strategy used in a population to identify the possible presence of an as-yet-undiagnosed disease in asymptomatic individuals. The majority of the screened population is supposed to be disease-free, and the points representing the individuals without a condition are therefore concentrated close to the origin of the coordinate system. As expected, the points representing the absence of impaired blood supply, foot deformity, neuropathic symptoms and previous or current foot ulceration were all located close to the origin. The decline of neurological functions (in particular vibrating sensation) with aging is well established (31, 32, 33). Accordingly, we found the position of points representing foot deformity and neuropathic symptoms close to the point of the age group over 70 years. Unfortunately, our protocol does not include the data on diabetes duration and metabolic control, which are known risk factors for chronic complications of diabetes (34). Some other apparently obvious associations have been demonstrated: the category point representing the history of foot ulceration was close to the points representing the loss of protective sensation, acute foot ulceration and fat pad atrophy. The points representing various foot deformities (claw toes, hallux valgus, toenail deformity) were close together and also close to abundant callus. Dry skin (indicating autonomic neuropathy) was close to neuropathic symptoms, and both were close to moderate loss of protective sensation (LPS1). The location of foot deformity points closer to female than male gender was unexpected, but could possibly be explained by the inappropriate footwear in women. 10.1515/sjph-2017-0009 Zdr Varst 2017; 56(1): 65-73 72 Foot screening in our centre is done by a nurse and evaluated by a physician –diabetologist. The risk category classification is a subjective decision based on the data from medical history and foot examination, and as such it might be prone to error. We conclude that it was done properly, since the points representing the four risk categories were distributed as expected: The risk group 1 point (normal sensation, no deformity) is among the points representing the absence of particular foot pathologies. The risk group 2 point (the loss of protective sensation, no deformity) is located away from the foot deformity points. The risk group 3 point (absent foot pulses, normal sensation, no deformity) is close to absent foot pulses and far away from neuropathy. The risk group 4 point (previous ulcer or amputation, Charcot, combination of neuropathy, ischemia and foot deformity) is close to those representing neuropathy, foot deformity, ulcer and amputation. The patients with impaired blood supply form a special group – possibly due to subjective symptoms, they seek help in earlier stages, and are therefore not discovered only by screening. It might be surprising that the point of group 3 (the group with pure ischaemia and no neuropathy) was not closer to the points representing absent foot pulses. A possible explanation is that the patients with poor blood supply may have other problems, a combination of disease changes, which classify them into the risk group 4. The post-screening decisions were also appropriate: the point representing footwear prescription was close to foot deformity, referral to foot clinic was close to acute foot ulceration, referral to angiologist close to absent foot pulses, and referral to the surgeon closer to acute foot ulceration than to absent pedal pulses. The main weakness of our study is that it included patients with different duration of diabetes. With longer duration of diabetes, the prevalence of neuropathy and peripheral arterial disease increase. Nevertheless, their influence on the development of foot ulceration mainly depends on their severity and not on duration. In human medicine, we are often faced with the situations where categorical (nominal and ordinal) variables are predominant. Even some laboratory results, although physical readings, are essentially of ordinal nature. For that reason, we believe that MCA can be useful in the analysis of medical data. CONFLICTS OF INTEREST The authors declare that no conflicts of interest exist. FUNDING The study had no funding. ETHICAL APPROVAL The data analysed in this study were collected at the University Medical Centre Ljubljana without information about the identity of individuals. The study was conducted in accordance with the code of Ethics of the World Medical Association (Declaration of Helsinki). REFERENCES 1. International Diabetes Federation. IDF diabetes atlas. 6th ed. Brussels: International Diabetes Federation, 2013. Available Feb 01, 2016 from: http://www.idf.org/diabetesatlas 2. Zaletel J. Sladkorna bolezen kroji mero današnjemu in prihodnjemu zdravju prebivalstva. Zdrav Vest 2012; 81: 681-3. 3. Nerat T, Kos M. Burden of type 2 diabetes from the healthcare payer perspective in Slovenia. Zdr Varst 2013; 52: 162-80. 4. Boulton AJ. The diabetic foot: from art to science; the 18th Camillo Golgi lecture. Diabetologia 2004; 47: 1343-53. 5. Vamos EP, Bottle A, Edmonds ME. et al. Changes in the incidence of lower extremity amputations in individuals with and without diabetes in England between 2004 and 2008. Diabetes Care 2010; 33: 2592-7. 6. Jeffcoate WJ, Harding KG. Diabetic foot ulcers. Lancet 2003; 361: 1545-51. 7. Urbancic-Rovan V. Causes of diabetic foot lesions. Lancet 2005; 366: 1675-6. 8. Maser RE, Steenkiste AR, Dorman JS. et al. Epidemiological correlates of diabetic neuropathy: report from Pittsburgh Epidemiology of Diabetes Complications Study. Diabetes 1989; 38: 1456-61. 9. Vinik AI, Park TS, Stansberry KB, Pittenger GL. Diabetic neuropathies. Diabetologia 2000; 43: 957-73. 10. Boulton AJ, Vinik AI, Arezzo JC. et al, on behalf of American Diabetes Association. Diabetic neuropathies: a statement by the American Diabetes Association. Diabetes Care 2005; 28: 956–62. 11. Murray HJ, Young MJ, Hollis S, Boulton AJ. The association between callus formation, high pressures and neuropathy in diabetic foot ulceration. Diabet Med 1996; 13: 979-82. 12. Pataky Z, Golay A, Faravel L. et al. The impact of callosities on the magnitude and duration of plantar pressure in patients with diabetes mellitus: a callus may cause 18,600 kilograms of excess plantar pressure per day. Diabetes Metab 2002; 28: 356-61. 13. Apelqvist J, Bakker K, van Houtum WH, Nabuurs-Franssen MH, Schaper NC, editors, International Working Group on the Diabetic Foot. Maastricht: International consensus on the diabetic foot, 1999. 14. Farlex Partner Medical Dictionary © Farlex 2012. Available Feb 01, 2016 from: http://medical-dictionary.thefreedictionary.com/ screening 15. Zaletel J. Sladkorna bolezen v letu 2015. Available Feb 01, 2016 from: http://www.nijz.si/sites/www.nijz.si/files/uploaded/jelka_ zaletel_sladkorna_bolezen_v_letu_2015.pdf 16. Rovan J, Urbančič-Rovan V, Slak M. Exploration of categorical screening procedure data by multiple correspondence analysis. In: Vidmar G, Stare J, editors. International Conference Applied Statistics 2007, September 23-26, Ribno (Bled), Slovenia. Ljubljana: Statistical Society of Slovenia, 2007: 95-6. 17. Greenacre MJ, Blasius J, editors. Multiple correspondence analysis and related methods. London: Chapman & Hall, 2006. 18. Greenacre MJ. Theory and applications of correspondence analysis. London: Academic Press, 1984. 19. Urbančič-Rovan V. [Personal survey]. Unpublished raw data, 2015. 20. Rovan J. Visualizing solutions in more than two dimensions. In: Greenacre MJ, Blasius J, editors. Correspondence analysis in the social sciences. London: Academic Press, 1994: 210-29. 21. Greenacre MJ. Multiple and joint correspondence analysis. In: Greenacre MJ, Blasius J, editors. Correspondence analysis in the social sciences. London: Academic Press, 1994: 141-61. 22. Benzécri JP. Sur le calcul des taux d’inertie dans l’analyse d’un questionnaire: addendum et erratum a [BIN.MULT]. Cahiers de L’analyse des Données 1979; 4: 377-8. 23. ECD. “Diabetes care”, in Health at a Glance 2015: OECD Indicators. Paris: OECD Publishing, 2015. 10.1515/sjph-2017-0009 Zdr Varst 2017; 56(1): 65-73 73