Displaying 1 result from an estimated 1 matches for "sameformula".
Did you mean:
rateformula
2010 May 11
1
nls() and nls2() behavior?
...ctive.xs <- function( b, m1, m2, m3, x, ydebug ) {
tgts.mean= c(rep( m1, 7), rep( m2, 3), rep( m3, 5));
yhat= tgts.mean + b*x
#debug
cat("Parameters: b=", b, " m123=", m1, m2, m3, " leading to ",
sum( (ydebug-yhat)^2 ), "\n")
return( yhat );
}
sameformula= (y ~ objective.xs(b, m1, m2, m3, x, y))
cat("\nNLS2 Function --- what do you do?:\n");
print(nls2( sameformula, data=d, control= nls.control(tol=1e-12),
algorithm="brute-force",
start=list( b=0.3, m1=0.5, m2=0.5, m3=0.5, trace=TRUE ) ) );
cat("\n---but with t...