Using perceptive computing in multiple sclerosis - the Short Maximum Speed Walk test
© Behrens et al.; licensee BioMed Central Ltd. 2014
Received: 11 December 2013
Accepted: 22 May 2014
Published: 27 May 2014
We investigated the applicability and feasibility of perceptive computing assisted gait analysis in multiple sclerosis (MS) patients using Microsoft Kinect™. To detect the maximum walking speed and the degree of spatial sway, we established a computerized and observer-independent measure, which we named Short Maximum Speed Walk (SMSW), and compared it to established clinical measures of gait disability in MS, namely the Expanded Disability Status Scale (EDSS) and the Timed 25-Foot Walk (T25FW).
Cross-sectional study of 22 MS patients (age mean ± SD 43 ± 9 years, 13 female) and 22 age and gender matched healthy control subjects (HC) (age 37 ± 11 years, 13 female). The disability level of each MS patient was graded using the EDSS (median 3.0, range 0.0-6.0). All subjects then performed the SMSW and the Timed 25-Foot Walk (T25FW). The SMSW comprised five gait parameters, which together assessed average walking speed and gait stability in different dimensions (left/right, up/down and 3D deviation).
SMSW average walking speed was slower in MS patients (1.6 ± 0.3 m/sec) than in HC (1.8 ± 0.4 m/sec) (p = 0.005) and correlated well with EDSS (Spearman’s Rho 0.676, p < 0.001). Furthermore, SMSW revealed higher left/right deviation in MS patients compared to HC. SMSW showed high recognition quality and retest-reliability (covariance 0.13 m/sec, ICC 0.965, p < 0.001). There was a significant correlation between SMSW average walking speed and T25FW (Pearson’s R = -0.447, p = 0.042).
Our data suggest that ambulation tests using Microsoft Kinect™ are feasible, well tolerated and can detect clinical gait disturbances in patients with MS. The retest-reliability was on par with the T25FW.
KeywordsComputerized gait assessment Gait impairment Multiple sclerosis Walking speed T25FW
Multiple sclerosis (MS) is a common chronic inflammatory and neurodegenerative disease that normally begins in young adulthood, typically affecting patient quality of life and leading to high rates of early retirement [1–4]. Impairment of gait and balance are major factors that restrict daily activity  and may occur as early as after a first clinical episode [6–8].
In current clinical practice, balance and gait impairment are quantified using a combination of clinical examination and patients’ reported maximum walking distance: The Kurtzke‘s Functional Systems and Expanded Disability Status Scale (EDSS)  are widely used both in clinical practice and for clinical trials. The EDSS provides a good overview of current neurologic status, gait impairment and mobility dysfunction, but it has limitations. It is a relatively subjective measure with high intra- and inter-rater variability and quantifying mild symptoms and symptom progression is challenging [10, 11]. The Timed 25-Foot Walk (T25FW) measures the time a patient takes to walk 25 feet at maximum speed. As part of the Multiple Sclerosis Functional Composite (MSFC), the T25FW is currently the most widely implemented method to objectively quantify gait disability in clinical MS trials and, to a lesser extent, clinical practice [12, 13]. Although an excellent method of quantifying overall gait disability, the T25FW measures only the time taken to walk a set distance . As alternatives to the T25FW, timed walking tests might improve reliability. In particular, the 2-minute walk test is proposed as an additional outcome parameter in clinical trials .
Recently, more objective motion capture systems (for an overview see ) have been proposed as tools for detecting the pattern of gait impairment more precisely than clinical examination using conventional tools (i.e. EDSS) [17, 18]. However, none of these systems has found its way into clinical routine yet, and are therefore rarely available in outpatient clinics and neurologic practices.
With an infrared light camera system, the Microsoft Kinect™, which was originally developed for video gaming, detects anatomical landmark positions in three dimensions (3D) . This function is based on a decision forest method that harnesses mass data of sensor-recorded skeletal joint movements . Using the programming interface made available via the Kinect software development kit (SDK) , it has in the past been employed to detect balance and motion of healthy subjects [22, 23] and stroke patients [24, 25]. Furthermore, some few studies have reported methods for gait parameter approximation using Kinect [26–28], however, none of these has been applied to patients with MS.
The objective of this study was to evaluate the feasibility of computerized versions of walking tests using Microsoft Kinect and to compare the resulting findings for MS patients with low to moderate neurological impairment to those for healthy controls. We established the Short Maximum Speed Walk test (SMSW) as a measure to analyse patients’ gait speed and degree of sway. We report feasibility, reliability and correlation of results with those of EDSS and the T25FW.
The local ethics committee of the Charité – Universitätsmedizin Berlin approved the study (EA1/225/12). It was conducted in accordance to the Declaration of Helsinki in its currently applicable version. All patients and healthy subjects gave written informed consent.
Twenty-two patients with diagnosed MS according to the current panel criteria  (age mean ± SD 43 ± 9 years, 13 female) and 22 age and sex matched healthy controls (HC) (age 37 ± 11 years, 13 female) were enrolled in the study. MS patients were recruited from on-going clinical trials at Charité’s neuroimmunology outpatient clinic. HC were recruited from volunteers. Blinding was not attempted, since the operators knew many subjects and visual contact was mandatory. MS patients were first clinically graded using EDSS (median 3.0, range 0.0 – 6.0) under supervision of a neurologist. The EDSS is calculated from a detailed neurological examination using functional system scores that assess the visual system function, brainstem function, pyramidal tract function, cerebellar function, sensory system function, bowel and bladder function, cerebral function and patient ambulation . Subjects additionally performed the T25FW component of the Multiple Sclerosis Functional Composite (MSFC) . T25FW results were not available from eight of the 22 HC, and one of the 22 MS patients. The first 8 HC had been tested before the T25FW was added to the study protocol and the MS patient could not complete the study in full due to time constraints.
Perceptive computing assisted motor assessment
Short Maximum Speed Walk test analysis
Overview of SMSW parameters
Average walking speed
Average speed during the measurement calculated as observed distance divided by observed time
Main outcome parameter; measures a subject’s average walking speed during the test
Standard deviation of speed between sequent frames
Measure of directional speed homogeneity
3D direction deviation
Mean square error from the Euclidean distance between the actual position and the position from a linear model approximation over the main z-vector
Combined expression of how much directional, up/down and left/right variability a subject’s gait shows
Left/right deviation between the actual position and the position as determined by a linear model approximation over the main z-vector
Parameter describes a subject’s left/right movement during the walking test.
Up/down deviation between the actual hip center position and the position from a linear model over the main z-vector.
Parameter describes a subject’s up/down movement during the walking test.
Total time the subject was tracked
Total distance the subject was tracked
Number of detected calibration jumps during the measurement
Mean distance of all calibration jumps detected during the measurement
Total number of frames collected
First, joint coordinates were normalized, correcting for sensor tilt and walking direction of the subject by computing a linear approximation of the hip-centre movement, translating the start of the movement to the point of origin and rotating the direction of movement onto the positive Z-axis.
As in any automated assessment tool, artefacts (in the form of tracking errors, or calibration jumps) were inevitable. Specifically, detection of the spatial location of the crucial hip-centre joint by the SDK system was clearly sometimes erroneous as determined by visual inspection. Furthermore, the SDK’s in-built error detection system (SDK Recognition Quality) was unable to identify these artefacts with sufficient accuracy, and instead falsely reported all joint detections as successful. To address this problem, we developed an alternative error correction technique. Here, we mathematically extracted the acceleration of the joint movement between two frames as recorded by the SDK system. We analysed the transition between the points in terms of a subject’s overall acceleration and averaged these values across the study cohort. We found that if aggregated acceleration over all coordinate axes exceeded a threshold of 0.6 m/frames2, the frames were likely to be an artefact. As such they were flagged as suspected tracking errors (“calibration jump”) and excluded from further analysis.
Using the calibration jump-corrected data, we computed the SMSW test parameters given in Table 1, which were then used for statistical analysis.
All data post-processing and analysis was performed using Matlab 2012A (Mathworks, Ismaning, Germany).
Analysis of measurement quality
Observed time (sec)
Observed distance (m)
Number of calibration jumps
Average calibration jump distance (m)
Measurement length (frames)
Analysis of immediate retest reliability
Immediate retest reliability
SMSW Average speed
SMSW Speed deviation
SMSW 3D direction deviation
SMSW Left/right deviation
SMSW Up/down deviation
Immediate retest reliability was analysed using intra-class correlation coefficients (ICC) with a one-sided, random model. Test results between MS patients and HC were compared using t-tests. When equal variances could be assumed (Levene’s test p > 0.05), p-values from independent sample Student’s t-tests are given. When equal variances could not be assumed (Levene’s test p < 0.05), results from the Welch’s t-test are given. The association between T25FW and SMSW measures was investigated using Pearson’s correlation analysis and a Bland-Altman analysis. Correlation with EDSS and functional system scores was investigated with Spearman’s Rho analysis due to the ordinal nature of EDSS data.
Statistical analysis was performed using SPSS version 21 (IBM, Armonk, NY, USA). A p < 0.05 was deemed significant. All tests should be understood as exploratory data analysis as no prior power calculation and subsequent corrections for multiple testing were applied.
Comparison between healthy controls and multiple sclerosis patients
Average speed (m/s)
Speed deviation (m/frame2)
3D direction deviation (m2)
Left/right deviation (cm)
Up/down deviation (cm)
Comparison between SMSW and T25FW
Association with EDSS and functional system scores
Correlation between test outcomes and EDSS
SMSW Average speed
SMSW Speed deviation
SMSW Derived T25FW
SMSW 3D direction deviation
SMSW Left/right deviation
SMSW Up/down deviation
In this cross-sectional study we investigated the applicability and feasibility of perceptive computing assisted motion analysis in MS patients using the Microsoft Kinect system. Our primary question was whether skeletal tracking data recorded by the device’s SDK program would yield reliable information for gait analysis. Using only the hip-centre joint, we established the test SMSW to analyse MS patients’ gait at maximum walking speed. Additionally, we analysed the SMSW left/right and up/down deviation of MS patients gait in comparison to HC.
Using the SMSW, we were able to detect differences in average walking speed between MS patients and HC. Several gait stability parameters also showed promising results: Especially 3D direction deviation and left/right deviation showed prominent differences between HC and MS patients. The lateral sway of MS patients (described by the parameter left/right deviation) detected by our test suggests that stability impairments are indeed detectable by perceptive computing systems, even over limited observation time and distance (although some types of gait variability may only become evident over longer distances) .
The linear distance covered by the recognition area was only 2.2 metres (see Figure 1 and Table 2), and consequently only a few steps of each subject were analysed. Despite this short walking and recording time, the overall detection quality of the target hip-centre joint of the SMSW over time was excellent. Calibration jumps (i.e. tracking errors by the SDK system) presented regularly (see Table 2), however, these could be filtered out by applying a custom-built filter developed during the analysis of this experiment. After filtering out the calibration jumps, the retest reliability of the SMSW average walking speed parameter was excellent and on par with T25FW. The gait stability parameters showed less retest reliability. This was likely due to the fact that most differences between MS patients and HC were only marginally significant. The reliability of the gait stability measurements was superior in MS patients, where a higher signal-to-noise ratio was expected due to weaker gait stability. Additionally, the detection limit was heavily influenced by the joint recognition quality of the sensor. In fact, one measure (3D direction deviation), originally designed to detect changes in walking stability, proved to be very sensitive to calibration jumps in walking direction, which primarily occurred at faster walking speed.
In combination, these data suggest that an assessment employing gait analysis at lower walking speed might improve detection of gait variability and stability in mobility impaired patients. This would not only lower detection errors in the form of calibration jumps but would also allow the system to record several step cycles for a more robust analysis. This is supported by a previous study investigating a Kinect-based system for gait measurements using 3D point cloud information instead of skeletal tracking data, which reported accurate estimation of gait parameters compared to a Vicon motion analysis system in healthy volunteers . A second study used skeletal data from all reported joints to analyse a postural control assessment battery in healthy subjects . Our study adds to these findings and shows that Kinect-based motion analysis is also feasible in MS patients and is able to reliably detect gait differences in comparison to healthy controls. Although the experimental setups and analysis approaches differ, all three studies show that inexpensive 3D perceptual computing can be reliably used to quantify motion impairment. Furthermore, a recent work of Stone and Skubic showed, Kinect can also be installed in the patient’s home as a means of frequently and conveniently assessing gait .
We expected that SMSW average speed and the T25FW would yield similar results. However, MS patients’ results sometimes differed strongly in the two tests. Although typical MS symptoms like reduced muscular strength, spasticity or balance deficiency are widely known to fluctuate , it is unlikely that random irregularities or gait variability are the key factors in this discrepancy, because retest reliability was high in both tests.
Instead, a different influence of functional impairments on both tests seems possible. SMSW and T25FW correlated well with overall EDSS degree of neurological disability. However, when drilling-down by comparing individual parameter results to the corresponding EDSS functional system scores (FS), we were surprised to find a difference between the SMSW and the T25FW: Whereas SMSW average walking speed showed a moderate correlation to cerebral function, T25FW results correlated with visual function, and SMSW walking speed showed correlation with the pyramidal scores of the EDSS.
Pyramidal impairment such as muscular weakness and spasticity clearly restricts the gait of many MS patients [34, 35]. However, as stated above, we found no correlation between the T25FW and the pyramidal function score of the EDSS. This is in discrepancy with previous findings, such as those of Phan-Ba et al. , and likely explained by the only mild pyramidal impairment of many patients in our study (41% patients with FS ≤ 2 ). In contrast to T25FW, SMSW correlated with pyramidal functional system scores. This might be explained by findings from previous studies where MS patients, performing the T25FW, accelerated walking speed after a relatively slow start [36, 37]. Because muscular rigidity is triggered by postural changes , such as starting from a static position to walk, we may have detected a starting weakness of spastic patients. Thus, a patient, who begins moving more slowly, might show a slower walking time at SMSW average walking speed compared to T25FW. In the T25FW a slower speed at start might level over the longer distance. In this regard SMSW average walking speed may be more sensitive to impaired initiation of movement in MS patients than the T25FW. We did not evaluate spasticity with a validated assessment, so we were not able to further pursue the question of whether the SMSW could be a marker for minor pyramidal weakness.
An exploratory pilot study of this type has some limitations. Only Kinect-derived skeletal tracking information generated by the SDK was used . Although not within the scope of the present study, it is conceivable that exclusive or hybrid inclusion of raw sensor data might yield additional or better results. Furthermore, we neither selected MS patients with specific and well-established motor symptoms or gait disability, nor did we compare SMSW test results to objective gait measures other than T25FW. The analysis between SMSW, T25FW and EDSS functional system scores should be treated as preliminary and the results should be interpreted with care. However, the demonstrated correlations fit well within the expected clinical framework and the results of previous studies. As such, our results can be expected to serve as a basis for confirmatory studies. Most importantly, the ability of SMSW measurements to detect changes longitudinally and the retest reliability of the test needs to be validated in a follow-up study.
Perceptive computing-based detection of ambulation speed via the joint hip-centre was feasible and reliable. SMSW average walking speed was a valid parameter as demonstrated by retest reliability results and the strong correlation with established clinical scores, such as the T25FW and EDSS. A notable difference was the slower walking speed of MS patients measured using SMSW (mean distance: 2.2 metres) compared to T25FW (mean distance: 7.6 metres). The gait stability parameters showed greater left/right deviation in MS patients compared to HC. The analysis approach presented here shows promise as an objective technique for detecting and assessing of gait pattern and as a simple and affordable tool in the clinician’s toolbox. Finally, this study provides evidence that further investigation of gait with periodic Kinect-measurements may give new insight into the disease progression of MS.
We thank Cynthia Kraut and Susan Pikol for excellent technical assistance. All primary data is write protected and archived under stewardship of the Charité – Universitätsmedizin Berlin’s office for Good Clinical Practice. It is available for verification upon request.
- Noseworthy JH, Lucchinetti C, Rodriguez M, Weinshenker BG: Multiple sclerosis. N Engl J Med 2000,343(13):938–952.PubMedView Article
- Patwardhan MB, Matchar DB, Samsa GP, McCrory DC, Williams RG, Li TT: Cost of multiple sclerosis by level of disability: a review of literature. Mult Scler 2005,11(2):232–239.PubMedView Article
- Krause I, Kern S, Horntrich A, Ziemssen T: Employment status in multiple sclerosis: impact of disease-specific and non-disease-specific factors. Mult Scler 2013,19(13):1792–1799.PubMedView Article
- Borisow N, Döring A, Pfueller CF, Paul F, Dörr J, Hellwig K: Expert recommendations to personalization of medical approaches in treatment of multiple sclerosis: an overview of family planning and pregnancy. EPMA J 2012,3(1):9.PubMed CentralPubMedView Article
- Findling O, Sellner J, Meier N, Allum JHJ, Vibert D, Lienert C, Mattle HP: Trunk sway in mildly disabled multiple sclerosis patients with and without balance impairment. Exp Brain Res 2011,213(4):363–370.PubMedView Article
- Kalron A, Achiron A, Dvir Z: Motor impairments at presentation of clinically isolated syndrome suggestive of multiple sclerosis: Characterization of different disease subtypes. NeuroRehabilitation 2012,31(2):147–155.PubMed
- Kalron A, Achiron A, Dvir Z: Muscular and gait abnormalities in persons with early onset multiple sclerosis. J Neurol Phys Ther 2011,35(4):164–169.PubMedView Article
- Martin CL, Phillips BA, Kilpatrick TJ, Butzkueven H, Tubridy N, McDonald E, Galea MP: Gait and balance impairment in early multiple sclerosis in the absence of clinical disability. Mult Scler 2006,12(5):620–628.PubMedView Article
- Kurtzke JF: Rating neurologic impairment in multiple sclerosis: an expanded disability status scale (EDSS). Neurology 1983,33(11):1444–1452.PubMedView Article
- Noseworthy JH, Vandervoort MK, Wong CJ, Ebers GC: Interrater variability with the Expanded Disability Status Scale (EDSS) and Functional Systems (FS) in a multiple sclerosis clinical trial. The Canadian Cooperation MS Study Group. Neurology 1990,40(6):971–975.PubMedView Article
- Goodkin DE: EDSS reliability. Neurology 1991,41(2 Part 1):332–332.PubMedView Article
- Fischer JS, Rudick RA, Cutter GR, Reingold SC: The Multiple Sclerosis Functional Composite Measure (MSFC): an integrated approach to MS clinical outcome assessment. National MS Society Clinical Outcomes Assessment Task Force. Mult Scler 1999,5(4):244–250.PubMedView Article
- Rosti-Otajärvi E, Hämäläinen P, Koivisto K, Hokkanen L: The reliability of the MSFC and its components. Acta Neurol Scand 2008,117(6):421–427.PubMedView Article
- Hobart J, Blight AR, Goodman A, Lynn F, Putzki N: Timed 25-foot walk: direct evidence that improving 20% or greater is clinically meaningful in MS. Neurology 2013,80(16):1509–1517.PubMedView Article
- Gijbels D, Dalgas U, Romberg A, de Groot V, Bethoux F, Vaney C, Gebara B, Medina CS, Maamâgi H, Rasova K, de Noordhout BM, Knuts K, Feys P: Which walking capacity tests to use in multiple sclerosis? A multicentre study providing the basis for a core set. Mult Scler 2012,18(3):364–371.PubMedView Article
- Pearson OR, Busse ME, van Deursen RWM, Wiles CM: Quantification of walking mobility in neurological disorders. QJM 2004,97(8):463–475.PubMedView Article
- Sosnoff JJ, Weikert M, Dlugonski D, Smith DC, Motl RW: Quantifying gait impairment in multiple sclerosis using GAITRite technology. Gait Posture 2011,34(1):145–147.PubMedView Article
- Sosnoff JJ, Sandroff BM, Motl RW: Quantifying gait abnormalities in persons with multiple sclerosis with minimal disability. Gait Posture 2012,36(1):154–156.PubMedView Article
- Criminisi A, Shotton J, Konukoglu E: Decision Forests for Classification, Regression, Density Estimation, Manifold Learning and Semi-Supervised Learning [Internet]. Microsoft Research 2011. Available from: http://research.microsoft.com/apps/pubs/default.aspx?id=155552
- Mendel J, Wu D: Perceptual Computing: Aiding People in Making Subjective Judgments. Hoboken, New Jersey: John Wiley & Sons, Inc.; 2010:339.View Article
- Kinect for Windows SDK from Microsoft Research [Internet]. [cited 2011 Oct 28]. Available from: http://www.microsoft.com/en-us/kinectforwindowsdev/Downloads.aspx
- Calderita LV, Bandera JP, Bustos P, Skiadopoulos A: Model-based reinforcement of kinect depth data for human motion capture applications. Sensors (Basel) 2013,13(7):8835–8855.View Article
- Clark RA, Pua Y-H, Fortin K, Ritchie C, Webster KE, Denehy L, Bryant AL: Validity of the Microsoft Kinect for assessment of postural control. Gait Posture 2012,36(3):372–377.PubMedView Article
- Lloréns R, Alcañiz M, Colomer C, Navarro MD: Balance recovery through virtual stepping exercises using Kinect skeleton tracking: a follow-up study with chronic stroke patients. Stud Health Technol Inform 2012, 181:108–112.PubMed
- Mobini A, Behzadipour S, Saadat Foumani M: Accuracy of Kinect’s skeleton tracking for upper body rehabilitation applications. Disabil Rehabil Assist Technol 2013. [http://dx.doi.org/10.3109/17483107.2013.805825]
- Khoshelham K, Elberink SO: Accuracy and resolution of kinect depth data for indoor mapping applications. Sensors (Basel) 2012,12(2):1437–1454.View Article
- Gabel M, Gilad-Bachrach R, Renshaw E, Schuster A: Full body gait analysis with Kinect. Conf Proc IEEE Eng Med Biol Soc. 2012, 2012:1964–1967.PubMed
- Stone EE, Skubic M: Capturing habitual, in-home gait parameter trends using an inexpensive depth camera. Conf Proc IEEE Eng Med Biol Soc. 2012, 2012:5106–5109.PubMed
- Polman CH, Reingold SC, Banwell B, Clanet M, Cohen JA, Filippi M, Fujihara K, Havrdova E, Hutchinson M, Kappos L, Lublin FD, Montalban X, O’Connor P, Sandberg-Wollheim M, Thompson AJ, Waubant E, Weinshenker B, Wolinsky JS: Diagnostic criteria for multiple sclerosis: 2010 revisions to the McDonald criteria. Ann Neurol 2011,69(2):292–302.PubMed CentralPubMedView Article
- Cutter GR, Baier ML, Rudick RA, Cookfair DL, Fischer JS, Petkau J, Syndulko V, Weinshenker BG, Antel JP, Confavreux C, Ellison GW, Lublin F, Miller AE, Rao SM, Reingold SM, Thompson A, Willoughby E: Development of a multiple sclerosis functional composite as a clinical trial outcome measure. Brain 1999,122(Pt 5):871–882.PubMedView Article
- Socie MJ, Sosnoff JJ: Gait variability and multiple sclerosis. Mult Scler Int. 2013, 2013:645197.PubMed CentralPubMed
- Stone EE, Skubic M: Passive in-home measurement of stride-to-stride gait variability comparing vision and Kinect sensing. Conf Proc IEEE Eng Med Biol Soc. 2011, 2011:6491–6494.PubMed
- Crenshaw SJ, Royer TD, Richards JG, Hudson DJ: Gait variability in people with multiple sclerosis. Mult Scler 2006,12(5):613–619.PubMedView Article
- Collongues N, Vermersch P: Multiple sclerosis spasticity: “state-of-the-art” questionnaire survey of specialized healthcare professionals. Expert Rev Neurother 2013,13(3 Suppl 1):21–25.PubMedView Article
- Donzé C, De Sèze J: Spasticity and everyday life in multiple sclerosis. Rev Neurol (Paris) 2012,168(3):S51-S56.View Article
- Phan-Ba R, Calay P, Grodent P, Delrue G, Lommers E, Delvaux V: Motor Fatigue Measurement by Distance-Induced Slow Down of Walking Speed in Multiple Sclerosis. PLoS One [Internet]. 2012 Apr 13 [cited 2013 Aug 22];7(4). Available from: http://dx.doi.org/10.1371/journal.pone.0034744
- Phan-Ba R, Calay P, Grodent P, Delrue G, Lommers E, Delvaux V, Moonen G, Nagels G, Belachew S: A corrected version of the Timed-25 Foot Walk Test with a dynamic start to capture the maximum ambulation speed in multiple sclerosis patients. NeuroRehabilitation 2012,30(4):261–266.PubMed
- Fleuren JF, Voerman GE, Snoek GJ, Nene AV, Rietman JS, Hermens HJ: Perception of lower limb spasticity in patients with spinal cord injury. Spinal Cord 2009,47(5):396–400.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.