Hypertabastic survival models
Hypertabastic survival models were introduced in 2007 by Mohammad Tabatabai, Zoran Bursac, David Williams, and Karan Singh. This distribution can be used to analyze time-to-event data in biomedical and public health areas and normally called survival analysis. In engineering, the time-to event analysis is referred to as reliability theory and in business and economics it is called duration analysis. Other fields may use different names for the same analysis. These survival models are applicable in many fields such as biomedical, behavioral science, social science, statistics, medicine, bioinformatics, medicalinformatics, data science especially in machine learning , computational biology, business economics, engineering, and commercial entities. They not only look at the time to event, but whether or not the event occurred. These time-to-event models can be applied in a variety of applications for instance, time after diagnosis of cancer until death, comparison of individualized treatment with standard care in cancer research, time until an individual defaults on loans, relapsed time for drug and smoking cessation, time until property sold after being put on the market, time until an individual upgrades to a new phone, time until job relocation, time until bones receive microscopic fractures when  undergoing different stress levels, time from marriage until divorce, time until infection due to catheter, and time from bridge completion until first repair. Hypertabastic cumulative distribution function. The Hypertabastic cumulative distribution function or simply the hypertabastic distribution function formula_1 is defined as the probability that random variable formula_2 will take a value less than or equal to formula_3. where formula_5 represents the hyperbolic secant function and formula_6, formula_7 are parameters. The parameters formula_6 and formula_7 are both positive with formula_5 and formula_11 as hyperbolic secant and hyperbolic cotangent respectively. The Hypertabastic probability density function is where formula_13 and formula_14 are  hyperbolic cosecant and hyperbolic tangent respectively and Hypertabastic survival function. The Hypertabastic survival function is defined as where formula_17 is the probability that waiting time exceeds formula_3. For formula_19, the Restricted Expected (mean) Survival Time of the random variable formula_2 is denoted by formula_21, and is defined as Hypertabastic hazard function. For the continuous random variable formula_2 representing time to event, the Hypertabastic hazard function formula_24, which represents the instantaneous failure rate at time formula_3 given survival up to time formula_3, is defined as The Hypertabastic hazard function has the flexibility to model varieties of hazard shapes. These different hazard  shapes could apply to different mechanisms for which the hazard functions may not agree with conventional models. The following is a list of possible shapes for the Hypertabastic hazard function: For formula_28, the Hypertabastic hazard function is monotonically decreasing indicating higher likelihood of failure at early times. For formula_29, the Hypertabastic hazard curve first increases with time until it reaches its maximum failure rate and thereafter the failure decreases with time (unimodal). For formula_30, the Hypertabastic hazard function initially increases with time, then it reaches its horizontal asymptote formula_6. For formula_32, the Hypertabastic hazard function first increases with time with  an upward concavity until it reaches its inflection point and subsequently continues to increase with a downward concavity. For formula_33, the Hypertabastic hazard function initially increases with an upward concavity until it reaches its point of inflection, thereafter becoming a linear asymptote with slope formula_6. For formula_35, the Hypertabastic hazard function increases with an upward concavity. The Hypertabastic cumulative hazard function is: Hypertabastic proportional hazards model. The hazard function formula_37 of the Hypertabastic proportional hazards model has the form where formula_39 is a p-dimensional vector of explanatory variables and formula_40 is a vector of unknown parameters. The combined effect  of explanatory variables formula_41 is a non-negative function of formula_42 with formula_43. The Hypertabastic survival function formula_44 for the proportional hazards model is defined as: and the Hypertabastic probability density function for the proportional hazard model is given by: Depending on the type of censoring, the maximum likelihood function technique along with an appropriate log-likelihood function may be used to estimate the model parameters. If the sample consists of right censored data and the model to use is Hypertabastic proportional hazards model, then, the proportional hazards log-likelihood function is: Hypertabastic accelerated failure time model. When the covariates  act multiplicatively on the time-scale, the model is called accelerated failure time model. The Hypertabastic survival function for the accelerated failure time model is given by: The Hypertabastic accelerated failure time model has a hazard function formula_37 of the form The Hypertabastic probability density function for the accelerated failure time model is: For the right censored data, the log-likelihood function for the Hypertabastic accelerated failure time model is given by where formula_53. A modified chi-squared type test, known as Nikulin-Rao-Robson statistic is used to test the goodness-of-fit for Hypertabastic accelerated failure time models  and its comparison with unimodal hazard rate functions. Simulation studies have shown that the Hypertabastic distribution can be used as an alternative to log-logistic and log-normal distribution because of its flexible shape of hazard functions. The Hypertabastic distribution is a competitor for statistical modeling when compared with Birnbaum-Saunders and inverse Gaussian distributions Likelihood functions for survival analysis. Consider a sample of survival times of n individuals formula_54 with associated p-dimensional covariate vectors formula_55 and an unknown parameter vector formula_56. Let formula_57 and formula_58 stand for the corresponding probability density function, cumulative distribution function, survival function and  hazard function respectively. and the log-likelihood formula_60 is or equivalently, and the log-likelihood function is or equivalently, where In the presence of left censored data, the likelihood function is and the corresponding log-likelihood function is where In the presence of interval censored data, the likelihood function is and the log-likelihood function is where formula_72 for all interval censored observations and If the intended sample consists of all types of censored data (right censored, left censored and interval censored), then its likelihood function takes the following form and its corresponding log-likelihood function is given by Applications  of hypertabastic survival models. Oil field quantities. The quantities of 49 locations of the same area of an oil field was examined to identify its underlying distribution. Using generalized chi-squared, the distribution of oil field quantities was represented by the Hyperbolastic distribution and compared with the lognormal (LN), log-logistic (LL), Birnbaum-Saunders (BS) and inverse Gaussian (IG) distributions. Remission duration for acute leukemia. The times of remission from clinical trial for acute leukemia of children study were used to analyze the remission duration of acute leukemia data for two groups of patients controlling for log of white blood  cell counts. The Hypertabastic accelerated failure time model was used to analyze the remission duration of acute leukemia patient. Brain tumor study of malignant glioma patients. A randomized clinical trial comparing two chemotherapy regimens for 447 individuals with malignant glioma. A total of 293 patients died within a five-year time period and the median survival time was about 11 months. The overall model fit, in comparison with other parametric distributions, was performed using the generalized chi-square test statistics and proportional hazards model. Analysis of breast cancer patients. The Hypertabastic proportional hazard model was used to analyze numerous breast  cancer data including the survival of breast cancer patients by exploring the role of a metastasis variable in combination with clinical and gene expression variables. Analysis of hypertensive patients. One hundred five Nigerian patients who were diagnosed with hypertension from January 2013 to July 2018 were included in this study, where death was the event of interest. Six parametric models such as; exponential, Weibull, lognormal, log-logistic, Gompertz and Hypertabastic distributions were fitted to the data using goodness of fit tests such as S.E., AIC, and BIC to determine the best fit model. The parametric models were considered because they  are all lifetime distributions. S.E., AIC, and BIC measures were used to compare these parametric models. Analysis of cortical bone fracture. Stress fractures in older individuals are very important due to the growing number of elderly. Fatigue tests on 23 female bone samples from three individuals were analyzed. Hypertabastic survival and hazard functions of the normalized stress level and age were developed using previously published bone fatigue stress data. The event of interest was the number of cycles until the bone gets microscopic fracture. Furthermore, Hypertabastic proportional hazard models were used to investigate tensile fatigue and cycle-to-fatigue for  cortical bone data. Analysis of kidney carcinoma patients. Using National Cancer Institute data from 1975 to 2016, the impact of histological subtypes on the survival probability of 134,150 kidney carcinoma patients were examined. The study variables were a race/ethnicity, age, sex, tumor grade, type of surgery, geographical location of patient and stage of disease. The Hypertabastic proportional hazards model was used to analyze the survival time of patients diagnosed with kidney carcinoma to explore the effect of histological subtypes on their survival probability and assess the relationship between the histological subtypes, tumor stage, tumor grade, and type of surgery . Kidney carcinoma SAS example code. Sample code in SAS: Proc nlp data=sasuser.KidneyCarcinoma tech=quanew cov=2 vardef=n pcov phes maxiter=250; run; Applications of hypertabastic survival models in bridge engineering. Although survival analysis tools and techniques have been widely used in medical and biomedical applications over the last few decades, their applications to engineering problems have been more sporadic and limited. The probabilistic assessment of service life of a wide variety of engineering systems, from small mechanical components to large bridge structures, can substantially benefit from the well-established survival analysis techniques. Modeling of time-to -event phenomena in engineering applications can be performed under the influence of numerical and categorical covariates using observational or test data. The "survival" of an engineering component or system is synonymous with the more commonly used term "reliability". The term "hazard rate" or "conditional failure rate" (defined as probability of survival per unit time assuming survival up to that time) is an important measure of the change in the rate of failure over time. In this context, failure is defined as reaching the target event in the time-to-event process. This could be defined as reaching a particular serviceability  condition state, localized/partial structural failure, or global/catastrophic failure applied the Hypertabastic parametric accelerated failure time survival model to develop probabilistic models of bridge deck service life for Wisconsin. Bridge decks are typically concrete slabs on which traffic rides as seen in the Marquette Interchange bridge. The authors used the National Bridge Inventory (NBI) dataset to obtain the needed data for their study. NBI records include discrete numerical ratings for bridge decks (and other bridge components) as well as other basic information such as Average Daily Traffic (ADT) and deck surface area (obtained by multiplying the provided bridge length  with bridge deck width). The numerical ratings range from 0 to 9 with 9 corresponding to brand new condition and 0 being complete failure. A deck condition rating of 5 was selected as the effective end of service life of bridge deck. The numerical covariates used were the ADT and deck surface area, while the categorical covariate was the superstructure material (structural steel or concrete). The hypertabastic Proportional Hazards and Accelerated Failure Time models are useful techniques in analyzing bridge-related structures due to its fleibility of hazard curves, which can be monotonically increasing or decreasing with upward or downward  concavity. It can also take the shape of a single mound curve. This flexibility in modeling various hazard shapes makes the model suitable for a wide variety of engineering problems. Tabatabai et al. extended the Hypertabastic bridge deck models developed for Wisconsin bridges to bridges in six northern US states and then to all 50 US states. The study of bridge decks in all 50 states indicated important differences in reliability of bridge decks in different states and regions. Stevens et al. further extended the use of Hypertabastic survival models to bridge superstructures. The covariates used were ADT, maximum bridge  span length and superstructure type. The survival function can be used to determine the expected life using the following equation (area under the entire survival curve): It is important to note that both the survival function and the expected life would change as the time passes by. The conditional survival function formula_77 is a function of time formula_3 and survival time formula_79 and is defined as Nabizadeh et al. used the Hypertabastic survival functions developed for Wisconsin to analyze conditional survival functions and conditional expected service lives formula_81: The conditional expected life would continue to increase as the survival time  formula_79 increases. Nabizadeh et al. term this additional expected life as "survival dividend.” An important mode of failure in bridge engineering is metal fatigue, which can result from repetitive applications of stress cycles to various details and connections in the structure. As the number of cycles formula_84 increase, the probability of fatigue failure increases. An important factor in fatigue life formula_85 is the stress range (Sr)(maximum minus minimum stress in a cycle). The probabilistic engineering fatigue problem can be treated as a "time"-to-event survival analysis problem if the number of cycles formula_84 is treated as a fictitious time variable  formula_87 This would facilitate the application of well-established survival analysis techniques to engineering fatigue problems and Tabatabai et al. The survival function formula_88, probability density function formula_89, hazard rate formula_90, and cumulative probability of failure formula_91 can then be defined as: The hypertabastic accelerated failure time model was used to analyze probabilistic fatigue life for various detailed categories in steel bridges.