load(file='../../BookCode/nlrwr/data/L.minor.rda')

# Plot the data to get an idea of what it looks like:
#
plot(L.minor$conc, L.minor$rate, xlab='conc', ylab='rate')
grid()

# Define the fitted function we will use:
#
MMfct = function(conc, K, Vm){
    Vm * conc/(K+conc)
}

# Guess at some initial values:
# 
K = 0.5*(min(L.minor$conc) + max(L.minor$conc))
Vm = max(L.minor$rate)

fit = nls(rate ~ MMfct(conc, K, Vm), data=L.minor, start=c(Vm=Vm, K=K))
print(summary(fit))