Fluvoxamine treatment response prediction in obsessive-compulsive disorder: association rule mining approach

Abstract

Background

Obsessive-compulsive disorder (OCD) is a debilitating psychiatric disorder characterized by intrusive thoughts or repetitive behaviors. Clinicians use serotonin reuptake inhibitors (SRIs) for OCD treatment, but 40%–60% of the patients do not respond to them adequately. Here, we described an association rule mining approach for treatment response prediction using an Iranian OCD data set.

Patients and methods

Three hundred and thirty OCD patients fulfilling DSM-5 criteria were initially included, but 151 subjects completed their pharmacotherapy which was defined as 12-week treatment with fluvoxamine (150–300 mg). Treatment response was considered as >35% reduction in the Yale-Brown Obsessive Compulsive Scale (Y-BOCS) score. Apriori algorithm was applied to the OCD data set for extraction of the association rules predicting response to fluvoxamine pharmacotherapy in OCD patients. We considered the association of each attribute with treatment response using interestingness measures and found important attributes that associated with treatment response.

Results

Results showed that low obsession and compulsion severities, family history of mental illness, illness duration less than 5 years, being married, and female were the most associated variables with responsiveness to fluvoxamine pharmacotherapy. Meanwhile, if an OCD patient reported a family history of mental illness and his/her illness duration was less than 5 years, he/she responded to 12-week fluvoxamine pharmacotherapy with the probability of 91%. We also found useful and applicable rules for resistant and refractory patients.

Conclusion

This is the first study where association rule mining approach was used to extract predicting rules for treatment response in OCD. Application of this method in personalized medicine may help clinicians in taking the right therapeutic decision.

Keywords:

• data mining
• apriori algorithm
• family history
• refractory patients

Introduction

Obsessive-compulsive disorder (OCD) is a heterogeneous neuropsychiatric condition that affects 3% of the population worldwide^1–3 and ~1.8% in Iran.⁴ Overall, the prevalence of OCD is slightly more in women than men with an earlier onset in the male patients.⁵

Evidence indicated that OCD patients compared to healthy subjects have functional impairment and poorer quality of life⁶ which is lower for samples with higher percentages of females and was negatively associated with the age.⁷ However, a significant number of OCD patients (40%–60%) fail to benefit sufficiently from pharmacologic treatments.^8–11 Different factors influence the effectiveness of OCD treatment suggesting personalized medicine as a solution. Research studies provide vast amounts of heterogeneous data that are too complex and voluminous to be processed and analyzed by traditional methods. Although the conventional methods have their application in medical data analyses, nonetheless, they have their limitations and we need inductive approaches such as data mining methods to discover knowledge with a small number of cause–effect mechanisms.¹²

Data mining approach is becoming increasingly popular for providing a deeper understanding of medical data, including disease pathogenesis and treatment leading to a new discovery from the medical data set that the conventional methods are unable to process due to their limitations.^13–16 Data mining algorithms are classified into two main categories: supervised and unsupervised learning.¹⁷ Supervised learning algorithms such as classification or regression predict the response values for a particular outcome. Unsupervised learning algorithms describe the form and hidden structure of data, using methods such as clustering and anomaly detection. One of the primary unsupervised algorithms is association rule mining that has been widely explored in different areas, including medicine. Association rule mining discovers associations and correlations among large sets of data and provides information in the form of if-then probabilistic statements.^15,18 Different algorithms have been proposed to generate association rules.^19–21 Apriori is a well-known primary association rule mining algorithm for extraction of frequent items in a set of transactions. It was first introduced by Agrawal et al.¹⁹ Sharma and Om used this algorithm for early diagnosis and treatment of oral cancer.²² Karabatak and Ince also proposed an expert system for detection of breast cancer based on association rules and neural network.¹⁴

