Displaying 2 results from an estimated 2 matches for "k_mod".
Did you mean:
s_mod
2011 Mar 15
1
Problem with nls.lm function of minpack.lm package.
...he 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, Thr, alp, n, L, ksat){
> Th_mod <- expression(Thr+(Ths-Thr)/(1+(alp*h)^(n-1/n))...
2006 Aug 26
1
problems with loop
...m[,1]))
{
if (tmp[k,1]==sim[n,1] && v[k]<=sim[n,2]){tmp[k,2]=1}
if(tmp[k,1]==sim[n,1] && v[k]>sim[n,2]) {tmp[k,2]=0}
}
}
return(tmp)
}
# The model to be fitted
Model_1=function(beta0_mod,mu_mod,k_mod,I_mod)
{
parms = c(beta0=beta0_mod, mu=mu_mod, k=k_mod)
my.atol = 1e-6
times1=c(0,0.002884615,0.003846154,0.005769231,0.009615385,0.01,0.019230769)
times2=c(0.028846154,0.038461538,0.057692308,
0.076923077,0.083333333,0.096153846)
times3=c(0.115384615,0.125,0.134615385,...