Next Article in Journal
Design and Metrological Analysis of a Backlit Vision System for Surface Roughness Measurements of Turned Parts
Previous Article in Journal
On Coverage of Critical Nodes in UAV-Assisted Emergency Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Approach to Clustering Accelerometer Data for Application in Passive Predictions of Changes in Depression Severity

1
Department of Psychiatry, University of Illinois at Chicago, Chicago, IL 60612, USA
2
Department of Biomedical Engineering, University of Illinois at Chicago, Chicago, IL 60612, USA
3
Department of Information and Decision Sciences, University of Illinois at Chicago, Chicago, IL 60612, USA
4
Department of Intelligence Computing, Hanyang University, Seoul 04763, Republic of Korea
5
Department of Computing, DePaul University, Chicago, IL 60604, USA
6
Department of Psychology, University of Illinois at Chicago, Chicago, IL 60612, USA
7
Department of Psychiatry, University of Utah, Salt Lake City, UT 84112, USA
*
Author to whom correspondence should be addressed.
Submission received: 9 December 2022 / Revised: 11 January 2023 / Accepted: 29 January 2023 / Published: 1 February 2023
(This article belongs to the Special Issue Smartphone Based Biosensing 2022)

Abstract

:
The treatment of mood disorders, which can become a lifelong process, varies widely in efficacy between individuals. Most options to monitor mood rely on subjective self-reports and clinical visits, which can be burdensome and may not portray an accurate representation of what the individual is experiencing. A passive method to monitor mood could be a useful tool for those with these disorders. Some previously proposed models utilized sensors from smartphones and wearables, such as the accelerometer. This study examined a novel approach of processing accelerometer data collected from smartphones only while participants of the open-science branch of the BiAffect study were typing. The data were modeled by von Mises-Fisher distributions and weighted networks to identify clusters relating to different typing positions unique for each participant. Longitudinal features were derived from the clustered data and used in machine learning models to predict clinically relevant changes in depression from clinical and typing measures. Model accuracy was approximately 95%, with 97% area under the ROC curve (AUC). The accelerometer features outperformed the vast majority of clinical and typing features, which suggested that this new approach to analyzing accelerometer data could contribute towards unobtrusive detection of changes in depression severity without the need for clinical input.

1. Introduction

Mood disorders, such as major depressive disorder and bipolar disorder, often require treatment that is challenging and sometimes a lifelong process [1,2,3]. Traditional methods to monitor mood and other symptoms of these disorders generally rely on self-reports and infrequent clinical visits, which can be time consuming and expensive for the individual and depict an inaccurate portrayal of daily experiences as a result of factors such as recall bias [4,5,6]. While the term mood can convey different meanings, such as momentary mood, which can fluctuate over the course of the day, clinicians instead use these methods to assess the individual’s mood in the context of a disordered state experienced by those with mental illness. As summarized by Hidalgo-Mazzei et al., researchers have begun to incorporate smart technologies into the development of novel methods to monitor mood disorders due in part to the ubiquity of smartphones and wearable devices in the recent years [7]. These devices, which are usually already integrated into the individual’s daily life, contain sensors that have been shown to be capable of unobtrusively detecting future changes in mood [8,9,10,11,12]. The use of technology to supplement traditional treatment approaches has been appealing due to the ability of extracting information on a more granular level than traditional approaches without the need for active input by the user [5]. The passive measures derived are independent of bias and might be a better reflection of everyday life [13].
Previous studies have analyzed a range of data obtained from smartphones and wearable devices, such as GPS location, phone and app usage patterns, voice and ambient noise, and motion sensor information, summarized in several reviews, such as those by Orsolini et al. and Victory et al. in 2020 [14,15]. Through harnessing these data, patterns in an individual’s life could be analyzed to evaluate potential changes related to mood without interrupting daily activities. Since varying moods can lead to differing activity levels [16,17], data recorded from the accelerometer in these devices have been used in combination with other features to passively monitor and predict mood, summarized by Highland and Zhou [18]. Many studies tended to examine metrics from the accelerometer related to the overall movement of the device (e.g., average displacement) as a proxy for the activity level of the individual [10,19,20,21,22,23]. Information about the orientation of the phone was disregarded, possibly due to the inability to identify the activity and location of the phone with respect to the individual without additional information. However, positional information during times of sedentary activities, such as identifying when the individual is laying down versus sitting or standing during different times of the day, could also contribute to passively tracking mood.
One way to obtain this information is by focusing on accelerometer patterns during smartphone keyboard typing. Smartphone typing, which can be performed while sedentary or active, is a common and frequent activity by most smartphone users. The cognitive processes involved while typing are thought to be influenced by mood, which has been supported in previous studies examining the relationship between typing dynamics and symptoms of mood disorders [20,21,24]. Identifying patterns in accelerometer signals during this specific activity could provide additional information related to individuals’ psychological wellbeing. Restricting recording to only during periods of active engagement with the phone in a known orientation allows us to better understand the diurnal patterns of individuals’ phone use and how diurnal changes are related to mood disorders, while also not draining the phone battery from constantly triggering the accelerometer.
In those with mood disorders, diurnal patterns have been found to be disrupted relative to healthy individuals [25,26]. Fluctuations in mood have been noted throughout the day, with mornings generally characterized by worse mood and overall improvement seen as the day progressed [25]. Continued disruptions in these diurnal mood fluctuation patterns may implicate overall deviations in mood [25]. We suspected that sustained deviations from norms in typing position over time might relate to changes in mood. We hypothesized that individuals tend to type on their smartphones in unique but specific orientations depending on their body position, which would be influenced by the time of day and week. Continuous alterations in these positions might be related to overall changes in depression severity. As the conventional clustering algorithms k-means [27], density-based spatial clustering of applications with noise (DBSCAN) [28], and Gaussian mixture models (GMM) were shown to be insufficient in clustering the accelerometer data to identify the preferred phone orientations for all participants, the objective of this study was to develop a novel approach of processing and analyzing accelerometer data longitudinally and to verify their utility to predict clinically relevant changes in depression severity.

2. Materials and Methods

2.1. The BiAffect iPhone Open Science Study

