 Research Article
 Open Access
 Open Peer Review
 Published:
Rough set theory based prognostic classification models for hospice referral
BMC Medical Informatics and Decision Makingvolume 15, Article number: 98 (2015)
Abstract
Background
This paper explores and evaluates the application of classical and dominancebased rough set theory (RST) for the development of datadriven prognostic classification models for hospice referral. In this work, rough set based models are compared with other datadriven methods with respect to two factors related to clinical credibility: accuracy and accessibility. Accessibility refers to the ability of the model to provide traceable, interpretable results and use data that is relevant and simple to collect.
Methods
We utilize retrospective data from 9,103 terminally ill patients to demonstrate the design and implementation RST based models to identify potential hospice candidates. The classical rough set approach (CRSA) provides methods for knowledge acquisition, founded on the relational indiscernibility of objects in a decision table, to describe required conditions for membership in a concept class. On the other hand, the dominancebased rough set approach (DRSA) analyzes information based on the monotonic relationships between condition attributes values and their assignment to the decision class. CRSA decision rules for sixmonth patient survival classification were induced using the MODLEM algorithm. Dominancebased decision rules were extracted using the VCDomLEM rule induction algorithm.
Results
The RSTbased classifiers are compared with other predictive and rule based decision modeling techniques, namely logistic regression, support vector machines, random forests and C4.5. The RSTbased classifiers demonstrate average AUC of 69.74 % with MODLEM and 71.73 % with VCDomLEM, while the compared methods achieve average AUC of 74.21 % for logistic regression, 73.52 % for support vector machines, 74.59 % for random forests, and 70.88 % for C4.5.
Conclusions
This paper contributes to the growing body of research in RSTbased prognostic models. RST and its extensions posses features that enhance the accessibility of clinical decision support models. While the nonrulebased methods—logistic regression, support vector machines and random forests—were found to achieve higher AUC, the performance differential may be outweighed by the benefits of the rulebased methods, particularly in the case of VCDomLEM. Developing prognostic models for hospice referrals is a challenging problem resulting in substandard performance for all of the evaluated classification methods.
Background
Hospice care reduces the emotional burden of illness on terminal patients by optimizing pain relief strategies [1] and provides a demonstrated, costeffective increase in the quality of endoflife care when compared to conventional programs [2]. This increase in quality of care elevates the quality of life of both patients and their families [3].
The advantages of hospice care are diminished for terminally ill patients who enter either prematurely or too late. In general, premature hospice referral represents a lost opportunity for the patient to receive potentially effective and lifeprolonging treatment. Conversely, late hospice referral is not desirable and negatively impacts both the quality of endoflife care and the quality of life of patients and their families [4, 5]. According to Medicare regulations, patient eligibility for hospice care is contingent upon a life expectancy of less than six months, as estimated by the attending physician and certified by the medical director of the hospice program [6]. Medicare claims data report that 14.9 % of hospice care patients lived for more than 180 days after enrollment, while 28.5 % were late referrals who died within 14 days [4, 6]. Accurate prognostication of life expectancy is crucial in endoflife care decisions and is consequently of vital importance for patients, their physicians and their families.
Prognostic models are an important instrument in prognostication as, in conjunction with direct physician observation, they increase the accuracy of prognostication when compared to physician observation alone [7]. However, a significant barrier to the widespread practical use of prognostic models is their perceived lack of clinical credibility [8].
The objective of this work is to explore and evaluate the application of rough set approaches in the development of datadriven prognostic models with respect to two criteria essential to clinical credibility: accuracy and accessibility. To this end, we will explore Rough Set Theory (RST) as it is applied to endoflife care and hospice referral decision support models. Additionally, we will compare the results of the RSTbased models with several widely known methods: logistic regression, support vector machines (SVM), C4.5, and random forests (RF).
This paper is organized as follows: The Motivation Section presents important features of clinically credible prognostic models and other characteristics of clinical data sets that motivate the use of RST. We then present an overview of the fundamental theory of rough sets for analyzing datasets (section Methods), followed with a similar overview of the theory of the Dominancebased Rough Set Approach (DRSA). In this section, we also discuss the use of decision rules in conjunction with the rough set approaches. The section Dataset description describes the dataset used for the demonstration of the proposed prognostic models. Section Experimental design presents the development of the prognostic models, followed by an overview of the performance evaluation methods used in this study. Finally, we report the results and conclusions, and discuss limitations and future directions of our work.
Motivation
The objective of a prognostic model is to determine relationships between covariates and a healthrelated outcome. In the case of life expectancy estimation, prognostic models improve the accuracy in critical clinical decisions and are shown to be superior to physicians’ prognostication alone [9]. Models for estimating the life expectancy of terminally ill patients include the use of statistical and probabilistic methods [10–18], artificial intelligence techniques such as neural networks and support vector machines (SVM) [19–21], decision trees [22, 23] and rough set methods [24, 25]. Survival models [6, 12, 14, 16, 18, 22, 23] focus on estimating the probability that a patient will survive a finite period of time. Classification models, based on methods such as neural networks, SVM and logistic regression [17, 19–21, 26], represent the survival outcome as a binary variable, predicting the status of a patient at a critical point in time (e.g. six months) by classifying the patient as surviving or not surviving the critical time frame. Classification models require the use of noncensored data where survival outcome is known for every patient in the dataset at the critical decision point in time.
A recent review [15] demonstrated that, despite the importance of accurate prognostication within the spectrum of medical care objectives, there is a lack of accessible and accurate prognostic models available to physicians in practice. To withstand clinical trials, and to meet the needs of physicians and patients, a prognostic model must have clinical credibility, meaning that the model must posses a high level of accuracy and accessibility for physicians to believe in the value of the model as a prognostic tool. That is, in addition to accurate prognostication, such a model should be traceable in its structure, meaning the “model’s structure should be apparent and its predictions should make sense to the doctors who will rely on them” [8]. Likewise, the model should provide interpretable results that facilitate explanation of the prognosis, the data required for the model must be relevant and simple to collect with high reliability, and physicians must be able to apply the modeling method correctly without violating the fundamental assumptions of the model.
Clinical datasets present unique challenges that must also be addressed when building datadriven prognostic models. Cios and Moore [27] argue that there are a number of features specific to medical data that result from the volume, heterogeneity and complexity of data that lack canonical form. Additionally, ethical, legal and societal concerns greatly affect the framework under which medical data may be used. The current US model encourages the use of deidentified, minimal risk medical data for research purposes, specifically data collected during routine treatment of patients. It is common for medical data collected in such a way to contain redundant, insignificant, incomplete or inconsistent data objects. Furthermore, the underlying conceptual structures of medicine are not easily formalized mathematically, as the medical field lacks the necessary constraints for the mathematical characterizations common to the physical sciences. As a result, many medical concepts are vaguely defined [28].
Rough Set Theory [29] is a mathematical tool for data analysis that has been used to address vagueness and inconsistencies present in datasets [30]. RST provides a systematic approach for analyzing data without implicit assumptions about relationships between covariates, an advantage that makes RST suitable for integration into medical applications [31]. The information extracted from the dataset by RST and its related methods can be represented in the form of “if–then” decision rules—an intuitive representation that offers significant advantage over “black box” modeling approaches [32] and that increases accessibility and thus clinical credibility.
In the medical field, applications of RST focus mainly on the diagnosis and prognostication of diseases, where it has been demonstrated that RST is useful for extracting medical prognostic rules from minimal information. Tsumoto [33] argues that the concepts of approximation established in RST reflect the characteristics of medical reasoning, explaining why RST performs well in the medical field. For example, RST can be used to highlight nonessential prognostic factors in a particular diagnosis, thus helping to avoid redundant, superfluous or costly tests [34–38]. Recently, methods that combine survival analysis techniques and RST have been used to generate prognostic rules that estimate the survival time of a patient [24, 25].
Methods
Classical rough set approach (CRSA)
Rough Set Theory, introduced by Pawlak in [29], provides methods for generalizing or reducing information so as to facilitate knowledge discovery by exploiting the relational indiscernibility of objects in an information table. Central to RST is the notion that an observed object has a certain amount of information associated with it. When considered in relation to a cohort of observed objects, this information is used to group similar objects into information granules. Together, the information provided by the set of observed objects can be generalized to describe the conditions required for membership in a concept class.
Notation
The methods of classical RST, hereafter referred to as the CRSA, act upon an information table of the form S=(U,A,V,f), where U is a nonempty finite set of objects, called the universe. A=C∪{d} is a set of attributes that describe a given object in U, comprised of a set C of condition attributes and an optional decision attribute d. When d is present, the information table is a decision table. The set of all values, V, contains the value sets V_{ a }, for every attribute a∈A. Given an object x∈U, f:U×A→V maps the condition attribute of object x to its associated value v=f(x,a)∈V_{ a }. A value attribute pair (a,v) for a given object is referred to as a descriptor.
Table 1 provides an example of a discretized decision table, where six prognostic factors, as the condition attributes, describe seven patients. The decision attribute, presence of coronary disease in the patient, is represented by the binary attribute d→{Yes,No}.
The objects in a decision table can be grouped according to their descriptors. For example, patients x_{5} and x_{6} have the same attribute values and are thus indiscernible from each other. In general, two objects x_{ i },x_{ j }∈U are indiscernible with respect to a set of condition attributes B⊆C if f(x_{ i },a)=f(x_{ j },a) ∀a∈B. This relation is called an indiscernibility relation, defined as R(B)={(x_{ i },x_{ j })∈U:∀a∈B,f(x_{ i },a)=f(x_{ j },a)}.
For example, the patients in Table 1 can be separated into four groups according to the indiscernibility relation R(C):X_{1}={x_{1}},X_{2}={x_{2}},X_{3}={x_{3},x_{4},x_{7}},X_{4}={x_{5},x_{6}}. These groups of objects are referred to as equivalence classes, or conditional classes for B⊆C. An equivalence class for the decision attribute is called a decision class or concept, and in this example there are two groups: Y_{ N o }={x_{1},x_{2},x_{3}} and Y_{ Y e s }={x_{4},x_{5},x_{6},x_{7}}. The equivalence class specified by the object x_{ i } with respect to R(B) is denoted as [ x_{ i }]_{ B }.
Set approximations
The goal of the CRSA is to provide a definition of a concept according to the values of the attributes of the equivalence classes that contain objects that are known instantiations of the concept. As such, in a consistent decision table, membership in a conditional class implies membership in a particular decision class. In Table 1, x∈X_{4} implies x∈Y_{ Y e s }. Membership in X_{3}, however, does not imply Y_{ Y e s } as x_{4},x_{7}∈Y_{ Y e s } but x_{3}∈Y_{ N o }. Thus Table 1 is inconsistent as f(x_{4},d)≠f(x_{3},d) and f(x_{7},d)≠f(x_{3},d).
To represent an inconsistent decision table, the CRSA establishes an upper and lower approximation for each decision class, Y. The lower approximation is comprised of all objects that definitely belong to Y, while the upper approximation includes all objects that possibly belong to Y. It can be said that an object x_{ i } definitely belongs to a concept Y if [ x_{ i }]_{ C }⊆Y and that x_{ i } possibly belongs to a concept Y if [ x_{ i }]_{ C }∩Y≠∅. Thus, the lower and upper approximations are defined as follows:
The boundary region, BND_{ B }(Y), contains those objects that possibly, but not certainly, belong to Y. Conversely, the set \(U  \overline {R}_{B} (Y)\) is the outside region containing those objects that certainly do not belong to Y. In our example, the lower and upper approximations for Y_{ Y e s } are \(\underline {R}_{C} (Y_{\mathit {Yes}}) = X_{4} = \{ x_{5}, x_{6} \}\) and \(\overline {R}_{C} (Y_{\mathit {Yes}}) = X_{4} \cup X_{3} = \{ x_{3}, x_{4}, x_{5}, x_{6}, x_{7} \}\), and the boundary region contains the objects BND_{ B }(Y_{ Y e s })={x_{3},x_{4},x_{7}}.
Let F={Y_{1},Y_{2},…,Y_{ n }} represent a classification, i.e. a set of decision classes. The quality of approximation of classification, γ_{ B }(F), with respect to attributes B, expresses the ratio of all objects covered by the lower approximation R̲_{ B }(F) = {R̲_{ B }(Y_{1}),R̲_{ B }(Y_{2}),…,R̲_{ B }(Y_{ n })} over all objects in U. The quality of approximation is expressed as:
Dominancebased rough set approach (DRSA)
Under the DRSA [39] the relations between objects are no longer made by the indiscernibility relation as described in the CRSA [29]. In its place, the DRSA introduces a new dominance relation that allows for ordinal attributes with preferenceordered domains wherein a monotonic relationship exists between the attribute and the decision classes. An example of such a relationship occurs when a “better” or “worse” value of an attribute leads to a “better” or “worse” decision class.
Notation
A decision table in the DRSA is expressed in the same way as the CRSA. To differentiate between attributes with and without a preferenceordered domain, those with a preference order are called criteria while those without are referred to as attributes, as in the CRSA.
In the DRSA the domain of criteria a∈A is completely preordered by the outranking relation ≽_{ a }, representing the preference order of the domain. The outranking relation is also applicable for comparing two objects such that for x_{ i },x_{ j }∈U, x_{ i }≽_{ a }x_{ j } means that x_{ i } is at least as good as (outranks) x_{ j } with respect to the criterion a∈A.
Commonly, the domain of a criteria a is a subset of real numbers, V_{ a }⊆R and the outranking relation is then a simple order “ ≥” on real numbers such that the following relation holds: x_{ i }≽_{ a }x_{ j }⇔f(x_{ i },a)≥f(x_{ j },a). This relation is straightforward for gaintype criteria (the more, the better), and can be easily reversed for costtype criteria (the less, the better).
Using Table 1 as an example, the decision criterion d is preferenceordered such that a positive diagnosis of coronary disease is assumed to be the “preferred” decision class. Criterion preference relations are then organized in the direction of the decision class; values which generally contribute to the incidence of coronary disease are preferred over those which indicate lower risk, much in the same way that a positive diagnosis indicates presence of coronary disease. For the criteria in Table 1, higher values are preferred to lower values—as in the case of Age, SystBP, and HDL—and “Yes” is preferred to “No”—as in the case of Smoker and Diabetic. No such preference relation exists for Gender; as such, it is considered an attribute.
Let T={1,…,n} represent increasing indexes corresponding to the order of preferences of the decision criterion d. Then, the decision table is partitioned into n classes Y_{ t }, t∈T, where each object x∈U is assigned to one and only one class Y_{ t }. The decision classes are preferenceordered according to the decision maker, i.e. for all r,s∈T such that for r>s the objects from class Y_{ r } are strictly preferred to the objects from class Y_{ s }.
For our example in Table 1, Y_{1}={x_{1},x_{2},x_{3}} corresponds to patients without a coronary disease and Y_{2}={x_{4},x_{5},x_{6},x_{7}} corresponds to the patients with a coronary disease. Therefore, each patient in Y_{2} is preferred to each patient in Y_{1}.
Set approximations
In the DRSA, the approximated sets are upwards and downwards unions of decision classes rather than individual decision classes as in the CRSA. Upward and downward unions of classes are defined as:
For any pair of objects (x_{ i },x_{ j })∈U, x_{ i } dominates x_{ j } with respect to a set of condition attributes P⊆C, denoted by x_{ i } D_{ P } x_{ j }, if the following conditions are satisfied simultaneously:
The dominance relation defines two sets called dominance cones, where for each x_{ i }∈U:
Considering the dominance cones, the lower and upper approximations of the union of decision classes are defined as follows. The lower approximation \(\underline {R}_{P} \left (Y_{t}^{\geq }\right)\) represents objects that certainly belong to \(Y_{t}^{\geq }\), such that there is no other object that dominates x and belongs to a decision class inferior to Y_{ t }. Similarly, the lower approximation \(\underline {R}_{P} \left (Y_{t}^{\leq }\right)\) represents objects that certainly belong to \(Y_{t}^{\leq }\), with no other object dominated by x and belonging to a decision class superior to Y_{ t }. The upper approximations represent objects that possibly belong to one of the upward or downward unions of decision classes.
Similar to the CRSA, the boundary regions are defined as:
Using our example decision table, Table 1, and considering the full set of condition attributes, it can be seen that x_{4} D_{ C } x_{3}, and furthermore \(D_{C}^{+} (x_{4}) = \{ x_{3}, x_{4}, x_{7} \}\), \(D_{C}^{} (x_{4}) = \{ x_{3}, x_{4}, x_{7} \}\). Considering the dominance cones for all patients, the lower and upper approximations of the union of decision classes are \(\underline {R}_{C} \left (Y_{2}^{\geq } \right) = \{ x_{5}, x_{6} \}\), \(\overline {R}_{C} \left (Y_{2}^{\geq }\right) = \{ x_{3}, x_{4}, x_{5}, x_{6}, x_{7} \}\), \(\underline {R}_{C} \left (Y_{1}^{\leq }\right) = \{ x_{1}, x_{2} \}\), \(\overline {R}_{C} \left (Y_{1}^{\leq } \right) = \{ x_{1}, x_{2}, x_{3}, x_{4}, x_{7} \}\) and the boundary regions are \(BND_{C} Y_{2}^{\geq } = BND_{C} Y_{1}^{\leq } = \{ x_{3}, x_{4}, x_{7} \}\).
For every subset of attributes P⊆C, the quality of approximation of the decision classes Y with respect to the attributes P, γ_{ P }(Y), is defined as the proportion among all objects in U of objects consistently defined with respect to the attributes P and the decision classes Y.
The variable consistency DRSA
The variable consistency DRSA (VCDRSA) allows the decision maker to relax the strictness of the dominance relation, thus accepting a limited number of inconsistent objects in the lower approximation, according to an object consistency level threshold, l∈(0,1]. In practice, by selecting this consistency level l, a patient x∈U becomes a member of the lower approximation of a given upward union if at least l∗100 % of the patients dominating x also belong to that decision class. By allowing inconsistencies, the VCDRSA avoids over fitting the training set and thus may be more effective in classifying new cases.
The lower approximations of the VCDRSAbased model are represented as follows:
Continuing with the example described in Table 1, setting l=0.6 moves the objects x_{4} and x_{7}, previously included in the upper approximation \(\overline {R}_{C} \left (Y_{2}^{\geq }\right)\), to the lower approximation of class \(Y_{2}^{\geq }\), i.e: \(\underline {R}_{C}^{0.6} \left (Y_{2}^{\geq } \right) = \left \{ x_{4}, x_{5}, x_{6}, x_{7} \right \}\). This follows from \(\frac { D_{C}^{+} (x_{i}) \cap Y_{t}^{\geq } }{  D_{C}^{+} (x_{i}) } = \frac {2}{3} \geq l\), for i=4,5,6,7.
Decision rules
There are a number of methods available for induction of decision rules from the lower or upper approximations of the decision classes [40–42] or from reducts extracted from the decision table [43]. Decision rules in this study were obtained using the MODLEM [40, 41] and VCDomLEM [42] algorithms for the induction of classical and dominancebased rough set rules, respectively. In both cases, decision rules are induced from approximations of decision classes. Both the MODLEM and VCDomLEM algorithms generate a minimal set of decision rules using a minimal number of rule conditions, thus the inclusion of MODLEM allows for an evaluation of the impact of accounting for the preference order information in the VCDRSA. Once decision rules have been induced, the collection of these rules can then be used to classify unseen objects—in the case of our example table, a new patient who may have cardiac disease.
A decision rule has the form if A then B, or A→B, where A is called the antecedent and B the consequent of the rule. The antecedent is a logical conjunction of descriptors and the consequent is the decision class or union of decision classes suggested by the rule.
Formally, in the CRSA, decision rules are generated from the lower or upper approximations. For example, for an approximation containing objects with descriptors r with respect to a set of condition attributes, B_{ r }⊆C, a decision rule is expressed as
where a_{ i }∈B_{ r } is an attribute found in the attribute set B_{ r }, and \(r_{a_{i}} \in V_{a_{i}}\) and Y_{ t } are the attribute values and a decision class, respectively, of the objects in the rulegenerating approximation. From our example in Table 1, a decision rule induced with MODLEM from the lower approximation \(\underline {R}_{\{Age, Smoker\}} (Y_{\mathit {Yes}}) = \{ x_{5}, x_{6} \}\) would be: if Age=H and Smoker=Yes then Coronary Disease=Yes.
In the DRSA, decision rules are induced from the lower approximations and the boundaries of the union of decision classes. From the lower approximations, two types of decision rules are considered. Decision rules generated from the Plower approximation of the upward union of decision classes \(Y_{t}^{\geq }\) are described by
where b_{ i }∈P are criteria, a_{ j }∈P are attributes, \(r_{b_{i}} \in V_{b_{i}}\) and \(r_{a_{j}} \in V_{a_{j}}\). From the example in Table 1, the Plower approximation of the upward union of the decision class, \(Y_{2}^{\geq } = \underline {R_{C}} Y_{2} = \{x_{5}, x_{6}\}\), leads to the following decision rule:

