Abstract

Background. Psoriasis is a complex skin disease and difficult to evaluate, and this study aimed to provide an objective and systematic approach for evaluating the efficacy of psoriasis. Methods. We sought to construct a Bayesian network from sixteen indicators in four aspects of psoriasis (skin lesion conditions, laboratory indexes, quality of life, and accompanying symptoms) and obtained weights of each index by combining the analytic hierarchy process with maximum entropy self-learning. Furthermore, we adopted stability analysis to calculate the minimum sample size of the system. The extended set pair analysis was utilized to evaluate the efficacy based on improved weights, which overcomes the limitation of set pair analysis (unable to evaluate the efficacy with uncertain grades and thresholds). Results. A total of 100 psoriasis vulgaris patients were included to evaluate the curative effect by the system. We obtained the weights of each index and the Euclidean distance for efficacy evaluation of 100 patients. The sensitivity analysis proved that the results had no significant change with the variation of single patient’s indexes, which indicated that our results were stable to assess the effectiveness. Conclusions. We provided an available method of comprehensive effective evaluation of various indicators of psoriasis and based on both subjective and objective weights.

1. Introduction

Psoriasis is a complex, chronic, immune-mediated inflammatory skin disorder that affects approximately 125 million people worldwide [1]. Moreover, complications are associated with increased exacerbations in subjects with psoriasis, including diabetes, metabolic syndrome, and chronic obstructive pulmonary disease [24]. As a refractory systemic disease, psoriasis has a great impact on human health, and some can even be life threatening [5, 6]. Due to the complexity and severity, it is currently difficult to evaluate the efficacy of psoriasis [7]. Although varied evaluation tools and models have been practiced to evaluate the efficacy of diseases, there are still some shortcomings including the inability to integrate multiple indicators to construct an evaluation system and the lack of impartiality in evaluating the weights of various indicators [810].

In the Bayesian network, a directed acyclic graph is constructed to intuitively reflect the potential relationship between factors, and a conditional probability distribution table is used to reflect the strength of association [11]. A Bayesian network can reflect the multifactor relationship of psoriasis, so we adopted it to evaluate the efficacy. Meanwhile, the importance of each factor affecting the evaluation of the curative effect is dissimilar, so a reasonable weight needs to be provided [12]. The Analytic Hierarchy Process (AHP) is based on expert experience and was relatively prejudiced [13]. Interestingly, maximum entropy can identify the probability distribution that is most consistent with the cost function and makes the fewest assumptions [14]. To obtain the comparatively actual weight, we combined AHP with maximum entropy to maximize the entropy of the evaluation network through self-learning [15]. On the basis of the evidence mentioned above, we evaluated the efficacy quantitatively through extended set pair analysis (ESPA) based on multifactorial network and reasonable weights [16].

In the current study, we aimed to develop a comprehensive efficacy evaluation system of psoriasis vulgaris, based on Bayesian maximum entropy self-learning and ESPA. A total of 100 patients were included by this system for efficacy evaluation. It is expected to provide a new approach for curative effectiveness evaluation of psoriasis and other complex diseases.

2. Methods

2.1. Patients

Adult patients (between 18 and 65 years of age) were eligible to participate if they satisfied the condition of both Western and traditional Chinese medicine (TCM) diagnosis standards for psoriasis vulgaris. The trial was performed in accordance with our previous study [17]. The study was approved by the Ethics Committee of the Yueyang Hospital of Integrated Traditional Chinese and Western Medicine (approval no. 2019-028). All participants provided written informed consent before entering the study.

2.2. Case Study Design

Eligible patients received oral TCM herbal medication tailored to the participant’s disease progression. Medication was administered twice every day during the intervention phase. After 8 weeks of treatments, clinical efficacy was assessed, and blood samples were collected for all eligible patients. The measurement of laboratory indexes was detected in accordance with our previous study [18], which is described in the supplementary materials.

2.3. Bayesian Network Construction

A Bayesian network consists of a directed acyclic graph (DAG) and a series of conditional probability tables (CPTs). The nodes represent random variables, and edges represent the conditional dependences among nodes in a DAG [19]. In our study, the nodes were the indexes related to psoriasis. The conditional probability of each node was obtained according to the relationship between the indexes. An overview of the study flow is shown in Figure 1.