The participants were a part of the open science branch of the BiAffect study and downloaded the BiAffect app from the Apple app store onto their personal iPhones without the requirement of being a part of a controlled study. Included participants comprised of a combination of those recruited for a controlled study, which utilized the BiAffect app, as well as citizen scientists, who downloaded the app of their own accord. Initially developed to predict changes in mood and cognition for those with bipolar disorder, the study aims to understand whether patterns passively detected through smartphone typing behaviors can be used to monitor mood disorder symptomatology, which could aid in symptom management without increasing the burden on the individual. The app provides users with a custom keyboard that records typing and accelerometer metadata, as well as active cognitive tasks, mood surveys, and rating scales for users to complete periodically. All data collected are de-identified. The data have been used in previous BiAffect studies that found relationships between smartphone keyboard typing patterns and mood disorders [12,20,21,22,29,30,31,32,33].
Specifically, this app recorded the category of keypresses (alphanumeric, backspace, punctuation, etc.) and timestamp while the person was typing using the customized keyboard but not the actual text. Accelerometer readings were also recorded at 10 Hz during typing sessions. In addition to recording typing events, participants were prompted weekly to report the Patient Health Questionaire 8 (PHQ without the suicidality item), which is a self-report of depression severity [34]. As usage of the BiAffect keyboard was entirely voluntary, participation was not consistent. For comparison, missing values were accounted for via two methods: (1) imputation and (2) filtering out of individuals with missing data. Generally, models using imputation performed a few percentage points higher than filtering in terms of accuracy and area under the ROC curve (AUC), which is consistent with previous results [12,33]. In total, there were 295 individuals available in this dataset when using imputation, but only 100 left when using filtering.

2.2. Accelerometer Processing

Accelerometer readings were normalized to gravity and recorded as x/y/z coordinates. Readings were filtered to only include coordinates with a magnitude between 0.95 and 1.05 m/s2. This filtering resulted in the inclusion of only typing sessions that occurred while the participant was sedentary (standing, sitting, etc.), which was empirically determined using test data obtained during internal testing in various sedentary and active activities. Data were grouped by week for each participant to account for within-week fluctuations in routine, and the von-Mises Fisher (vMF) distribution was calculated for each group (spherical_kde package, version 0.1.0) [35]. This type of distribution was chosen due to the spherical nature of the accelerometer readings, which resembled the unit sphere following normalization and filtering. The vMF distributions were sampled at 1000 equidistant points across a unit sphere to determine the densities of the distribution at a resolution that captured the shape of the distribution while also being sparse enough to not be too computationally expensive [36].

2.3. Clustering

To cluster the accelerometer readings per week, the number of clusters in each distribution was first determined through identifying the number of local maxima in the vMF distributions. We reasoned that the peaks in the distribution corresponded to the locations where the majority of the accelerometer points, as well as, subsequently, the number of clusters, lied. Local maxima were defined as points sampled from the vMF distribution with an associated density larger than the 8 neighboring sampled points, with neighbors defined by the number of directly surrounding equidistant sampled points on the vMF distribution. The local maxima with an associated density below a set threshold were attributed to noise and discarded. The threshold was set according to the second bin value of the vMF distribution’s histogram with fifty bins, since the first bin of the histogram contained the portions of the vMF distribution with minimal to no accelerometer readings. The cluster centers were labeled according to the coordinates of the local maxima on the sampled unit sphere.
The common methods to cluster the accelerometer data we chose to compare to our network graph-based method were spherical k-means, DBSCAN, and GMM. The Scikit-Learn package in Python (version 0.21.3) [37] was used to cluster the accelerometer data for DBSCAN and GMM, and the modification of the Scikit-Learn k-means function by the spherecluster package (version 0.1.7) [38] was used for the spherical k-means method. The number of clusters was determined for the spherical k-means and GMM methods using the vMF distribution-based method described above, and the distance metric used for the DBSCAN method was cosine distance due to the spherical nature of the data. All other parameters were left as default.
For our network graph-based method of clustering (shown in Figure 1), the accelerometer readings were assigned to a cluster using graph distance. First, an adjacency matrix was constructed for the equidistant sphere points sampled from the vMF distribution. The edges were weighted using an average of the density sampled from the neighboring points (vMFi,j), shown by Equation (1).
w ij = e ( vMF i + vMF j ) 2
This weighting was used in place of the distance between the neighboring points, since the points were all equally spaced across the unit sphere. Using the average density between two nodes to weight the edge created a graph in which nodes located in a high-density region of the vMF distribution were close together, while nodes located in low density regions were farther apart.
A network graph was constructed using the weighted adjacency matrix (networkx package, version 2.6.3) [39], and the sampled sphere points were assigned to a cluster using Dijkstra’s shortest path algorithm (illustrated in Figure 2). All accelerometer points were matched to the closest sampled sphere point using a nearest neighbor algorithm (scipy.spatial.cKDTree package, version 1.3.1) [40] and assigned the corresponding cluster label.
All processing was conducted in Python, version 3.7.4 [41], using the pandas package (version 1.2.0) [42,43] and NumPy package (version 1.17.2) [44]. Plots were constructed using the matplotlib package, version 3.5.3 [45].

2.4. Modeling

Since the accelerometer was selectively recorded only during periods of typing activity, the readings possible were reduced due to the generally limited phone orientations feasible while typing. With that restriction, information about the orientation of the phone was extracted from processed accelerometer readings for longitudinal analysis. This information was possible to infer from the accelerometer alone, since the data were filtered to only include readings of no accelerations aside from that with respect to gravity. By observing the projection of gravity onto the three axes, the relative orientations of the phones were deduced. Orientation and time-based variables calculated from the clustering of weekly accelerometer readings (Table 1) were used, along with typing features, which captured typing speed and variability, as well as clinical factors developed by Bennett et al. [12], to predict clinically relevant changes in PHQ score (difference of 4 or more) [46] using random forest, gradient boosting, and deep learning neural networks methods. Since the majority of the participants included in the dataset did not experience a clinically relevant change in depression, the imbalances were adjusted using synthetic minority oversampling technique (SMOTE) [47]. Feature rankings were determined by a filter-based feature selection method using a random forest model in the python package Scikit-Learn, based on information gain [48]. Odds ratios were calculated in Excel.