If Gender = M and Age ≥ H and HDL ≥ L and Diabetic ≥ Yes and Smoker = Yes, then Coronary Disease = Yes.
Decision rules generated from the Plower approximation of the downward union of classes \(Y_{t}^{\leq }\) are described by
where b_{ i }∈P are criteria, a_{ j }∈P are attributes, \(r_{b_{i}} \in V_{b_{i}}\) and \(r_{a_{j}} \in V_{a_{j}}\). From the example in Table 1, the Plower approximation of the downward union of classes \(Y_{1}^{\leq } = \underline {R_{C}} Y_{1} = \{x_{1}, x_{2}\}\), leads to the following decision rules:

If Gender = F and Age ≤ H and SystBP ≤ M and HDL ≤ L and Diabetic ≤ No and Smoker ≤ Yes, then Coronary Disease = No

If Gender = M and Age ≤ H and SystBP ≤ M and HDL ≤ L and Diabetic ≤ No and Smoker ≤ Yes, then Coronary Disease = No
The boundaries \(BND_{P} Y_{t}^{\geq }\) and \(BND_{P} Y_{t}^{\leq }\) generate the following rules
where b_{ i },b_{ j }∈P are criteria, a_{ k }∈P are attributes, \(r_{b_{i}} \in V_{b_{i}}\), \(r_{b_{j}} \in V_{b_{j}}\) and \(r_{a_{k}} \in V_{a_{k}}\) (note i and j are not necessarily different). From the example in Table 1, the boundary decision classes \(BND Y_{2}^{\geq } = BND Y_{1}^{\leq } = \{x_{3}, x_{4}, x_{7}\}\), leads to the following decision rule:

