Computational and Mathematical Methods in Medicine

Computational and Mathematical Methods in Medicine / 2013 / Article

Research Article | Open Access

Volume 2013 |Article ID 904860 |

Benjamin W. Y. Lo, R. Loch Macdonald, Andrew Baker, Mitchell A. H. Levine, "Clinical Outcome Prediction in Aneurysmal Subarachnoid Hemorrhage Using Bayesian Neural Networks with Fuzzy Logic Inferences", Computational and Mathematical Methods in Medicine, vol. 2013, Article ID 904860, 10 pages, 2013.

Clinical Outcome Prediction in Aneurysmal Subarachnoid Hemorrhage Using Bayesian Neural Networks with Fuzzy Logic Inferences

Academic Editor: Nestor V. Torres
Received24 Jan 2013
Accepted23 Mar 2013
Published10 Apr 2013


Objective. The novel clinical prediction approach of Bayesian neural networks with fuzzy logic inferences is created and applied to derive prognostic decision rules in cerebral aneurysmal subarachnoid hemorrhage (aSAH). Methods. The approach of Bayesian neural networks with fuzzy logic inferences was applied to data from five trials of Tirilazad for aneurysmal subarachnoid hemorrhage (3551 patients). Results. Bayesian meta-analyses of observational studies on aSAH prognostic factors gave generalizable posterior distributions of population mean log odd ratios (ORs). Similar trends were noted in Bayesian and linear regression ORs. Significant outcome predictors include normal motor response, cerebral infarction, history of myocardial infarction, cerebral edema, history of diabetes mellitus, fever on day 8, prior subarachnoid hemorrhage, admission angiographic vasospasm, neurological grade, intraventricular hemorrhage, ruptured aneurysm size, history of hypertension, vasospasm day, age and mean arterial pressure. Heteroscedasticity was present in the nontransformed dataset. Artificial neural networks found nonlinear relationships with 11 hidden variables in 1 layer, using the multilayer perceptron model. Fuzzy logic decision rules (centroid defuzzification technique) denoted cut-off points for poor prognosis at greater than 2.5 clusters. Discussion. This aSAH prognostic system makes use of existing knowledge, recognizes unknown areas, incorporates one's clinical reasoning, and compensates for uncertainty in prognostication.

1. Introduction

Advances in biostatistics and computing in the past several decades have led to creation of different types of clinical outcome prediction models. Three of these include artificial neural networks, fuzzy logic and bayesian analysis [13]. These techniques complement classical or frequentist approaches, such as regression analysis.

Artificial neural networks mimic biological neural systems. In biological systems, incoming dendrites collect signals which are fed to the neuron. A signal summation is then sent as a spike of electrical current along an axon, with resultant discharge at the synapse, connecting it to other neurons. Examples of biological neural networks include the human brain and the human retina. Analogous to the biological system, artificial neural networks are made up of a group of input variables which converge on a number of nodes. Nodes are grouped in layers, with interconnection links among themselves. Hidden or latent variables can exist in one or two layers. After processing from different activation functions, output signals are then sent onto output nodes in the network. artificial neural networks assume all or none logic, that is, subjects are regarded as having or not having a diagnosis. Nodes in the neural nets are connected with each other via connection links. Each of these links has an associated weight and activation function. Neural networks are intelligent systems that can learn and change behaviour by themselves as they gain experience. In addition, they also take into account unobservable variables that the researcher is not aware of while designing the neural net.

Assuming a basic artificial neural network with inputs, hidden or latent variables (in 1 layer), and an outputs, using a multilayer perceptron model as illustrated in Figure 1, output is equal to the summation of input to hidden layer, as well as hidden layer to output layer.

If the activation function is the nonlinear hyperbolic tangent function, then, where hidden unit = hyperbolic tangent function × (bias term + sum of (weights from input unit to hidden unit   × input unit)), and, where output = = bias term + sum of (weights on connection from hidden unit to output unit   ×  hidden unit).

Bayesian analysis allows the researcher to make use of existing states of knowledge before incorporation of new data. Simply put, it reflects the fact that knowledge is cumulative. Here, existing knowledge is expressed in the form of distributions (such as the normal bell-shaped distribution). The “prior” distribution is then combined with its likelihood of occurrence, forming a posterior probability. The end result (or posterior probability) represents a revised or updated belief after taking new data into account. If there is a lack of existing knowledge on the subject of interest, the researcher can still rely on Bayesian techniques. In this case, vague or uninformed prior probabilities are used.

