Displaying 1 result from an estimated 1 matches for "pka1".
Did you mean:
pka
2010 Feb 14
0
gnls not optimizing values
...model to data using the gnls function,
it doesn't return optimized values of the model parameters. Instead
it just returns the exact same values I used as initial guesses, but
with standard errors calculated. Example is as follows:
two_site_mHb <- gnls(y ~ (CS_AH2 + CS_AH1*10^(n1*(x-pKa1)) + CS_A_*10^
((n1+1)*x-n1*pKa1-pKa2)) / (1 + 10^(n1*(x-pKa1)) + 10^((n1+1)*x-
n1*pKa1-pKa2)), start=list(CS_AH2=y[1], CS_A_=y[length(y)], CS_AH1=y
[length(y)/2], pKa1=3.3, pKa2=5.9, n1=1.5))
and the results are:
Generalized nonlinear least squares fit
Model: y ~ (CS_AH2 + CS_AH1 * 10^(n1 *...