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

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

Re: [R] Uncertainty propagation

2010-09-26 Thread Ben Bolker
Maayt m.lupker at 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

[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

Re: [R] Uncertainty propagation

2010-09-25 Thread Ben Bolker
Maayt m.lupker at 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) *