Predictive Apriori,²¹ motivated by Apriori, maximizes the accuracy of an association rule in contrast to the confidence related to focus of Apriori. Apriori ranks the rules based on confidence only, whereas predictive Apriori considers both confidence and support in the rule ranking. Nahar et al consider three rule generation algorithms – Apriori, predictive Apriori, and Tertius²⁰ – to extract cancer prevention factors.²³ Their results showed that Apriori is the most useful association rule mining algorithm to be used in the discovery of prevention factors. In another study by this group, they used association rule mining approach to identify factors that contributed to heart disease considering sex.²⁴ Chen et al¹² used association rule mining to detect possible side effects due to the exposure time to drugs during pregnancy, which resulted in the discovery of interesting information.

There are limited studies in the field of psychiatry that use association rule mining approach for data extraction. Shen et al applied association rule mining to the Taiwan National Health Insurance Research Database to explore associations among comorbidities of borderline personality disorder.²⁵ In another study, Schweikert et al applied Combinatorial Fusion Analysis (CFA) and association rule mining to explore the relationship between autism prevalence and lead mercury concentration, which led to a deeper understanding of its pathogenesis.²⁶ In the current study, we applied data mining algorithms to an OCD treatment response data set of Iranian patients to explore the most significant factors contributing to the treatment response. At first, we converted our data set variables which were chosen from previous studies^9,27–31 and based on experts’ opinions to an appropriate format for application of Apriori algorithm. Next, we evaluated the association of each attribute with treatment response using support, confidence, and lift measures and found the most critical attributes that were associated with treatment response in OCD for pharmacotherapy with fluvoxamine. Finally, we used Apriori algorithm to explore the most associated rules which predict treatment response in OCD.

Patients and methods

The preprocessing and analyses of all data were conducted using Matlab version 2014a (The Math Works Inc.).

Subjects and treatment protocol

Three hundred and thirty outpatients with an Iranian origin and meeting DSM-4 criteria for OCD were recruited from Imam Hossein hospital (Tehran, Iran) between 2014 and 2017. OCD diagnosis was established by a psychiatrist, and its subtypes and severity were investigated by an expert psychologist through the structured interview using Yale-Brown Obsessive Compulsive Scale (Y-BOCS) checklist and severity scale.³² The sociodemographic data were collected through a questionnaire consisting of the patient’s full name, address, age at assessment, the age of onset, gender, marital status, educational level, occupation, chief complaint, illness duration, history of substance usage, and the family history of psychiatric disorders, specifically OCD. Participants were in the age range of 18–65 years, having OCD symptoms for more than 1 year and with total Y-BOCS severity score more than 9 (suggested cutoff point for Iranian patients³³); a total of 17 patients out of 151 had total Y-BOCS scores 9–16. The exclusion criteria were having a history of psychotic disorders or intellectual disability; reporting severe neurological pathology; history of substance usage; diagnosis with other DSM-4 Axis l disorders except for depression, anxiety, or tic disorder; and under other pharmacotherapy or cognitive behavioral therapy (CBT).

Pharmacotherapy was defined as a 12-week treatment with fluvoxamine (150–300 mg), which has been described previously.²⁸ Briefly, the escalating fluvoxamine daily dose was initiated from 25 mg/day and increased up to 100 mg/day in the third week. Fluvoxamine daily dose for the next 3 weeks was 150 mg/day, and after the sixth week, patients were visited by the psychiatrist and received maximum tolerated dose of the fluvoxamine for the next 6 weeks. No concomitant therapy was allowed during the whole treatment period, either pharmacological or non-pharmacological. Y-BOCS severity was evaluated at the beginning of the treatment period and after 12 weeks, and based on the reduction in patients’ Y-BOCS score, the patients were divided into two groups: group A (responders) included patients who showed >35% reduction in Y-BOCS scores after treatment with fluvoxamine, and group B (non-responders) was composed of patients who exhibited <35% reduction in Y-BOCS scores.³⁴ We included another group (refractory patients) consisting of patients who experienced various selective serotonin reuptake inhibitor (SSRI) trials at maximum tolerated dose during their illness period but did not respond to them adequately.³⁴