If Age ≥ M and SystBP ≥ M and HDL ≥ H and Diabetic ≥ No and Smoker ≥ No and Age ≤ M and SystBP ≤ M and HDL ≤ H and Diabetic ≤ No and Smoker ≤ No and Gender = F, then Coronary Disease = (No, Yes)
The MODLEM and the VCDomLEM algorithms utilize a heuristic strategy called sequential covering [44] to iteratively construct a minimal set of minimal decision rules. The sequential covering strategy successively constructs a set of decision rules for each upward and downward union of decision classes in a training set by selecting, at each iteration, the “best” decision rule, after which the training objects described by the rule conditions are removed. Subsequent iterations again select the best decision rule and remove the covered objects until reaching a stopping criteria or until all of the objects in the unions of decision classes are described by a rule in the rule set.
To ensure minimality, antecedent descriptors, called elementary conditions, of each rule are checked at each iteration and redundant elementary conditions are removed. Similarly, redundant rules are removed from the final rule set.
In both algorithms, decision rules are grown by consecutively adding the best available elementary condition to the rule. CRSA elementary conditions are evaluated in the MODLEM algorithm in terms of either the class entropy measure [45] or Laplacian accuracy [46]; the former was used in this study. MODLEM does not restrict elementary conditions to those attributes not currently in the rule; as such, multiple elementary conditions may contain the same attribute. Therefore, a decision rule induced by MODLEM may contain antecedents in which attribute values are described as belonging to a range or a set of values or as being greater or less than a particular value.
Dominancebased elementary conditions are evaluated according to a rule consistency measure. VCDomLEM provides three such measures; the rule consistency measure used in this study is μ, as described in [47]. For the sake of clarity, Y_{ t } shall be used to represent an individual decision class in the CRSA or alternatively an upward or downward union of decision classes, \(Y^{\geq }_{t}\) or \(Y^{\leq }_{t}\), with respect to the DRSA. The consistency, μ, of a proposed rule, \(r_{Y_{t}}\), suggesting assignment to Y_{ t } is defined as
Here \(\left [\Phi (r_{Y_{t}})\right ]\) indicates the set of objects described by the elementary conditions in \(r_{Y_{t}}\). The elementary condition, ec, that is selected for inclusion is that which leads to the highest rule consistency measure \(\mu (r_{Y_{t}} \cup ec)\) when combined with the current set of elementary conditions in the proposed rule. In the event of a tie, the elementary condition providing greatest coverage of the new rule is selected, by \(\left  \left [ \Phi (r_{Y_{t}} \cup ec) \right ] \cap \underbar {Y}_{t} \right \). The rule consistency measure, μ, was also implemented in MODLEM to relax consistency requirements and to allow more general rules to be induced. For further details on the MODLEM and VCDomLEM algorithms, the reader is referred to [40–42, 47].
To classify an unseen object, a standard voting process [43] is used to allow all rules to participate in the decision process, arriving at a patient classification by majority vote. Each rule is characterized by two support metrics. The left hand side (LHS) support is the number of patients in the table whose attributes match the antecedent, i.e: [Φ(r)], while the right hand side (RHS) support indicates the number of patients matching both the antecedent and the consequent of the rule, i.e: [Φ(r)]∩Y_{ t }. For a new, unseen patient, any rule whose antecedent descriptors match the patient descriptors “fires” by contributing as votes the RHS support for each decision class. For example, drawing up the example Table 1, the decision rule If Age = H and Smoker = Yes, then Coronary Disease = Yes has LHS = 2 since its antecedent matches patient x_{5} and x_{6} and RHS = 2 since its antecedent and consequent match the same patients. A new patient matching the antecedent of this rule will receive two votes for decision class Yes and zero votes for decision class No.
Once all rules have “voted”, the number of votes for each decision class is normalized against the total number of LHS support for all fired rules. The resultant ratio of RHS to LHS support is considered a frequencybased estimate of the probability that the patient belongs to the given decision class.
A final classification is therefore determined according to a threshold value, τ∈[ 0,1]. A patient is classified as not surviving six months if the estimated probability of death in six months is greater than τ. In the event of an estimated probability equal to τ, or in the absence of any fired rules (no rule matches the patient profile), classification is not possible and the patient is labeled undefined. For example, if the threshold value is set as 0.5 and the voting process yields an estimated probability of 70 %, then the patient is classified as not surviving the six month period.
Dataset description
SUPPORT dataset
The dataset used in this study is the SUPPORT (Study to Understand Prognoses and Preferences for Outcomes and Risks of Treatments) prognostic model dataset [48], a study of 9,105 terminally ill patients. SUPPORT enrolled patients, 18 years or older, who met specific criteria for one of nine serious illnesses, who survived more than 48 hours but were not discharged within 72 hours. Patients were followed such that survival and functional status were known for 180 days after entry. The result of the SUPPORT study is a prognostic model for 180day survival estimation of seriously ill hospitalized adults based on cubic splines and a Cox regression model. Given the inclusion criteria (described in full in Appendix 1 of [12]), the dataset is ideal for the present research in regards to clinical applicability, completeness of data, and comparability of results.
We consider as condition attributes the variables used in the SUPPORT prognostic model equation [12] to ensure consistency. The SUPPORT variables include ten physiologic variables in addition to the diagnosis groups, age, number of days in the hospital before entering the study, presence of cancer, and neurologic function as recorded in the SUPPORT data. Attribute names, descriptions and value ranges are listed in Table 2.
The median survival time for the patients in the study is 223 days.
Figure 1 shows the patients KaplanMeier survival curve with respect to number of days until death. The SUPPORT study inclusion criteria was designed to include patients with 50 % risk of death at 180 days; as seen in Table 2 death prior to 180 days was observed in approximately 47 % of patients.
General observations regarding the influence of condition attributes can be made by analyzing their relation in the proportion of patients surviving the six month period. For example, the KaplanMeier survival curve in Fig. 2 shows that a significant portion (75 %) of patients with coma or multiorgan system failure with malignancy (MOSF w/ malig) do not survive longer than 180 days, but patients with congestive heart failure (CHF) or chronic obstructive pulmonary disease (COPD) tend to live longer than 180 days.
Data preprocessing
In its published form, the SUPPORT dataset contains 9,105 cases. Missing physiological attribute values are filled in with a standard fillin value representing a normal physiological response, as provided by the SUPPORT authors in [48]. It is also worth noting that in the SUPPORT study, where neurologic functioning of the patient is recorded in the SUPPORT coma score (scoma), a patient for whom it was not possible to establish a Glasgow coma score was given a scoma value of zero. After missing data imputation, two cases have missing values in physiological attributes not addressed in the SUPPORT data set. The two incomplete cases were removed and the remaining 9,103 cases were considered in the development of the prognostic models.
Discretization
Discretization is the process by which appropriate categorical ranges are found for variables with a continuous value range. There are a number of methods available for unsupervised discretization that operate without input from the decision maker and are based only on the information available in the data table. In this work, however, discretization was primarily performed using the Acute Physiology and Chronic Health Evaluation (APACHE) III scoring system [11], a clinically accepted scoring system designed to estimate the risk of death in ICU patients. In this sense, the use of the APACHE III scoring system represents a researchvalidated, clinically appropriate, expert discretization scheme. This choice is founded on the proposition that expert discretization via APACHE III will result in medically and contextually relevant classification rules and data collection requirements, thus increasing the accessibility of the proposed prognostic model and ensures directly comparable rule sets for all evaluated rulebased methods.
APACHE III scores are designed to increase monotonically with respect to risk of death and thus provide the necessary preference relations for the DRSA. APACHE III scores for any given variable are close to zero for normal or only slightly abnormal values of that variable and increase according to increased severity of disease. For example, normal pulse rates of 50–99 bpm are given a score of 0, while elevated and lowered levels, 100–109 and 40–49 bpm respectively, are both given a score of 5. Thus, higher APACHE III scores are preferred to lower scores, as the higher scores indicate greater severity of disease and therefore greater risk of death within six months (considered the positive diagnosis). Discretization is not a requirement of any of the methods used in this study, however the APACHE III scores provide the monotonic preference relations for the DRSA and are used for the all of the rulebased methods.
For the rulebased methods considered in this study, the nine physiologic variables and the age variable were transformed to their representative APACHE III scores. The remaining physiologic variables not included in APACHE III—neurologic function, scoma, and blood gasses, pafi—were discretized using clinically accepted categorizations [49, 50]. The variable hday was discretized using the boolean reasoning algorithm [43]. Table 3 shows the categories defined in this process. Higher values of each of these variables are preferred to lower values.
Experimental design
This section provides details on the implementation and performance evaluation procedures for the comparison of the classification methods used in this study. The following two sections, describe the RSA and comparative methods respectively, the software used for their implementation and the selection of appropriate parameters for each of the methods. Finally, the methods for performance evaluation are discussed.
The general schema of the experimental design is as follows: after selecting appropriate parameters for each of the methods, 5fold cross validation was used to divide the data into training and testing sets. Methods with decision rule outputs were trained and tested on the discretized data set to demonstrate expected performance of a clinically credible rule set. Methods without decision rule outputs were trained on the raw, nondiscretized, data set. For these methods, designed to be applied to continuous variables, discretization does not improve clinical credibility and would likely hinder performance [51, 52].
Rough set rule induction and classification
MODLEM algorithm for CRSA decision rules
CRSA decision rules were obtained using the MODLEM algorithm as described in [40] and [41], implemented by the authors in the R programming language [53]. Decision rules were generated from the lower approximations with a rule consistency level μ≥m. The rule syntax follows the presentation in section Decision rules.
VCDomLEM algorithm for VCDRSA decision rules
Dominancebased rules were obtained using the VCDRSA as described in section The variable consistency DRSA and the VCDomLEM algorithm as implemented in jMAF [54]. VCDomLEM decision rules were generated from the lower approximation of each decision class, with an object consistency level threshold l. The syntax of the VCDRSA decision rules is as shown in section Decision rules. Only decision rules with rule consistency measure μ greater than the rule consistency threshold l are included in the classification model. Note that the rule consistency threshold and the object consistency threshold are equal and set at l.
Parameter selection
In order to select the most appropriate models for comparison, the performance of the rough set based models was evaluated for varying levels of rule consistency, m and l, for the CRSA and VCDRSA respectively. Classifier performance at a particular value of m or l is datasetdependent; however, in general, values close to one provide rule sets that are more conservative in describing the training set objects, while values closer to zero provide rule sets that are more general. Thus, to find the appropriate balance between strict, descriptive models that are prone to overfitting and overly general models that provide little useful information, the RSA models were evaluated at m,l=0.1,0.2,0.4,0.6,0.8,1.0.
Comparative methods
To evaluate the performance of the RSAbased prognostic models, logistic regression, SVM, and RF were applied to the nondiscretized SUPPORT dataset. To ensure directly comparable rule sets, C4.5 was applied to the discretized SUPPORT dataset. Each of these methodologies was applied using the software package Weka 3.6.9 [55], within which appropriate parameters were selected for SVM, C4.5 and RF using GridSearch with 10fold cross validation settings. Logistic regression was selected for its popularity in classification models using noncensored data and in clinical settings [18, 56].
Support vector machines, originally presented in [57], find separating boundaries between decision classes after input vectors are nonlinearly mapped into a high dimensional feature space. Support vector machines have been investigated in survival analysis applications [58] as they—similar to the RSAbased methods—automatically incorporate nonlinearities and do not make a priori assumptions about factor interactions. SVMbased models are known to perform well at classification tasks, however they do not provide clinicianinterpretable justification for their results [59]. Support vector machines were selected to evaluate whether the increased accessibility of the RSAbased methods involves a tradeoff in accuracy.
C4.5 is a well known algorithm for generating a decision tree using information entropy to select the best splitting criteria at each node [60]. A decision tree built by C4.5 can be expressed as a set of ifthen decision rules, thus providing a comparative decision rule based method. Decision trees were obtained using the Weka J48 implementation [60] of the C4.5 algorithm.
Random forests is a popular ensemble classification method based on decision trees [61]. The random forests algorithm builds an ensemble of decision trees, where each tree is built on bootstrap samples of training data with a randomly selected subset of factors.
Performance evaluation methods
The performance of the models was tested by measuring the discriminatory power of both the m and lconsistent decision rules sets when applied to the reserved testing data. For our notation, a classification of d.6months=Yes is referred to as a positive classification, and d.6months=No is negative. Sensitivity is defined as the fraction of patients who did not survive six months and are correctly classified by the model, or the fraction of true positive classifications of all test patients who did not survive six months. Conversely, specificity is defined as the fraction of patients who did survive six months and were correctly classified by the model, or the fraction of true negatives of all test patients who did survive six months.
The overall accuracy of the classification models is reported in terms of area under the receiver operating characteristic (ROC) curve, or AUC (area under the curve). The ROC curve graphs the sensitivity of the classifier, or the true positive rate, versus 1 − specificity, the false positive rate, as the threshold probability, τ, for positive classification is varied from 0 to 1. The best overall classification performance is realized when AUC is equal to 1, while an AUC of 0.5 indicates a classifier performance no better than random selection. Best separation between decision classes is realized at the threshold corresponding to the point on the ROC curve closest to the point (0,1).
In order to select the most appropriate MODLEM and VCDomLEMbased models for comparison, two performance issues related to the generated rule set were considered: coverage and AUC of the model. The coverage of the classification model is defined as the percentage of testing set patients for whom a classification is possible. Additionally, to evaluate the number of rules that would fire for an unseen patient, we collected information on the number of rules matching each test case patient for the evaluated levels of m and l.
Cohen’s Kappa coefficient was computed for both the selected RSAbased models and the comparative models [62]. Cohen’s Kappa coefficient is designed to measure the agreement between two classification methods, but it is commonly used to measure model performance by comparing a classifier with a random allocation of patients among the decision classes. A value of zero indicates classification accuracy equivalent to chance (zero disagreement).
Performance of the prognostic models was evaluated using a 5fold cross validation procedure [63] wherein training and testing sets are repeatedly selected. Cross validation is a well known method that provides a reasonable estimate of the generalization error of a prediction model. In 5fold cross validation, the entire dataset is randomly divided into five subsets, or folds, and then each fold (20 % of the dataset) is used once as a testing set, with the remaining folds (80 %) used for training.
Results
This section presents the results obtained using MODLEM, VCDomLEM, logistic regression, SVM, C4.5 and RF models for sixmonth life expectancy prognostication of terminally ill patients. The results are analyzed and compared.
In order to select appropriate m and l values for MODLEM and VCDomLEMbased models, respectively, the performance of these models was evaluated first. AUC and coverage for each evaluated m and l level are shown in Table 4. Figures 3 and 4 display the number of rules that fire for each patient in the five testing folds for each m and l value. Based on these results, m=l=0.6 was chosen as the rule consistency parameter for both algorithms for further evaluation with the comparative methods.
The quality of approximation is 0.9244 for the CRSA, 0.3110 for the DRSA and finally 0.9014 for the VCDRSA where the object consistency parameter l=0.6.
Table 5 describes the number of rules and the number of descriptors in each rule for the two rough set approachbased classifiers at the selected consistency level of 0.6. The average number of MODLEM decision rules in the five rule sets generated by cross validation is 773 rules, with mean and maximum length of 3.65 and 8 descriptors, respectively. In Fig. 3, it can be seen that at rule consistency levels of m=0.2 and m=0.1, the number of rules fired per patient is always 2. This is because the rule set is generated by only two attributes and each rule contains only one attribute in the antecedent. The VCDomLEM decision rules are on average slightly longer, with mean and maximum length of 6.85 and 13 elementary conditions, respectively. The mean total number of VCDomLEM rules is 1,095 rules.
For SVM, the gamma (γ) and cost parameter (C) were evaluated between 10^{−1} and 10^{5} at increments of 10^{−1}; final selected parameters were γ=0.1 and C=100. For RF, the number of trees was explored between 10 and 1,000 trees at intervals of 10; the optimal number of trees thus obtained was 500. The maximum number of attributes selected at each bootstrap iteration was also explored in the range of 1 to 15 attributes, with best performance observed when the number of attributes was limited to 1. In the case of C4.5, the confidence factor used for pruning was evaluated between 0.1 and 0.9 with increments of 0.1 and 0.5 was selected. The minimum number of instances per leaf for the C4.5 decision tree was explored in steps of 1 between 1 and 100, with best performance achieved with a minimum of 40 instances per leaf. The pruned C4.5 trees contained an average of 74 nodes over the 5 cross validation folds.
The performance of all of the evaluated classification models is shown in Table 6, where Cohen’s kappa coefficient [62] and AUC are reported for each classifier, averaged over the 5 cross validation folds. Highest average kappa coefficient was achieved by RF with \(\bar {\kappa } = 0.37\). Second highest average kappa coefficient was achieved by VCDomLEM, logistic regression and SVM at \(\bar {\kappa } = 0.35\). The MODLEMand C4.5 classifiers achieved \(\bar {\kappa } = 0.32\) and 0.31, respectively. Average sensitivity and specificity for each of the models are also shown in Table 6. For each model and cross validation fold configuration, the sensitivity and specificity were recorded at the threshold at which both values are simultaneously maximized. This threshold is equivalent to the point on the ROC plot closest to the upper left corner and represents the point of maximum accuracy of the model.
Discussion
All of the methodologies show fair classification accuracy given that Kappa coefficients are in the range of 0.20 to 0.40 [64]. The results presented in Table 6 show that all of the methods have similar AUC with the best performing algorithm being RF (AUC = 0.7459) and the worst being MODLEM (AUC = 0.6974). The best performing method among the decision and rulebased methods was VCDomLEM with an average AUC of 0.7173.
With respect to MODLEM and VCDomLEM, m and l are clearly critical values in determining model performance. Together, Table 4 and Figs. 3 and 4 demonstrate that selecting m=l=0.6 balances the accuracy and coverage achieved by the rough set based classifiers against the amount of inconsistency allowed in each. In the case of MODLEM, m=0.6 is associated with highest AUC and acceptable coverage. However, in the case of VCDomLEM, reducing l below 0.6 provides only marginal benefits in terms of AUC and coverage but greatly increases the amount of inconsistency allowed in the generated rules.
The quality of approximation for the CRSA classifier is 0.9244. The difference between the quality of approximation in the CRSA and the DRSA is the inclusion of the preferenceordering information, determined by the APACHE III scores. In the case of the DRSA, a strict application of this information in determining the lower approximation leads to few patients in the lower approximations, thus reducing the overall quality of approximation. Consequently, decision rules generated from this approximation are too specific and less suitable for generalizing to the classification of new cases. It is therefore reasonable to relax the conditions for assignment of objects to lower approximations. Thus, using the VCDRSA and setting the object consistency parameter l=0.6, results in an improved quality of approximation of 0.9014.
All of the rule or decisiontreebased methods demonstrated somewhat reduced performance when compared with the nonrule based classifiers. The worstperforming rulebased method, MODLEM, had an AUC 0.049 below the bestperforming method, RF (0.6974, MODLEM, vs. 0.7459, RF). In contrast, VCDomLEM demonstrated average AUC much closer to that of RF, with an average difference of only 0.029 (0.7173, VCDomLem, vs. 0.7459, RF). In practice, this relatively small difference in performance is likely to be outweighed by the accessibility of the rulebased format of VCDomLEM, while such benefits would be less justified in the case of MODLEM.
Interpretation and usability of decision rules
Clinical credibility in prognostic models depends in part on the ease with which physicians and patients can understand and interpret the results of the models, in addition to the accuracy of the information they provide. The RSAbased prognostic models present the physician with a list of matched decision rules, offering significant advantages by increasing both the traceability of the model and the amount of information included in its results. This advantage is further increased in the case of VCDomLEM, where dominancebased decision rules permit greater information density per rule by including attribute value ranges in each rule.
Table 7 contains the decision rules that fire for an example patient selected from the SUPPORT data set. This patient was 41 years old with a primary diagnosis of coma. The patient displayed moderate head injury on the Glasgow Coma Scale, elevated levels of creatinine (1.60 mg/dL) and respiratory rate (26 bpm), normal levels of sodium (133 mEq/L), low white blood cell count (1.90 cells/nL) and mean blood pressure of 107 bpm. Both the MODLEM and VCDomLEM classifiers correctly predict that the patient will not survive six months (the patient in fact survived only 4 days).
The VCDomLEM classifier predicts d.6months = Yes with an associated score of 80 %, based on the two rules (Rules 5 and 6). As can be seen in Table 7, Rule 5 isolates the combination of Coma and elevated creatinine and sodium levels as a key predictor of sixmonth survival. In the case of Rule 5, 51 patients in the training set have similar conditions as the example patient, of which 47 did not survive six months. On the other hand, Rule 6 somewhat counterbalances this prediction, pointing to 8 young patients with moderate coma who have been in the hospital less than 44 days, of whom all 8 survived six months.
The MODLEM classifier provides a less specific prediction, classifying the example patient as not surviving six months with an associated score of 55 %. Upon further investigation, the rules matching the example patient (Rules 1–4) are more general than the rules provided by the VCDomLEM classifier. Rules 1–3 provide general rules that point to the age, level of head trauma and primary diagnosis of the patient. Considering only these three rules, the associated score would be d.6months = Yes with a score of 54 %, but this score is revised slightly by Rule 4 further in favor of d.6months = Yes. Rule 4 isolates normal average heart beat, high respiratory rate and low (and also very high) white blood cell counts.
For all of the classifiers, a final prediction and associated score are presented by the classifier. However, only in the case of MODLEM and VCDomLEM is the prediction further supported by the set of rules from which said prediction derived. Thus, the gestalt survival expectation is presented without loss of contradictory information, providing the physician with both the prognostication as well as supporting and contradicting information. In contrast, while a decision tree obtained using C4.5 can be represented as a set of rules, only a single rule representing the matching terminal node is returned to the physician. Among the rulebased methods, those rules derived from the VCDRSA naturally include attribute value ranges for which the rule is valid, succinctly providing even more information to the physician and further increasing the utility of the life expectancy prediction. In a clinical setting, this set of rules serves to support clinical decisions for future treatment or palliative care strategies as well as to support the explanation of these decisions to the involved patient and their family.
Decision tree models offer the additional benefit of visually representing the entire model in a single structure, and given their hierarchical structure can be used to guide the decision process of the physician [65]. Decision trees models are most useful when built with the input of domain experts as pruning visually complex decision trees must balance tradeoffs between accuracy and simplicity [66]. Many methods exist for the visualization of decision trees and the performance of visuallytuned decision trees may be comparable to more complex versions of the same model [67].
A further benefit of the rulebased methods is that the rules clearly indicate the patient characteristics most relevant to their survival expectation. This increases the transparency and interpretability of the classification process, strengthening the accessibility, and hence credibility, of the model. Additionally, the decision rules do not individually involve all of the condition attributes. This offers the advantage of providing potentially acceptable results should a particular prognostic factor be difficult or too costly to ascertain for a patient [34].
This is in stark contrast to SVM, neural networks, and other blackbox methods where very little insight is available to a decision maker as to how an outcome was predicted. While similar performance in terms of accuracy was seen for all of the classification models, the rulebased models naturally express results in terms of a set of decision rules, a benefit that is not present in logistic regression, RF, or the mentioned blackbox methods. As an ensemble method, the RF method functionally reduces to a blackbox style model, despite its use of decision trees.
Decision analysis for hospice referral
Consider the costs—economic, emotional and physical—associated with the decision to enter hospice care. These costs are justified for patients who either enter hospice care at the appropriate time or for those who do not enter hospice care when they could benefit from curative treatment. These cases represent true positive and true negative classifications. A higher emotional and physical cost is born by patients sent to hospice care but who ultimately survive six months—a false positive. The highest cost of all, emotionally, economically and physically is born by the patient and his or her family when costly treatment is prolonged for a patient who should have been referred to a hospice care program—a false negative. In this last case, some or all of the benefits of hospice care would be lost while the stresses and economic burden of aggressive treatment are endured.
In this light, the threshold parameter, τ (described in section Decision rules), can be seen as a representation of the patient and family’s preference for hospice care treatment and their risk tolerance for a mistaken referral. The threshold parameter relates sensitivity to specificity and stipulates the required level of certainty for a positive classification. A higher threshold value requires a higher probability of not surviving six months for the classification of a patient as a hospice candidate, decreasing the sensitivity and increasing specificity (indicating a preference for continued treatment). Conversely, a lower threshold value increases sensitivity while reducing specificity, indicating a preference for avoiding the costly mistake of unnecessary treatment.
As this threshold value is a subjective matter and varies between physicians, patients and family members, one suggested approach [68] involves the measurement of the amount of regret the decision maker would have should an incorrect decision be made. As medical decisions must take into account the preferences of those ultimately affected by the decision, this application of regret theory allows for the formal treatment of those preferences by calculating the threshold value as a function of the measured anticipated regret.
Conclusions
This paper contributes to the growing body of research in RST—and its extensions—as a prognostic modeling framework and highlights the strengths of this approach in terms of accessibility. The nonrulebased methods—RF, logistic regression, and SVM—were found to more accurately predict death within six months, however the benefits of the rulebased methods may outweigh the performance differential, particularly in the case of VCDomLEM where this difference was small. The intuitive structure of the rough set approaches, built on indiscernibility and dominance relations and expressed in terms of ifthen decision rules, offers both more insight into the modeling process and more opportunity for the knowledge extraction process to incorporate the personal preferences of those making and being affected by the decision.
The performance of the classifiers presented in this study is good but suboptimal, indicative of a challenging problem in need of further research. The increased performance achieved by the variable consistency approach suggests a dataset of highly diverse patients. Future research will explore methods to improve the overall classifier performance and address this diversity by building localized models for patient subgroups using rough sets concepts to group patients with similar differentiating characteristics.
A recent study developed a sixmonth survival prognostic model primarily based on the Medicare Health Outcomes Survey responses of communitydwelling elderly patients [69]. This model, named the PatientReported Outcome Mortality Prediction Tool (PROMPT), achieved comparable AUC using only basic medical information, indicating that the performance of classification models for sixmonth survival is still a major issue for the targeted domain of hospice referral recommendation.
An important limitation of this study is that patientspecific disease progression over time is not considered, in part due to the static nature of the data set used. Future research must address the temporal aspect of disease progression, a consideration often missing in other prognostic models for hospice referral. The progression of a terminal illness is often highly nonlinear by nature and generally does not present as a steady decline over time but rather as periods of relative stability marked by turning points of acute decline. A prognostic model that takes into account this temporal aspect may possibly provide both more accurate life expectancy prognoses and more useful information for palliative care planning.
Availability of supporting data
The data set supporting the results of this article is publicly available at http://biostat.mc.vanderbilt.edu/wiki/Main/SupportDesc.
References
 1