3. Results

3.1. Clustering

The accelerometer data collected was clustered to identify the predominant orientations each participant held their phone while typing. The clustering methods k-means, DBSCAN, and GMM were tested on the accelerometer data itself and generally performed well, but were found to not reliably cluster the points in the expected way that resembled the associated vMF distribution, as shown through two users’ data in Figure 3. For the majority of the participants data, DBSCAN labeled all of the accelerometer points as belonging to one cluster and so was not deemed a feasible method to cluster the data. The spherical k-means algorithm clustered the data evenly between the identified cluster centers and did not account for the density of points when labeling clusters. As seen in Figure 3A, the spherical k-means method split grouped points into separate clusters in a linear fashion. Additionally, unconnected points were sometimes grouped into one cluster (Figure 3B). GMM performed somewhat similarly overall to spherical k-means but tended to group all sparsely located accelerometer points as belonging to one cluster, rather than grouping those points to the nearest cluster of dense accelerometer points (Figure 3A).
Since the common clustering algorithms did not consistently cluster the accelerometer points across different participants and weekly groupings, we developed a new approach to label the points by cluster. Compared to the conventional clustering methods, clustering of the accelerometer data performed using the respective vMF distribution and network graph was found to be the most effective in correctly identifying the number of clusters and labeling the accelerometer readings to the appropriate cluster, as determined through visual inspection (Figure 3). By creating a network graph that resembled the vMF distribution of the accelerometer points, the clusters were mapped to the sampled points from the vMF distribution using the local maxima and weighted graph as a guide. Cluster labels for the accelerometer points themselves were then transferred from the sampled points of the vMF distribution using the nearest neighbor algorithm. This method accommodated irregular and inconsistent cluster shapes between participants’ data by clustering based on the individual network graph. Unlike with the DBSCAN algorithm, separate clusters were identified within regions of overall higher density of accelerometer points within the distribution. Moreover, clusters were allowed to encompass the entire area of higher density, regardless of cluster shape and distance from the cluster center.
The clusters Identified represented the phone orientations that each participant used each week. For example, in Figure 4, one participant held their phone in multiple orientations throughout the course of a week, ranging from deviations of upright (gray, beige, and green clusters) to facing upwards (orange cluster) to horizontally (teal, yellow, and blue clusters).

3.2. Changes in Phone Orientation over Time

The labeled accelerometer data was used to identify consistency or changes in typing orientation over time that were unique to each individual. The predominant cluster label per hour was identified and plotted to analyze shifts in typing orientation over time, shown in Figure 5. Some participants had consistent typing orientations, shown by the zero distance from previous cluster on the plot, while other participants shifted typing orientation regularly throughout the day and weeks, which was visualized by frequent changes in the distance from previous cluster on the plot. Larger distances between consecutive clusters corresponded to a more drastic phone orientation change, while smaller distances corresponded to only slight changes in phone orientation. In Figure 5, participant A shifted typing orientation frequently for the majority of the time using the BiAffect keyboard, but spent one week typing only in one orientation. Participant B, on the other hand, predominantly typed in one orientation over the course of a month.
Accelerometer features were calculated from the clustered data to determine their efficacy in predicting clinically relevant changes in PHQ score when combined with typing and clinical features. Model accuracy ranged from 94 to 95.5% (2% standard deviation), with 96 to 98% AUC, similar to that reported in [33]. Importantly, the accelerometer features calculated from the cluster labels contributed greatly to model predictions, even performing better than some otherwise highly ranked clinical features previously reported in the literature [12,49], as can be seen in Table 2.
To understand the directionality of how these PHQ changes were affected by the features, odds ratios were also calculated. The odds ratios for number of clusters and cluster transitions were 1.26 and 1.2, respectively, which suggested that participants with a larger number of clusters and cluster transitions per week had a higher probability of a clinically relevant change in PHQ.

4. Discussion