From 330 OCD patients initially included in the study, 151 persons completed their pharmacotherapy. Others were excluded because of refusal to come for follow-up (N=108), not taking medication properly (N=38), discontinuity in pharmacotherapy due to complications, including allergy to the drug (N=8), meeting exclusion criteria during pharmacotherapy (N=6), and other reasons (N=19). Since we used the previously published data set to investigate the association rule data mining approach, the CONSORT diagram which summarizes the flow of participants through different stages of the trial was presented previously.²⁸ The excluded and included patients showed no significant differences in mean baseline scores considering their demographic (sex and marital status) and clinical characteristics (history of psychiatric disorders, age at assessment, age of onset, obsession severity, compulsion severity, and total severity) as well as symptom dimensions from the Y-BOCS symptom checklist (data not shown, all P>0.1).

Sixty-four persons out of the 151 included patients were female, 66% were married, 58% of patients were unemployed, and 76% of them did not have academic degrees. Sixteen percentage of patients reported no history of mental illness in their family, 33% reported a history of mental illness, and 51% were from the families with a history of OCD.

Patients’ age at onset (mean±SD) was 23±10.4 years and age at assessment (mean±SD) was 34.2±10.2 years. The mean of Y-BOCS scores at intake (mean±SD) was 10.84±4.7 for obsession, 9.3±5.39 for compulsion, and 20.1±8.8 for the total severities. This research was performed in accordance with the latest version of Declaration of Helsinki and with the approval of the research ethics committee of Neuroscience Research Center, Shahid Beheshti University of Medical Sciences (IR.SBMU.PHNS.REC.1397.14). Written consent for participation was obtained from participants after being informed about the nature of the research.

Variables and preprocessing

The analyzed OCD data set contains sociodemographic and clinical variables that were considered as independent variables (predictors) and treatment responses (responder, non-responders, refractory) which were considered as dependent variables.

Twenty-two variables were identified as predictors, including sex, marital status, family history, initial Y-BOCS obsession severity, initial Y-BOCS compulsion severity, Y-BOCS obsession subtypes, Y-BOCS compulsion subtypes, avoidance, depression score, and age of onset. These variables have to be converted to an appropriate format (true or false values) for applying associative rule mining algorithm. For example, sex translates into two columns: one is for being male and another for being female or marital status saved to one column named as married. Finally, we had a data set containing 37 columns for each sample, 34 columns for predictors, and three columns for treatment response status (Table 1). The time complexity of Apriori algorithm is exponential and takes too much time to extract associative rules,³⁵ so feature selection may help to overcome these problems. Some studies used association rules to find the most critical features and excluded redundant ones to reduce data set dimensions. Karabatak et al¹⁴ developed an expert system for the detection of breast cancer based on the association rules and neural network. At first, they detected large item sets for each studied class. Next, they removed redundant features and used essential items to classify their data set. They applied this method in another study for diagnosis of erythemato-squamous diseases.³⁶ Using their strategy, we found out that being male, no history of psychiatric disorders and counting appeared in none of the frequent item sets, so we excluded them from our further analysis.

Table 1 Features used in the association rule mining algorithm

No	Feature	No	Feature	No	Feature	No	Feature
1	Female	2	Male	3	Married	4	Early onset (onset <18)
5	Late onset (onset ≥18)	6	No history	7	Family history of mental illness	8	Family history of OCD
9	Aggressive obsession	10	Contamination obsession	11	Sexual obsession	12	Hoarding obsession
13	Religious obsession	14	Symmetry obsession	15	Somatic obsession	16	Cleaning compulsion
17	Checking compulsion	18	Repeating compulsion	19	Counting compulsion	20	Ordering compulsion
21	Obsession severity ≤10	22	Obsession severity >10	23	Compulsion severity ≤10	24	Compulsion severity >10
25	Total severity ≤18	26	Total severity >18	27	Presence of insight	28	No insight
29	Presence of avoidance	30	Lack of avoidance	31	Gaf ≤3	32	Gaf >3
33	Illness duration ≤5 years	34	Illness duration >5 years	35	Responder	36	Non-responder
37	Refractory

Abbreviations: Gaf, global assessment of functioning; OCD, obsessive-compulsive disorder.

Interestingness measures