Kane RL, Klein SJ, Bernstein L, Rothenberg R, Wales J. Hospice role in alleviating the emotional stress of terminal patients and their families. Medical Care. 1985; 23(3):189–97.
 2
Bulkin W, Lukashok H. Rx for dying: The case for hospice. N Engl J Med. 1988; 318(6):376–8.
 3
Dawson NJ. Need satisfaction in terminal care settings. Soc Sci Med. 1991; 32(1):83–7.
 4
Christakis NA. Timing of referral of terminally ill patients to an outpatient hospice. J Gen Intern Med. 1994; 9(6):314–20.
 5
Teno JM, Shu JE, Casarett D, Spence C, Rhodes R, Connor S. Timing of referral to hospice and quality of care: length of stay and bereaved family members’ perceptions of the timing of hospice referral. J Pain Symptom Manag. 2007; 34(2):120–5.
 6
Christakis NA, Escarce JJ. Survival of medicare patients after enrollment in hospice programs. N Engl J Med. 1996; 335(3):172–8.
 7
Lee KL, Pryor DB, Harrell FE, Califf RM, Behar VS, Floyd WL, et al. Predicting outcome in coronary disease statistical models versus expert clinicians. Am J Med. 1986; 80(4):553–60.
 8
Wyatt JC, Altman DG. Commentary: Prognostic models: clinically useful or quickly forgotten?BMJ. 1995; 311(7019):1539–41.
 9
