computeStdDevMonteCarlo

Return to directory of Todd's CASA extensions

For a linear relationship, this function computes Y and its uncertainty given a value of X and the uncertainties on the fitted slope and intercept that relate these two variables. This is a wrapper for the method of the same name in the linfit class.

Usage:

au.computeStdDevMonteCarlo(slope, slopeSigma, intercept, interceptSigma, x, trials=10000)
  • trials: the higher the number, the more accurate is the result (at the cost of time)

Example:

Find the uncertainties in the expected 12C/13C isotope ratio using the curve fit of Milam et al. 2005 for an object at a galactocentric radius of 4.79 kpc. There is also a convenient wrapper for this: au.coGradient(4.79).
CASA <2>: au.computeStdDevMonteCarlo(5.41, 1.07, 19.03, 7.9, 4.79)
44.844182 +- 9.520469 (scaled MAD = 9.568575)

-- ToddHunter - 2013-08-27
Topic revision: r2 - 2013-08-28, ToddHunter
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding NRAO Public Wiki? Send feedback