R/survival.R

In dschnelldavis/survivalcalculator: Bankruptcy Survival Calculator

#' Survival
#' 
#' Bankruptcy survival function
#' 
#' @export
#' @param input. Required.
survival <- function(inputformula, inputvalues) {
	useformula <- as.formula(inputformula)
	usevalues <- as.data.frame(inputvalues)
	mylogit <- glm(useformula, data=successdata, family=binomial)
	newdata1 <- with(successdata, usevalues)
	newdata2 <- predict(mylogit, newdata=newdata1, type="response", se.fit=TRUE)
	efit = newdata2$fit
	lfit = efit - (newdata2$se.fit * 1.96)
	hfit = efit + (newdata2$se.fit * 1.96)
	list(efit=efit, hfit=hfit, lfit=lfit)
}