Displaying 1 result from an estimated 1 matches for "ssbilin".
Did you mean:
sibiling
2007 Apr 20
1
nlme trouble
...mCall[["x"]], LHS, data)
x <- xy[["x"]]
x0 <- (min(x)+max(x))/2 + nu*((log(a, exp(1)) - log(b, exp(1))))/(a+b)
#log defaults to base e
value = c(nu,a,b,c,x0)
names(value)=mCall[c("nu","a","b","c","x0")]
value
}
SSBiLin=selfStart(~-nu*log(exp(-a*(x-x0)/nu)+exp(b*(x-x0))/nu, exp(1))+c,
initial=BiLinInit0, parameters= c("nu","a","b","c","x0")) #log
defaults to base e, but I used exp(1) just to be sure (tried with just
log(x) as well)
ge = read.table(&qu...