source('read_neutropenia_data.R')

DF = read_neutropenia_data()

## Fit a first logistic regression model:
##
resp = cbind(DF$Neutropenia_cases, DF$Treatment_courses)
dose = DF$Dose

## logit is linear in dose:
m1 = glm(resp ~ dose, family=binomial)
print(summary(m1))

## logit is nonlinear in dose:
m2 = glm(resp ~ dose + I(dose^2), family=binomial)
print(summary(m2))

## logit is linear in dose with G-CSF:
m3 = glm(resp ~ dose + DF$Use_G_CSF, family=binomial)
print(summary(m3))

## Find the safe does:
##
betas = coefficients(m1)
safe_dose = (log(1/5) - betas[1])/betas[2]
print(sprintf('safe_dose= %f', safe_dose) )