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.13, c=0.50,
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 can
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
3 matches
Mail list logo