In the Bayesian approach to artificial neural networks [4], the goal is to find the predictive distribution for target values in the new test case/model, given inputs for that case and inputs/targets in training cases. Here, represents the probability of data according to a particular model. It is an integral, representing the summation of all possible parameter values weighted by the strength of belief (as assigned by the researcher) in these parameter values, or The probability of a new case/model given existing cases and associated parameters (, incorporating both weights and biases) is expressed as In general, Bayesian neural networks can be expressed as: Posterior probability density is proportional to product of prior probability density and its associated likelihood. Likelihood, as explained above, is the product of probabilities of data given parameters (weights and biases) as Fuzzy logic is an extension of neural nets, but with the distinct advantage that it can assume functions with any value from 0 to 1, accounting for the entire spectrum of certainty of diagnoses and spectrum of severity of diseases studied. It registers a mild case of a certain disease, as it recognizes grey zones in diagnoses. Another strength of fuzzy logic lies in its explicit knowledge representation; that is, it allows the clinician to explicitly state its inputs, control actions, and outputs. The clinician can also clarify, or defuzzify, the entire process by carrying out crisp control actions, such as adding a cut-off level for prognoses and diagnoses, and trigger thresholds for treatment. By doing so, all actions (verification) in fuzzy logic are accounted for, and optimization is achievable. Thus, fuzzy logic systems are efficient ones. On the other hand, some regard this strength of fuzzy logic as a limitation. Unlike neural nets, where learning is done by the nets themselves as they gain experience from datasets, fuzzy logic systems cannot train themselves. The designer has to derive all action commands (in the form of “if-then” rules) manually, which is labour intensive with large datasets. When the degrees of various problems are combined in an equation, however, the resulting calculation may represent a more accurate prediction of the real data and individual patient.

Bayesian neural networks with fuzzy logic inferences can be represented as follows: where defuzz can be max-min, centroid, left of mean, right of mean, or another defuzzification crisp control action rule.

2. Objectives and Relevance

In this paper, we aim to use advanced biostatistical methods to create clinical prognostic decision rules in aneurysmal subarachnoid hemorrhage derived from a large aneurysmal subarachnoid hemorrhage database, which can be tailored to a specific patient population. We explore novel methods that account for existing states of knowledge (Bayesian meta-analysis and regression), complex nonlinear relationships between independent, latent, and dependent variables (artificial neural networks), and grey zones in prognoses (fuzzy Logic decision rules). In fact, the combination of such techniques can represent a novel health research method for clinical outcome prediction applicable to many diseases in medicine.

3. Methods

Intracranial aneurysmal subarachnoid hemorrhage affects about 45, 000 individuals in North America annually. Aneurysmal subarachnoid hemorrhage is associated with a mortality rate of at least 45% in the first 30 days following rupture [5, 6]. Apart from the primary neurological injury from the aneurysmal rupture itself, other secondary injury processes can further worsen an individual’s neurological condition and eventual clinical outcome. These processes include both neurological processes (such as delayed stroke, rebleeding, brain swelling, vasospasm induced strokes, seizures, and hydrocephalus) and systemic medical complications (such as myocardial infarction, fever, pulmonary edema). and Taken together, these processes can lead to long-term disability. Types of disability include physical, neurocognitive, and psychological impairment [7].

Tirilazad aneurysmal subarachnoid hemorrhage database is used to illustrate prognostic decision principles derived from a combination of techniques from multiple linear regression, artificial neural networks, fuzzy Logic and bayesian analysis. Tirilazad is a 21-aminosteroid compound produced by Pharmacia & Upjohn, Kalamazoo, MI, USA, originally investigated by the University of Virginia Health Sciences Center, as a free radical scavenger for potential treatment of cerebral vasospasm. This medication was investigated in five randomized clinical trials [812] involving patients with aneurysmal subarachnoid hemorrhage between 1990 and 1997 in 162 centers from 21 countries across North America, Europe, Australia, New Zealand, and South Africa. Tirilazad was found to have no effect on clinical outcome in patients with aneurysmal subarachnoid hemorrhage. The resultant database from these five studies contains 3550 patients, with its primary outcomes being Glasgow outcome score at 3 months and death from any cause. Glasgow outcome score is a 5-point neurological scale with the following designations: 5: good recovery—normal life activities despite minor deficits, 4: moderate disability—disabled but independent, 3: severe disability—conscious but disabled, 2: persistent vegetative state—unresponsive and speechless, and 1: death [13].

