OP: You can do this pretty easily in R. Just write a function
f(CL, V, F, ka) = \alpha E_\text{iv}(CL, V) + (1-\alpha) E_\text{po}(CL, V, F, ka)
where E_\text{iv} is a function returning the sum of squared errors for the IV model given parameters CL, V, and E_\text{po} is the same thing for the...