Next Article in Journal
Information Flow Network of International Exchange Rates and Influence of Currencies
Previous Article in Journal
Encoding Individual Source Sequences for the Wiretap Channel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interactive System for Similarity-Based Inspection and Assessment of the Well-Being of mHealth Users

1
Knowledge Management and Discovery Lab, Otto-von-Guericke University, 39106 Magdeburg, Germany
2
Institute of Clinical Epidemiology and Biometry, University of Würzburg, 97078 Würzburg, Germany
3
Institute of Databases and Information Systems, Ulm University, 89081 Ulm, Germany
4
Institute DigiHealth, Neu-Ulm University of Applied Sciences, 89231 Neu-Ulm, Germany
5
WSAudiology, Sivantos GmbH, 91058 Erlangen, Germany
6
Department of Psychiatry and Psychotherapy, University of Regensburg, 93053 Regensburg, Germany
*
Authors to whom correspondence should be addressed.
Submission received: 26 October 2021 / Revised: 13 December 2021 / Accepted: 13 December 2021 / Published: 17 December 2021
(This article belongs to the Special Issue Challenges of Health Data Analytics)

Abstract

:
Recent digitization technologies empower mHealth users to conveniently record their Ecological Momentary Assessments (EMA) through web applications, smartphones, and wearable devices. These recordings can help clinicians understand how the users’ condition changes, but appropriate learning and visualization mechanisms are required for this purpose. We propose a web-based visual analytics tool, which processes clinical data as well as EMAs that were recorded through a mHealth application. The goals we pursue are (1) to predict the condition of the user in the near and the far future, while also identifying the clinical data that mostly contribute to EMA predictions, (2) to identify users with outlier EMA, and (3) to show to what extent the EMAs of a user are in line with or diverge from those users similar to him/her. We report our findings based on a pilot study on patient empowerment, involving tinnitus patients who recorded EMAs with the mHealth app TinnitusTips. To validate our method, we also derived synthetic data from the same pilot study. Based on this setting, results for different use cases are reported.

1. Introduction

Ecological Momentary Assessments (EMAs) are utilized to capture the immediate behavioral experience for a medical phenomenon. The behavioral experience can range from the present moment to a few minutes earlier, or can be a recollection of events that occurred at earlier, longer time frames. Presently, EMAs are mainly recorded with help of mobile technology, namely, by digital devices that notify the users multiple times for a period of days or weeks, so that they record current or recent medical states, behaviors, or environmental conditions [1]. By using this approach, the period of recall can be reduced to hours or minutes. Hence, EMAs allow observing the natural set of behaviors and moods of the participants [2]. EMAs are used for behavioral monitoring and tracking the progression of bipolar disorders [3,4], for studying tinnitus distress [5], and for monitoring mood in major depressive disorders [6]. Importantly, EMAs reduce the recall bias [7,8,9], and can thus be exploited by physicians for reliable patient monitoring and decision support. Comprehensible Artificial Intelligence (cAI) [10] is a transition framework that encompasses multiple disciplines such as AI, Human-Computer Interaction (HCI), and End User explanations, along with the combinations of techniques and approaches such as visual analytics, interactive ML, and dialog systems. Visual analytics act as an intersection between AI and end-user explanations to provide rich visualizations, which are helpful for humans to understand, interpret and further improve the trust of a developed system. A Clinical Decision Support System (CDSS) involves the provision of effective assistance to clinicians during the process of patient treatment and diagnosis [11]. It was utilized for effective and interactive communication between a physician and patients through alerts that are provided during self-monitoring [12]. The ability for the patients to make clinical decisions are also utilized in CDSS [13]. Recently, Intellicare [14] platforms attempt to provide remote therapies for anyone at any point in their mental health journey through smartphone-based applications to reduce stress, depression, and anxiety. By including the parts of cAI transition framework, the existing systems can be enhanced.
Note that considerations of mHealth data in the context of well-being are not fundamentally new. Studies and concepts can be found that deal with the combined perspective of well-being and mHealth. In [15], for example, a technical concept is discussed, in which mHealth monitoring is related to well-being monitoring in general, concretely by a comparison with an online social network scenario. The results shown by [16], in turn, reveal that emotional bonding with mHealth apps can be related to the well-being of the users. Many other recent and related works can be found in this context [17,18,19]. However, the decision to investigate the similarity among the users to assess their well-being is not considered so far, to the best of the authors’ knowledge. In addition, to accomplish the similarity inspection based on advanced visualization techniques is not being pursued in the same way by other works.
In this work, we contribute to CDSS through a medical analytics interactive tool that facilitates the inspection of a user’s EMA and clinical data. It juxtaposes timestamped EMA recordings of users with related clinical data, and predicts future user recordings with respect to symptoms of interest. In particular, we investigate the following research questions:
RQ1- 
How to predict a user’s EMA in the near and far future, on the basis of similarities to other users?
RQ2- 
How to identify and show outlierness in users’ EMAs?
RQ3- 
How to show similarities and differences in the EMAs of users who are similar in their clinical data?
The core idea of our approach is the exploitation of similarity among users in their static clinical data as well as in their dynamic, timestamped EMAs. We build upon similarity for prediction in time series (see RQ1), upon identification of users with outlier behavior (see RQ2), and upon visual and quantitative juxtaposition among users (see RQ3). For the validation of our work, we use the data of a pilot study on the role of mHealth tools for patient empowerment. The study involved 72 tinnitus patients, who recorded EMAs with the mHealth tool TinnitusTipps over an 8-weeks period between 2018 and 2019.
Our contributions can be summarized as follows:
  • We demonstrate the users neighborhood comparisons over data (i.e., both static, dynamic, and timestamped EMA) and utilize them for predicting user’s EMA recordings and show that users neighborhood are indeed useful in making the ahead predictions.
  • We introduce a voting-based outlier detection methodology to identify users who behave differently in their interaction with the app and also introduce tailored interactive visualizations that can be inspected.
  • We introduce a medical analytics tool with the introduction of tailored interactive visualizations to demonstrate the nearest neighboring user’s behaviors recorded through the app, and a visualization to show the ahead predictions for a study based on the identified nearest neighbors by constructing pathways.
The remainder of this paper is structured as follows. In Section 2 we provide the necessary information about the mobile health app and data of the pilot study on TinnitusTipps. In Section 3, Section 5 and Section 6 we detail the methods introduced for answering our RQ’s. In Section 7 we report on the results of our analysis. In Section 8 we elaborate on the obtained findings and discuss improvements and limitations.

2. Materials

The tinnitus study was approved by the Ethical Review Board of the University Clinic Regensburg. The ethical approval number is 17-544-101. All study participants provided informed consent.

2.1. The TinnitusTipps Mobile Health App

