[SciPy-user] Announce: ErrorVal.py 1.0 for error-bounds calculations

Gary Ruben gazzar at email.com
Sun May 9 23:12:58 CDT 2004


I finally got around to putting my errorbar calculation module on my website.
<http://users.bigpond.net.au/gazzar/python.html>
This should be of interest to Numeric Python users who deal with calculation of error-bounds on experimental data.

Synopsis:

A module providing a Python number type class and helper functions to ease
the task of computing error bounds on experimental data values.

This module defines an abstract data type, Err, which carries a central/prime
value and upper and lower error bounds.
Helper functions are provided to allow Python Numeric rank-1 arrays of Err objects
to be constructed with ease and to apply Numeric Ufuncs.
Written under Python 2.3.3 and Numeric 23.0

Example of usage:

upperPressure = ArrayOfErr([909., 802., 677., 585., 560., 548.], 1.0)
lowerPressure = ArrayOfErr([144., 246., 378., 469., 493., 505.], 1.0)
pressureDiff = upperPressure - lowerPressure
V_RStandard = ArrayOfErr([2.016, 2.016, 2.020, 2.017, 2.021, 2.019], 0.001)
R = 100.2                                   # standard resistor value [Ohm]
I = V_RStandard / R                         # current [A]
V_RAB = ArrayOfErr([6.167, 6.168, 6.170, 6.160, (6.153, 0.02), (5.894, 0.01)], 0.002)
R_AB = V_RAB / I
logR_AB = ApplyUfuncToErr(R_AB, log10)      # This means log10(R_AB)
print PrimeVals(logR_AB)
print MinVals(logR_AB)
print MaxVals(logR_AB)


Enjoy,
Gary Ruben

-- 
___________________________________________________________
Sign-up for Ads Free at Mail.com
http://promo.mail.com/adsfreejump.htm



More information about the SciPy-user mailing list