An association rule has two parts: an antecedent (if) and a consequent (then) in the form X → Y. This rule means that when X occurs, then Y also happens with a certain probability.¹⁹ Usually, an overwhelming number of rules are created in association rule mining, and evaluating this amount of rules by a user is tedious and time-consuming. Interestingness measures are beneficial and play an essential role in the detection of strong and interesting rules from the set of discovered rules. These indexes are intended for selecting and ranking rules according to their potential interest to the user. Confidence and support are two traditional measures of rule interestingness that reflect the certainty and usefulness of a rule, respectively.³⁷ The support is defined as a number of transactions in “data set” that include all items in the antecedent (if) and consequent (then) parts of the rule. The support is sometimes expressed as the percentage of the total number of records in the database (Equation 1$Support(X\Rightarrow Y)=P(X\cap Y)$(1) ). Item sets whose frequencies are greater than a predefined minimum support are called frequent item sets. Confidence is the conditional probability that a randomly selected transaction will consist of all the items in the consequent given that the transaction includes all the items in the antecedent (Equation 2$\begin{array}{l}Confidence(X\Rightarrow Y)=P(Y|X)=\\ (Support(X\cap Y))/(Support(X))\end{array}$(2) ). Confidence is the accuracy of the rule and is used in Apriori to rank the rules. The task of association rule mining is generation off all association rules from a set of transactions that have a support value greater than or equal to minimum support (Minsup) and confidence value greater than or equal to minimum confidence (Minconf). The values of Minsup and Minconf are set by the users and are parameters of the Apriori algorithm for association rule generation.

Rules with the high support and high confidence values are not necessarily impressive. Most of the rules, resulting from applying association rules mining to a medical data set, are irrelevant and do not have clinical importance.³⁸ So another parameter such as lift is used as a metric for evaluation of medical significance and reliability of association rules. Lift is nothing but the ratio of confidence to expect confidence (Equation 3$Lift(X\Rightarrow Y)=(P(X\cap Y))/(P(X)\ast P(Y))$(3) ). A lift ratio larger than 1.0 implies that the relationship between the antecedent and the consequent is more significant than it would be expected and makes those rules potentially useful for predicting the consequent in the future. The larger the lift ratio, the more significant is the association. The support, confidence, and lift values were calculated using the following formulas: $Support(X\Rightarrow Y)=P(X\cap Y)$(1) $\begin{array}{l}Confidence(X\Rightarrow Y)=P(Y|X)=\\ (Support(X\cap Y))/(Support(X))\end{array}$(2) $Lift(X\Rightarrow Y)=(P(X\cap Y))/(P(X)\ast P(Y))$(3) where X and Y are two non-overlapping subsets of items. P(X) denotes the probability that a randomly selected transaction from data set will include all the items in X. X ∩ Y refers to all transactions that contain both X and Y.

Applying Apriori algorithm to the OCD data set for the detection of association rules

The Apriori algorithm is a well-known association rule mining algorithm that works well on complete data sets with discrete values. Since we need to find the associations between attributes and treatment response, we considered OCD data set attributes as the antecedent parts of A → B rules and the treatment response class (responder, non-responder, and refractory) as the consequent part.

The applied association rule mining algorithm consisted of two sections. In the first step, frequent patterns (item sets with the supports greater than predefined Minsup) were generated. In the second step, frequent patterns’ confidences were calculated and those with confidence greater than Minconf were selected as the final rules. The pseudo-code of the used algorithm is presented in Table 2.

Table 2 Pseudo code for Apriori algorithm that is suitable for supervised learning

Function rules=Apriori(T, minsup, minconf, classitem)

% T is the transaction data

% class item is the item that corresponds to class label

F2=Add class item to each element of F1

K=3

Repeat

); {generate candidate k-item sets}

Ck=prune(Ck); {prune candidate item sets}

If t contains c

Count(c)=count(c)+1

end for

end for

Until

//rule generation section

For each candidate item set f

A=f – class item % all items except class item, antecedent part of the rule

B=class item % consequent part of the rule

If confidence(A→B)>minconf

Add A→B to rules

end if

end for

return rules