Understanding patterns in a person’s phone orientation while typing could uncover information about changes in their depression severity. While many studies have primarily examined the movement of the phone or wearable device compiled over all activities in analyses using the accelerometer [10,19,20,21,22,23], we sought to determine how the orientation of the phone specifically during typing related to changes in depression severity instead. By limiting the recording times to only during smartphone typing, we were exploring the utility of passive tracking during activities that individuals generally already do on a daily basis while also reducing the toll on the smartphone’s battery. In this study, we developed a novel method of processing accelerometer data to discover how personalized features related to phone orientation can predict clinically relevant changes in depression.
To identify individualized phone orientation tendencies that each participant preferred, clustering of the accelerometer data was performed on a weekly basis to account for day-to-day fluctuations in participants’ schedules due to work and other activities. After first testing the efficacy of conventional clustering methods, we discovered that these algorithms did not consistently provide adequate identification of the different clusters present for all of the participants. As shown in Figure 3, spherical k-means, DBSCAN, and GMM were applied to the accelerometer data and compared to one another.
Although DBSCAN did not require prior knowledge about the number of clusters in the data unlike the other methods used, the algorithm labeled clusters based on the density of the points. Since the accelerometer points for the vast majority of the participants were not well separated between clusters, all points were labeled as belonging to one large cluster, and no distinctions were made between areas of higher and lower densities within the distribution, shown in Figure 3. This method might have worked better if there were clear distinctions between the phone orientations used by the participants, but the overall spread of the points prevented any clear separations.
The spherical k-means algorithm performed well for many of the participants’ data and could reasonable replicate the vMF distributions overall. However, the method works by partitioning the points into clusters, such that each point is labeled as belonging to the nearest cluster center, independent of its location in the distribution [27]. The clusters formed when using this method follow a circular shape, which was not accurate for every instance in our data. Using this method, points that appeared to lie on the edges of a larger cluster might instead be assigned to another cluster solely due to the distance to each cluster center instead of taking into account the shape of the distribution. Moreover, points that were physically separated were sometimes grouped together due to their distances from the cluster centers (Figure 3B).
GMM was an improvement from the spherical k-means algorithm since the algorithm could handle non-circular cluster shapes, but did not seem to handle the sparse points well likely due to noise in the data, as seen in Figure 3A. The distribution of the accelerometer data did not follow a normal distribution, which might have contributed to the subpar performance of the algorithm on the participants’ data overall.
Since the conventional clustering methods tested did not appropriately identify and label the accelerometer data into clusters for all participants, a new approach was developed, outlined in Figure 1 and Figure 2. At the step of clustering the data, the points had already been modeled by vMF distributions in order to determine the number of clusters in each distribution. We then used this representation of the data to create a customized mapping of each distribution for the clustering to preserve their unique characteristics. The clustering using this method was able to accommodate irregular and inconsistent cluster shapes across data from different participants, as well as take into account the varying densities in the vMF distributions when forming the clusters, shown in Figure 3. These clusters represented the different phone orientations participants used each week and depending on the data could suggest several different corresponding body positions, ranging from standing or sitting upright to lounging or laying down, as shown through one participant’s data in Figure 4.
Accelerometer features were then designed as a proxy for body positioning while typing to provide more information about the individual’s environment in models to predict changes in depression severity (Table 1). It is well known that sleep and behavior can be disrupted during periods of depression [50,51,52]. Therefore, we investigated patterns derived from the clustered accelerometer data over time to reveal information about the state of the individual. Reinersten et al. outlined several studies describing that changes in activity can be indicative of changes in depression [16], so we speculated that individuals who typed inconsistent to their usual typing patterns might show signs of a change in depressed mood relative to their previous state. As seen in the comparison between two participants in Figure 5, the patterns of movement between clusters throughout days and weeks ranged in consistency between participants and within participants over time. One participant typed in one orientation over the course of a month aside from a few instances generally in the mornings and nights in which the shift in phone orientation was drastic due to the large distance between consecutive clusters. On the other hand, the other participant had regular shifts in phone orientation throughout the day for many weeks, with some being minor shifts in phone orientation (small distances between consecutive clusters) and others being major shifts in phone orientation (larger distances between consecutive clusters), which suggested that their body position changed from an upright position to variations of lounging or laying down multiple times throughout the day on a regular basis. The shift to consistent phone orientation in week 4 suggested a change in the participant’s behavior during that time, which might have been the result of a change in depression severity. We derived features to capture this information in order to further investigate the relationship between participants’ chosen phone orientation while typing and changes in depression severity. We suspected that how often and the degree to which phone orientations (and therefore body positions) changed over time would be related to participants’ depression severity, with major shifts being more indicative of a change.
We constructed models that examined typing dynamics and clinical measures to predict clinically relevant changes in depression severity in order to evaluate our newly-derived features. Model accuracy was around 95%, and importantly, we observed that the features derived from clustered accelerometer data were very highly ranked in feature importance (Table 2). This ranking suggested that these features, which captured information about the participants’ position while typing, were just as, if not more, important than demographic and clinical information in predicting whether a participant would have a clinically relevant change in depression severity the following week. Moreover, odds ratios, which inform on the direction of change, suggested that the higher number of clusters (i.e., number of typing orientations per week) and number of cluster transitions (i.e., number of times the participant changed between typing positions per week) resulted in an increased likelihood that the participant would experience a change in depression severity the following week. This finding suggested that participants who chose multiple body positions while typing on their smartphone and changed often between them were more likely to experience a fluctuation in their mood, consistent with previous studies evaluating the relationship between psychomotor disturbances, circadian rhythm disruptions, and depression [53,54].
Conventionally, assessments of mood to evaluate treatment efficacy and symptom management rely heavily on a person’s ability to accurately recall their experiences leading up to clinical visits, which can be sparse or difficult to access [6]. These recollections have been previously noted in the literature to often be subject to recall bias [5,55], making regular and accurate monitoring of symptoms difficult. Our approach to processing and analyzing accelerometer data collected during smartphone typing can help to better understand the signals present in this modality and uncover their relationship to mood. Pending further work, models using this approach could be more advantageous in tracking changes in mood due to the greater reliance on passive measures recorded on personal smartphones already used by the majority of the population. The increased granularity of the recorded data and input that does not necessitate access to medical professionals could benefit many individuals by providing objective supplemental information during clinical visits and potentially serve to identify early signs of changes in mood that otherwise might be recognized too late [15,56]. Additional work to investigate this implementation needs to be conducted, though, to determine the extent of beneficial effects of feedback about depression severity, as well as to identify any potentially harmful side effects due to incorrect predictions of changes in depression severity or other factors [57].
This analysis, however, does not come without limitations. First, as the data analyzed belongs to the open-science branch of the BiAffect study, the demographic and clinical information submitted by the participants is not verified by a psychologist and so might not be as accurate as if the data was obtained through a controlled study. Even though the adherence to using the BiAffect keyboard was less consistent as a result of the nature of the participation, the ease to participate facilitated the recruitment and volume of participants.
Furthermore, for the analysis, we elected to filter the data to only include accelerometer readings when the phone was not accelerating independent of gravity, which likely excluded typing sessions while participants were walking or otherwise moving. This exclusion might have impacted the results, but since we did not have any information from other sensors, we would have been unable to deduce the rationale for the accelerations and incorporate them into the models. Future directions could explore the inclusion of this data and the impact on predictions of depression severity.

5. Conclusions

In recent years, focus has been placed on more effective and efficient methods to monitor treatment and progression of mood disorder symptoms in order to relieve the burden on individuals. Passive and unobtrusive measures obtained from smartphones and wearable devices already incorporated into most people’s everyday lives have become targets for models to predict future changes in symptomatology. The accelerometer, which is found in most smart devices, has become a popular choice due in part to the ease of information that can be gathered passively without much infringement upon the privacy of the individual compared to other possible information gathered from these devices. This study developed a novel approach of processing accelerometer data that has the potential to augment predictions of changes in depression severity with less dependence on clinical input by the individual.

Author Contributions

Conceptualization, M.K.R., T.T., F.H., A.P.D., E.N., S.A.L., O.A. and A.D.L.; Data curation, M.K.R.; Formal analysis, M.K.R., C.C.B., E.B. and D.K.; Funding acquisition, O.A. and A.D.L.; Investigation, M.K.R.; Methodology, M.K.R., T.T., C.C.B., E.B., D.K. and A.D.L.; Software, A.D.L.; Supervision, T.T. and A.D.L.; Validation, M.K.R. and T.T.; Visualization, M.K.R.; Writing—original draft, M.K.R.; Writing—review and editing, M.K.R., T.T., C.C.B., F.H., A.P.D., E.N., S.A.L., O.A. and A.D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partially funded by Mood Challenge for Research kit and NIMH, grant number 1R01MH120168.