2.4. Calculation of the Initial Weight (IW) by the AHP

We invited experts to score according to the significance of the indexes (Supplementary Table S1). Then, the IWs of indexes were calculated by the AHP according to the following steps [20]. The scores were regarded by the expert, and the IWs were back in calculation if they failed to pass the consistency test.

2.5. Calculation of Maximum Entropy and the Final Weight (FW) by Self-Learning

In the process of self-learning, the maximum entropy was taken as the output condition by the gradient descent [21]. The IWs were input into the system to self-learn. Finally, the FW of the indexes was calculated according to the Bayesian network.

2.6. Efficacy Evaluation by ESPA

The efficacy of patients was evaluated by ESPA as described in the previous studies [16]. For the data of included patients in the index , the greatest and smallest value of the data represented the upper threshold and lower threshold , respectively. The assuming arbitrary value belongs to . In this study, if the smaller values mean a better level of the effect, the connection degree (CD) of the patient l with respect to the index k could be calculated by equation (1). Conversely, if the greater values of data mean a lower level of efficacy, then the CD could be calculated by equation (2) [22].where denotes CD of the patient l with respect to the index k, indicates the uncertainty coefficient of discrepancy and its value range is , and connotes the contradictory coefficient with the value defined as .

After that, the similarity degree (SD) was proposed to reflect a couple of patients’ similarities, and the range of SD was . Then, the ideal patient denoted a patient with the optimal response and was availed to make a comparison with the patients who are evaluated for the efficacy. As set pair analysis (SPA) defines [23], the CD of an ideal patient was calculated by the following equation:where denotes the CD of an ideal patient.

Moreover, the Euclidean distance (ED) was set to evaluate the SD of each patients, and the specific calculation method was mentioned in Yan’s study [16]. The comprehensive evaluation was confirmed, and the weights were also considered in the calculation of the ED [3]. The formula is given as follows:where denotes the ED between the assessed patient and the ideal patient, connotes the CD of the assessed patient in index , represents the SD of the CD for patient and the ideal patient, and index is the FW of index .

2.7. Stability and Sensitivity Analysis

The stability analysis was practiced to test the rationality of the sample size. When the sample size changed, the corresponding entropy and self-learning times were determined by the Bayesian maximum entropy self-learning model. The sample size was set from 5 to 100, and the interval was 5. This process was repeated 20 times.

As inputs of some indexes are uncertain, the results were affected by the uncertainties. Consequently, a sensitivity analysis was performed to check the consistency. Each patient was set as the variable sample in turn, and data in each index were set as the variable inputs. Assuming an error of ±10% in the inputs determined, that is to say, the range of input values was between 90% and 110% of the reference values [24]. Herein, the interval of input values was set as 1%. Then, ED with different input values were made by ESPA.

3. Results

3.1. Results of the Index Weights

A total of 100 psoriasis vulgaris patients were included from outpatients of 20 clinical centers. The characteristics of the patients and 16 indicators are shown in Table 1. The details are provided in Supplementary Table S2.

To assess the efficacy of psoriasis vulgaris, we sought to construct the Bayesian network with four layers (Figure 2). The top layer was efficacy evaluation, that is, the target value of the network. The intermediate layers corresponded to the attribute values for the upper layer, and the target values for the next layer. The second layer consisted of four aspects of psoriasis vulgaris, including skin lesion conditions, laboratory indexes, quality of life, and accompanying symptoms. The third layer contained psoriasis area and severity index (PASI), body surface area (BSA) [25], squamous cell carcinoma antigen (SCC-Ag), tumor necrosis factor-α (TNF-α), interleukin (IL) and complement levels [26], Dermatology Life Quality Index (DLQI), Self-Rating Anxiety Scale (SAS) and Self-Rating Depression Scale (SDS) [27], the scales of Xerostomia Questionnaire (XQ), Cleveland Clinic Score (CCS), and Pittsburgh Sleep Quality Index (PSQI) [2831]. Of these, IL comprised of IL-10, IL-17, IL-22, and IL-23, and the complement involved complement 3 (C3) and complement 4 (C4) [26]. Each node of the bottom layer was the attribute value of assessment, which was as well the outcome measure of this study.