Dawes RM, Faust D, Meehl PE. Clinical versus actuarial judgment. Science. 1989; 243(4899):1668–74.
 10
Beck JR, Pauker SG, Gottlieb JE, Klein K, Kassirer JP. A convenient approximation of life expectancy (the “deale”) ii: use in medical decisionmaking. Am J Med. 1982; 73(6):889–97.
 11
Knaus WA, Wagner DP, Draper EA, Zimmerman JE, Bergner M, Bastos PG, et al. The APACHE III prognostic system. Risk prediction of hospital mortality for critically ill hospitalized adults. Chest. 1991; 100(6):1619–36.
 12
Knaus WA, Harrell FE, Lynn J, Goldman L, Phillips RS, Connors AF, et al. The SUPPORT prognostic model: Objective estimates of survival for seriously ill hospitalized adults. Ann Intern Med. 1995; 122(3):191–203.
 13
Christakis NA. Predicting patient survival before and after hospice enrollment. Hosp J. 1998; 13:71–88.
 14
Gripp S, Moeller S, Bölke E, Schmitt G, Matuschek C, Asgari So. Survival prediction in terminally ill cancer patients by clinical estimates, laboratory tests, and selfrated anxiety and depression. J Clin Oncol. 2007; 25(22):3313–20.
 15