Centers followed strict treatment protocols, and variables had fewer than 5% missing entries. Tirilazad was administered in the intravenous form from day 3 to day 10 after subarachnoid hemorrhage onset. Only one percent of patients were lost to followup.

Patients in each treatment group were managed in a similar manner. Over 85% underwent surgical clipping, with 50% operated within the first 48 hours. Baseline demographics in both treatment and control arms were balanced in terms of gender, age, number of preexisting medical conditions (including hypertension, myocardial infarction, and angina), mean time to treatment, mean admission systolic blood pressure, admission neurological grade, ruptured aneurysm location, and admission amount of blood. These potential confounders were accounted for in statistical analysis. Proportions of patients experiencing vasospasm were similar in different treatment groups, as were the percentages of patients experiencing both neurological and systemic disabilities, including cerebral hemorrhage, cerebral ischemia, second stroke, rebleeding, hydrocephalus, sepsis, pulmonary embolus, brain herniation, pneumonia, and renal insufficiency.

A Cochrane systematic review [14] on the five trials on Tirilazad found no substantial heterogeneity among the trials, and there were no significant differences in adverse events between the treatment group and placebo group. The Tirilazad database represents the largest currently available aneurysmal subarachnoid hemorrhage clinical trial database worldwide.

Derivation of prognostic decision rules comprised of the following investigations: multivariable linear regression (IBM SPSS version 19.0 (Armonk, NY, USA)), artificial neural networks (IBM SPSS), Bayesian regression with uninformed and informed prior likelihoods (WinBUGS 1.4.3), neural networks with Bayesian regularization (MATLAB R2012a (Natick, MA, USA); version R 2.15.2), and creation of fuzzy logic decision rules, with clinical case application.

4. Results & Discussions

4.1. Multivariable Linear Regression Analysis

Frequentist linear regression was created using IBM SPSS. Included in the analysis were predictor variables from the Tirilazad database, without recoding, renaming, reclassification, or data transformation. Treatment variables were not included in this analysis. Predictor variables included were age, gender, neurological grade, intraventricular hemorrhage, subarachnoid hemorrhage thickness, time to treatment, clinical vasospasm, mean arterial pressure, aneurysm location, prior anticoagulation, eye opening, normal motor response, normal speech, admission angiographic vasospasm, intracerebral hemorrhage, hydrocephalus, prior subarachnoid hemorrhage, history of hypertension, history of myocardial infarction, history of angina, history of migraines, history of diabetes mellitus, antiepileptic use, ruptured aneurysm location, cerebral edema, pulmonary edema, vasospasm day, baseline temperature, fever on day 8, and cerebral infarction. Dependent variable is the patient’s clinical outcome (Glasgow outcome score) at 3 months.

Table 1 lists the statistically significant predictor () variables for neurologic outcome which include normal motor response, cerebral infarction, history of myocardial infarction, cerebral edema, history of diabetes mellitus, fever on day 8, prior subarachnoid hemorrhage, admission angiographic vasospasm, neurological grade, intraventricular hemorrhage, ruptured aneurysm size, history of hypertension, vasospasm day, age, and mean arterial pressure.

Independent variable value coefficientsCollinearity diagnostics
(95% confidence interval)(tolerance, VIF)

Normal motor response<0.001−0.329 (−0.496, −0.161)(0.27, 3.71)
Cerebral infarction<0.0010.790 (0.695, 0.885)(0.86, 1.16)
History of myocardial infarction0.0090.386 (0.097, 0.675)(0.92, 1.09)
Cerebral edema<0.0010.322 (0.190, 0.453)(0.96, 1.05)
History of diabetes mellitus0.0280.239 (0.026, 0.452)(0.98, 1.03)
Day-8 fever<0.0010.231 (0.150, 0.311)(0.93, 1.08)
Prior subarachnoid hemorrhage0.0040.197 (0.063, 0.332)(0.98, 1.02)
Admission angiographic vasospasm0.0150.175 (0.035, 0.315)(0.93, 1.08)
Neurological grade<0.0010.167 (0.093, 0.242)(0.16, 6.43)
Intraventricular hemorrhage0.0010.142 (0.056, 0.229)(0.80, 1.25)
Ruptured aneurysm size0.0010.130 (0.053, 0.206)(0.97, 1.03)
History of hypertension0.0090.119 (0.030, 0.208)(0.85, 1.18)
Vasospasm day0.050.112 (0.001, 0.225)(0.20, 5.11)
Age<0.0010.018 (0.015, 0.021)(0.86, 1.17)
Mean arterial pressure0.0120.003 (0.001, 0.006)(0.91, 1.10)

