Displaying 3 results from an estimated 3 matches for "laveggio".
2009 Feb 10
2
plotting the result of a nonlinear regression
..., na.action=na.omit, lower=list(a=0, b=10), upper=list(a=100, b=100))
plotfit(r.PTG.V)
I tried to use the function plotfit on the result of the following for loop but I got an error message:
for (i in 1:length(formList))
{
resultList[[i]] <- nls(formList[[i]], data=subset(dati, Fiume=="Laveggio"), start=startList7[[i]],
nls.control(maxiter=1000, warnOnly=TRUE), algorithm='port', na.action=na.omit,lower=lowerList7[[i]],
upper=upperList7[[i]])
}
plotfit(resultList[[1]]
"Error in diff(as.numeric(y[ord])) :
(list) object cannot be coerced to type 'double...
2009 Feb 03
3
non linear regression with nls
...5: start=list(a=10, b=10, d=1), lower=list(a=0, b=0, d=-50), upper=list(a=1000, b=1000, d=50)
so far i manage to do non linear regression one at a time that is, using one function for one river) using nls(), for
example:
r.NT.lav<-nls(NT.N~ fz1(Portata, a,b), data=subset(dati,river.name=="Laveggio"), start=list(a=10, b=10), nls.control
(maxiter=200), algorithm='port', trace=TRUE, na.action=na.omit, lower=list(a=0, b=0), upper=list(a=1000, b=1000))
and then I get the results with summary() and the graph using plotfit()
but this will get extremly long since I have 12 rivers to an...
2009 Feb 10
3
summary of a list
...fz3(Portata,a, b, d, e),
PTG.P ~ fz4(Portata, a, b), PTG.P ~ fz5(Portata, a, b, d), PO4.P ~ fz1(Portata, a, b),
PO4.P ~ fz2(Portata, a, b), ...
And the loop I use is:
resultList <- list()
for (i in 1:length(formList))
{
resultList[[i]] <- nls(formList[[i]], data=subset(dati, Fiume=="Laveggio"), start=startList[[i]],
nls.control(maxiter=1000, warnOnly=TRUE), algorithm='port', na.action=na.omit,lower=lowerList[[i]],
upper=upperList[[i]])
}
When the computation ends I get 5 warning messages (one of false convergence, 4 of singular convergence:
In nls(formList[[i]],...