Glare P, Sinclair C, Downing M, Stone P, Maltoni M, Vigano A. Predicting survival in patients with advanced disease. Eur J Cancer. 2008; 44(8):1146–56.
 16
Hyodo I, Morita T, Adachi I, Shima Y, Yoshizawa A, Hiraga K. Development of a predicting tool for survival of terminally ill cancer patients. Jpn J Clin Oncol. 2010; 40(5):442–8.
 17
Han PKJ, Lee M, Reeve BB, Mariotto AB, Wang Z, Hays RD, et al. Development of a prognostic model for sixmonth mortality in older adults with declining health. J Pain Symptom Manag. 2012; 43(3):527–39.
 18
Hosmer DW, Lemeshow S, Inc N. Applied Survival Analysis: Regression Modeling of Time to Event Data. New Jersey, USA: John Wiley & Sons, Inc; 2008, p. 404.
 19
Dybowski R, Gant V, Weller P, Chang R. Prediction of outcome in critically ill patients using artificial neural network synthesised by genetic algorithm. The Lancet. 1996; 347(9009):1146–50.
 20
Jayasurya K, Fung G, Yu S, DehingOberije C, De Ruysscher D, Hope A, et al. Comparison of bayesian network and support vector machine models for twoyear survival prediction in lung cancer patients treated with radiotherapy. Med Phys. 2010; 37(4):1401–7.
 21