VIF: variance inflation factor.

All significant prognostic variables, with the exception of normal motor response, point to poorer prognosis. Presence of normal motor response at presentation signifies more favourable outcome. Its collinearity diagnostics reveal close correlation with other prognostic variables, namely, neurological grade.

Vasospasm day is closely correlated with other prognostic variables, namely, clinical and angiographic vasospasm.

Closer examination of the tirilazad patient dataset, using nontransformed data points, reveals that heteroscedasticity is present in the model. As heteroscedasticity is present, we cannot be confident that the strength of prediction of the linear regression equation from this multiple linear regression model is equally strong across all levels of the included independent variables.

Therefore, artificial neural networks were used to explore presence of complex nonlinear relationships and latent variables inherent in the database.

4.2. Bayesian Analysis

Beta coefficients of prognostic variables generated with Bayesian regression (WinBUGS version 1.4.3) using uninformed priors are similar in magnitude to those generated from frequentist multiple regression analysis, as demonstrated in Table 2.

Node Mean SD MC error 2.5%Median97.5%StartSample

b.MOTOR 500095001
b.CVA 500095001
b.BSWELL 500095001
b.MI 500095001
b.DM 500095001
b.ADMITVSP 500095001
b.PREVSAH 500095001
b.ANSIZE 500095001
b.IVH 500095001
b.HTN 500095001
b.VSPDAY 500095001
b.AGE 500095001
b.MAP 500095001
b.D8TEMP 500095001

Bayesian meta-analysis synthesizes research evidence from multiple independent studies. It has a hierarchical structure, placing one layer of sampling above another, with each study population having an observed odds ratio estimated from its sample of subjects. Bayesian meta-analysis models created using WinBUGS version 1.4.3 for this study gave generalizable posterior distributions of consensus odds ratios with representative medians, standard deviations, and 95% credible intervals, for the prognostic variables age [1521], neurological grade [1620, 22, 23], and aneurysmal size [17, 19, 2327]. Results of hierarchical meta-analysis are shown in Table 3.

Prognostic cariableOR

Neurological grade2.170.401.672.093.13
Aneurysm size1.

These values are very useful clinically by themselves, as they can be applied to patient prognostication. They also represent the informed priors for Bayesian regression models. Because of this database’s very large sample size (), beta coefficients generated with Bayesian regression (uninformed priors) are similar in magnitude to those generated from Bayesian regression (informed priors).

4.3. Artificial Neural Networks

Artificial neural networks (ANN), created using IBM SPSS, for prognosis in aneurysmal subarachnoid hemorrhage using the Tirilazad database revealed the following features:Model type: Multilayer Perceptron (MLP)Number of layers: 3Layer 1: 30 input variablesActivation function: hyperbolic tangentLayer 2: hidden layer, 11 hidden or latent variablesLayer 3: 5 output nodesTraining: 60% of sample sizeTesting: 40% of sample sizeTraining algorithm: gradient descentModel sensitivity: area under ROC curve = 0.85

This ANN model recognizes the presence of complex nonlinear relationships between variables, not identified by multiple linear regression.

Variables, in order of magnitude of normalized importance in Table 4, include age, second stroke, myocardial infarction, temperature, mean arterial pressure, neurological grade, ruptured aneurysm size, diabetes mellitus, angina, subarachnoid clot thickness, lung edema, admission angiographic vasospasm, previous subarachnoid hemorrhage, vasospasm day, cerebral edema, vasospasm during treatment, aneurysm location, time to treatment, normal motor response, intracerebral hematoma, normal speech, day 8 temperature, gender, eye opening, migraine history, intraventricular hemorrhage, hypertensive history, anticoagulant use, seizures, and hydrocephalus.

Artificial neural networks
independent variable
Type of prognostic factorImportance

Age Demographic0.111
Second strokeNeurologic0.081
Myocardial infarctionSystemic0.075
Mean arterial pressureSystemic0.054
Neurological gradeNeurologic0.048
Ruptured aneurysm sizeNeurologic0.039
Diabetes mellitusSystemic0.037
SAH clot thicknessNeurologic0.033
Lung edemaSystemic0.032
Admission angiographic vasospasmNeurologic0.029
Previous subarachnoid hemorrhageNeurologic0.028
Vasospasm dayNeurologic0.028
Cerebral edemaNeurologic0.028
Vasospasm during treatmentNeurologic0.027
Aneurysm locationNeurologic0.025
Time to treatmentDemographic0.025
Normal motor responseNeurologic0.024
Intracerebral hematomaNeurologic0.022
Normal speechNeurologic0.021
Day-8 temperatureSystemic0.021
Eye openingNeurologic0.018
Migraine historyNeurologic0.015
Intraventricular hemorrhageNeurologic0.015
Hypertensive historySystemic0.014
Anticoagulant useSystemic0.014