Choosing a single cutoff point for Minsup and Minconf may cause some problems in conducting Apriori algorithm. If we used high Minsup threshold, Apriori algorithm would not detect association rules for the low-frequency classes, while with the low Minsup, Apriori algorithm took hours to complete and resulted in a lot of rules that most of them were redundant. To deal with this challenge, multiple minimum supports for different classes were assigned. Moreover, minimum support for each class was set by considering the proportion of its distribution over the whole data set.^39,40

Results

Association of OCD data set attributes with the treatment response classes

Table 3 summarizes the most important attributes that associated with different treatment response classes (rules with lift <1 and support <0.08 were excluded and those left were sorted based on their confidence). We find rules in the form of A → B such that A is a single attribute and B is the treatment response class (responder, non-responder, and refractory). As the results showed, no single attribute had very high confidence that could be used alone for the detection of treatment response class.

Table 3 Most important attributes associated with treatment response (support >0.08, lift >1; sorted based on confidence)

Responder		Non-responder		Refractory
Feature	Confidence	Feature	Confidence	Feature	Confidence
Obsession severity ≤10	0.72	Sexual obsession	0.38	Obsession severity >10	0.23
Family history of mental illness	0.72	Illness duration >5 years	0.32	Family history of OCD	0.20
Illness duration ≤5 years	0.68	Being male	0.30	Early onset	0.20
Compulsion severity ≤10	0.67	Compulsion severity >10	0.30	Presence of avoidance	0.18
Married	0.65	Contamination obsession	0.30
Being female	0.64

Abbreviation: OCD, obsessive-compulsive disorder.

Results showed that low obsession severity (confidence=0.72), low compulsion severity (confidence=0.67), family history of mental illness (confidence=0.72), illness duration less than 5 years (confidence=0.68), being married (confidence=0.64), and being female (confidence=0.65) have strong associations with treatment responsiveness. For example, if an OCD patient has low obsession severity, he/she will be a responder to fluvoxamine pharmacotherapy with a probability equal to 0.72.

The most critical attributes associated with non-responsiveness were the sexual obsession, illness duration more than 5 years, being male, high compulsion severity, and having contamination obsession. High obsession severity, family history of OCD, early-onset OCD, and the presence of avoidance also showed association with treatment refractory class. However, the confidence values for associated attributes with non-responder and refractory classes were low (0.38–0.18).

Association rules for OCD treatment response

The most important rules associated with each class that were extracted by Apriori algorithm are summarized in Table 4. For responder class, selected rules were those with confidence levels above 70%, support levels above 15%, and lifts greater than 1. The most associated rule with response to fluvoxamine pharmacotherapy (the family history of mental illness, aggressive obsessions, and illness duration less than 5 years → responder) implies that these three attributes are good predictors of being a responder. If someone has these characteristics simultaneously, he/she will be a responder with the probability of 92%. Comparing this rule with the following rules in Table 4 and results of the single attribute association analysis with treatment response class (Table 3) showed that a family history of mental illness and illness duration less than 5 years are the most important predictors of responder class. Late-onset OCD, being female, and married as well as low obsession and compulsion severities are other attributes that appeared in the subsequent rules, showing their importance in treatment responsiveness.

Table 4 Best discovered rules for different treatment responses

Class	Rule number	Antecedent	Supp	Conf
Responder	1	Family history of mental illness and Aggressive and Illness duration ≤5	0.16	0.92
	2	Family history of mental illness and Illness duration ≤5	0.19	0.91
	3	Late onset and Family history of mental illness and Illness duration ≤5	0.16	0.89
	4	Female and Married and Late onset and Compulsion severity ≤10	0.16	0.85
	5	Married and Obsession severity ≤10 and Compulsion severity ≤10	0.19	0.82
	6	Married and Hoarding and Obsession severity ≤10	0.16	0.82
Non-responder	7	Contamination and Cleaning and Repeating and Illness duration >5	0.08	0.73
	8	Contamination and Obsession severity >10 and Total severity >18 and Illness duration >5	0.08	0.69
	9	Cleaning and Illness duration >5	0.08	0.69
	10	Cleaning and checking and Obsession severity >10 and Compulsion severity >10	0.07	0.68
	11	Cleaning and checking and Obsession severity >10 and Compulsion severity >10 and Total severity >18	0.07	0.68