Clermont G, Angus DC, DiRusso SM, Griffin M, LindeZwirble WT. Predicting hospital mortality for patients in the intensive care unit: a comparison of artificial neural networks with logistic regression models. Crit Care Med. 2001; 29(2):291–296.
 22
Segal MR. Features of treestructured survival analysis. Epidemiology. 1997; 8(4):344–6.
 23
Hwang SS, Scott CB, Chang VT, Cogswell J, Srinivas S, Kasimis B. Prediction of survival for advanced cancer patients by recursive partitioning analysis: role of karnofsky performance status, quality of life and symptom distress. Cancer Invest. 2004; 22(5):678–87.
 24
Bazan J, Osmólski A, Skowron A, Ślçezak D, Szczuka M, Wroblewski J. Rough set approach to the survival analysis In: Alpigini J, Peters J, Skowron A, Zhong N, editors. Rough Sets and Current Trends in Computing. Lecture Notes in Computer Science, vol. 2475. Malvern, PA, USA: Springer: 2002. p. 522–529.
 25
Pattaraintakorn P, Cercone N, Naruedomkul K. Hybrid rough sets intelligent system architecture for survival analysis In: Peters JF, Skowron A, Marek VW, Orłowska E, Słowinski R, Ziarko W, editors. Transactions on Rough Sets VII. Berlin Heidelberg: Springer: 2007. p. 206–24.
 26
Delen D, Walker G, Kadam A. Predicting breast cancer survivability: a comparison of three data mining methods. Artif Intell Med. 2005; 34(2):113–27.
 27
Cios KJ, Moore GW. Uniqueness of medical data mining. Artif Intell Med. 2002; 26(12):1–24.
 28
Simons P. Critical notice of timothy williamson, nagueness. Int J Philos Stud. 1996; 4:321–7.
 29
Pawlak Z. Rough Sets: Theoretical Aspects of Reasoning About Data. Norwell, MA: Springer; 1992.
 30
Pawlak Z. Vagueness a rough set view In: Mycielski J, Rozenberg G, Salomaa A, editors. Structures in Logic and Computer Science, vol. 1261. Springer Berlin Heidelberg: 1997. p. 106–117, doi:10.1007/3540632468_7.
 31
Pattaraintakorn P, Cercone N. Integrating rough set theory and medical applications. Appl Math Lett. 2008; 21(4):400–3.
 32
Hart A, Wyatt J. Evaluating blackboxes as medical decision aids: issues arising from a study of neural networks. Inform Health Soc Care. 1990; 15(3):229–36.
 33
