Re: [R] function in nls argument

2008-05-09 Thread elnano
Find the data (data_nls.lm_moyano.txt) here: ftp://ftp.bgc-jena.mpg.de/pub/outgoing/fmoyano Katharine Mullen wrote: > > Thanks for the details - it sounds like a bug. You can either send me the > data in an email off-list or make it available on-line somewhere, so that > I and other people ca

Re: [R] function in nls argument

2008-05-09 Thread elnano
Thank you Katharine. I am certain nprint is affecting my solution. Let me know how I can send the data (~300Kb). The script I used it: ST1 <- ST04 SM1 <- SM08 SR1 <- SRch2 ST <- ST1 [!is.na(SR1)] SM <- SM1 [!is.na(SR1)] SR <- SR1 [!is.na(SR1)] q <- 0.90 p <- c("a"=-0.003, "b"=0.1

Re: [R] function in nls argument

2008-05-08 Thread elnano
I've basically solved the problem using the nls.lm function from the minpack.lm (thanks Katharine) with some modifications for ignoring residuals above a given percentile. This is to avoid the strong influence of points which push my modeled vs. measured values away from the 1:1 line. I based it o