Refractory	12	Female and Early onset and Aggressive and Obsession severity >10	0.04	0.83
	13	Female and Married and Early onset and Checking and Gaf ≤3	0.04	0.71
	14	Family history of OCD and Obsession severity >10 and Presence of avoidance and Gaf >3 and Illness duration ≤5	0.04	0.71
	15	Female and Married and Early onset and aggressive and checking and Gaf ≤3	0.03	0.71
	16	Female and Married and Early onset and checking and Gaf ≤3 and Illness duration >5	0.03	0.71

Abbreviations: Conf, confidence; Gaf, global assessment of functioning; OCD, obsessive-compulsive disorder; Supp, support.

As previously mentioned, we set different Minsup and Minconf for each class. So for the non-responder class, Minsup and Minconf were considered as 0.07 and 0.6, respectively, and the most important rules are mentioned in Table 4. As it is shown in Table 4, the number of attributes that appeared in the antecedent part of the rules is more than the attributes in the responder class rules. For instance, if someone has five attributes, including contamination obsession, cleaning, repeating compulsions, and illness duration more than 5 years, he/she will not respond to treatment with the probability of 73%. High obsession and compulsion severities were the other attributes that appeared in the rules of non-responder class.

For refractory class, we set Minsup =0.03 and Minconf =0.6, and extracted rules were sorted based on their confidence (Table 4). The most important rule was detected as being female, early onset, aggressive obsession, and obsession severity more than 10. This rule implies that if a woman has an early-onset OCD with great obsession severity and aggressive obsession, she will not respond to most of the treatment strategies adequately with the probability of 83%. Early-onset OCD is an attribute that appeared in almost all the rules in this class, implicating the importance of this attribute for the refractory class.

Discussion

This research has evaluated the efficacy of association rule mining for the first time to predict treatment response in OCD for pharmacotherapy with fluvoxamine using real data. In this study, we applied Apriori algorithm to the OCD pharmacotherapy data set and extracted the most associated features and the most important rules for each treatment response class based on the highest confidence level. Since the results of association rule data mining are in the form of if-then statements, their interpretability is high and easy to use in personalized medicine.

Personalized medicine centers on the assumption that each person has a unique physiologic characteristic which plays a significant role in both disease susceptibility and response to specific treatment. So the significant goals of personalized medicine are the prediction of an individual’s vulnerability to developing an illness, achieve an accurate diagnosis, and optimize the most efficient response to treatment with the least side effects.⁴¹

Personalized medicine has vital importance in psychiatric disorders, since they are multifactorial diseases with a variable clinical course and considerable inter-individual variability in patient responses to therapies. About a third of patients respond to antidepressants, so psychiatrists use their clinical skill and experience to adapt the medical protocol for each patient.⁴² However, in many cases, patients do not respond adequately to the medication or they report considerable side effects that result in a change in pharmacotherapy protocol. So, a series of drug alterations increase treatment duration as well as social and economic costs. Detection of predictors of treatment response may help psychiatrists to prescribe the most effective treatment with the least side effects considering individual characteristics.⁴³

