Displaying 1 result from an estimated 1 matches for "th_mod".
Did you mean:
chmod
2011 Mar 15
1
Problem with nls.lm function of minpack.lm package.
...>
> I have another question?
>
> It is possible to use two models with the same parameters (alp, n) for fitting Th~h and k~h
> simultaneously with following models
> I guess, but I dont know how I can do that
>
> f <- function(h, Th, k, Ths, Thr, alp, n, L, ksat){
> Th_mod <- expression(Thr+(Ths-Thr)/(1+(alp*h)^(n-1/n)))
> eval (Th_mod)
> k_mod <-
> expression(ksat+(((1+(alp*h)^n)-((alp*h)^(n-1)))^2)/((1+(alp*h)^n)^((n-1/n)*(L+2))))
> eval (k_mod)
> }
> #fonction d'aide pour le gradient analytique
> j <- function(h, Th, K, Ths, Th...