Institutional Review Board Statement

The study involving human participants was reviewed and approved by University of Illinois at Chicago Institutional Review Board (IRB approval number 2016-1261).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to privacy concerns.

Acknowledgments

We thank the citizen scientists who elected to participate in and contribute their data to the BiAffect study.

Conflicts of Interest

Alex Leow is a cofounder of KeyWise AI, has served as a consultant for Otsuka US, and is currently on the medical board of Buoy Health. Olusola Ajilore is a cofounder of KeyWise AI and is on the advisory boards of Embodied Labs, Blueprint Health, and Sage Therapeutics. The other authors have no conflict to report.

References

  1. Frank, E.; Cassano, G.B.; Rucci, P.; Thompson, W.K.; Kraemer, H.C.; Fagiolini, A.; Maggi, L.; Kupfer, D.J.; Shear, M.K.; Houck, P.R.; et al. Predictors and Moderators of Time to Remission of Major Depression with Interpersonal Psychotherapy and SSRI Pharmacotherapy. Psychol. Med. 2011, 41, 151–162. [Google Scholar] [CrossRef] [PubMed]
  2. Alang, S.; McAlpine, D. Treatment Modalities and Perceived Effectiveness of Treatment among Adults with Depression. Health Serv. Insights 2020, 13, 1178632920918288. [Google Scholar] [CrossRef]
  3. McIntyre, R.S.; Berk, M.; Brietzke, E.; Goldstein, B.I.; López-Jaramillo, C.; Kessing, L.V.; Malhi, G.S.; Nierenberg, A.A.; Rosenblat, J.D.; Majeed, A.; et al. Bipolar Disorders. Lancet 2020, 396, 1841–1856. [Google Scholar]
  4. Patel, A. The Cost of Mood Disorders. Psychiatry 2009, 8, 76–80. [Google Scholar] [CrossRef]
  5. Hidalgo-Mazzei, D.; Young, A.H. Psychiatry Foretold. Aust. N. Z. J. Psychiatry 2019, 53, 365–366. [Google Scholar] [CrossRef] [PubMed]
  6. Sheikh, M.; Qassem, M.; Kyriacou, P.A. Wearable, Environmental, and Smartphone-Based Passive Sensing for Mental Health Monitoring. Front. Digit. Health 2021, 3, 662811. [Google Scholar] [CrossRef] [PubMed]
  7. Hidalgo-Mazzei, D.; Llach, C.; Vieta, E. MHealth in Affective Disorders: Hype or Hope? A Focused Narrative Review. Int. Clin. Psychopharmacol. 2020, 35, 61–68. [Google Scholar] [CrossRef] [PubMed]
  8. Canzian, L.; Musolesi, M. Trajectories of Depression: Unobtrusive Monitoring of Depressive States by Means of Smartphone Mobility Traces Analysis. In Proceedings of the 2015 ACM International Joint Conference on Pervasive and Ubiquitous Computing, Osaka, Japan, 7–11 September 2015; Association for Computing Machinery: New York, NY, USA, 2015; pp. 1293–1304. [Google Scholar]
  9. Saeb, S.; Lattie, E.G.; Schueller, S.M.; Kording, K.P.; Mohr, D.C. The Relationship between Mobile Phone Location Sensor Data and Depressive Symptom Severity. PeerJ 2016, 4, e2537. [Google Scholar] [CrossRef] [PubMed]
  10. Beiwinkel, T.; Kindermann, S.; Maier, A.; Kerl, C.; Moock, J.; Barbian, G.; Rössler, W. Using Smartphones to Monitor Bipolar Disorder Symptoms: A Pilot Study. JMIR Ment. Health 2016, 3, e4560. [Google Scholar] [CrossRef]
  11. Cho, C.-H.; Lee, T.; Kim, M.-G.; In, H.P.; Kim, L.; Lee, H.-J. Mood Prediction of Patients With Mood Disorders by Machine Learning Using Passive Digital Phenotypes Based on the Circadian Rhythm: Prospective Observational Cohort Study. J. Med. Internet Res. 2019, 21, e11029. [Google Scholar] [CrossRef] [PubMed]
  12. Bennett, C.C.; Ross, M.K.; Baek, E.; Kim, D.; Leow, A.D. Predicting Clinically Relevant Changes in Bipolar Disorder Outside the Clinic Walls Based on Pervasive Technology Interactions via Smartphone Typing Dynamics. Pervasive Mob. Comput. 2022, 83, 101598. [Google Scholar] [CrossRef]
  13. Insel, T.R. Digital Phenotyping: A Global Tool for Psychiatry. World Psychiatry 2018, 17, 276–277. [Google Scholar] [CrossRef] [PubMed]
  14. Orsolini, L.; Fiorani, M.; Volpe, U. Digital Phenotyping in Bipolar Disorder: Which Integration with Clinical Endophenotypes and Biomarkers? Int. J. Mol. Sci. 2020, 21, 7684. [Google Scholar] [CrossRef] [PubMed]
  15. Victory, A.; Letkiewicz, A.; Cochran, A.L. Digital Solutions for Shaping Mood and Behavior among Individuals with Mood Disorders. Curr. Opin. Syst. Biol. 2020, 21, 25–31. [Google Scholar] [CrossRef]
  16. Reinertsen, E.; Clifford, G.D. A Review of Physiological and Behavioral Monitoring with Digital Sensors for Neuropsychiatric Illnesses. Physiol. Meas. 2018, 39, 05TR01. [Google Scholar] [CrossRef]
  17. Vancampfort, D.; Firth, J.; Schuch, F.B.; Rosenbaum, S.; Mugisha, J.; Hallgren, M.; Probst, M.; Ward, P.B.; Gaughran, F.; De Hert, M.; et al. Sedentary Behavior and Physical Activity Levels in People with Schizophrenia, Bipolar Disorder and Major Depressive Disorder: A Global Systematic Review and Meta-Analysis. World Psychiatry 2017, 16, 308–315. [Google Scholar] [CrossRef]
  18. Highland, D.; Zhou, G. A Review of Detection Techniques for Depression and Bipolar Disorder. Smart Health 2022, 24, 100282. [Google Scholar] [CrossRef]
  19. Bardram, J.E.; Frost, M.; Szántó, K.; Faurholt-Jepsen, M.; Vinberg, M.; Kessing, L.V. Designing Mobile Health Technology for Bipolar Disorder: A Field Trial of the Monarca System. In Proceedings of the SIGCHI Conference on Human Factors in Computing Systems, Paris, France, 27 April–2 May 2013; Association for Computing Machinery: New York, NY, USA, 2013; pp. 2627–2636. [Google Scholar]
  20. Cao, B.; Zheng, L.; Zhang, C.; Yu, P.S.; Piscitello, A.; Zulueta, J.; Ajilore, O.; Ryan, K.; Leow, A.D. DeepMood: Modeling Mobile Phone Typing Dynamics for Mood Detection. In Proceedings of the 23rd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, Halifax, NS, Canada, 13–17 August 2017; Association for Computing Machinery: Halifax, NS, Canada, 2017; pp. 747–755. [Google Scholar]
  21. Zulueta, J.; Piscitello, A.; Rasic, M.; Easter, R.; Babu, P.; Langenecker, S.A.; McInnis, M.; Ajilore, O.; Nelson, P.C.; Ryan, K.; et al. Predicting Mood Disturbance Severity with Mobile Phone Keystroke Metadata: A BiAffect Digital Phenotyping Study. J. Med. Internet Res. 2018, 20. [Google Scholar] [CrossRef]
  22. Huang, H.; Cao, B.; Yu, P.S.; Wang, C.; Leow, A.D. DpMood: Exploiting Local and Periodic Typing Dynamics for Personalized Mood Prediction. In Proceedings of the 2018 IEEE International Conference on Data Mining (ICDM), Singapore, 17–20 November 2018; pp. 157–166. [Google Scholar]
  23. Spathis, D.; Servia-Rodriguez, S.; Farrahi, K.; Mascolo, C.; Rentfrow, J. Passive Mobile Sensing and Psychological Traits for Large Scale Mood Prediction. In Proceedings of the 13th EAI International Conference on Pervasive Computing Technologies for Healthcare, Trento, Italy, 20–23 May 2019; Association for Computing Machinery: New York, NY, USA, 2019; pp. 272–281. [Google Scholar]
  24. Mastoras, R.-E.; Iakovakis, D.; Hadjidimitriou, S.; Charisis, V.; Kassie, S.; Alsaadi, T.; Khandoker, A.; Hadjileontiadis, L.J. Touchscreen Typing Pattern Analysis for Remote Detection of the Depressive Tendency. Sci. Rep. 2019, 9, 13414. [Google Scholar] [CrossRef] [PubMed]
  25. Peeters, F.; Berkhof, J.; Delespaul, P.; Rottenberg, J.; Nicolson, N.A. Diurnal Mood Variation in Major Depressive Disorder. Emotion 2006, 6, 383–391. [Google Scholar] [CrossRef]
  26. Bechtel, W. Circadian Rhythms and Mood Disorders: Are the Phenomena and Mechanisms Causally Related? Front. Psychiatry 2015, 6, 118. [Google Scholar] [CrossRef] [Green Version]
  27. Hartigan, J.A.; Wong, M.A. Algorithm AS 136: A K-Means Clustering Algorithm. J. R. Stat. Soc. Ser. C Appl. Stat. 1979, 28, 100–108. [Google Scholar] [CrossRef]
  28. Ester, M.; Kriegel, H.-P.; Sander, J.; Xu, X. A Density-Based Algorithm for Discovering Clusters in Large Spatial Databases with Noise. In Proceedings of the Second International Conference on Knowledge Discovery and Data Mining, Portland, Oregon, USA, 2 August 1996; AAAI Press: Portland, OR, USA, 1996; pp. 226–231. [Google Scholar]
  29. Stange, J.P.; Zulueta, J.; Langenecker, S.A.; Ryan, K.A.; Piscitello, A.; Duffecy, J.; McInnis, M.G.; Nelson, P.; Ajilore, O.; Leow, A. Let Your Fingers Do the Talking: Passive Typing Instability Predicts Future Mood Outcomes. Bipolar Disord. 2018, 20, 285–288. [Google Scholar] [CrossRef] [PubMed]
  30. Vesel, C.; Rashidisabet, H.; Zulueta, J.; Stange, J.P.; Duffecy, J.; Hussain, F.; Piscitello, A.; Bark, J.; Langenecker, S.A.; Young, S.; et al. Effects of Mood and Aging on Keystroke Dynamics Metadata and Their Diurnal Patterns in a Large Open-Science Sample: A BiAffect IOS Study. J. Am. Med. Inform. Assoc. 2020, 27, 1007–1018. [Google Scholar] [CrossRef]
  31. Zulueta, J.; Demos, A.P.; Vesel, C.; Ross, M.; Piscitello, A.; Hussain, F.; Langenecker, S.A.; McInnis, M.; Nelson, P.; Ryan, K.; et al. The Effects of Bipolar Disorder Risk on a Mobile Phone Keystroke Dynamics Based Biomarker of Brain Age. Front. Psychiatry 2021, 12, 2284. [Google Scholar] [CrossRef] [PubMed]
  32. Hussain, F.; Stange, J.P.; Langenecker, S.A.; McInnis, M.G.; Zulueta, J.; Piscitello, A.; Ross, M.K.; Demos, A.P.; Vesel, C.; Rashidisabet, H.; et al. Passive Sensing of Affective and Cognitive Functioning in Mood Disorders by Analyzing Keystroke Kinematics and Speech Dynamics. In Digital Phenotyping and Mobile Sensing: New Developments in Psychoinformatics; Montag, C., Baumeister, H., Eds.; Studies in Neuroscience, Psychology and Behavioral Economics; Springer International Publishing: Cham, Switzerland, 2023; pp. 229–258. ISBN 978-3-030-98546-2. [Google Scholar]
  33. Bennett, C.C.; Ross, M.K.; Baek, E.; Kim, D.; Leow, A.D. Smartphone Accelerometer Data as a Proxy for Clinical Data in Modeling of Bipolar Disorder Symptom Trajectory. NPJ Digit. Med. 2022, 5, 1–10. [Google Scholar] [CrossRef] [PubMed]
  34. Löwe, B.; Unützer, J.; Callahan, C.M.; Perkins, A.J.; Kroenke, K. Monitoring Depression Treatment Outcomes with the Patient Health Questionnaire-9. Med. Care 2004, 42, 1194–1201. [Google Scholar] [PubMed]
  35. Handley, W. Spherical_kde: Spherical Kernel Density Estimation; 2020. Available online: https://github.com/williamjameshandley/spherical_kde. (accessed on 9 December 2022).
  36. Deserno, M. How to Generate Equidistributed Points on the Surface of a Sphere. Polym. Ed 2004, 99, 2. [Google Scholar]
  37. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-Learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  38. Laska, J. Clustering on the Unit Hypersphere in Scikit-Learn; 2018. Available online: https://github.com/jasonlaska/spherecluster. (accessed on 9 December 2022).
  39. Hagberg, A.A.; Schult, D.A.; Swart, P.J. Exploring Network Structure, Dynamics, and Function Using NetworkX. In Proceedings of the 7th Python in Science Conference, Pasadena, CA, USA, 19–24 August 2008; Varoquaux, G., Vaught, T., Millman, J., Eds.; pp. 11–15. [Google Scholar]
  40. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [PubMed]
  41. Van Rossum, G.; Drake, F.L. Python 3 Reference Manual; CreateSpace: Scotts Valley, CA, USA, 2009; ISBN 1-4414-1269-7. [Google Scholar]
  42. McKinney, W. Data Structures for Statistical Computing in Python. In Proceedings of the 9th Python in Science Conference, Austin, TX, USA, 28 June–3 July 2010; pp. 56–61. [Google Scholar]
  43. Reback, J.; Jbrockmendel; McKinney, W.; den Bossche, J.V.; Augspurger, T.; Cloud, P.; Hawkins, S.; Gfyoung; Roeschke, M.; Sinhrks; et al. Pandas-Dev/Pandas: Pandas 1.3.3. Zenodo, 2021. Available online: https://zenodo.org/record/5501881#.Y9oAFHZByUk (accessed on 28 January 2023).
  44. Harris, C.R.; Millman, K.J.; van der Walt, S.J.; Gommers, R.; Virtanen, P.; Cournapeau, D.; Wieser, E.; Taylor, J.; Berg, S.; Smith, N.J.; et al. Array Programming with NumPy. Nature 2020, 585, 357–362. [Google Scholar] [CrossRef] [PubMed]
  45. Hunter, J.D. Matplotlib: A 2D Graphics Environment. Comput. Sci. Eng. 2007, 9, 90–95. [Google Scholar] [CrossRef]
  46. Kroenke, K.; Spitzer, R.L. The PHQ-9: A New Depression Diagnostic and Severity Measure. Psychiatr. Ann. 2002, 32, 509–515. [Google Scholar] [CrossRef] [Green Version]
  47. Chawla, N.V.; Bowyer, K.W.; Hall, L.O.; Kegelmeyer, W.P. SMOTE: Synthetic Minority Over-Sampling Technique. J. Artif. Intell. Res. 2002, 16, 321–357. [Google Scholar] [CrossRef]
  48. Bolón-Canedo, V.; Sánchez-Maroño, N.; Alonso-Betanzos, A.; Benítez, J.M.; Herrera, F. A Review of Microarray Datasets and Applied Feature Selection Methods. Inf. Sci. 2014, 282, 111–135. [Google Scholar] [CrossRef]
  49. de Jonge, P.J.; Wardenaar, K.J.; Lim, C.C.W.; Aguilar-Gaxiola, S.; Alonso, J.; Andrade, L.H.; Bunting, B.; Chatterji, S.; Ciutan, M.; Gureje, O.; et al. The Cross-National Structure of Mental Disorders: Results from the World Mental Health Surveys. Psychol. Med. 2018, 48, 2073–2084. [Google Scholar] [CrossRef]
  50. Benca, R.M.; Okawa, M.; Uchiyama, M.; Ozaki, S.; Nakajima, T.; Shibui, K.; Obermeyer, W.H. Sleep and Mood Disorders. Sleep Med. Rev. 1997, 1, 45–56. [Google Scholar] [CrossRef] [PubMed]
  51. Peterson, M.J.; Benca, R.M. Sleep in Mood Disorders. Sleep Med. Clin. 2008, 3, 231–249. [Google Scholar] [CrossRef]
  52. Schuch, F.; Vancampfort, D.; Firth, J.; Rosenbaum, S.; Ward, P.; Reichert, T.; Bagatini, N.C.; Bgeginski, R.; Stubbs, B. Physical Activity and Sedentary Behavior in People with Major Depressive Disorder: A Systematic Review and Meta-Analysis. J. Affect. Disord. 2017, 210, 139–150. [Google Scholar] [CrossRef]
  53. Scott, J. Clinical Parameters of Circadian Rhythms in Affective Disorders. Eur. Neuropsychopharmacol. 2011, 21, S671–S675. [Google Scholar] [CrossRef]
  54. Smit, A.C.; Snippe, E. Real-Time Monitoring of Increases in Restlessness to Assess Idiographic Risk of Recurrence of Depressive Symptoms. Psychol. Med. 2022, 1–10. [Google Scholar]
  55. Ben-Zeev, D.; Young, M.A.; Madsen, J.W. Retrospective Recall of Affect in Clinically Depressed Individuals and Controls. Cogn. Emot. 2009, 23, 1021–1040. [Google Scholar] [CrossRef]
  56. Hidalgo-Mazzei, D.; Young, A.H.; Vieta, E.; Colom, F. Behavioural Biomarkers and Mobile Mental Health: A New Paradigm. Int. J. Bipolar Disord. 2018, 6, 9. [Google Scholar] [CrossRef] [PubMed]
  57. Maher, N.A.; Senders, J.T.; Hulsbergen, A.F.C.; Lamba, N.; Parker, M.; Onnela, J.-P.; Bredenoord, A.L.; Smith, T.R.; Broekman, M.L.D. Passive Data Collection and Use in Healthcare: A Systematic Review of Ethical Issues. Int. J. Med. Inform. 2019, 129, 242–247. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Diagram of steps to cluster accelerometer data.