The first-line medications for OCD are antidepressants such as SSRIs; however, ~50% of OCD patients show a poor response.⁴⁴ Genetic factors as well as environmental (eg, smoking, foods, and beverages) and personal (eg, age, gender, ethnicity, and illnesses) variables affect the response to treatment.⁴⁵ Various studies detected different predictors for poor treatment response, including living without a partner,⁴⁶ sexual subtype,⁴⁷ higher OCD severity at intake,^48,49 poor insight into obsessions, symmetry/hoarding and contamination/washing dimension, and the presence of specific personality disorders.⁵⁰ Moreover, anterior cingulate cortex, orbitofrontal cortex, caudate, putamen, and thalamus volume changes and their altered activity were reported as predictors of poor response.⁵⁰ There are several pharmacogenetic reports related to the importance of genes involving serotonergic and glutamatergic signaling pathways, but their results are inconsistent owing to the small sample size, lack of standardized methodology, and rigorous statistical analyses.⁴⁴ In the case of non-responsiveness, there are other strategies to augment the treatment response ranging from adding CBT to pharmacotherapy augmentations, such as antipsychotics, a combination of another serotonin reuptake inhibitors (SRIs), anticonvulsants, or even the use of neurosurgical procedures.⁵¹ For example, the presence of anger associated with sexual obsessions and contamination may require the use of mood stabilizers as augmentation.⁵² OCD patients may also benefit from cognitive (CT) and behavioral therapy (BT) separately or in combinations (CBT) which have been proven to be efficacious in multiple studies.^53,54 A recent meta-analysis comparing the efficacy of cognitive and behavioral therapies separately and in combination (CBT) in OCD treatment found no predictors for treatment response, but they indicated that improved patients had lower baseline depression and stronger beliefs about responsibility/threat and importance/control of thoughts. They showed that higher baseline scores on depression adversely affected outcomes for BT but not CT or CBT, and higher education and lower OCD severity were associated with better response to CT only.⁵⁵

Data mining methods are useful tools in personalized medicine for the discovery of new knowledge, appropriate to diagnosis and treatment. In the current study, Apriori algorithm that is a well-known association rule mining algorithm was used to extract the most essential rules for the prediction of treatment response in OCD patients. Results showed that low obsession and compulsion severities, family history of mental illness, illness duration less than 5 years, being married, and being female are the most associated variables with responsiveness to fluvoxamine pharmacotherapy. Our findings are consistent with previous studies.^{9,29,30,46,49} If an OCD patient reports a family history of mental illness and his/her illness duration is less than 5 years, he/she will respond to 12-week fluvoxamine pharmacotherapy with a probability of 91%. Although other clinical variables appeared in the associated rules, including late-onset OCD, being female and married as well as low obsession and compulsion severities are good predictors of responsiveness to treatment with fluvoxamine. Extracted rules revealed that if an OCD patient has contamination obsession, cleaning and repeating compulsions, and illness duration more than 5 years, he/she will not respond to fluvoxamine pharmacotherapy with the probability of 73%. As the psychiatric disorders are multifactorial diseases, it is highly unlikely that a single variable predicts treatment response in mental disorders, so considering these rules could be very useful for making clinical decisions.

To the best of our knowledge, this is the first study that applied association rule mining to detect rules associated with treatment response. This strategy has the advantage over traditional analytical approaches because of its higher accuracy. Moreover, the formats of results are easy to use by clinicians. The other strength of this study is prescribing a single drug, fluvoxamine, which avoided having confounded factors related to the different SSRI treatments. However, similar to any other research, this study also has some limitations. First, the small sample size which makes it necessary to apply this method to similar data sets with more samples. Second, although comparison of completers with non-completers considering their mean baseline characteristics showed no significant differences between these two groups, a large number of patients who did not complete their pharmacotherapy make it difficult to generalize these rules to the Iranian OCD patients. Third, treatment response was evaluated after 12-week pharmacotherapy, while assessment of treatment response’s outcome in more extended period may detect stronger predictors of fluvoxamine responsiveness. Fourth, previous pharmacological or non-pharmacological treatments may affect response to fluvoxamine treatment. Fifth, baseline systematic assessments were not performed routinely for the included subjects. Altogether, it is suggested that this approach should be applied in a similar study to achieve more accurate rules across different populations.

Conclusion

The application of data mining methods, especially association rule data mining coupled with personalized medicine, helps to extract feasible rules considering individual clinical features in disease diagnosis and treatment.

Author contributions

All authors contributed toward data analysis, drafting and critically revising the paper, gave final approval of the version to be published and agree to be accountable for all aspects of the work.

Acknowledgments

This research was supported by a grant from Shahid Beheshti University of Medical Sciences, Iran (grant No 14553-6). We would like to thank patients and their families for participating in this study and Imam Hossain staff for their cooperation with our research team. We also gratefully thank NRC and Tarbiat Modarres graduate students and staff for participating in this study as a control group.

Disclosure

The authors report no conflicts of interest in this work.