Re: [R] Uncertainty propagation

2010-09-27 Thread Maayt
Thanks for the help, I start to get reasonable errors on the model... I finally turned to the simpler lm() fitting. As my data from which I fit has only 8 points in each case, I guess it does not make much sense to downweight outliers and use rlm() in this case. -- View this message in contex

Re: [R] Uncertainty propagation

2010-09-26 Thread Ben Bolker
Maayt hotmail.com> writes: > I linearized my power relations en fitted them with a linear rlm() function. > When I re-sample my pairs from a bivariate normal distribution for my power > law what transformation do I need to apply a transformation to my covariance > (vcov) matrix to get back from m

Re: [R] Uncertainty propagation

2010-09-26 Thread Maayt
Thanks a lot for the help, I linearized my power relations en fitted them with a linear rlm() function. When I re-sample my pairs from a bivariate normal distribution for my power law what transformation do I need to apply a transformation to my covariance (vcov) matrix to get back from my linear

Re: [R] Uncertainty propagation

2010-09-25 Thread Ben Bolker
Maayt hotmail.com> writes: [snip] > My first intention was to use a kind of monte carlo routine and run the > model many times by changing the power law parameters. These power laws were > obtained by fitting data points under R. I thus have std error associated to > them: alpha (±da) * WaterHigh

[R] Uncertainty propagation

2010-09-25 Thread Maayt
I have a small model running under R. This is basically running various power-law relations on a variable (in this case water level in a river) changing spatially and through time. I'd like to include some kind of error propagation to this. My first intention was to use a kind of monte carlo routi