search for: nnls2

...= 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 this b starting value, this is so much better (?!):\n")...