[R] Problem fitting a non-linear regression model with nls

2010-01-13 Thread Nathalie Yauschew-Raguenes
Hi, I'm trying to make a regression of the form : formula - y ~ Asym_inf + Asym_sup * ( (1 / (1 + (n1 * (exp( (tmid1-x) / scal1) )^(1/n1) ) ) ) - (1 / (1 + (n2 * (exp( (tmid2-x) / scal2) )^(1/n2) ) ) ) ) which is a sum of the generalized logistic model proposed by richards. with data such

Re: [R] Problem fitting a non-linear regression model with nls

2010-01-13 Thread Gabor Grothendieck
You could try the brute force of nls2 package; however, note that you have 8 parameters and only 16 points so you might look for a more parsimonious model. Plotting it it seems somewhat gaussian in shape so: mod - nls(y ~ a * dnorm(x, b, c), start = c(a = mean(y)/dnorm(0, 0, sd(x)), b = mean(x),

Re: [R] Problem fitting a non-linear regression model with nls

2010-01-13 Thread Bert Gunter
My question is how could I estimate those initial values so that the nls fitting works. You can't. Your parameters are almost certainly nonidentifiable (which is what Gabor told you more gracefully). Just because you believe in a complex (often mechanistic) nonlinear model and have some data

Re: [R] Problem fitting a non-linear regression model with nls

2010-01-13 Thread Nathalie Yauschew-Raguenes
Actually, the data that I used are measurements of plant growth during an entire year.It is usual to model the growth with logistic models. I have already tried the simple logistic model (which works). But the problem is that with this model the inflexion point occurs half-way up or down the