Each abovementioned index contributed variously to the efficacy evaluation. Therefore, we confirmed the IW by the AHP according to expert score (Supplementary Table S3). Then, we performed consistency checks (Supplementary Table S4). The results passed the consistency check, so the weights were available for the efficacy evaluation, whereas the AHP required repeated artificial modification of the judgment matrix, and the evaluation of each index with different experts tended to be various [13]. Therefore, AHP was combined with the Bayesian network and self-learning to obtain more unbiased and accurate results. The Bayesian maximum entropy self-learning weights were obtained when the entropy values were the maximum. Furthermore, the FW was calculated according to the hierarchical relationships of the Bayesian network (Table 2).

To further demonstrate the availability of our method, we compared the weights from two methods (AHP and AHP with maximum entropy self-learning). The results revealed that the weight orderings of recognized indexes of psoriasis (PASI and BSA) were same, which confirmed our method was reliable to an extent. Interestingly, the weight orderings of complement (C3 and C4), quality of life (DLQI, SAS, and SDS), and accompanying symptoms (PSQI, XQ, and CCS) ascended after maximum entropy self-learning, which was consistent with the latest research that psoriasis is a systemic disease [14, 7]. Unexpectedly, the weights of IL-17, IL-22, and IL-23 were decreased that were correlated to the interaction between the interleukin family, or the specificity of interleukins in psoriasis is actually not high.

3.2. Stability Analysis

Considering the influence of sample size, stability analysis was carried out. The weights and entropy gradually altered and finally tended to be stable in the course of self-learning (Figures 3 and 4(a)). It indicated the results are relatively stable when the sample size is 100. On this basis, we explored the relationship among sample size, weights, self-learning times, and entropy values. When sample size input was increased, the entropy altered. In addition, the error bars of the entropy and weights were shorter gradually, indicating that the larger the sample size is, the more stable the entropy and weights are (Figures 4(b) and 5). Besides, the self-learning times showed a linear increase (r = 0.987) with the growth of the sample size (Figure 4(c)), revealing that the computing cost was enhanced with increasing sample size. Next, we showed the relations between the three. The vertical error bars and horizontal error bars varied with the numbers of self-learning increasing. With the increase of sample size, the entropy was more stable, and the numbers of self-learning increased, which reminded the three-way interaction (Figure 4(d)). Above, the stability analysis demonstrated that entropy increases from starting levels and becomes gradually stable with the growth of sample size and self-learning times. Approximately 50 patients were needed to make the model stable.

3.3. Results of Efficacy Evaluation

We included 100 patients for efficacy evaluation based on the results of stability analysis.

Table 3 shows means and standard deviations (Stds) of all patients with respect to each index (complete results are in the Supplementary Table S5). The SD reflected the similarities between the patient and the ideal patient. When the value of SD was closer to 1, it indicated that a single index of the patient being tested was more similar to that of an ideal patient.

Aggregating SD of each index to obtain ED: the ED from small to large corresponded to the efficacy from superior to inferior (Table 4). The results showed that the patients with preferable curative effect were L1 (0.3196), L18 (0.3713), L17 (0.3722), L37 (0.3906), and L21 (0.4126), whereas those with unfavorable effect were L100 (0.7904), L90 (0.7749), L61 (0.7656), L82 (0.7469), and L99 (0.7400). Interestingly, ESPA was feasible for efficacy evaluation, overcoming the major limitation of uncertain grades and thresholds.

3.4. Sensitivity Analysis

As the inputs of some indexes of efficacy evaluation are uncertain, efficacy evaluation results were affected by their uncertainties. Thus, a sensitivity analysis was performed to check the consistency of the obtained efficacy evaluation ranking. If the change of a patient’s data affected the thresholds of the indexes, the ED of other patients changed accordingly. Sensitivity analysis was performed on 100 patients, where 20 patients are randomly shown in Figure 6. The other patients’ ED had no noticeable changes when a patient’s indicators were changed. This result illustrated that the ESPA efficacy evaluation model is relatively stable.

4. Discussion

Psoriasis is an immune-mediated chronic inflammatory skin disease with a high incidence. Till now, the question of how to comprehensively evaluate the treatment for psoriasis is a major clinical issue [32]. In this study, we first constructed a novel evaluation system of psoriasis by adopting Bayesian maximum entropy self-learning and ESPA.