Figure 1. Diagram of steps to cluster accelerometer data.
Sensors 23 01585 g001
Figure 2. Each of the 1000 equidistant points (that jointly discretized the unit sphere) is assigned to one of the local maxima by following a path of increasing kernal density (i.e., a gradient ascent procedure on the kernal density function), as illustrated using the points marked by an X. This is algorithmically implemented by forming a weighted graph followed by computing the shortest path length between points.
Figure 2. Each of the 1000 equidistant points (that jointly discretized the unit sphere) is assigned to one of the local maxima by following a path of increasing kernal density (i.e., a gradient ascent procedure on the kernal density function), as illustrated using the points marked by an X. This is algorithmically implemented by forming a weighted graph followed by computing the shortest path length between points.
Sensors 23 01585 g002
Figure 3. Comparison of conventional clustering methods tested (DBSCAN, spherical k-means, GMM) to our method for two users’ data (A,B) over a week.
Figure 3. Comparison of conventional clustering methods tested (DBSCAN, spherical k-means, GMM) to our method for two users’ data (A,B) over a week.
Sensors 23 01585 g003
Figure 4. One participant’s accelerometer data over a week labeled by cluster and plotted in the (A) xz axes, (B) xy axes, and (C) yz axes.
Figure 4. One participant’s accelerometer data over a week labeled by cluster and plotted in the (A) xz axes, (B) xy axes, and (C) yz axes.
Sensors 23 01585 g004
Figure 5. Two participants’ plots of the haversine distance between consecutive cluster centers per hour over the weeks using the BiAffect keyboard.
Figure 5. Two participants’ plots of the haversine distance between consecutive cluster centers per hour over the weeks using the BiAffect keyboard.
Sensors 23 01585 g005
Table 1. Accelerometer features calculated using clustered accelerometer data.
Table 1. Accelerometer features calculated using clustered accelerometer data.
Feature NameFeature Description
n_clustersNumber of distinct phone typing orientations per week
total_distance_between_clustersSum of haversine distances traveled between session cluster labels per week
avg_n_clusters_perSessionAverage number of cluster labels of x/y/z readings within each session per week
avg_n_transitions_perSessionAverage number of changes between consecutive cluster labels of x/y/z readings within each session per week
n_cluster_transitionsNumber of changes between consecutive session cluster labels per week
median_X, median_Y, median_ZMedian x/y/z reading of session’s cluster center per week
sum_X_motion, sum_Y_motion, sum_Z_motionSum of differences between consecutive cluster center x/y/z readings of session cluster labels per week
X_motion_sd, Y_motion_sd, Z_motion_sdStandard deviation of differences between consecutive cluster center x/y/z readings of session cluster labels per week
arc_sumThree-dimensional rotational motion per week (calculated based on session cluster center x/y/z accelerometer readings)
Table 2. Feature rankings based on information gain filter-based method (using a random forest model) [48] with accelerometer features calculated from cluster labels in bold.
Table 2. Feature rankings based on information gain filter-based method (using a random forest model) [48] with accelerometer features calculated from cluster labels in bold.
Feature Ranking
median_Y
median_X
Age
n_clusters
median_Z
BD
n_cluster_transitions
sum_Z_motion
BD_binary
MDQdiag
phoneSize
Gender
Anxiety
medianPressDur
PTSD
sum_Y_motion
sum_X_motion
Depression
OCD
percent_upright_night
medianDistCenter
arc_sum
avg_n_transitions_perSession
ADHD
count_X_horizontal
n_XYZ
avg_n_clusters_perSession
X_motion_sd
NoneOfTheseDiag
total_distance_between_clusters
distToCenterPrevRatioAA
autocorrectRate_wkSD
medianIKD
Avg90PercentileAA
autocorrectRate
AvgVarAB
SubstanceAddictionDisorder
backspaceRate
AvgVarAA
Z_motion_sd
medIKD_wkSD
Avg_nBackspace
percent_upright_afternoon
percent_upright_morning
bkspRate_wkSD
Avg_nAutocorrect
AvgVarBB
percent_upright
SeasonalAffectiveDisorder
nKeypresses
Y_motion_sd
percent_upright_evening
AvgMedAA
Avg_medPressDuration
AvgMAD_AA
Diag_PreferNotAnswer
AvgMedAB
Avg_nAlphanum
Schizophrenia
AvgMedBB
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ross, M.K.; Tulabandhula, T.; Bennett, C.C.; Baek, E.; Kim, D.; Hussain, F.; Demos, A.P.; Ning, E.; Langenecker, S.A.; Ajilore, O.; et al. A Novel Approach to Clustering Accelerometer Data for Application in Passive Predictions of Changes in Depression Severity. Sensors 2023, 23, 1585. https://0-doi-org.brum.beds.ac.uk/10.3390/s23031585

AMA Style

Ross MK, Tulabandhula T, Bennett CC, Baek E, Kim D, Hussain F, Demos AP, Ning E, Langenecker SA, Ajilore O, et al. A Novel Approach to Clustering Accelerometer Data for Application in Passive Predictions of Changes in Depression Severity. Sensors. 2023; 23(3):1585. https://0-doi-org.brum.beds.ac.uk/10.3390/s23031585

Chicago/Turabian Style

Ross, Mindy K., Theja Tulabandhula, Casey C. Bennett, EuGene Baek, Dohyeon Kim, Faraz Hussain, Alexander P. Demos, Emma Ning, Scott A. Langenecker, Olusola Ajilore, and et al. 2023. "A Novel Approach to Clustering Accelerometer Data for Application in Passive Predictions of Changes in Depression Severity" Sensors 23, no. 3: 1585. https://0-doi-org.brum.beds.ac.uk/10.3390/s23031585

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop