need suggestion about building formual
Have you considered writing a function to do the complex math and
then calling nls referring to that function? Consider the following
(not tried):
expg <- function(a., x, G0, R, T){
exp((a.[1]+(a.[2]*x)^(1/2)-G0)/(R*T))
}
mynls<-nls(formula=y~expg(a.=c(a, b), x=x, G0=G0, R=R, T=T),
data=mydata,...)
If "nls" stops prematurely, I then write another function, "SSE" to
compute the sum of squares of deviations from y and then ask "optim" to
minimize "SSE", using "hessian=TRUE". If you try this and have trouble
making it work, please send another post.
spencer graves
Simple wrote:
Thanks for your kind respond. Although the answer didn't solve my question clearly,maybe I still not understand the art of R. I'm sorry that I had not talked the problem clearly, maybe a example with more detail will be suitable as suggested in the the posting guide. In function fitting program, such as Sigmaplot, a fitting formula, can be write in separate form: G=a+(b*x)^(1/2) k=exp((G-G0)/(R*T)) fit k to y of course,in R's nls, can write as: mynls<-nls(formula=y~exp((a+(b*x)^(1/2)-G0)/(R*T)),data=mydata,...) In this example, the formula is simple and acceptable. However, when the formula is more complexity,writing all in one formula,the readability will be damaged.So I'm looking for a way to write simple and readable code in this situation. Spencer Graves wrote:
I'm not certain what you are asking.
You can build expressions in R as character strings and then execute
them. Example:
expr <- paste("two <-", 1, "+", 1)
eval(parse(text=expr))
two
If this does not answer your question, PLEASE do read the posting
guide, "www.R-project.org/posting-guide.html". It can help increase the
chances of a quick and useful reply.
spencer graves
Simple wrote:
hi, I'm an newbie for R,I want do some fitting in R. I wander if it is possible to write a few of equations but only one formual when fitting Currently,My problem is,in R, is there methods combination a few equations into one formual? For example, y=f1(k); k=f2(t); t=f3(x); although it is certain that the can be equations turn into one formual as y~f(x),but write such a complexity string make me painful. I have searched the web and found out there were only examples with one formual.any suggestion? I hope that I have omit something.
------------------------------------------------------------------------
______________________________________________ R-help at stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA spencer.graves at pdf.com www.pdf.com <http://www.pdf.com> Tel: 408-938-4420 Fax: 408-280-7915