The most obvious finding to emerge from the analysis was that PASI and BSA, as the most common evaluation tools of psoriasis, were the core indicators with the highest weights either before or after self-learning, which is consistent with expert knowledge and clinical experience. Guidelines of European Association of Dermatology and venereal Diseases recommended that PASI score and BSA should be the first choice of an objective indicator [33]. Even though the United States guidelines indicated that PASI score is cumbersome, it is still employed as a criterion for assessing the severity of psoriasis patients [34]. One previous study showed that compared with Patients Global Assessment, BSA and PASI have higher weights, which is consistent with our result [35].

One interesting finding was that the weights of C3 and C4 increased significantly in the self-learning system. Previous studies showed that the levels of C3 and C4 in patients with psoriasis were significantly higher, and C3 could be considered as a reliable marker of cardiac metabolic risk in psoriasis [36, 37]. Animal experiments have demonstrated that psoriasiform dermatitis was significantly alleviated and a marked reduction in C3 has been observed throughout when the S100A9 gene was deleted, in an imiquimod-induced mouse psoriasis model [38]. These studies indicate that there is a correlation between psoriatic lesions and C3 level, and the complement factor should be monitored as a considerable indicator of evaluation in various therapeutic interventions. Combining the AHP and maximum entropy criterion, we acquired more subjective and accurate results than the AHP on the weights of C3 and C4. However, C3 and C4 have not been integrated into the existing system, and the current study will fill this gap.

Besides, accompanying symptoms and quality of life received higher weights after self-learning. Previous research suggested that DLQI should be recognized as a major indicator and has even assisted in deciding patient-specific treatment strategies [27, 39]. In addition, the current study showed that the weight of SAS is higher than that of the AHP. One study has revealed that the psoriasis patients are more prone to anxiety than normal people, even though there is no obvious correlation between anxiety and the severity of psoriasis lesion, suggesting that doctors should not ignore the anxiety level of patients with mild psoriasis [40].

What is curious about this result is that the weight of interleukin was significantly lower than that by the AHP, and there were few differences among IL-10, IL-17, IL-22, and IL-23. In contrast, these indicators were considered as essential indicators in one of our earlier studies [26]. On the one hand, the reason may be related to the interactions among these indicators, and the evidence of the complex interplay has been proven. Previous studies reported that IL-23 stimulates the process of IL-17 secretion by Th17 cells, while the secretion of IL-22 needs to be encouraged by Th17 [41]. Moreover, enhancing IL-22 could lead to inhibit the production of IL-17 or stimulate the production of IL-10 [42], indicating that there may be a mechanism of interaction among IL-17, IL-23, and IL-22, which needs to be further unraveled. Another possible explanation for this is that it likely relates to specificity of interleukin. The occurrence and progression of psoriasis is the result of the combined force of a series of inflammatory cytokines and signaling pathways, including IL-6, IL-17, IL-23, IL-27, TNF-α, and INF-γ, but only IL-17 is the core mediator directly involved in the inflammation and disease progression of psoriasis. However, it has been proved that the IL-17 family has been found to play a causative role in tumor progression and autoimmune disease [43]. Therefore, this phenomenon may be correlated with the complex downstream signaling pathways resulting in less specificity.

In order to guarantee the accuracy of unsupervised self-learning, we verified the stability and sensitivity of the system. The large sample sizes increased statistical power that makes the results reliable. With the increase of the sample size, the entropy values were gradually stable, and the numbers of self-learning were also gradually increasing. Moreover, we conducted sensitivity analysis arbitrarily changing the clinical data of patients one by one and found that EPSA results are almost unaffected by the sample size, which indicated that our results were stable.

Compared with the AHP, the present system contains the following characteristic and advantages. The main strength of this study is that we applied a Bayesian network to construct a comprehensive measurement of treatment effectiveness among a variety of symptom indicators. Secondly, the judgment matrix needs to repeatedly modify in the weight’s calculation process by the AHP, and the scores assessed by a variety of experts have considerable variation by subjective factors. Therefore, on the basis of the known weights and self-learning method, the optimal weights were output according to the maximum entropy criterion, which has higher objectivity [44]. In addition, SPA or the set pair cloud model all relied on hierarchical classification so that the symptom grades and thresholds were needed during calculation of the CD. In the study, ESPA was an improvement over the SPA, which is applicable at uncertain grades and thresholds.