Interrelationships between input nodes (predictor variables), hidden variables (11 of them in one hidden layer), and output nodes (Glasgow outcome score) are illustrated in Figure 2. Possible latent (unobservable) variables, not measured by investigators, include the following: disrupted cerebral autoregulation contributing to both ischemia and cerebral edema after subarachnoid hemorrhage,biochemical markers of brain injury predisposing to cortical depression,cellular markers demonstrating physiologic dysfunction (such as mitochondrial dysfunction as reflected by imbalance between oxygen supply and consumption),genetic factors affecting outcome (such as inheritance of genes making patients more prone to microthrombotic events in the cerebral microvasculature disrupting cerebral blood flow),multiple drug-drug interactions, including drug hypersensitivities, especially in the elderly aneurysmal subarachnoid hemorrhage patient with multiple preexisting comorbidities,variables known to affect clinical outcome that were not captured in the Tirilazad database, including smoking and alcohol consumption,multiorgan system dysfunctions and their influence on neurologic outcome (including interactions between the central nervous system and cardiovascular, respiratory, renal, immune/hematologic, gastrointestinal (including hepatic and splenic), endocrinologic, and metabolic homeostasis).

4.4. Neural Networks with Bayesian Regularization

Using MATLAB version R2012a (Neural Network Toolbox), the following network was created:Number of layers: 3Layer 1: 30 input nodes (hyperbolic tangent activation)Layer 2: 11 hidden nodes (linear activation)Layer 3: 5 output nodesTraining: 60% of sample size (2131 of 3551)Testing: 40% of sample size (1420 of 3551)Training algorithm: neural networks with bayesian regularization (“training” algorithm)

MATLAB generated the following covariance matrix from layer 2 to layer 3, representing weights, associated with their bayesian ranges, between hidden nodes and output nodes: [; ; ; ; ; ; ; ; ; ; ].

Due to Tirilazad database’s large sample size (), the posterior distributions of weights from neural networks with Bayesian regularization are very narrow. Hence, choosing the mean of each distribution is the most appropriate for the calculation of each variable’s normalized importance.

4.5. Fuzzy Logic Decision Rules

Risk factors (according to their adjusted synaptic weights from artificial neural networks and neural networks with Bayesian regularization models) can be classified into three clusters, (1) demographic (3 factors—age, time to treatment, and gender), (2) systemic (9 factors—myocardial infarction, temperature, mean arterial pressure, diabetes mellitus, angina, pulmonary edema, fever on day 8, hypertensive history, anticoagulant use), and (3) neurologic (18 factors—second stroke, neurological grade, ruptured aneurysm size, subarachnoid clot thickness, admission angiographic vasospasm, previous subarachnoid hemorrhage, vasospasm day, cerebral edema, vasospasm during treatment, aneurysm location, normal motor response, intracerebral hematoma, normal speech, eye opening, migraine history, intraventricular hemorrhage, seizures, and hydrocephalus).

4.6. Fuzzy Logic Rules

Individualized fuzzy logic rules can be derived based on one’s experience. In this case, our three linguistic variables are (1) demographic risk factor cluster, (2) systemic risk factor cluster, and (3) neurologic risk factor cluster.

Fuzzification begins with assigning each linguistic variable a range of membership functions. When all members of each cluster are present, then maximum membership function of 1 is reached for that particular linguistic variable.

Fuzzy inferences then proceed with derivation of “if-then” rules that define system behaviour. In our case (Figure 3), if cluster one (demographic risk factor cluster) is fulfilled, then, one has low suspicion for poor neurologic outcome. If members of both cluster one (demographic risk factor cluster) and cluster two (systemic risk factor cluster) are present, then one has raised suspicion for poor neurologic outcome. One has high suspicion for poor neurologic outcome if some members of cluster one (demographic risk factor cluster), cluster two (systemic risk factor cluster), and cluster three (neurologic risk factor cluster) are present.

