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
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
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
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
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) *
5 matches
Mail list logo