Despite these promising results, limitations remain. Firstly, the current study only included patients with psoriasis vulgaris, so that one should be cautious when generalizing to other types of psoriasis and the indicators may need to be adjusted. Secondly, a larger sample for model testing may be needed before practical application. Thirdly, considering the large numbers of indicators for psoriasis, other indicators could be included to evaluate the efficacy of psoriasis.

5. Conclusions

Given diversified indicators of psoriasis, we proposed a system that considers multiple indicators and uses more objective weights to evaluate the efficacy, while overcoming the application limitations of SPA. In addition, 100 patients were included to confirm the stability and effectiveness of the system through stability analysis and sensitivity analysis. Theoretically, the system can be applied to other diseases, which requires further research for clarification.

Abbreviations

AHP:Analytic hierarchy process
ESPA:Extended set pair analysis
TCM:Traditional Chinese medicine
DAG:Directed acyclic graph
CPTs:Probability tables
IW:Initial weight
FW:Final weight
CD:Connection degree
SD:Similarity degree
SPA:Set pair analysis
ED:Euclidean distance
PASI:Psoriasis area and severity index
BSA:Body surface area
SCC-Ag:Squamous cell carcinoma antigen
TNF-a:Tumor necrosis factor-α
IL:Interleukin
DLQI:Dermatology Life Quality Index
SAS:Self-Rating Anxiety Scale
SDS:Self-Rating Depression Scale
XQ:Xerostomia Questionnaire
CCS:Cleveland Clinic Score
PSQI:Pittsburgh Sleep Quality Index
C3/4:Complement 3/4
Stds:Standard deviations.

Data Availability

All of the data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

The data from preclinical studies [17] have been approved by the Ethics Committee of the Yueyang Hospital of Integrated Traditional Chinese and Western Medicine, protocol 2019-028.

All clinical study participants provided written informed consent before the study.

Disclosure

Le Kuai, Xiao-ya Fei, and Jing-si Jiang are the co-first authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Dr and Bin Li had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. Le Kuai and Xiao-ya Fei conceptualized the study;Xin Li curated data; Le Kuai and Xin Li conducted formal analysis; Bin Li, Le Kuai, and Xin Li acquired funding; Jing-si Jiang, Ying Zhang, Ying Luo, and Yi Ru conducted investigation; Jing-si Jiang and Yi Ru formulated the methodology; Wei Li and Jian-kun Song were involved in project administration; Wei Li, Jian-kun Song, and Shuang-yi Yin obtained resources; Xiao-ya Fei and Shuang-yi Yin were responsible for software; Jing-si Jiang, Shuang-yi Yin, Jian-kun Song, and Bin Li supervised the work; Xin Li and Ying Luo performed validation; Shuang-yi Yin performed visualization; Xiao-ya Fei, Jing-si Jiang, and Ying Zhang prepared the original draft; Le Kuai, Xiao-ya Fei, Jing-si Jiang, Xin Li, Ying Zhang, Ying Luo, Yi Ru, Jian-kun Song, Wei Li, Shuang-yi Yin, and Bin Li reviewed and edited the manuscript.

Acknowledgments

This work was supported by a grant from the National Key Research and Development Program of China (No. 2018YFC1705301); the National Natural Science Foundation of China (Nos. 81671959, 81874470, 81973860, 81904214, 82004235, and 82074427); the Xinglin Young Scholar, Shanghai University of Traditional Chinese Medicine (No. RY411.33.10); the School-level Postgraduate Innovation Training Program, Shanghai University of Traditional Chinese Medicine (No. JY611.02.03.83); the Joint Funds of the NSFC and Henan Province (No. U1704171); the Shanghai Pujiang Talent Program (No. 2020PJD067); the Shanghai Development Office of TCM (No. ZY(2018-2020)-FWTX-4010); and Dermatology Department of Traditional Chinese Medicine, Clinical Key Specialty Construction Project of Shanghai (No. shslczdzk05001). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Supplementary Materials

Supplementary Table S1: experts’ judgment on importance of each index. Supplementary Table S2: the data of each index of 100 patients. Supplementary Table S3: expert scores of the indexes. Supplementary Table S4: results of the weights calculated by the AHP and the consistency test. Supplementary Table S5: the SD of 100 patients. (Supplementary Materials)