The TinnitusTipps app (https://tinnitustipps.lenoxug.de/, accessed on 15 December 2021) was developed by computer scientists, psychologists, and the Sivantos GmbH (which is a company specialized in hearing aids). TinnitusTipps is based on two considerations. Firstly, we wanted to utilize and exploit the experiences and benefits of the TrackYourTinnitus (TYT) mobile app. TYT combines Ecological Momentary Assessments and Mobile Crowdsensing (EMA-MCS) [20] to perform Digital Phenotyping [21] for users affected by tinnitus. Digital Phenotyping, in turn, is based on the idea to perform Ecological Momentary Assessments by the use of digital technology like smartphones. Secondly, we wanted to improve TYT and therefore added feedback features. In TYT, we only gathered EMA data from participating users, but no feedback was sent back to them. Compared to TYT, TinnitusTipps contains following new features: daily tinnitus tipps, data visualizations, and direct feedback from a healthcare provider.

2.2. Data of the CHRODIS+ Pilot Study on Tinnitus

From TinnitusTipps we collected following types of data:
  • Registration data with the Tinnitus Sample Case History Questionnaire (TSCHQ, 35 items) [22] and the Tinnitus Hearing Questionnaire (THQ, 6 questions) (cf. Appendix A Table A1) completed by all app users.
  • EMA data: with an EMA questionnaire (8 items) that captures within-day fluctuations of tinnitus loudness, distress, of hearing ability, stress and further aspects of mood and health condition (cf. Appendix A Table A2). As it can be seen in Appendix A Table A2, this questionnaire consists of 6 short questions to be answered with a numerical value between 0 and 100%, where larger values indicate worse condition.
    Users are notified by the smartphone multiple times within a day, and randomly, to record the current EMA.
For our analyses, we considered only users who have recorded both types of data, i.e., we excluded 19 users for which we did not have the registration data or who did not record any EMA. The total number of EMA recordings of the 53 users was 7255.

3. Methods

Our workflow consists of methods for data conversion (1), neighborhood creation (2), and generation of question combinations as features (3), which are used for all three RQs. The backend and frontend components of our workflow are shown in Figure 1 and explained hereafter.
On Table 1 we summarize the notation we used throughout this paper.

3.1. Data Conversion

The collected user data were stored in the JavaScript Object Notation (JSON) format, from which they were converted into Comma Separated Values (CSV) format. The processed data encompass two types of user answers: the static data from the two registration questionnaires and the time series data from the EMA questionnaire. The numeric values of the EMA variables were rescaled to be between 0 and 1. We used the time series on loudness (s02) and tinnitus distress (s03) for prediction, and the time series of the items s02, s03, …, s08 for visualizations.

3.2. Time Series Alignment

The users start recording their observations over different time periods and after a few days, there are possibilities to stop. To deal with such gaps and variations in the user’s time series recordings, a time series alignment approach is taken into account. The approach involves the aggregating the time series at a daily-level using either mean, maximum or minimum. The first observation within a given month is initialized to 0, and is assigned to the p r e v _ e n c o d e d _ v a l variable. For the next observation within the given month, the p r e v i o u s _ e n c o d e d _ v a l is incremented and the updated value is reassigned to p r e v i o u s _ e n c o d e d _ v a l . For the first observations for the next month, the last day index value of p r e v i o u s _ e n c o d e d _ v a l is obtained and incremented accordingly until the last recording for a user. For the ahead predictions, minimum, maximum and mean aggregations are investigated.

3.3. Creation of Neighborhood

We created two types of neighborhoods—using static data and using time series data, as described hereafter.

3.3.1. Neighborhoods on Static Data

Let U = u 1 , u 2 , , u n represent a set of users, with n being the number of users and features of a user u n is denoted by s t a t i c F = s 1 , s 2 , , s r , with | s t a t i c F | is the total number of features in the static data. The response attributes in s t a t i c F contain the user’s static information (age, background, and comorbidities) and their past medical conditions (presence/absence of neck pain, etc.).
For each user, we compute the set of k-nearest neighbors subject to a threshold, namely, the average of the distances over all pairs of users, so that only closer neighbors contribute to the similarity.
The Heterogeneous Euclidean Overlap Metric (HEOM) is used to measure the similarity. The definitions hereafter come from [23,24], but we use a different notation.
Let sr(x), sr(z) be the feature values of s r s t a t i c F , then the overlap metric distance for the users x , z U is defined as:
d O M ( x , z ) = i = 1 n f O M s r ( x ) , s r ( z )
with, f O M s r ( x ) , s r ( z ) = 0 if s r ( x ) = s r ( z ) 1 otherwise
When the attributes value of sr(x), sr(z) are continuous, then we define the range difference between the two attribute values as follows:
r n _ d i f f s r ( x , z ) = | s r ( x ) s r ( z ) | r a n g e s r
where, r a n g e s r = m a x ( s r ) m i n ( s r ) is used to normalize the attributes.
From the overlap metric and range difference, we derive HEOM as:
d H E O M ( x , z ) = i = 1 n f H E O M s r ( x ) , s r ( z ) 2
with, f H E O M s r ( x ) , s r ( z ) = f O M s r ( x ) , s r ( z ) if s r is nominal s r ( x ) s r ( z ) if s r is continuous 1 otherwise

3.3.2. Neighborhoods on the Dynamic Data of the EMA Time Series—One per EMA Item

Next to the static information, the users also record an ordered sequence of observations that constitute the EMA recordings. Let U = u 1 , u 2 , , u n represent the set of users with the EMA recordings. We denote the sequence of observations of an u n U as u n = o 1 , o 2 , , o i , , o m , where each observation contains a set of features along with a timestamp o i = { t i , F i }, where t i represents the timestamp of the user u n , F i represent the feature space of EMA recordings which is denoted as F i = f 1 , f 2 , , f j ,for j = 1 , 2 , , | F i | ; | F i | = 8, and m represents the total number of EMA recordings recorded by u n over the days denoting the length of the time series. The multiple observed values recorded by the user within a given day for the considered EMA variable are averaged for the neighborhood computations.
We compute the similarity between a user x and user z for each EMA item separately, after aligning their time series at the same (nominal) day 0:
  • ‘Day matching’ between x and z: number of days from day 0 onward and until the last day of the shortest between the two time series. For example, if x has EMA observations for 30 days and z for 60 days, then the matching is on the first 30 days, with day k = 0 , , 29 of x matched to the corresponding day of z.
  • For each of the matched days for x and z, a euclidean distance is computed. A counter is maintained to capture the number of days when both the users have reported their observations, denoted by, c o u n t s and the distance obtained for each day is summed up and is denoted as total _ dist .
  • Finally, a fraction of the t o t a l _ d i s t c o u n t s is returned as the similarity between the user x and z.

3.4. Creation of Feature Combinations

The similarity between users and, accordingly, their neighborhoods over static data, can be computed using the whole set of registration items or a subset of them. To assess the effect of different aspects of a user’s recordings, we construct the following overlapping ‘subspaces’, i.e., combinations of features from the set of registration items:
C1:
user background and tinnitus complaints information (items tschq02-04); tinnitus historical information (items tschq05-18), including the initial onset of tinnitus, loudness, tinnitus awareness, and different past treatments
C2:
experienced effects of tinnitus (tschq19-25), and questions on hearing quality/loss (hq01, hq03)
C3:
further conditions, such as neck pain, dizziness etc. (tschq28-35), as well as the items on hearing quality/hearing loss (hq02, hq04)
C4:
all the TSCHQ and HQ items
Along with these combinations, the numerical value loudness (tschq12) is also included, so that there are no ties in the computation of user similarity.
By using these combinations, a similarity is built for the specified neighborhood size (k); i.e., for each of the user and the obtained similar users are utilized for the prediction of tinnitus distress by computing the k nearest neighbors over the registration data.

4. RQ1: How to Predict a User’s EMA in the Near and Far Future, on the Basis of Similarity to Other Users?

A prediction of a class label for a test user u, for an EMA recording f j F at a timepoint t i is performed through obtaining the nearest neighbor users of u and by utilizing the nearest neighbors EMA recordings until t i . The work focuses on the prediction of a numeric class label only.

4.1. Ahead Prediction of the Target Variable

Let T be the set of test users, let u be a test user, and let t i be the timepoint at which we want to predict the values of l 1 future EMA recordings for u. We first compute the set of nearest neighbors of u for a given k. We denote it as N N k ( u ) . Then, we consider the following options:
  • Weighted average in a user’s neighborhood: At timepoint t i , the predicted value of an EMA recording is then the average of the values of this EMA recording over the users in N N k ( u ) . This computation can be done for l 1 timepoints at once.
  • Linear regression in a user’s neighborhood: for a test user u T , a linear regressor is built on the time series observations of each of the nearest neighbor users u j N N k ( u ) until the timepoint and averaging the slope ( m s l o p e , u j ) and the intercept: ( m i n t e r c e p t , u j ) parameters of each of the nearest neighbor user to create a combined model as introduced in [25].

4.2. Evaluation for RQ1

We investigate how the value of k, the choice of subspace (from the registration data), the kind of neighborhood (static vs dynamic) and the aggregation function over the EMA within a day (min, mean, max) affect the quality of the prediction. For the evaluation, we consider Root Mean Squared Error (RMSE) over a set T test users. In particular:
  • For a test user u i T , the nearest neighbors N N k ( u i ) are obtained (either for static or dynamic). Once after obtaining the nearest neighbors of u i , their respective time series are aligned, and wherever multiple observations are present for a given day, they are aggregated using either the minimum(), mean(), or maximum() functions.
  • Next, based upon ahead prediction timepoint (in a day), future EMA recordings (l), and the nearest neighbors time series information; ahead predictions are obtained based on the proposed methodologies as per Section 4.1.
  • For a test user, the error for a feature f i F is measured using Root Mean Squared Error (RMSE) against the true value of that feature as follows:
    R M S E ( f i , f i ^ ) = 1 l i = 1 n f i f i ^ 2
    where f i ^ denotes the predicted values obtained for a feature f i F , and l denotes the number of ahead days to predict.
Let the above three steps of evaluation be denoted as e v a l . The same process is carried out for all the test users and compute the final quality value for a given neighborhood size k as
a v g R M S E ( T , k ) = 1 | T | i = 1 | T | e v a l ( u i )
and the achieved RMSE for all the test users is averaged for the provided neighborhood size (k).
To assess the changes in ahead-prediction from earlier days to days in the distant future, a comparison between near predictions (earlier days) and far predictions (far-future days) are also done based on the concepts of [25]. The neighborhood size from 1 to 30 is chosen in this work.

5. RQ2: How to Identify and Show Outlierness in User’s EMA

An outlier instance is one whose properties deviate from most of the other instances. We highlight our workflow next, and detail the concept of visualizing the tendency of outlierness in Section 6.2.

Methodology

Our approach improves the rank-based outlier detection approach [26] using the voting-based strategy proposed in [27]. Instead of working over a single similarity measure and manually defining the neighborhood, a voting-based strategy with the capabilities to work over the tailored similarity measure is employed. The algorithm runs the detection over a set of k neighborhoods, over each of the defined subspace combinations and votes the outlying data object at each of the neighborhood iterations based upon the outlier score as per Appendix C Equation (A1) to pick out the top-n outliers, which are highest voted over the entire set of the defined neighborhoods. According to our proposed methodology, if the z-score value is greater than or equal to 2.0 , then it is treated as an outlier since the algorithm runs over multiple values of k and the scores are combined as overall average scores across the set of the neighborhoods. Additionally, the proposed voting strategy also utilizes the Local Outlier Factor (LOF) algorithm in the workflow to indicate a local outlier. The descriptions of both the algorithms are provided in the Appendix B and Appendix C. In the comprehensive survey of outlier detection techniques for data streams [28], the authors explain the concept of incremental LOF [29], in which the past data objects together with the local outlier factors are maintained at k. Close to their approach, by utilizing both the variants of the outlier detection techniques, different types of outliers can be obtained. The outlier scores are separately obtained for both the detection algorithms, and there may be no common outlier data objects amongst the two algorithms.
The detection and the voting mechanism involves to detect outliers at each k and subspace by both LOF and RBDA by maintaining a global outlier data store that keeps the count of number of times the data object have the tendency of outlierness, and the outlier scores obtained at each k and subspace. By the intersection of both the data stores, finally the common outliers are produced.

6. RQ3: How to Show Similarities and Differences in the EMA of Users Who Are Similar in Their Clinical Data?

For the visualization of similarities and differences, we propose a system as shown in the activity diagram Figure 2, which combines the mechanisms for predictions and for outlierness computation. The figure shows the interactions between the user and the domain expert (physician or clinician) with the system.

6.1. Similarity Overview User Panel

An overview of the proposed visual interface is shown in the Figure 3. The design of the visual interface is based upon the visualization guidelines introduced in [30], which is to firstly provide an overview, allowing to zoom to required regions in the visualization, filter, and allowing to provide details on demand. The panel is designed based on iterative discussions and some inspirations from the existing visualization literature. The designed interface is divided into six parts: (1) denotes the selection panel of the selected study user with 1(a) representing the user-identifier; 1(b) denotes the selection of the criteria of similarity; 1(c) is a selection drop-down to choose the question combinations; (2) represents the study user time-series observations displayed over a line chart with capabilities of interactions; (3) shows the similar users to the study user highlighted in orange; (4) represents the comparison of the time series recordings grouped by sessions of the day for the selected user and the nearest neighbor; 4(a) denotes the selection of various attributes to dynamically render the comparison plot; (5) represents the comparisons of the static questions between the study user and the selected nearest neighbor; and (6) represents a slider for allowing to interactively change the number of the nearest neighboring users and updating the visualizations accordingly. Furthermore, we also depict the user panel for EMA (Refer Appendix D Figure A1).

6.2. Nearest Neighbor Visualization

The visualization consists of two steps. At first, through the analysis steps, the required similar user data and information are obtained. Secondly, by using the obtained similar users and the information, color encoding rules are created as represented in Figure 4b and visualized as a graph.
Construction of the visualization: To visualize the nearest neighbors to the study user and to represent the user characteristics, a node-link diagram is utilized as per [31], and the users are placed in the best fit by using the force-directed scheme. However, we add another aspect, namely the representation of each user node with multiple visualization elements.
The Figure 4 shows the nearest neighbor visualization: (1) the study user (orange) and nearest neighbors (blue), (2) the closest neighbor with the edge in bold. The tailored visualization of a single user node is depicted in Figure 5, with the following interactive visualization elements:
Outer Circle: when clicked, the study user is highlighted in orange and circles depicted in blue show the nearest neighbors to the study user.
Inner Circle: the color of the area inside the inner circle is determined by the rules derived from the answers of the user, scaled from 0 to 100 percent. Based on the values, blue color bands ranging from light to dark are used to indicate high to low scale values, represented as shown in Figure 4b. The circle radius is drawn dynamically based on the numerical valued answer of the users. Firstly, a default value for the radius is defined as r and based on the answer value, the radius is dynamically updated by adding the answer values to the radius denoted by r : r + a n s w e r _ v a l .
Outlier Score Indicator: As outlier scores, we use the z-scores provided by the Rank Based Detection Algorithm (RBDA). To visualize the scores within the nearest neighbor visualization, we introduce a rule-based color coding with values ranging from blue to red, with the shades of red indicating the suspected outliers. The color encoding legend is shown in Figure 4c, representing the color-coding methodology.
Circle Border: the boundary of the circle is colored with purple when the users answered YES to the question “TSCHQ26-Do you think you have a hearing problem?”.

Interactions

  • Zoom into each of the user’s nodes to understand the level of the loudness (TSCHQ12) shown in the inner circle and the presence or absence of the hearing problem (TSCHQ26) answer shown in the circle boundary.
  • Hover in the center of user node to obtain the user clinical data. (Refer: Figure A2 in Appendix D for more information). By such information easily, a domain expert can identify the differences amongst the obtained similar users.
  • Click on one of the nearest neighbors to obtain comparison plots to compare the time series observations as sessions, and the registration properties with the study user.
  • Additionally, the distance information is provided within the nearest neighbor visualization on click, which generates a table (Refer Use Case Example: Figure A4) calculated as the difference in distance of each of the nearest neighbors from the first nearest neighbor. A cell within the generated table is colored pink when the obtained difference in distance value of a nearest neighbor is greater than a threshold ( τ ) defined as the mean ( μ ) + 1 * standard_deviation ( σ ) of the distance differences amongst all the nearest neighbors. This provides the information of the dissimilar users in the computed neighborhood.
  • A slider is provided to the domain expert to dynamically increase or decrease the number of nearest neighbors for a study user.

6.3. Comparison Visualization

When a click interaction is made from the nearest neighbor visualization, two plots are obtained. The Figure 3(4) shows the comparison between the EMA recordings of the study user and the interacted nearest neighbor user. The plot of study user is colored orange, and the nearest neighbor in blue follows the same coloring standard. The users record their observations multiple times in a day, containing the date and timestamp. From the timestamp, the hour recording value can be utilized to group the sequence of observations into Early Morning (4:00–8:00), Morning (08:00–12:00), Afternoon (12:00–16:00), Evening (16:00–20:00), Night (20:00–00:00) and Late Night (00:00–4:00). By utilizing this information, a box plot is created for each of the user recorded observations and interactions through the tool is provided. These plots are useful in comparing two users who are similar in their static properties may also be similar in their EMA recordings at parts of the day. The Figure 3(5) shows an important plot for the comparison of the static properties of the study user and the interacted nearest neighbor. For representation, the user’s static data are ordinal encoded and presented in the form of a heatmap.

6.4. Predictive Visualization over the Line Plot of the Study User

Let the study user trajectory for a variable under consideration be represented by the set of observations o 1 , o 2 , , o i , , o n , where n is the number of observations and each of the observation o i contain { t i , v i } , where t i is the aligned time index called as d a y _ s e s s i o n _ i n d e x and v i is the value of the EMA variable (s01–s08) under consideration. This trajectory is then visualized over the line plot, with the x-axis denoting the time aligned index and the y-axis showing the variable under consideration. This work concentrates on visualizing one of the attributes of the Ecological Momentary Assessment (EMA), which is tinnitus distress (s03); recordings are shown in the Figure 6. To visualize the predictions for the ahead timepoints, a pathway prediction visualization is introduced.
Construction of the pathways: Firstly, the ahead predictions for the reference timepoint are obtained, and then they are dynamically redrawn over the line plot to show the pathways.

Interactions

  • Click over each of the data points on the plot and a prediction is made for the next timepoints over a selected data point, and the pathway is visualized in Figure 6(2).
  • Hover functionality; compare the predicted with the actual values.
  • Brush (Figure 6(1)) a region of timepoint; perform a prediction in that region.

7. Results

7.1. RQ1: Experimental Setup

The experimental setup involves creating a split of the users into 80%(train) and 20% (test) to investigate the effectiveness of the proposed method for the prediction of the target variable tinnitus distress-(s03) from the Ecological Momentary Assessment (EMA) data. A neighborhood is created first by using the static registration data over the defined subspaces, and secondly, by using the Ecological Momentary Assessment (EMA) variable Loudness-(s02) for the neighborhood creation amongst users as per the proposed methodology in Section 3.3.2. After the neighborhood creation, a train and test split through randomization is created. For static registration similarity, the registration data are used for the neighborhood creation, and the obtained nearest neighbor time series from the EMA data is used for the ahead predictions of the test users. In the case of loudness-based similarity, both the neighborhood creation and the prediction are over EMA data.
The user’s multiple recordings within the day are aggregated as mean (mean_day), maximum (max_day), and minimum (min_day) observations for the day. This work addresses such a scenario to investigate which of the aggregation function leads to the best prediction of tinnitus distress-(s03).

7.2. RQ1: Comparison across Combinations

At first, all the questions were considered for the evaluation. The distance weighting was not useful in the ahead prediction, and hence, was not considered further in our analysis. By utilizing the combination C1–C3, and constructing a neighborhood, it was observed that the combination C3, mean_day aggregation is best in the ahead prediction of tinnitus distress-s03 as shown in the Figure 7. Additionally, for larger neighborhoods, there is not much difference in the error rate of the subspace combinations.
Table 2 provides the average for k ranging from 3 < = k < = 30 for each mean_day, min_day, and max_day observations of the neighbors assisting in the predictions. The lowest average RMSE value obtained across the neighborhood range is highlighted in bold, along with the standard deviation. Combination C3 has the lower RMSE with C3 performing better for both mean_day and min_day with mean_day having low RMSE values.

7.3. RQ1: AvgRMSE for Near and Far Predictions

The previous analysis showed that mean_day and combination C3 performs best with respect to the ahead prediction of tinnitus distress. To assess if the ahead predictions remain the same even for the further days for the test users, execution is performed following the same evaluation process by choosing the starting day at far (80%) of the test user’s time series recordings, over the k values, and the subspace combination as C3. A comparison is made over the early and far predictions. The Figure 8 denotes the near and far comparisons of both mean_day aggregations of nearest neighbor time series recordings. It is observed that predictions at earlier days are better than at far days. A possible explanation is that for the far off days there are data from fewer users, thus more noise.

7.4. RQ1: Comparison over the Similarity Constructed through Loudness-s02

To compare if the loudness values recorded by the users can achieve better similarity, and in-turn improves the predictions of tinnitus distress-(s03), a neighborhood is constructed, and ahead predictions are evaluated for the earlier days as seen from the Section 7.3. The mean_day aggregations of the nearest neighbors are better (Refer: Figure S1 of the Supplementary Material). K values between 5 and 11 have lower RMSE and the error rate increases gradually at the higher neighborhood sizes.

7.5. Important Findings of RQ1

As seen from previous results that the RMSE decreases as k increases but at the higher k values, there is not much change to the RMSE. So, the best value for the k chosen is at the point where the error rate stops changing or increases (elbow method). Based upon this, it can be seen that k = 9 is suitable by considering static registration-based similarity over the combination C3, and mean_day. By considering the loudness-(s02) based similarity, k = 11 is suitable with mean_day aggregations of nearest neighbors. Next, two results, one considering the neighborhood over static similarity and the other considering loudness-based similarity, are shown with the obtained values and following the evaluation procedure as per Section 4.2.
By fixing the value of k to 9, for the combination-C3 over the static registration similarity, for the loudness-based similarity k to 11, and choosing the mean_day observations of the nearest neighbors, 3 ahead predictions are reported by selecting multiple bounds (20%, 30%, 50%) as the start day of the test user’s time series recordings. The RMSE obtained for the ahead predictions at each of the bounds is averaged and results as an average RMSE over multiple bounds. A discretization process is applied for the target variable tinnitus distress-(s03) to assess if the test users in the neighbor ranking positions move similarly to know if the test users have low tinnitus distress values, or high tinnitus distress. The discretization process involves converting the numerical tinnitus distress values into categories. Here, the categories created are [0–0.4] as Low, [0.4–0.7] as Moderate, and [0.7–1.0] as High as per [25]. Some important findings of this analysis are that the test users with the lower RMSE values tend to have most of the observations recorded in the Low category, and the RMSE is higher when the observations are in the High categories. This provides an understanding that the users who have their tinnitus distress in the range of Low to Moderate categories are similar together as others. However, it was also observed that there were few exceptions for certain users whose most observations are recorded in Moderate and High categories. This was further verified by combining (train+test) recordings and was seen that the probability of the observations belonging to the Low category is highest with a value of 0.497 providing an insight that better predictions are obtained when observations fall into the Low category. We tabulate the frequencies and prior probabilities of tinnitus distress against each of the categories in the Table 3.

7.6. RQ2: Results of Outlier Detection Mechanisms

7.6.1. Results over the Created Synthetic Data

In this work, the synthesized data are created using PyOD [32] library, where the outlier instances are created by the uniform distribution and normal instances from the Gaussian distribution. The created data is evaluated with Area under the ROC-curve (AUC) [33]. Two scenarios are highlighted over the synthesized data to assess the ability of the Ranked Based Detection Algorithm (RBDA) in comparison with LOF (Local Outlier Factor) to understand the differences in the detection rate as the neighborhood (k) increases.
Scenario-1: In this scenario, 100 instances are chosen as a global detection with number of features as 35, k is varied from 10 to 60, and 10 outliers are ingested. The chosen distance measure is Euclidean for both RBDA and LOF. From the Figure 9(Left), both algorithms are able to detect outliers, while LOF detects them at earlier neighborhood size and do not change with increasing the neighborhood size. The RBDA detects outliers at full extent at the higher neighborhood values.
Scenario-2: In this scenario, the number of instances are chosen as 100, containing 4 clusters of data points containing 15 features, the cluster have different size and density. The k value is varied from 10 to 60, and 10 outliers are ingested. The chosen distance measure is Euclidean for both RBDA and LOF. From the Figure 9(right), at initial neighborhood sizes, LOF is able to detect the outliers well. However, for neighborhood with larger size k, the detection rate decreases and ranking based detection performs better. This is because of the clusters are of varying densities and the LOF is not be able to identify them for large neighborhood sizes.
From this analysis, we make sure that both of these detection mechanisms can be included in our proposed workflow.

7.6.2. By Using the Tinnitus Data

From the proposed voting-based methodology as per Section 5, differently behaving users are identified from both the introduced similarity methods. The outlier scores from the Ranked Based Outlier Detection obtained from the workflow are used to indicate differently behaving users in the visualization tool (Refer: Section 7.7). By using the same discretization process as per Section 7.5, most of the user’s belonged to either high loudness+high tinnitus distress or low loudness+low tinnitus distress groups. However, there are few users who had high loudness+low tinnitus distress. A plot showing the total number of observations of loudness and tinnitus distress of the users are shown in Supplementary Material (Refer: Figures S2 and S3).

7.7. RQ3: Validation of the Visual Interface through Usage Scenarios

To validate the effectiveness of the introduced visualizations, various usage scenarios are introduced. For some use cases, simulated data are utilized.

7.7.1. Simulated Users Data Creation for Validation

The simulated users are created by utilizing the same data from the pilot study. We basically create three different type of users: (1) Normal users: these users are created by obtaining a random value for a specified feature from the existing static data and for a given set of a range (in days) generating a random set of timestamped values for creating the EMA recordings; (2) outlying user: the top-n outlying users detected as differently behaving in the interaction of the app from the introduced outlier detection methodology for both static and EMA recordings are used for the study; (3) twin user: the twin user is just the exact copy of one of the selected user from the existing static data and having a randomly generated EMA recordings for a given range.

7.7.2. Outlier User Scenario

Study of a user whose static properties deviate from the rest of the users, also termed as an outlying user. For this scenario, user-20 is studied from the visualization tool. From the tool, combination C4 is chosen for the analysis so that all features can be observed. The same procedure is followed for other subspace combinations as well. The Figure 10 shows the 5 nearest neighbors through the visualization. The outlier score is indicated as a coloring encoding box. More details of the visual explanations are provided in Figure 4. Here, the user-20 is an outlier identified by the drawn dark red rectangle, indicating the higher z-score values obtained from the voting-based method (Refer Appendix D Figure A5 for the zoomed view).

7.7.3. Twin User Scenario

Study of a user whose static properties are exactly the same as the one of the other user in the underlying dataset. The main focus of this scenario is the understanding and assessment of the constructed similarity in identifying the closest neighbor as the user with exactly the same static properties. From the Figure 11, for the simulated user-444154, the closest neighbor is 154. The simulated user is the copy of user-154 shown by the last 3 digits of the user-id of the simulated user. This indicates that the proposed user similarity in this work is able to identify it correctly. The distance value will be 0 in this case and from the visualization perspective, an offset of distance is added between the study user and closest neighbor to avoid overlap as shown by the dotted marker line in Figure 11. The distance information is also provided (Refer Appendix D Figure A4).

8. Discussion

The previous Section 7 highlighted some important results of the research questions concerning the ability to predict user’s EMA based on similarity, identification of individuals whose properties are different, and the overview of key usage scenarios addressed in the visualization tool.
Regarding the RQ1, from the Figure 7 it confirms that the subspaces do help in the ahead predictions of the user’s EMA recordings as opposed to considering all the questions. Similar behavior was also observed in [34] where the predictive power of the model was assessed through clustering. However, considering only the user’s EMA recordings to construct similarity has a lower error rate as compared to considering the registration data. Therefore, it can be inferred that a similarity measure considering both forms of data together must be explored. As seen from Figure 8, the prediction of user’s EMA for the far future is indeed difficult, which is in agreement with the findings of [25].
Regarding the RQ2, we showed that there exist users whose EMA observations deviate from most of the other users. For example, it is already known that there are users with high tinnitus loudness and low distress [35]. These users behave differently in their interaction with the app, and closer inspection are of interest to the domain expert through interacting with the tool. The predictive mechanism introduced in [34] can be enhanced by including our outlier detection methodology within the workflow.
The proposed medical analytics interactive tool introduces two innovative visualizations regarding RQ3. As seen from the user-scenarios in the Section 7.7, we were able to correctly detect the nearest neighboring user’s based on their daily life behaviors recorded through the app. As seen from the outlier scenario, the visual encoding of the color code provides useful information to the physician about a user’s outlierness. The proposed ahead prediction visualization can help in reconsidering the patient treatment or monitoring policy.
Finally, the results are briefly reflected in a broader scope. In general, we have learned in recent years that mHealth can deliver interesting results for medical purposes [36]. However, the use of mHealth has revealed too many flavors without comparable construction principles like presented by [37]. Consequently, we have to investigate comparable and powerful tools to approximate ourselves to tangible medical results based on the use of mHealth. The idea to compare users might open one more perspective to learn more about mHealth behaviour and users. In this light, with this work, we may enable both computer scientists to create new algorithms better and empower medical experts to more quickly reveal user differences. As tinnitus is a very heterogeneous phenomenon, works like this may help to better demystify this heterogeneity.
From a practitioner point of view, the proposed method could be useful to identify those patients that need special care or will need additional special care any time soon. Especially in the case of sparse resources, it would allow focusing on those patients that need more or specific care to prevent more severe clinical conditions (and associated costs). By preventing higher costs and severe conditions, a commercial solution could gain market share over others that were not able to predict future data based on currently available patient information.
Threats to validity: The experimental design for RQ2 did not contain class labels to verify whether a user is behaving differently or not. To truly assess RQ2, an experiment setup must be made over a labeled dataset and verified. The underlying data did not have user answers for Mini-TQ questionnaire, and we firmly believe that making this mandatory in the app can help to better assess the user’s EMA recordings [34]. In the proposed interactive tool, our ahead prediction visualization is concerned with only one of the EMA variables. Including predictions for the other EMA variables will be of greater benefit to the physician. The use cases explored in our interactive tool were developed with limited input from very few physicians, and a wider user study with more experts would give better insights into whether the use cases are exhaustive and sufficient from the point of view of a physician.
Future actions: As part of the future work, we are looking forward to exploring and exploiting similarity methods that can work on both static and dynamic data, and perform a comparison experiment on our introduced and existing approaches. In terms of the visualization, we look to explore new approaches of visualization to compare EMA recordings as event sequences.

9. Conclusions

To conclude, we showed from the tool that it is possible to predict a user’s future EMA based on the data of similar users, with the ability to identify outlying individuals. Our findings on outlierness indicate how important it is to closely monitor the users’ EMAs as a medical practitioner, who can assess whether action is needed as EMA values change. Our findings on the evolution of users with similar registration data over time show that registration data alone are not adequate to assess similarity during EMA recordings. Apart from the quality of predictions, we also see that the results of such an interactive system to explore the usefulness of neighborhoods and the outlierness of an individual can help bring scientists in two fields come closer together. The medical practitioner can use such a system to better intuit factors that make users different, which researchers in computer science can use to develop better algorithms to discover these differences and visualize them in a way that non-experts can intuitively understand.
On a technical front, we discuss some opportunities for extensions on both the prediction of EMAs in the near and far future, and also the possibility to simultaneously exploit static and dynamic data during the neighborhood discovery process. The results of these methods can be compared against the current proposed method, which can serve as a baseline. The design of the visualization can be adapted to this special case with minimal effort. We also make an initial attempt with the facilitation for collecting tinnitus static data (Refer: Figure A3). We further look to improvise regarding automatic collection, information processing and visualize within the tool itself.
More work is needed to understand whether there are subpopulations of users who have similar registration data and evolve similarly. Our findings on near and far future prediction of EMA indicate that after the first days of interaction with the mHealth app, the involvement of the users may decrease. Hence, mHealth app developers and medical practitioners should think of incentives to stimulate regular interaction (for example, gamification).

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/e23121695/s1. Please refer to the supplementary material Section S1 for the results on constructing similarity using Loudness-s02, and Section S2 for comparison between reference users with differences in their tinnitus loudness and distress.

Author Contributions

S.P. and M.S. were involved in the conceptualization, methodology, and writing—original draft preparation and review, with the support of B.L. and W.S. on medical questions, R.P., R.K. and J.S. for mHealth technology questions, R.P. and R.H. for mHealth app design questions, and W.S. for the mHealth app experiment questions; V.U. was involved in methodology and provided support on prediction questions; R.P., J.S. and R.K. provided the necessary software. M.S. provided supervision and project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially funded by the the CHRODIS PLUS Joint Action, which has received funding from the European Union, in the framework of the Health Programme (2014–2020), Grant Agreement 761307 “Implementing good practices for chronic diseases”. The development of the TinnitusTipps mHealth app was partially financed by Sivantos GmbH.

Institutional Review Board Statement

Ethical approval was granted by the Ethical review board of the University of Regensburg. The approval number is 17-544-101.

Informed Consent Statement

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

Data Availability Statement

The data used in this study are not available publicly. The authors may be contacted in case of data requests.

Acknowledgments

This work is partially inspired by the European Union’s Horizon 2020 Research and Innovation Programme, Grant Agreement 848261 “Unification of treatments and Interventions for Tinnitus patients” (UNITI).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
CDSSClinical Decision Support System
CSVComma Separated Values
EHRElectronic Health Records
EMRElectronic Medical Records
PSMPatient Similarity Metric
EMAEcological Momentary Assessment
HEOMHeterogeneous Euclidean Overlap Metric
ICUIntensive Care Unit
JSONJavascript Object Notation
kNNk Nearest Neighbors
LOFLocal Outlier Factor
MIMICMultiparameter Intelligent Monitoring In Intensive Care
OMOverlap Metric
RBDARank Based Detection Algorithm
RMSERoot Mean Squared Error
TYTTrackYourTinnitus
TSCHQTinnitus Sample Case History Questionnaire

Appendix A. Data Description

The data description of the static registration data are based on the Tinnitus Sample Case History Questionnaire (TSCHQ) (https://www.tinnitusresearch.net/images/files/migrated/consensusdocuments/en/TINNITUS_SAMPLE_CASE_HISTORY_QUESTIONNAIRE.pdf, Accessed on 15 December 2021), the data description of hearing questionnaire (HQ) is highlighted in Appendix A Table A1 and the data description of the EMA data (tinnitus diary) containing the daily recordings of the users is depicted in Appendix A Table A2.
Hearing Questionnaire (HQ)
Table A1. Data Description of HQ.
Table A1. Data Description of HQ.
Variable NameQuestionData Type
hq01How often do you have hearing problems?Single Choice
hq02How much do your hearing problems affect you?Single Choice
hq03Did an ENT physician or acoustician diagnose you with a hearing loss?Boolean-(YES/NO)
hq04Do you wear hearing aids?Single Choice
hq05If yes, which hearing aids do you wear?Free Text
hq06If yes, how often do you wear hearing aids?Single Choice
Tinnitus Diary (EMA)
Each of the time the user answers these questions, a timestamp is also associated with it, making it a time series data. These questions are asked multiple times in a day at random intervals.
Table A2. Data Description of Tinnitus diary.
Table A2. Data Description of Tinnitus diary.
Variable NameQuestionData Type
s01Did you perceive the tinnitus right now?Boolean-(Yes/No)
s02How loud is the tinnitus right now?Numeric-[0–100]
s03How stressful is the tinnitus right now?Numeric-[0–100]
s04How well are you hearing right now?Numeric-[0–100]
s05How much are you limited by your hearing ability right now?Numeric-[0–100]
s06Do you feel stressed right now?Numeric-[0–100]
s07How exhausted do you feel right now?Numeric-[0–100]
s08Are you wearing a hearing aid right now?Boolean-(Yes/No)

Appendix B. Description of LOF

The provided descriptions are based upon the works of [33,38]. A local outlier are those data objects that are at a distance apart as compared to other data objects within a cluster. Below are the steps taken to perform the outlier detection to compute LOF score:
  • For each of the points, p D where D represents the underlying dataset. The k-nearest neighbors q N N k ( p ) are obtained and in case of a tie in the distance then more than the provided k neighbors are used.
  • By using the obtained k-nearest neighbors N N k , a local density estimation for the data object is calculating using Local reachability density (LRD) given as:
    L R D k ( p ) = 1 / q N N k ( p ) d k ( p , q ) N N k ( p )
    where d k ( p , q ) denotes the reachability distance. More details are provided in [38].
  • By comparing the LRD with the data object and the obtained k-nearest neighbors, the LOF score is obtained as:
    L O F ( p ) = q N N k ( p ) L R D k ( q ) L R D k ( p ) N N k ( p )

Appendix C. Description of RBDA

Based upon the shortcomings of LOF to detect outliers in the varying densities, [26] proposed a new approach for detecting outliers based upon the reverse neighborhood ranking. The idea of the outlier detection is over the mutual proximity between the data object and the obtained nearest neighbors: Let p D be a data object with D as the underlying dataset. Then, q N N k ( p ) are a set of nearest neighbors of p. For a q that is in the k-nearest neighborhood of p and close to each other, then a question to ask is about closeness of p to q. The data objects p and q are close with higher confidence when they are close to each other and are not anomalous data objects. According to [26], for a data object, p D where D represent the underlying dataset and q N N k ( p ) are the nearest neighbors of p. Then the degree of outlierness of p can be defined as:
O k ( p ) = q N N k ( p ) r q ( p ) N N k ( p )
where r q ( p ) indicate the rank of the data object p amongst its neighbors q for the provided k value denoting the reverse neighborhood ranking. The higher value of the O k ( p ) denote that p is an outlier.
The final score involves normalizing the O k ( p ) . A z-score normalization is applied and the value greater than or equal to 2.5 , then, p is an outlier.

Appendix D. Other Visualizations and Facilities

User panel view when EMA is selected as the similarity criteria
When the similarity criteria as EMA is chosen, the loudness-based similarity is executed, and as shown from the user panel the question combination is disabled. The rest of the components remain the same.
Figure A1. User Panel-EMA. (1) denotes the selection panel of the selected study user, with 1(a) representing the user-identifier; 1(b) denotes the selection of the criteria of similarity. Here, EMA; 1(c) is a selection drop-down to choose the question combinations, when EMA it is disabled for selection; (2) represents the study user time-series observations displayed over a line chart with capabilities of interactions; (3) shows the similar users where the study user is highlighted in orange and the nearest neighbors in blue; (4) represents the comparison of the time series recordings grouped by sessions of the day for the selected user and an interacted nearest neighbor; 4(a) denotes the selection of various attributes of EMA which generates rendering of the comparison plot.
Figure A1. User Panel-EMA. (1) denotes the selection panel of the selected study user, with 1(a) representing the user-identifier; 1(b) denotes the selection of the criteria of similarity. Here, EMA; 1(c) is a selection drop-down to choose the question combinations, when EMA it is disabled for selection; (2) represents the study user time-series observations displayed over a line chart with capabilities of interactions; (3) shows the similar users where the study user is highlighted in orange and the nearest neighbors in blue; (4) represents the comparison of the time series recordings grouped by sessions of the day for the selected user and an interacted nearest neighbor; 4(a) denotes the selection of various attributes of EMA which generates rendering of the comparison plot.
Entropy 23 01695 g0a1
Hover on the user node to obtain clinical information
Figure A2. User clinical information on hover.
Figure A2. User clinical information on hover.
Entropy 23 01695 g0a2
A facility to collect the tinnitus static information data (TSCHQ) is also provided in the proposed visualization tool
Note: This feature is provided to support future expansion, and was not the source of the data used in the current analysis. All data that was part of this study was downloaded from the database which collected users’ mHealth questionnaire responses.
Figure A3. Support for future expansion through a data collection capability.
Figure A3. Support for future expansion through a data collection capability.
Entropy 23 01695 g0a3
Figure A4. Distance information of the simulated user.
Figure A4. Distance information of the simulated user.
Entropy 23 01695 g0a4
Figure A5. Zoom into single node.
Figure A5. Zoom into single node.
Entropy 23 01695 g0a5

References

  1. Pryss, R. Mobile crowdsensing services for tinnitus assessment and patient feedback. In Proceedings of the 2017 IEEE International Conference on AI & Mobile Services (AIMS), Honolulu, HI, USA, 25–30 June 2017; pp. 22–29. [Google Scholar]
  2. Stone, A.A.; Shiffman, S. Ecological momentary assessment (EMA) in behavioral medicine. Ann. Behav. Med. 1994, 16, 199–202. [Google Scholar] [CrossRef]
  3. Ortiz, A.; Grof, P. Electronic monitoring of self-reported mood: The return of the subjective? Int. J. Bipolar Disord. 2016, 4, 28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Porras-Segovia, A.; Molina-Madueño, R.M.; Berrouiguet, S.; López-Castroman, J.; Barrigón, M.L.; Pérez-Rodríguez, M.S.; Marco, J.H.; Díaz-Oliván, I.; de León, S.; Courtet, P.; et al. Smartphone-based ecological momentary assessment (EMA) in psychiatric patients and student controls: A real-world feasibility study. J. Affect. Disord. 2020, 274, 733–741. [Google Scholar] [CrossRef]
  5. Wilson, M.B.; Kallogjeri, D.; Joplin, C.N.; Gorman, M.D.; Krings, J.G.; Lenze, E.J.; Nicklaus, J.E.; Spitznagel, E.E.; Piccirillo, J.F. Ecological momentary assessment of tinnitus using smartphone technology: A pilot study. Otolaryngology 2015, 152, 897–903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Colombo, D.; Palacios, A.G.; Alvarez, J.F.; Patané, A.; Semonella, M.; Cipresso, P.; Kwiatkowska, M.; Riva, G.; Botella, C. Current state and future directions of technology-based ecological momentary assessments and interventions for major depressive disorder: Protocol for a systematic review. Syst. Rev. 2018, 7, 1–26. [Google Scholar] [CrossRef] [PubMed]
  7. Schlee, W.; Pryss, R.C.; Probst, T.; Schobel, J.; Bachmeier, A.; Reichert, M.; Langguth, B. Measuring the moment-to-moment variability of Tinnitus: The TrackYourTinnitus smart phone app. Front. Aging Neurosci. 2016, 8, 1–8. [Google Scholar] [CrossRef] [Green Version]
  8. Shiffman, S.; Stone, A.A.; Hufford, M.R. Ecological Momentary Assessment. Annu. Rev. Clin. Psychol. 2008, 4, 1–32. [Google Scholar] [CrossRef] [PubMed]
  9. Pryss, R. Prospective crowdsensing versus retrospective ratings of tinnitus variability and tinnitus–stress associations based on the TrackYourTinnitus mobile platform. Int. J. Data Sci. Anal. 2019, 8, 327–338. [Google Scholar] [CrossRef] [Green Version]
  10. Bruckert, S.; Finzel, B.; Schmid, U. The Next Generation of Medical Decision Support: A Roadmap Toward Transparent Expert Companions. Front. Artif. Intell. 2020, 3, 75. [Google Scholar] [CrossRef] [PubMed]
  11. Alther, M.; Reddy, C.K. Clinical decision support systems. Healthc. Data Anal. 2015, 625–656. [Google Scholar] [CrossRef] [Green Version]
  12. Caballero-Ruiz, E.; García-Sáez, G.; Rigla, M.; Villaplana, M.; Pons, B.; Hernando, M.E. A web-based clinical decision support system for gestational diabetes: Automatic diet prescription and detection of insulin needs. Int. J. Med. Inform. 2017, 102, 35–49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Feuillâtre, H.; Auffret, V.; Castro, M.; Le Breton, H.; Garreau, M.; Haigron, P. Study of similarity measures for case-based reasoning in transcatheter aortic valve implantation. Comput. Cardiol. 2017, 44, 1–4. [Google Scholar] [CrossRef]
  14. Graham, A.K.; Greene, C.J.; Kwasny, M.J.; Kaiser, S.M.; Lieponis, P.; Powell, T.; Mohr, D.C. Coached Mobile App Platform for the Treatment of Depression and Anxiety Among Primary Care Patients. JAMA Psychiatry 2020, 77, 906–914. [Google Scholar] [CrossRef] [PubMed]
  15. Khorakhun, C.; Bhatti, S.N. Wellbeing as a proxy for a mHealth study. In Proceedings of the 2014 IEEE International Conference on Bioinformatics and Biomedicine (BIBM), Belfast, UK, 2–5 November 2014; pp. 32–39. [Google Scholar]
  16. Li, J.; Zhang, C.; Li, X.; Zhang, C. Patients’ emotional bonding with MHealth apps: An attachment perspective on patients’ use of MHealth applications. Int. J. Inf. Manag. 2020, 51, 102054. [Google Scholar] [CrossRef]
  17. Lin, Y.; Tudor-Sfetea, C.; Siddiqui, S.; Sherwani, Y.; Ahmed, M.; Eisingerich, A.B. Effective behavioral changes through a digital mHealth app: Exploring the impact of hedonic well-being, psychological empowerment and inspiration. JMIR mHealth uHealth 2018, 6, e10024. [Google Scholar] [CrossRef] [PubMed]
  18. De Korte, E.M.; Wiezer, N.; Janssen, J.H.; Vink, P.; Kraaij, W. Evaluating an mHealth app for health and well-being at work: Mixed-method qualitative study. JMIR mHealth uHealth 2018, 6, e6335. [Google Scholar] [CrossRef]
  19. Aboelmaged, M.; Hashem, G.; Mouakket, S. Predicting subjective well-being among mHealth users: A readiness–value model. Int. J. Inf. Manag. 2021, 56, 102247. [Google Scholar] [CrossRef]
  20. Kraft, R.; Schlee, W.; Stach, M.; Reichert, M.; Langguth, B.; Baumeister, H.; Probst, T.; Hannemann, R.; Pryss, R. Combining mobile crowdsensing and ecological momentary assessments in the healthcare domain. Front. Neurosci. 2020, 14, 164. [Google Scholar] [CrossRef]
  21. Insel, T.R. Digital phenotyping: Technology for a new science of behavior. JAMA 2017, 318, 1215–1216. [Google Scholar] [CrossRef]
  22. Langguth, B.; Goodey, R.; Azevedo, A.; Bjorne, A.; Cacace, A.; Crocetti, A.; Del Bo, L.; De Ridder, D.; Diges, I.; Elbert, T.; et al. Consensus for tinnitus patient assessment and treatment outcome measurement: Tinnitus Research Initiative meeting, Regensburg, July 2006. Prog. Brain Res. 2007, 166, 525–536. [Google Scholar] [CrossRef] [Green Version]
  23. Wilson, D.R.; Martinez, T.R. Improved heterogeneous distance functions. J. Artif. Intell. Res. 1997, 6, 1–34. [Google Scholar] [CrossRef]
  24. Hielscher, T.; Spiliopoulou, M.; Völzke, H.; Kühn, J.P. Using participant similarity for the classification of epidemiological data on hepatic steatosis. Proc. IEEE Symp. Comput. Based Med. Syst. 2014, 1–7. [Google Scholar] [CrossRef]
  25. Unnikrishnan, V.; Beyer, C.; Matuszyk, P.; Niemann, U.; Pryss, R.; Schlee, W.; Ntoutsi, E.; Spiliopoulou, M. Entity-level stream classification: Exploiting entity similarity to label the future observations referring to an entity. Int. J. Data Sci. Anal. 2020, 9, 1–15. [Google Scholar] [CrossRef]
  26. Huang, H.; Mehrotra, K.; Mohan, C.K. Rank-based outlier detection. J. Stat. Comput. Simul. 2013, 83, 518–531. [Google Scholar] [CrossRef] [Green Version]
  27. Dogan, A.; Birant, D. A Two-Level Approach based on Integration of Bagging and Voting for Outlier Detection. J. Data Inf. Sci. 2020, 5, 111–135. [Google Scholar] [CrossRef]
  28. Salehi, M.; Rashidi, L. A Survey on Anomaly detection in Evolving Data: [with Application to Forest Fire Risk Prediction]. ACM sigKDD Explor. Newsl. 2018, 20, 13–23. [Google Scholar] [CrossRef]
  29. Pokrajac, D.; Lazarevic, A.; Latecki, L.J. Incremental local outlier detection for data streams. In Proceedings of the 2007 IEEE Symposium on Computational Intelligence and Data Mining, Honolulu, HI, USA, 1 March–5 April 2007; pp. 504–515. [Google Scholar]
  30. Shneiderman, B. The eyes have it: A task by data type taxonomy for information visualizations. In Proceedings of the 1996 IEEE Symposium on Visual Languages, Boulder, CO, USA, 3–9 September 1996; pp. 336–343. [Google Scholar] [CrossRef] [Green Version]
  31. Badam, S.K.; Zhao, J.; Sen, S.; Elmqvist, N.; Ebert, D. TimeFork: Interactive prediction of time series. Conf. Hum. Factors Comput. Syst. Proc. 2016, 5409–5420. [Google Scholar] [CrossRef]
  32. Zhao, Y.; Nasrullah, Z.; Li, Z. Pyod: A python toolbox for scalable outlier detection. arXiv 2019, arXiv:1901.01588. [Google Scholar]
  33. Goldstein, M.; Uchida, S. A comparative evaluation of unsupervised anomaly detection algorithms for multivariate data. PLoS ONE 2016, 11, e0152173. [Google Scholar] [CrossRef] [Green Version]
  34. Unnikrishnan, V.; Schleicher, M.; Shah, Y.; Jamaludeen, N.; Pryss, R.; Schobel, J.; Kraft, R.; Schlee, W.; Spiliopoulou, M. The Effect of Non-Personalised Tips on the Continued Use of Self-Monitoring mHealth Applications. Brain Sci. 2020, 10, 924. [Google Scholar] [CrossRef]
  35. Hiller, W.; Goebel, G. When tinnitus loudness and annoyance are discrepant: Audiological characteristics and psychological profile. Audiol. Neurotol. 2007, 12, 391–400. [Google Scholar] [CrossRef] [PubMed]
  36. Perez, M.V.; Mahaffey, K.W.; Hedlin, H.; Rumsfeld, J.S.; Garcia, A.; Ferris, T.; Balasubramanian, V.; Russo, A.M.; Rajmane, A.; Cheung, L.; et al. Large-scale assessment of a smartwatch to identify atrial fibrillation. New Engl. J. Med. 2019, 381, 1909–1917. [Google Scholar] [CrossRef] [PubMed]
  37. Agarwal, S.; LeFevre, A.E.; Lee, J.; L’engle, K.; Mehl, G.; Sinha, C.; Labrique, A. Guidelines for reporting of health interventions using mobile phones: Mobile health (mHealth) evidence reporting and assessment (mERA) checklist. BMJ 2016, 352, i1174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Breunig, M.M.; Kriegel, H.P.; Ng, R.T.; Sander, J. LOF: Identifying density-based local outliers. In Proceedings of the 2000 ACM SIGMOD International Conference on Management of Data, Dallas, TX, USA, 16–18 May 2000; Volume 29, pp. 93–104. [Google Scholar]
Figure 1. Overview of methods and components of our approach.
Figure 1. Overview of methods and components of our approach.
Entropy 23 01695 g001
Figure 2. Activity Diagram of proposed visualization system. Main interactions with users (left), system (center), and expert (right) are represented. Various functionalities and actions are shown. Gray boxes indicate actions not included as part of this work.
Figure 2. Activity Diagram of proposed visualization system. Main interactions with users (left), system (center), and expert (right) are represented. Various functionalities and actions are shown. Gray boxes indicate actions not included as part of this work.
Entropy 23 01695 g002
Figure 3. Similarity Overview User Panel.
Figure 3. Similarity Overview User Panel.
Entropy 23 01695 g003
Figure 4. Nearest Neighbor Visualization.
Figure 4. Nearest Neighbor Visualization.
Entropy 23 01695 g004
Figure 5. Tailored visualization for each user.
Figure 5. Tailored visualization for each user.
Entropy 23 01695 g005
Figure 6. Visualization to depict predictions.
Figure 6. Visualization to depict predictions.
Entropy 23 01695 g006
Figure 7. Prediction comparison over mean_day observations of users.
Figure 7. Prediction comparison over mean_day observations of users.
Entropy 23 01695 g007
Figure 8. Near and Far future predictions.
Figure 8. Near and Far future predictions.
Entropy 23 01695 g008
Figure 9. Simulated data: different scenarios for outlier detection.
Figure 9. Simulated data: different scenarios for outlier detection.
Entropy 23 01695 g009
Figure 10. 5 nearest neighbors for outlier user-20.
Figure 10. 5 nearest neighbors for outlier user-20.
Entropy 23 01695 g010
Figure 11. Nearest neighbors of simulated user-444154.
Figure 11. Nearest neighbors of simulated user-444154.
Entropy 23 01695 g011
Table 1. Overview of notation used in following sections.
Table 1. Overview of notation used in following sections.
NotationDescription
ua user of the application
x and zWhenever a sentence refers to more than one user, these users are denoted them as denote as x , z
Uthe set of all users who use the mHealth application, which is the data source of this study
Dthe (sub)set of users used for training the models; D U
Tthe (sub)set of test users, who are not used during model training, but are used to assess model quality; T U
kneighborhood size
N N k ( u ) the k nearest neighbors of u
staticF the feature space of the registration data, with s r being a feature, for r = 1 , , | staticF |
Fthe feature space of EMA recordings, with f j being a feature, for j = 1 , | F | ; | F | = 8
t i a timepoint
f i A feature belonging to EMA feature space F
l 1 the timepoints in the future for which we make predictions, i.e., we predict o i + l at t i for user u
f i ^ A predicted value for feature f i F
Table 2. Average of WA—weighted average; LR—linear regression over combinations.
Table 2. Average of WA—weighted average; LR—linear regression over combinations.
Feature CombinationsMax_DAYMean_DAYMin_DAY
WALRWALRWALR
C10.234070.23320.211050.20960.204810.20414
±0.00656±0.0059±0.0084±0.00948±0.00857±0.01003
C20.21560.21880.194030.195390.192520.1973
±0.00808±0.0048±0.00859±0.00529±0.00899±0.00454
C30.20770.20990.188820.191910.189480.19619
±0.00925±0.00923±0.00745±0.007320±0.0116±0.0122
Table 3. Frequencies and prior probabilities of target variable tinnitus distress-(s03).
Table 3. Frequencies and prior probabilities of target variable tinnitus distress-(s03).
CategoryNum UsersTotal Number of Observations (s03)Probability
Low4519110.497
Moderate4811550.30
High377730.201
Total-38391.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Prakash, S.; Unnikrishnan, V.; Pryss, R.; Kraft, R.; Schobel, J.; Hannemann, R.; Langguth, B.; Schlee, W.; Spiliopoulou, M. Interactive System for Similarity-Based Inspection and Assessment of the Well-Being of mHealth Users. Entropy 2021, 23, 1695. https://0-doi-org.brum.beds.ac.uk/10.3390/e23121695

AMA Style

Prakash S, Unnikrishnan V, Pryss R, Kraft R, Schobel J, Hannemann R, Langguth B, Schlee W, Spiliopoulou M. Interactive System for Similarity-Based Inspection and Assessment of the Well-Being of mHealth Users. Entropy. 2021; 23(12):1695. https://0-doi-org.brum.beds.ac.uk/10.3390/e23121695

Chicago/Turabian Style

Prakash, Subash, Vishnu Unnikrishnan, Rüdiger Pryss, Robin Kraft, Johannes Schobel, Ronny Hannemann, Berthold Langguth, Winfried Schlee, and Myra Spiliopoulou. 2021. "Interactive System for Similarity-Based Inspection and Assessment of the Well-Being of mHealth Users" Entropy 23, no. 12: 1695. https://0-doi-org.brum.beds.ac.uk/10.3390/e23121695

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