Next, defuzzification step translates the linguistic variable results into the crisp control action of denoting high likelihood for poor outcome. The centroid rule is applied to designation for poor prognostication, whereby a patient fulfills risk factors from 2.5 clusters.

As an example, an elderly patient (demographic cluster), with a number of medical comorbidities (examples from systemic cluster, such as coronary artery disease, hypertensive, and diabetic), who experiences a number of neurological complications after treatment (examples from neurologic cluster, such as second stroke, cerebral ischemia, and seizures) is predicted to have a poor long-term neurologic outcome (poor three-month Glasgow outcome score).

5. Limitations

The techniques of bayesian neural networks with fuzzy logic inferences were applied to the Tirilazad database. We note that case mix in this database were patients who underwent surgical clipping of cerebral aneurysms. Since the conduct of the Tirilazad trials, there are advancements in both medical management and surgical treatment of cerebral aneurysms. These include improved neurocritical care of aneurysmal subarachnoid hemorrhage patients and aneurysmal coiling. In addition, the Tirilazad database did not include important prognostic variables such as smoking, alcohol consumption, rebleeding, and infection. In order to overcome these limitations, ongoing efforts are now underway to combine a number of aneurysmal subarachnoid hemorrhage databases worldwide in the multinational Subarachnoid Hemorrhage International Trialists (SAHIT) collaboration. Important prognostic variables as well as aneurysmal coiling patients will be included in this database. In addition, Bayesian neural networks with fuzzy logic inferences will be applied to the SAHIT database.

6. Conclusions

Complex relationships exist among heterogeneous groups of prognostic factors. The accuracy of clinical outcome prediction depends on clarification of these relationships. General linear models have been used frequently for decades. These popular techniques produce interpretable coefficients for explanatory variables and are easily estimated using commercially available statistical programs. In real life, however, data points rarely fit perfectly linear relationships. Greater deviations from linearity point to the need for exploratory analyses with complex nonlinear systems.

Typical artificial neural networks can fit training data with high precision and detect nonlinear relationships among predictor variables, with the overall aim of predicting yet to be seen observations. Neural networks also incorporate latent (unobserved variables) in one or two hidden layers. Interrelationships between independent, latent, and outcome variables are assigned synaptic weights, or connection strengths. A weighted average of these connection strengths gives a variable’s normalized importance, or the percentage contribution of each predictor variable to the overall clinical outcome, taking into account the error between predicted and actual values. The sum of all relative importance values of input variables (representing influences of predictor variables on clinical outcome in relation to the rest of the independent variables) equals 100 percent. Small sample sizes can affect model building, making it more difficult to distinguish between true signal and noise. Over- and underfitting can occur in these cases, which, in turn, can affect model generalization. Neural networks with Bayesian regularization technique have been devised to overcome the above problem, whereby weights are assigned probability density distributions, incorporating Bayesian statistics to estimate weight uncertainty, or the relative degree of belief in the different values for synaptic weights.

Typical artificial neural networks give fixed structures, whereas Bayesian neural networks give flexible structures. Bayesian neural networks are ones that assign probability distributions to all elements of the network, including inputs, hidden nodes, outputs, as well as their associated weights. Prior likelihoods for probability distributions can be generated using Bayesian meta-analysis. Bayesian regularization terms are included, which prevent model over- and underfitting. In addition, Bayesian Neural Networks give generalizable posterior distributions without compromising nonlinearity properties. Bayesian Neural Networks are trained by sampling from joint posterior likelihoods of network structure and weights by Monte Carlo sampling methods. This training avoids the problem of convergence at local minima. Posterior distributions of weights can be used to evaluate uncertainties of predictions of trained networks and can also be used to assess network sensitivities. The larger the sample size, the smaller the Bayesian posterior probability distribution ranges. If data points fall into a linear relationship, typical artificial neural networks and Bayesian neural networks can detect this relationship. Hence, linear regression can be viewed as a special case of neural networks.

Typical artificial neural network and Bayesian neural network error can be due to network weight uncertainty (model uncertainty due to imperfections in data, nonoptimum network structure, and nonoptimum learning algorithms), and error from remaining sources, including intrinsic noise that includes random error due to measurement noise, and error due to finite resolution of observation system. If posterior distributions of weights are very narrow in relation to noise distribution, then the width of distribution of networks outputs can be influenced by noise. On the other hand, if the posterior distributions of weights are larger than noise distribution, then, the width of network outputs is dominated by distribution of network weights.

