the presently preferred and alternative embodiments of the invention , including the best mode for practicing the invention known at this time , are now described in detail in connection with the accompanying drawings . it is , however , expressly noted that the present invention is not limited to these embodiments , but rather the intention is that modifications that are apparent to the person skilled in the art and equivalents thereof are also included . in part , the methods and apparatus disclosed herein relate to a bayesian belief network ( hereinafter “ bbn ”) method ( called pharminator ), for calculating the posterior probability that a specific nce will succeed or fail . the success or failure determination can be based on at least on of : ( 1 ) prior data regarding success rates for nces of the same therapeutic class and source , and ( 2 ) the nce &# 39 ; s therapeutic indices , in vitro | in vivo proof of concept data , and proof of concept data in humans from phase i and early phase ii studies . the main distinction between pharminator and previously described methods is that pharminator focuses on evaluating a specific nce and determining if it warrants further consideration in light of certain parameters . other methods have taken more of a population - based analysis approach , yielding valuable data on overall success rates , but not really addressing the needs of drug developers concerned with the termination decision for a single , specific nce . additionally , in various embodiments pharminator includes a graphic user interface representation that facilitates executing the bayesian network approaches and nce data manipulation techniques described herein . other bayesian approaches described in the literature differ from pharminator with respect to the domain to which bayes theorem is applied . published bayesian approaches to pharmacovigilance compare the benefits of bayesian statistics over frequentist approaches and focus on the utilization of bayesian statistics for the analysis of clinical trial data which is in turn used to define “ stopping boundaries ” ( explored in detail below ). bayesian theory has also been used to facilitate drug development - related tasks such as determining clinical trial sample size and designing clinical trials . yet another proposed use for bayes theorem is as an alternate approach to utilizing population pharmacokinetic data to predict toxicity in ongoing clinical trials . these are all clearly different tasks from the aim of this method , which is neither concerned with the long - standing bayesian - frequentist debate , nor with the utilization of bayes theorem to analyze clinical trial results . pharminator is specific to individual nces rather than individual patients or individual studies , and is much more broad in scope in that it attempts to predict safety , efficacy and nce clinical success for a specific nce in question . currently , a review of the literature for decision analytic approaches to pharmacovigilance yields several publications of interest . berry et al adopted a bayesian decision - theoretic approach to determine “ stopping boundaries ” for the development of an nce . their approach utilizes accumulating information on the nce &# 39 ; s performance to determine at which point the clinical trial &# 39 ; s evidence of efficacy is sufficiently negative that the trial should be stopped . the authors argue that if prior data are positive , then one should be willing to tolerate somewhat more negative results in the current clinical trial than if previous evidence is also negative . given that this manuscript was published in 1988 , without the benefit of hindsight of the past 15 years &# 39 ; 70 - 90 % nce failure rate , the authors &# 39 ; argument is representative of a dangerous and costly assumption . namely , that the nce &# 39 ; s clinical trial data can be ignored if it is negative in the context of positive prior data . a counter - argument could demand that the posterior probability distribution should be relied upon to reflect the updated belief that incorporates prior data and the nce &# 39 ; s most current evidence , and that if this posterior probability distribution reveals poor performance , serious consideration should be given to terminating the nce &# 39 ; s development . another difference between berry &# 39 ; s approach and pharminator is that berry &# 39 ; s approach focuses on efficacy in isolation . pharminator utilizes a bayesian belief network to relate safety and efficacy as independent variables , conditional on the common parent ( root ) node , clinical success . the root node &# 39 ; s prior probability distribution is constructed based upon extensive data on nce failure rates stratified by therapeutic class and nce source . berry et al also state the prior distribution is subjective where the subject is the pharmaceutical company . this forces one to predict the posterior probability for a pharmaceutical company rather than for a specific nce . spiegelhalter et al make a cogent argument for the superiority of bayesian over frequentist models for the analysis of clinical trial data . the authors promote the use of bayesian statistics to analyze the outcome of specific ongoing clinical trials . their argument is not relevant to this method . pharminator utilizes the nce &# 39 ; s characteristics within the framework of a bayes network to predict the outcome for a given nce ; clearly a different use of bayes theorem compared to approach advocated by spiegelhalter . similarly to spiegelhalter , johns and andersen describe the utility of predictive probabilities for interim analyses of phase ii and phase iii clinical trials . pharminator focuses on earlier phase decisions so as to avoid costly phase ii and phase iii trials . pallay describes the use of bayes theorem for economically oriented futility analyses of ongoing phase ii clinical trials . this is vastly different from the conditional dependencies incorporated into pharminator , which relates clinical success , efficacy , safety , therapeutic indices and proof of concept data to update prior beliefs pertaining to the nce &# 39 ; s therapeutic class and source . taken together , the persistently high rate of nce failures is the strongest evidence that these previously published and widely adopted approaches do not appear to enable drug developers to be sufficiently accurate in their pharmacovigilance decisions . the present method and apparatus will facilitate improving the efficiency of development of nces . the product of this method is an application ( pharminator ) to be used by the drug development team at the phase i | early phase iia time point for a given nce that has already passed the ind screening process . the user is prompted to answer several key questions about the nce ( detailed below ). the answers provide information regarding which prior probabilities and conditional probability tables are to be used in the method , as well as the observed data upon which prediction will be made . the output is a numeric and graphical ( binomial distribution plot ) report of the prior and posterior probability distributions for clinical success , safety and efficacy . a bayesian network ( often referred to as a bayesian belief network ( bbn ) or a ‘ bayes net ’) is defined as a directed acyclic graph encoding assumptions of conditional independence , with stochastic variables represented as nodes within the network , and inter - variable dependencies represented as inter - nodal links . in addition to a graphical model , a bbn also requires the definition of certain parameters for probabilistic inference purposes . therefore , it is necessary to specify the conditional probability distribution for each node . for distributions with a binary outcome ( i . e . 2 states ), the conditional probability distribution can be represented as a 2 × 2 conditional probability table ( cpt ). these tables specify the probabilities that the node is in state ( 0 , 1 ) given that its parent is in state ( 0 , 1 ). the cpt for the top ( root ) node , which has no parent node , is that root node &# 39 ; s prior probability distribution . assuming conditional independence , and utilizing the chain rule of probability , the joint probability for a network consisting of a root node , r , that has n child nodes , c 1 , c 2 , . . . c n can be calculated : p ( r , c 1 , c 2 , . . . c n )= p ( r )* p ( c 1 | r )* p ( c 2 r ) * . . . p ( c n r ) fig1 is a schematic illustration of a 3 - layer baysian belief network known to the prior art in which the root node ( r ) has 2 child nodes ( c ) and each child node has 2 child nodes ( i . e . that are “ grandchildren ” ( g ) to the root node ), the joint probability for the network can be calculated as in formula 1 below : formula 1 : calculation of the joint probability over a bayesian belief network p ⁡ ( r , c 1 , c 2 , g 1 ⁢ a , g 1 ⁢ b , g 2 ⁢ a , g 2 ⁢ b ) = ⁢ p ⁡ ( r ) * p ⁡ ( c 1 ⁢ ❘ ⁢ r ) * p ⁡ ( c 2 ⁢ ❘ ⁢ r ) * p ⁡ ( g 1 ⁢ a ⁢ ❘ c 1 ) * ⁢ p ⁡ ( g 1 ⁢ b | c 1 ) * p ⁡ ( g 2 ⁢ a ⁢ ❘ ⁢ c 2 ) * p ⁡ ( g 2 ⁢ b ⁢ ❘ ⁢ c 2 ) = ∏ i = 0 n for ⁢ ⁢ ⁢ n ⁢ ⁢ nodes ⁢ ⁢ p ⁡ ( node i ⁢ ❘ ⁢ parent ) the inner - layer nodes for this network are referred to as hidden nodes , and the lowest layer nodes as leaf nodes . calculating the bayesian posterior probability distribution of the root node given that the leaf nodes are in a specified state is a ratio of the sum of joint distributions . for the example network described above , calculating the probability that the root node is in state ‘ f ” ( false ), given that all 4 leaf nodes are in state ‘ t ” ( true ) is achieved as follows : p ⁡ ( r = f ❘ c 1 , c 2 , g 1 ⁢ a = ⁢ t , g 1 ⁢ b = t , g 2 ⁢ a = t , g 2 ⁢ b = t ) ⁢ = ∑ v ∈ { t , f } ⁢ ⁢ ∑ w ∈ { t , f } ⁢ p ⁡ ( r = f , c 1 = v , c 2 = w , g 1 ⁢ a = t , g 1 ⁢ b = t , g 2 ⁢ a = t , g 2 ⁢ b = t ) ⁢ ∑ x ∈ { t , f } ⁢ ⁢ ∑ v ∈ { t , f } ⁢ ⁢ ∑ w ∈ { t , f } ⁢ p ⁡ ( r = x , c 1 = v , c 2 = w , g 1 ⁢ a = t , g 1 ⁢ b = t , g 2 ⁢ a = t , g 2 ⁢ b = t ) each of the joint probability distributions in formula 2 can be calculated utilizing formula 1 . in this way , bayes theorem can be applied to a given bbn . at the core of the pharminator algorithm is a bbn encompassing relevant clinical variables in pharmacovigilance . fig2 a is an illustrative embodiment of a bbn of clinical variables believed to be relevant to clinical success for an nce , with implicit assumptions of conditional independence . the structure of this network is designed to include only those variables believed to be most critical to predicting nce failure , based on the literature , the author &# 39 ; s training , and consultation with drug development experts . the downward direction of the arrows encodes the nce &# 39 ; s deep or hidden knowledge within “ clinical success ” that is manifested as “ safety ” and “ efficacy ”. the knowledge embedded within each of these two hidden nodes is in turn manifested as therapeutic indices or proof - of - concept signals , respectively . this may seem counter - intuitive , since it may seem logical to believe that safety and efficacy cause clinical success ( or failure ) rather than represent a manifestation of an unknown degree of clinical success . however , the goal of pharminator is to predict the clinical success inherent in the nce , rather than to identify causes of clinical success . the nce has an inherent true degree of “ safety ” and “ efficacy ”. a major goal of clinical trials is to determine what these true values are by studying samples and by assuming that the safety or efficacy in the studied samples are accurate estimates of the nce &# 39 ; s true safety and efficacy . the same explanation can be used to justify the links going from safety to the therapeutic index nodes , and from efficacy to the signal nodes . it is this representation that allows pharminator to predict what the nce &# 39 ; s inherent clinical success is based upon the observed therapeutic indices and proof - of - concept signal data . as this is a bbn , the assumption of conditional independence is required . therefore , ti_vital is assumed to be independent of ti_disease ( see fig2 a and section b 2 for definitions ) conditional on the common parent node , safety . the same is true for the in vitro , in vivo and human signal nodes , and their common parent , efficacy . fig2 a is designed to be a “ best guess ” representation of those variables deemed most important in predicting clinical success . other embodiments of pharminator will encode a bbn with slightly different leaf nodes , based upon accumulated data on nce failure specific to therapeutic classes . also absent from the bbn is a representation of idiosyncratic severe adverse events which are unpredictable , by definition . other embodiments of pharminator include pharmacogenomic markers of adverse events and drug resistance , as depicted in fig2 b . fig2 b is an illustrative embodiment of a bbn of clinical variables believed to be relevant to clinical success for an nce , with pharmacogenomics added as a leaf variable under efficacy and safety . as used herein , the term “ clinical success ” refers to an nce that is still on the market 1 year after nda approval . as used herein , the term “ efficacy ” refers to an nce that produces a “ sufficient ” degree of change in a surrogate or true marker , compared to control ( placebo or current gold standard therapy ). as used herein , the term “ sufficient ” degree of change depends on ( 1 ) the clinical indication and ( 2 ) the development phase ( during phase ii , the signal need not be statistically significant , while phase iii studies must show statistical significance in at least two separate trials ). certain indications may require only modest effect from an nce in order to be successful ( e . g . acute , relatively benign disorders ), while others may require extreme effects ( e . g . life - saving therapies ). as used herein , the term “ life - saving ” refers to an nce the disease for which it is indicated is fatal and , there are no alternate life - extending therapies . the pharminator utilizes the life - saving status of the nce to determine the influence that the safety data will have on the calculation of the posterior probability of clinical success . the assumption is that a higher degree of toxicity is tolerated for an nce that is truly life - saving , as defined above , thereby making the probability of clinical success largely dependent on efficacy . for example , an nce that truly extends life expectancy but causes acute renal failure may still have a high probability of success because it is assumed that the initiation of dialysis is preferable to death . this assumption is open to argument from the point of view of quality of life issues , since truly curative therapies for lethal disorders are rare , however assuming that supportive therapies are preferable to death is reasonable . as used herein , the term “ marker ” refers to an indicator of response to therapy . a surrogate marker is a marker that is not directly or primarily involved in the pathogenesis of the disease , whereas a true marker is primarily integral to the disease mechanism . an example of a surrogate marker is the cd 4 count in hiv . a true marker for hiv is viral load . as used herein , the term “ acquired nce ” refers to an nce that a pharmaceutical company has licensed - in from another company , such as a biotechnology firm , or that has been acquired from some other source . as used herein , the term “ self - originated nce ” refers tp an nce for which the initial pre - clinical development occurred within the same pharmaceutical company that will assume responsibility for conducting clinical trials . the prior probability of success differs significantly between acquired and self - originated nces . not surprisingly , nces that have undergone initial clinical testing abroad ( and demonstrate potential effectiveness in humans ) are more likely to succeed . pharminator gives the user the option of selecting whether the prior probability should be optimistic or pessimistic . as described in detail in “ prior probabilities and conditional probability tables ”, the nce &# 39 ; s intended therapeutic class affects the selection of the nce &# 39 ; s prior probability of clinical success . prior data on nce success rates are stratified by therapeutic source , and include the total number of nces within each therapeutic class , the fraction of the total nces that have failed , and the fraction of the total nces that are still under development , for a total of 671 nces spanning ind filing dates from 1981 to 1992 . however , it should be noted that if an nce is sufficiently safe and effective , the nce &# 39 ; s posterior probability of success will be high , regardless of the prior probability . the tufts center for the study of drug development &# 39 ; s ( tcsdd ) published reports provide the prior probability of failure based on the current failure rate , as well as the probability of failure assuming that all nces still under development are successful ( i . e . a more optimistic prior probability of failure ). when the user sets pharminator &# 39 ; s prior bias to “ pessimistic ” ( the default setting ), the former prior probability is used — i . e . the prior probability of failure based on the current failure rate . when the user sets pharminator &# 39 ; s prior bias to “ optimistic ”, the latter prior probability is used — i . e . the prior probability of failure based on the assumption that all nces still under development will not fail . as used herein , the term “ safety ,” when referring to an nce , is essentially synonymous with toxicity . every nce that is not an inert placebo has some degree of “ toxicity ” in that even the desired effects of an nce become toxic if a large enough dose is given . an nce &# 39 ; s “ safety ” is therefore defined as a degree of toxicity that is an acceptable balance against the benefit to the patient . ultimately , an nce can only be deemed safe once it has undergone phase iv post - marketing surveillance . prior to phase iv , insufficient numbers of patients have received the nce such that rare but severe idiosyncratic reactions would not likely be detected . as used herein , the term “ therapeutic class ” refers to the organ system affected by the disease process for which the nce is indicated . this definition of therapeutic class is utilized rather than the more traditional chemical class because prior probabilities of success are known for a total of 671 nces , stratified by therapeutic class , and stratifying by chemical class would fractionate the data beyond use with too many categories and too few nces in each category . the therapeutic classes included in pharminator are : analgesic / anesthetic , antimicrobial , antineoplastic , cardiovascular , central nervous system ( cns ), endocrine , gastrointestinal ( gi ), immunologic , respiratory , and miscellaneous . clearly these are less specific categories than those used by clinicians ( e . g . calcium channel blockers , ace inhibitors , mono - amine oxidase inhibitors etc .) however , as stated above , there do not appear to be sufficient data to allow for a more specific stratification without fractionating the data beyond utility . as used herein , the term “ therapeutic index ” ( ti ) refers to the ratio of the nce dose that produces an undesired effect to the nce dose that produces the desired effect in a proportion of the study population . the numerator is the td x ( toxic dose in x % of the population ) and the denominator is the ec y ( effective dose in y % of the population ). each nce has several therapeutic indices , depending on the number of specific adverse events ( e . g . the ti for hepatotoxicity is different from the ti for nephrotoxicity ), the number of specific desired effects ( e . g . ace inhibitors reduce blood pressure and reduce proteinuria ), and depending on the definition of the proportion of the study population ( i . e . the values of x and y ). a larger ti represents a generally safer nce . a smaller ti will be either too unsafe to be used clinically , or will require very close therapeutic drug monitoring in order to ensure safety ( e . g . digoxin ). the pharminator is designed to be inherently pessimistic , given the high rate of nce failures to date , and the extreme costs associated with these failures . therefore , the current implementation of pharminator requires input for two specific tis : the lowest ti for an undesired effect on any vital organ ( brain , heart , lungs , liver , kidney , exocrine pancreas , bone marrow ), and the lowest ti for an undesired effect on any organ or system that is already adversely affected by the disease / system for which the nce is indicated . an example of the latter is retinal toxicity caused by an nce indicated for the treatment of diabetes mellitus . in other embodiments of the pharminator , additional or alternate ti variables are added to the method . one aspect of the invention recognizes the reasons why therapeutic class and nce source are included in the method as “ prior probability modifiers ” and not as stochastic variables ( nodes ) in fig2 a . as discussed , “ network structure and rationale ”, safety and efficacy do not “ cause ” clinical success in the pharminator method . the safety and efficacy nodes ( and indeed all child nodes in the method ) are in fact manifestations of the inherent degree of clinical success of the nce . contrary to this , therapeutic class and nce source have a direct impact on nce clinical success . additionally , while the nce &# 39 ; s true safety , therapeutic indices , efficacy and proof - of - concept signals are not known , the nce &# 39 ; s intended therapeutic class and source are known with certainty . it is therefore nonsensical to represent therapeutic class and nce source as stochastic variables . the approaches disclosed herein make use of this distinction . however , fig3 is provided below as an adjunct to fig2 a , in order to demonstrate the dependencies between therapeutic class , nce source and clinical success , and to demonstrate that the noisy - or assumption can be utilized to calculate the prior probability of clinical success from the prior data on nce failure rates stratified by therapeutic class and nce source . fig3 is a schematic illustration of how therapeutic class and nce source relate to clinical success ( see also , fig2 a ). referring to fig3 , the direction of the arrows indicates that therapeutic class and nce source have a causal effect on clinical success . determining p ( clinical success | therapeutic class , nce source ) is therefore not a bayesian posterior probability and the noisy - or assumption is applicable . tcsdd publications are the most extensive , accessible , and reliable source of the prior probability of nce success ( and failure ). dimasi recently analyzed the causes of failure and reported the success rates for 671 nces for which inds were filed between 1981 and 1992 . in his report , he provides “ current and maximum possible success rates ” stratified by therapeutic class for 503 self - originated nces . the “ current success rate ” is the fraction of the number of nces ( in that class ) that have been successful over all nces in that class . this is in fact a pessimistic prior because the implicit assumption is that all open nces ( i . e . nces still in development ) will fail . the “ maximum possible success rate ” is the success rate assuming that “ all open nces will eventually be approved ”— an optimistic assumption . dimasi also provides probabilities of nce success stratified by nce source . pharminator asks the user to indicate the nce &# 39 ; s therapeutic class , nce source , as well as the user &# 39 ; s desired “ prior bias ”, which may be either pessimistic or optimistic . the prior bias determines which therapeutic class prior probability is utilized : if the user selects “ pessimistic ” ( the default setting ), the current success rate is used to calculate the nce &# 39 ; s prior probability of success . conversely , if the user selects “ optimistic ”, the maximum possible success rate is used . although dimasi &# 39 ; s analysis is quite informative , he did not sub - stratify by therapeutic class and nce source combinations . in order to allow pharminator to choose a prior probability that most accurately reflects the nce &# 39 ; s therapeutic class and its source , the algorithm utilizes the noisy - or assumption . paraphrasing szolovits , the noisy - or assumption states that the probability that some set of variables causes an outcome equals the probability that at least one of the variables does so . the probability of interest is p ( clinical success | therapeutic class , nce source ). for an explanation as to why this is not a posterior probability distribution , see section b2 and fig3 . given the noisy - or assumption , the probability of interest can be calculated : therefore , p ( clinical success | therapeutic class , source )= 1 −[( 1 − p ( clinical success | therapeutic class ))*( 1 − p ( clinical success | source ))]≈ the prior probability , p ( clinical success ) for the nce in question and p ( clinical failure )= 1 − p ( clinical success ) formula 3 : noisy - or pharminator utilizes the prior bias selected by the user to determine which prior probability of success to utilize for the selected therapeutic class , then uses this probability along with the probability of success for the selected nce source to calculate p ( clinical success | therapeutic class , source ), given the noisy - or assumption as in formula 3 above . this calculated probability is utilized as the prior probability of clinical success for the nce in question . fig4 is an illustrative flowchart of one embodiment of the invention wherein the prior probability of clinical success is determined . as described in section b1 , a bbn requires cpts for each node in order to be utilized for probabilistic inference . although the data provided in the tcsdd reports is valuable , the format of those reports is not directly applicable to the construction of bbn conditional probability tables . the main limitation of pharminator is the absence of appropriate conditional probability data and the need to make certain assumptions to allow utilization of the data from the tcsdd published reports . other embodiments of pharminator focus on the task of obtaining appropriate data to populate the cpts . these assumptions will be tested by sensitivity analyses once characteristics and outcomes for specific successful and failed nces become available . in the absence of such data , the methods by which the cpts are currently constructed are described in this section . dimasi analyzed the causes of failure for 348 nces that were withdrawn from development . it should be noted that nces that proceeded through all clinical trial phases but failed to achieve nda approval are not included in dimasi &# 39 ; s analysis . as well , dimasi stratified the causes of failure by “ primary ” cause , thereby not disclosing any degree of overlap — i . e . nces that failed primarily due to one reason , but may have also failed for another reason ( e . g . an nce that failed because it was not safe , but was also not very effective ). his analysis demonstrated that of a total of 348 nces that were terminated , the primary reason for termination was efficacy in 121 , safety in 72 , economics in 109 , and “ other ” in 46 . since pharminator is concerned only with safety and efficacy , the probability that safety is the primary cause of failure is 72 /( 72 + 121 )= 0 . 37 , and the probability that efficacy is the primary cause of failure is 121 /( 72 + 121 )= 0 . 63 . assuming that the proportions of causes of failure are consistent across the withdrawn drugs , the cpt probabilities , p ( safety = f | clinical success = f ) and p ( efficacy = f | clinical success = f ), can be calculated by an “ overlap ” function , as follows : p ( safety = f | clinical success = f )= total number of primary safety failures +( total number of primary safety failures * proportion of primary efficacy failures )= 72 +( 72 * 121 / 193 )= 0 . 606 formula 4 : “ overlap ” function p ( efficacy = f | clinical success = f )= total number of primary efficacy failures +( total number of primary efficacy failures * proportion of primary safety failures )= 121 +( 121 * 72 / 193 )= 0 . 860 these values ( and their respective complement values ) occupy the first rows of their respective cpts . while the overlap function permits estimation of the first row of each of the safety and efficacy cpts ( i . e . p ( node = f | parent = f ) and p ( node = t | parent = f )), currently , there are no adequate , available data for the second rows of the safety and efficacy cpts ( p ( node = f | parent = t ) and p ( node = t | parent = t )). for now , these values are currently set as pessimistic estimates . for the efficacy cpt , the probability p ( efficacy = f | clinical success = t ), i . e . the probability that an nce is not efficacious given that it is clinically successful , is logically estimated to be very low . until data for sensitivity analyses become available , this value is set at 0 . 01 , and its complement , p ( efficacy = t | clinical success = t ) is therefore 1 - 0 . 01 = 0 . 99 . for the safety cpt , the probability p ( safety = f | clinical success = t ), i . e . the probability that an nce is not safe given that it is clinically successful is also estimated to be low , but likely not as low as for p ( efficacy = f | clinical success = t ). until data for sensitivity analyses become available , this value is set at 0 . 05 for nces that are not life - saving , and for nces that are life - saving , this value is set at 0 . 5 . this difference is to reduce the influence that tis have on the posterior probability of clinical success for life - saving nces . fig5 is a schematic illustration of the construction of the hidden nodes &# 39 ; cpts . fig5 is a schematic illustration of an algorithm demonstrating one embodiment of the invention wherein the overlap function ( formula 4 ) and the life - saving preference are utilized to determine the cpts for the safety and efficacy nodes . just as for the hidden nodes ( safety and efficacy ), there are no easily accessible data on therapeutic indices and proof of concept signal data for nces that have failed . however , devising models that approximate these relationships is a somewhat less arduous task than for safety and efficacy . with respect to the relationship between ti and safety , it is known that ti is directly proportional to the degree of safety because a larger ti simply means more prescribing “ room ” between the effective dose and the toxic dose . most prescription medications have tis that are in the 8 - 10 range . given that ti is a ratio , and that the lowest rational value for a ti is 1 , the relationship between ti and safety can be approximated by a logistic sigmoid model ( formula 5 , fig6 ): the slope and intercept of this model were selected to reflect what is believed to be an accurate approximation of the relationship between ti and p ( ti | safety ). fig6 a is a graph showing the logistic sigmoid functions that are used to approximate the p ( ti | safety = t ) cpt values from the ti values . fig6 b is a graph showing the logistic sigmoid functions that are used to approximate the p ( ti | safety = f ) cpt values from the ti values . a similar assumption is made for the proof - of - concept signal data in that the quantity of the signal is proportional to the degree of efficacy . in contrast to the sigmoid relationship between ti and safety , the relationship between proof - of - concept signal and efficacy is assumed to be a simple linear function ( y = mx + b ). however , the signal data must first be transformed to a standardized measure so that different ranges and scales will not influence the interpretation of the signal . for this purpose a modified signal - to - noise ratio is used ( formula 6 ), requiring the user to enter the mean and variance ( s . d . 2 ) for the control and experimental groups , for each experimental environment ( in vitro , in vivo — highest - order species , human ), as well as whether each signal is a true or surrogate marker . this formula provides a variance - corrected measure of the degree of signal as a value between 0 and 1 . the resultant signal - to - noise ratio value is utilized by the linear models to estimate p ( signal | efficacy ) for in vitro , in vivo and human signals , stratified by true and surrogate markers . the slopes ( and intercepts ) of the linear models are adjusted to reflect differences between in vitro , in vivo and human signals , and between true and surrogate markers ( fig7 ). specifically , the slope of the function is proportionate to the environment order ( in vitro & lt ; in vivo & lt ; human ), and the marker ( surrogate & lt ; true ). fig7 is a group of linear functions demonstrating how the modified signal : noise ratio value is utilized to approximate p ( signal | efficacy ) for a series of randomly - generated nce and control means & amp ; variances . the linear functions are stratified by experimental environment ( in vitro , in vivo , human ), and type of marker ( surrogate vs . true marker ). the slopes and intercepts are adjusted to reflect what is believed to be an adequate approximation . fig8 is an illustrative embodiment of an algorithm for constructing the leaf node cpts . the default state for all leaf nodes is “ true ”, and it is the leaf node cpts that change in response to the ti and signal data entered by the user . this represents one major departure from bbn methodology : pharminator &# 39 ; s use of the input data to determine the specific cpt to be used for a given set of fixed leaf node states . the acquisition of data on nces that have failed will facilitate the modification of the linear and sigmoid functions rather than direct changes to the leaf node cpts . justification for this approach is that this can actually work to the advantage of each drug development institution that utilizes pharminator . many pharmaceutical companies develop medications in a small number of therapeutic classes and therefore have nce failure data that is highly specific to that pharmaceutical company &# 39 ; s future development projects . therefore the use of these data to modify pharminator &# 39 ; s cpt functions will result in a company - specific implementation of pharminator , the predictive accuracy of which will be directly proportionate to the specific company &# 39 ; s development history and prior investments in nce failures ( i . e . accuracy proportionate to their losses ). smaller companies with little or no development history will not have the ability to implement a company - specific implementation , but will benefit from the prior knowledge of the entire industry , excluding confidential and privileged information from other companies . data from a larger pharmaceutical company will remain exclusive to that specific pharmaceutical company unless that company agrees to allow pharminator to utilize their data for the benefit of the entire industry , always maintaining confidentiality regarding specific nces that have failed . pharminator requires the following input from the user : nce “ demographics ” including nce name , therapeutic class , source , and life - saving status ; the user &# 39 ; s prior bias preference ( pessimistic vs . optimistic ; default is pessimistic ); signal data for each of in vitro , in vivo ( highest - order species ) and human ; the minimum ti_vital ; and the minimum ti_disease . the signal data are entered for the maximal dose given , regardless of toxicity . therefore , toxicity ( safety ) is not taken into account for signal data because safety and efficacy are assumed to be independent variables , conditional on the common parent , clinical success . the signal nodes &# 39 ; inputs include nce mean and variance , control mean and variance , and the type of marker : true or surrogate . the current implementation of pharminator requires all of these values to be entered in order for the posterior probability to be calculated . other embodiments of pharminator perform probability inference on partial nets ( i . e . nets that are missing one or more leaf node variables ). with this information , pharminator selects the appropriate prior probability of clinical success , and calculates the posterior probability distribution for clinical success , safety , and efficacy . the hidden node “ prior ” probabilities are required in order to calculate hidden node posterior probabilities . these “ prior ” probabilities are calculated from the hidden and root nodes ° cpts : p ( hidden node )= σ [ p ( hidden node | parent node )* p ( parent node )] formula 7 : calculating a hidden node &# 39 ; s “ prior ” probability the prior and posterior probability distributions are displayed graphically as binomial distributions . the “ n ” for the clinical success prior probability distribution (“ prior n ”) is the total number of nces from which the prior data were attained . the “ n ” for the clinical success posterior probability distribution ( post n ) is ( prior n + 1 ). this is likewise for the prior n and post n for the safety and efficacy probability distributions ( see section “ implementation and examples ” for a pictorial demonstration ). conglomerating the discussion of the previous sections , including fig4 and 8 and formula 2 results in the algorithm shown in fig9 . fig9 is a schematic illustration of one embodiment of the invention wherein an overview algorithm is used in conjunction with the invention , combining the components as described . pharminator is implemented in java 1 . 4 , using apple projectbuilder v2 . 1 , on apple os x . 2 jaguar . an object - oriented , model - view approach was utilized to structure to the program . the accuracy of the bbn was validated against bayesware discoverer ® ( http :// bayesware . com ). the methods and systems described herein can be performed in software on general purpose computers , servers , or other processors , with appropriate magnetic , optical or other storage that is part of the computer or server or connected thereto , such as with a bus . the processes can also be carried out in whole or in part in a combination of hardware and software , such as with application specific integrated circuits . the software can be stored in one or more computers , servers , or other appropriate devices , and can also be kept on a removable storage media , such as a magnetic or optical disks . curonc is a fictional anti - neoplastic agent devised solely for the purpose of illustrating some key features of pharminator . curonc is self - originated in the usa and meets the definition of “ life - saving ”. the signal inputs , ti inputs , and clinical success probability distribution plots are shown in fig1 . the safety and efficacy probability distribution plots are shown in fig1 . the effect of changing the life - saving option to “ not life - saving ” is shown in fig1 . the effect of changing the prior bias to optimistic is shown in fig1 . overall , the probability distributions generated by pharminator suggest that curonc has a high probability of efficaciousness ( 0 . 7872 ), but is also very likely to have significant toxicity ( p ( safety = t )= 0 . 0645 ). therefore , if curonc is indeed “ life - saving ”, it has a probability of clinical success of 0 . 4951 with little overlap between the clinical success prior ( 0 . 2304 ) and posterior probability distributions . however , if curonc is not truly life - saving , its probability of clinical success is 0 . 2016 ( less than the prior probability ) when prior bias is pessimistic , and at best ( prior bias = optimistic ), the probability of clinical success is 0 . 3359 , which is still less than the prior probability . given these results , development of curonc should be continued only if it is deemed to be truly life - saving . fig1 is a screen shot of one embodiment according to the invention illustrating the prior and posterior probability distributions for clinical success for the fictional antineoplastic agent , curonc . fig1 is a screen shot of one embodiment according to the invention illustrating the prior and posterior probability distributions for safety and efficacy for the fictional antineoplactic agent , cur onc . fig1 is a screen shot of one embodiment according to the invention illustrating the effect of selecting the “ not life saving ” the posterior probability distribution for clinical success for the fictional antineoplastic agent , curonc . fig1 is a screen shot of one embodiment according to the invention illustrating the effect of setting the prior bias to “ optimistic ” on the prior and posterior probability distribution for clinical success for the fictional antineoplastic agent , curonc . recombinant human activated protein c ( rhapc ) is a relatively novel agent that is known for its anti - coagulant , pro - fibrinolytic , and anti - inflammatory properties . eli lilly ™ research laboratories has developed ly203638 ( rhapc ) as a novel therapy for sepsis ( clinical investigator &# 39 ; s brochure kindly provided by dr . robert rubin ). in general , this example is limited in that several unpublished pre - clinical efficacy studies are listed in the clinical investigator &# 39 ; s brochure , but no data are accessible . the most relevant in vitro study was used . this in vitro study was performed prior to the go / no - go decision time point . bajzar et al reported dose - dependent lysis times , but did not include any measures of variability . therefore , in vitro variance is set to 0 for both the nce and control ( the in vitro variance entries are actually set to 0 . 000001 because the program &# 39 ; s current implementation will not calculate posterior probabilities if any value is 0 . this minor problem will be resolved with future implementations ). published in vivo studies performed prior to the go / no - go decision time point evaluated pre - clinical efficacy in primates , canines , guinea pigs , and rats . the primate data are used for this example because primates are the highest - order species studied . early phase ii study data in humans were provided in the clinical investigator &# 39 ; s brochure . this example is an approximation based upon the accessible information only . note that the true outcome marker is successful treatment of sepsis . the phase ii endpoints reported are therefore all surrogate markers : organ failure - free days , number of transfusion requirements , icu -, hospital -, and ventilator - free days , and 28 - day all - cause mortality . for the purpose of this example , organ failure - free days ( shock ) was chosen as a good sepsis - specific surrogate marker in that multi - organ failure and sepsis - related morbidity are very tightly correlated . no specific data on therapeutic indices could be found either in the clinical investigator &# 39 ; s brochure or in the literature from the go / no - go decision time point . however , the clinical investigator &# 39 ; s brochure contains data from phase i studies at doses ranging from 12 - 48 μg / kg / hour suggesting that the ti is at least 4 ( 48 / 12 ). toxicology studies in primates demonstrated that the “ no - observed - adverse - effect level ” was 2 mg / m 2 / hour with toxic effects observed at a dose of 8 mg / m 2 / hour . taken together , these data suggest that the ti is approximately 4 . in the absence of more accurate ti data , this value is used in the example . ly203638 is classified as a cardiovascular agent ( since there are no prior data for hematologic agents and it is not antimicrobial ). the limitations in acquiring appropriate data for ly203638 underscores the requirement to have unfettered access to the nce &# 39 ; s data in order to optimize pharminator &# 39 ; s predictive accuracy . fig1 is a screen shot of one embodiment according to the invention illustrating the prior and posterior probability distributions for clinical success for ly203638 demonstrating that it has a very low probability of clinical success of 0 . 0521 , much lower than the prior probability of 0 . 246 . fig1 is a screen shot of one embodiment according to the invention illustrating the prior and posterior probability distributions for safety and efficacy for ly203638 ( rhapc ) demonstrating that the probabilities of safety and efficacy are both very low ( 0 . 0211 and 0 . 0667 , respectively ). even when the prior bias is set to optimistic ( fig1 ), the probability of clinical success is only 0 . 0526 , also much lower than the optimistic prior probability of 0 . 2916 . fig1 is a screen shot of one embodiment according to the invention illustrating the effect of setting the prior bias to “ optimistic ” on the prior and posterior probability distribution for clinical success for ly203638 ( rhapc ). therefore , based only on data available prior to later phase ii studies : even as a life - saving nce , and when assuming an optimistic prior probability of success , ly203638 has a very low probability of clinical success based only on data available prior to phase iii studies . of interest , after ly203638 received nda approval , subsequent post - approval studies raised several concerns about ly203638 &# 39 ; s safety and efficacy , calling for phase iv studies to be performed . it should be appreciated that various aspects of the claimed invention are directed to portions of the systems described , the methods and the processes of the pharminator embodiments disclosed herein . further , the terms and expressions employed herein are used as terms of description and not of limitation , and there is no intention , in the use of such terms and expressions , of excluding any equivalents of the features shown and described or portions thereof , but it is recognized that various modifications are possible within the scope of the invention claimed . accordingly , what is desired to be secured by letters patent is the invention as defined and differentiated in the following claims , including all equivalents .