Tsumoto S. Modelling medical diagnostic rules based on rough sets In: Polkowski L, Skowron A, editors. Rough Sets and Current Trends in Computing. Lecture Notes in Computer Science, vol. 1424. Berlin Heidelberg: Springer: 1998. p. 475–82, doi:10.1007/3540691154_65.
 34
Komorowski J, Øhrn A. Modelling prognostic power of cardiac tests using rough sets. Artif Intell Med. 1999; 15(2):167–91.
 35
Paszek P, WakuliczDeja A. Applying rough set theory to medical diagnosing In: Kryszkiewicz M, Peters J, Rybinski H, Skowron A, editors. Rough Sets and Intelligent Systems Paradigms. Lecture Notes in Computer Science, vol. 4585. Berlin Heidelberg: Springer: 2007. p. 427–35.
 36
Ningler M, Stockmanns G, Schneider G, Kochs HD, Kochs E. Adapted variable precision rough set approach for EEG analysis. Artif Intell Med. 2009; 47(3):239–61.
 37
LongJun H, Lipin D, CaiYing Z. Prognosis system for lung cancer based on rough set theory. In: Proceedings of the 2010 Third International Conference on Information and Computing, vol. 4. Washington, DC, USA: IEEE Computer Society: 2010. p. 7–10.
 38
Son CS, Kim YN, Kim HS, Park HS, Kim MS. Decisionmaking model for early diagnosis of congestive heart failure using rough set and decision tree approaches. J Biomed Inform. 2012; 5(45):999–1008.
 39
Greco S, Matarazzo B, Slowinski R. Rough sets theory for multicriteria decision analysis. Eur J Oper Res. 2001; 129(1):1–47.
 40
Stefanowski J. The rough set based rule induction technique for classification problems. In: Proceedings of 6th European Congress on Intelligent Techniques and Soft Computing EUFIT 98. Aachen: IEEE Computer Society: 1998. p. 109–113.
 41
Stefanowski J. On combined classifiers, rule induction and rough sets In: Peters J, Skowron A, Düntsch I, GrzymałaBusse J, Orłowska E, Polkowski L, editors. Transactions on Rough Sets VI, Lecture Notes in Computer Science, vol. 4374. Springer: 2007. p. 329–50.
 42
Blaszczyński J, Slowiński R, Szelag M. Probabilistic rough set approaches to ordinal classification with monotonicity constraints In: Hüllermeier E, Kruse R, Hoffmann F, editors. Computational Intelligence for KnowledgeBased Systems Design. Lecture Notes in Computer Science, vol. 6178. Berlin Heidelberg: Springer: 2010. p. 99–108.
 43
Bazan JG, Nguyen HS, Nguyen SH, Synak P, Wróblewski J. Rough set algorithms in classification problem In: Polkowski L, Tsumoto S, Lin T, editors. Rough Set Methods and Applications. Studies in Fuzziness and Soft Computing, vol. 56. Heidelberg, Germany: PhysicaVerlag GmbH: 2000. p. 49–88.
 44
Han J. Data Mining: Concepts and Techniques. San Francisco, CA, USA: Morgan Kaufmann Publishers Inc.; 2005.
 45
GrzymalaBusse JW, Stefanowski J. Three discretization methods for rule induction. Int J Intell Syst. 2001; 16(1):29–38.
 46
Clark P, Boswell R. Rule induction with cn2: Some recent improvements In: Kodratoff Y, editor. Machine Learning – EWSL91. Lecture Notes in Computer Science, vol. 482. Berlin Heidelberg: Springer: 1991. p. 151–63.
 47
Blaszczynski J, Slowinski R, Szelag M. Sequential covering rule induction algorithm for variable consistency rough set approaches. Inf Sci. 2011; 181(5):987–1002.
 48
Harrell FE. SUPPORT Datasets. 2010. http://biostat.mc.vanderbilt.edu/wiki/Main/SupportDesc. Accessed 17 Nov 2015.
 49
Stein SC. Minor head injury: 13 is an unlucky number. The Journal of Trauma and Acute Care Surgery. 2001; 50(4):759–60.
 50
Martin L. Reviews, notes, and listings: pulmonary medicine: All you really need to know to interpret arterial blood gases. Ann Intern Med. 1993; 8(118):656.
 51
Cohen J. The cost of dichotomization. Appl Psychol Meas. 1983; 7(3):249–53.
 52
Van Belle G. Statistical Rules of Thumb vol. 699. Hoboken, NJ, USA: John Wiley & Sons; 2011.
 53
Team RC. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. http://www.Rproject.org/. Accessed 17 Nov 2015.
 54
Błaszczyński J, Greco S, Matarazzo B, Słowiński R, Szelag M. jmaf  dominancebased rough set data analysis framework In: Skowron A, Suraj Z, editors. Rough Sets and Intelligent Systems  Professor Zdzisław Pawlak in Memoriam. Intelligent Systems Reference Library, vol. 42. Springer: 2013. p. 185–209.
 55
Chang CC, Lin CJ. LIBSVM: A library for support vector machines. ACM Trans Intell Syst Technol. 2011; 2(3):1–27.
 56
Bender R, Grouven U. Ordinal logistic regression in medical research. J R Coll Physicians Lond. 1997; 31(5):546–51.
 57
Cortes C, Vapnik V. Supportvector networks. Mach Learn. 1995; 20(3):273–97.
 58
Van Belle V, Pelckmans K, Van Huffel S, Suykens JA. Support vector methods for survival analysis: a comparison between ranking and regression approaches. Artificial Intelligence In Medicine. 2011; 53(2):107–18.
 59
Barakat N, Bradley AP. Rule extraction from support vector machines: a review. Neurocomputing. 2010; 74(13):178–90.
 60
Quinlan JR. C4. 5: Programs for Machine Learning vol. 1. San Francisco, CA, USA: Morgan Kaufmann Publishers Inc.; 1993.
 61
Breiman L. Random forests. Mach Learn. 2001; 45(1):5–32.
 62
Cohen J. A coefficient of agreement for nominal scales. Educ Psychol Meas. 1960; 20(1):37–46.
 63
Kohavi R. A study of crossvalidation and bootstrap for accuracy estimation and model selection. In: Proceedings of the 14th International Joint Conference on Artificial Intelligence  Volume 2. IJCAI’95. San Francisco, CA, USA: Morgan Kaufmann Publishers Inc.: 1995. p. 1137–1143.
 64
BenDavid A. Comparison of classification accuracy using Cohen’s weighted kappa. Expert Systems with Applications. 2008; 34(2):825–32.
 65
Podgorelec V, Kokol P, Stiglic B, Rozman I. Decision trees: an overview and their use in medicine. J Med Syst. 2002; 26(5):445–463.
 66
Bohanec M, Bratko I. Trading accuracy for simplicity in decision trees. Mach Learn. 1994; 15(3):223–50.
 67
Stiglic G, Kocbek S, Pernek I, Kokol P. Comprehensive decision tree models in bioinformatics. PLoS ONE. 2012; 7(3):33812.
 68
Tsalatsanis A, Hoz I, Andrew V, Djulbegovic B, Hozo I, Vickers A, Djulbegovic B. A regret theory approach to decision curve analysis: a novel method for eliciting decision makers’ preferences and decisionmaking. BMC Med Inform Decis Mak. 2010; 10(51):51.
 69
Han PKJ, Lee M, Reeve BB, Mariotto AB, Wang Z, Hays RD, et al. Development of a prognostic model for sixmonth mortality in older adults with declining health. J Pain Symptom Manag. 2011; 43(3):527–39.
Acknowledgements
This work was supported by the Department of Army under grant # W81 XWH0920175.
The authors thank Dr. Jurek Blaszczynski from the Laboratory of Intelligent Decision Support Systems at the Poznan University of Technology, Poznan, Poland, for his collaboration in providing valuable information about the Dominance Based Rough Set Approach.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
EG and GA are the principal authors; EG performed the data analysis and implemented the different algorithms. EG and GA wrote the manuscript. GA and AY participated in writing/editing and critically revising the manuscript with substantial contributions. AY, AT, BD and LB contributed to the interpretation of the data and revised the manuscript for important intellectual content. All authors read and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Rough set theory
 Dominancebased rough set approach
 Hospice referral
 Prognostic models