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
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
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
3 matches
Mail list logo