[R] gaussian error propagation
Martin Maechler
maechler at stat.math.ethz.ch
Tue Oct 19 10:11:51 CEST 2004
>>>>> "Bernhard" == Bernhard Reiter <ockham at gmx.net>
>>>>> on Tue, 19 Oct 2004 06:54:59 +0200 writes:
Bernhard> Dear list, I'm intending to use R (1.9.0 installed
Bernhard> on my linux machine) for my applied experimental
Bernhard> physics course, requiring it to do rather basic
Bernhard> statistical things (like calculation of mean,
Bernhard> std-deviation and some linear fitting). There is,
Bernhard> though, a task for which I haven't found a
Bernhard> solution for after some days of search: in
Bernhard> experiments involving multiple (supposedly)
Bernhard> independent measurements I am required to
Bernhard> calculate the total error according to Gaussian
Bernhard> error propagation law. [Sqrt (Sum of (squares of
Bernhard> (partials(mean) times corresponding std-dev) ) )]
Bernhard> After failing to find an appropriate function,
Bernhard> I've started trying to build my own (something
Bernhard> that takes a (mathematical) function, does the
Bernhard> partial derivations, takes the corresponding mean
Bernhard> values and sd's (i.e. from a data frame) as
Bernhard> arguments and does the remaining necessary
Bernhard> calculations. However, even my most basic
Bernhard> experiments with D() and derive() failed, though I
Bernhard> stuck to the derive() help page and the FAQ, but
Bernhard> (x and y unspecified,) typing
>> D(expression(x^2), "x")
Bernhard> yields 0
not for me, see below
Have you accidentally redefined D() to be something else?
> (dd <- D(expression(x^2), "x"))
2 * x
> str(dd)
language 2 * x
> x <- 1:3
> eval(dd)
[1] 2 4 6
>
Bernhard> And so does derive (that is, for the
Bernhard> gradients; this time x and y set to integer
Bernhard> values, as R complains about failing objects
Bernhard> otherwise), when, i.e., I do
>> eval(derive(expression(x*y), c("x","y")))
Bernhard> I'd be very glad if anybody could help me with
Bernhard> this rather strange behaviour and tell me if/what
Bernhard> I'm doing wrong and if maybe there's an even
Bernhard> simpler way to perform gaussian error propagation.
Bernhard> Greetings, Ockham
Bernhard> ______________________________________________
Bernhard> R-help at stat.math.ethz.ch mailing list
Bernhard> https://stat.ethz.ch/mailman/listinfo/r-help
Bernhard> PLEASE do read the posting guide!
Bernhard> http://www.R-project.org/posting-guide.html
More information about the R-help
mailing list