Fuzzy logic approach to Bayesian neural networks allows the clinician to explore where within the Bayesian range the nonlinear relationship for a particular case is most likely to exist. Results generated from Bayesian neural networks with fuzzy logic inferences will, then, slightly differ from case to case, accounting for the special characteristics of that certain case.

Fuzzy logic inferences should be applied at the end of Bayesian neural network formulation. In other words, one should allow the Bayesian neural network learning machine to do its own learning before applying fuzzy logic rules, so that all probability distributions are explored.

Bayesian neural networks with fuzzy logic inference can be conceptually interpreted as follows. Based on one’s own experience (summation of existing parameters weighted by strength of belief in what happened beforehand), one can predict (based on one’s assigned strength of belief) where along a spectrum of probabilities of the unknown quantity a value will end up. If it falls outside the spectrum in real life, then, one has to check whether there are still unknown elements influencing the outcome variable in question.

Conflict of Interests

The authors of the paper do not have any financial relationship with pharmaceutical company, Pharmacia & Upjohn, Kalamazoo, MI, USA, manufacturer of the trial agent, Tirilazad.


  1. S. Haykin, Neural Networks: A Comprehensive Foundation, Macmillan, New York, NY, USA, 1994.
  2. K. Tanaka, An Introduction To Fuzzy Logic For Practical Applications, Springer, New York, NY, USA, 1997.
  3. D. Spiegelhalter, K. Abrams, and J. Myles, Bayesian Approaches To Clinical Trials and Health-Care Evaluation, Wiley, Chichester, UK, 2002.
  4. R. M. Neal, Bayesian Learning For Neural Networks, Springer, New York, NY, USA, 1996.
  5. R. L. Macdonald, M. D. Cusimano, N. Etminan et al., “Subarachnoid hemorrhage international trialists data repository (SAHIT),” World Neurosurgery, 2013. View at: Publisher Site | Google Scholar
  6. B. N. Jaja, M. D. Cusimano, R. L. Macdonald et al., “Clinical prediction models for aneurysmal subarachnoid hemorrhage: a systematic review,” Neurocritical Care, vol. 18, no. 1, pp. 143–153, 2013. View at: Google Scholar
  7. B. Lo and R. L. Macdonald, “Health related quality of life after aneurysmal subarachnoid hemorrhage: new data from a large series in Germany,” World Neurosurgery, 2012. View at: Publisher Site | Google Scholar
  8. E. C. Haley Jr., N. F. Kassell, W. M. Alves et al., “Phase II trial of tirilazad in aneurysmal subarachnoid hemorrhage. A report of the Cooperative Aneurysm Study,” Journal of Neurosurgery, vol. 82, no. 5, pp. 786–790, 1995. View at: Google Scholar
  9. N. F. Kassell, E. C. Haley, C. Apperson-Hansen, and W. M. Alves, “Randomized, double-blind, vehicle-controlled trial of tirilazad mesylate in patients with aneurysmal subarachnoid hemorrhage: a cooperative study in Europe, Australia, and New Zealand,” Journal of Neurosurgery, vol. 84, no. 2, pp. 221–228, 1996. View at: Google Scholar
  10. E. C. Haley Jr., N. F. Kassell, C. Apperson-Hansen, M. H. Maile, and W. M. Alves, “A randomized, double-blind, vehicle-controlled trial of tirilazad mesylate in patients with aneurysmal subarachnoid hemorrhage: a cooperative study in North America,” Journal of Neurosurgery, vol. 86, no. 3, pp. 467–474, 1997. View at: Google Scholar
  11. G. Lanzino and N. F. Kassell, “Double-blind, randomized, vehicle-controlled study of high-dose tirilazad mesylate in women with aneurysmal subarachnoid hemorrhage. Part II. A cooperative study in North America,” Journal of Neurosurgery, vol. 90, no. 6, pp. 1018–1024, 1999. View at: Google Scholar
  12. G. Lanzino, N. F. Kassell, N. W. C. Dorsch et al., “Double-blind, randomized, vehicle-controlled study of high-dose tirilazad mesylate in women with aneurysmal subarachnoid hemorrhage. Part I. A cooperative study in Europe, Australia, New Zealand, and South Africa,” Journal of Neurosurgery, vol. 90, no. 6, pp. 1011–1017, 1999. View at: Google Scholar
  13. B. Jennett and M. Bond, “Assessment of outcome after severe brain damage. A practical scale,” The Lancet, vol. 1, no. 7905, pp. 480–484, 1975. View at: Google Scholar
  14. L. C. Wang, M. Liu, B. Wu, and S. H. Zhang, “Tirilazad for aneurysmal subarachnoid haemorrhage,” Cochrane Database of Systematic Reviews, no. 2, Article ID CD006778, 2010. View at: Publisher Site | Google Scholar
  15. M. Holling, A. Jeibmann, J. Gerss et al., “Prognostic value of histopathological findings in aneurysmal subarachnoid hemorrhage: clinical article,” Journal of Neurosurgery, vol. 110, no. 3, pp. 487–491, 2009. View at: Publisher Site | Google Scholar
  16. M. J. McGirt, G. F. Woodworth, A. Mohammed, K. D. Than, R. J. Tamargo, and R. E. Clatterbuck, “Persistent perioperative hyperglycemia as an independent predictor of poor outcome after aneurysmal subarachnoid hemorrhage,” Journal of Neurosurgery, vol. 107, no. 6, pp. 1080–1085, 2007. View at: Publisher Site | Google Scholar
  17. C. S. Ogilvy, A. C. Cheung, A. P. Mitha, B. L. Hoh, and B. S. Carter, “Outcomes for surgical and endovascular management of intracranial aneurysmsusing a comprehensive grading system,” Neurosurgery, vol. 59, no. 5, pp. 1037–1043, 2006. View at: Publisher Site | Google Scholar
  18. A. A. Rabinstein, J. A. Friedman, D. A. Nichols et al., “Predictors of outcome after endovascular treatment of cerebral vasospasm,” American Journal of Neuroradiology, vol. 25, no. 10, pp. 1778–1782, 2004. View at: Google Scholar
  19. A. J. Rosengart, K. E. Schultheiss, J. Tolentino, and R. L. Macdonald, “Prognostic factors for outcome in patients with aneurysmal subarachnoid hemorrhage,” Stroke, vol. 38, no. 8, pp. 2315–2321, 2007. View at: Publisher Site | Google Scholar
  20. R. E. Temes, E. Tessitore, J. M. Schmidt et al., “Left ventricular dysfunction and cerebral infarction from vasospasm after subarachnoid hemorrhage,” Neurocritical Care, vol. 13, no. 3, pp. 359–365, 2010. View at: Publisher Site | Google Scholar
  21. B. E. Zacharia, A. F. Ducruet, Z. L. Hickman et al., “Renal dysfunction as an independent predictor of outcome after aneurysmal subarachnoid hemorrhage: a single-center cohort study,” Stroke, vol. 40, no. 7, pp. 2375–2381, 2009. View at: Publisher Site | Google Scholar
  22. J. M. Kahn, E. C. Caldwell, S. Deem, D. W. Newell, S. R. Heckbert, and G. D. Rubenfeld, “Acute lung injury in patients with subarachnoid hemorrhage: incidence, risk factors, and outcome,” Critical Care Medicine, vol. 34, no. 1, pp. 196–202, 2006. View at: Publisher Site | Google Scholar
  23. E. J. Roos, G. J. E. Rinkel, B. K. Velthuis, and A. Algra, “The relation between aneurysm size and outcome in patients with subarachnoid hemorrhage,” Neurology, vol. 54, no. 12, pp. 2334–2336, 2000. View at: Google Scholar
  24. K. A. Hanafy, R. Morgan Stuart, L. Fernandez et al., “Cerebral inflammatory response and predictors of admission clinical grade after aneurysmal subarachnoid hemorrhage,” Journal of Clinical Neuroscience, vol. 17, no. 1, pp. 22–25, 2010. View at: Publisher Site | Google Scholar
  25. A. H. Kramer, I. Mikolaenko, N. Deis et al., “Intraventricular hemorrhage volume predicts poor outcomes but not delayed ischemic neurological deficits among patients with ruptured cerebral aneurysms,” Neurosurgery, vol. 67, no. 4, pp. 1044–1053, 2010. View at: Publisher Site | Google Scholar
  26. S. Juvela, “Prehemorrhage risk factors for fatal intracranial aneurysm rupture,” Stroke, vol. 34, no. 8, pp. 1852–1857, 2003. View at: Publisher Site | Google Scholar
  27. C. S. Ogilvy and B. S. Carter, “A proposed comprehensive grading system to predict outcome for surgical management of intracranial aneurysms,” Neurosurgery, vol. 42, no. 5, pp. 959–970, 1998. View at: Publisher Site | Google Scholar

Copyright © 2013 Benjamin W. Y. Lo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.