Comments (1)
I looked into this more carefully today. Shinsuke was completely right: the fractional Stokes errors were not being computed properly, and were not used properly through the rest of the program. I thought about it a bit, and there was also a somewhat deeper problem: the equation in the fractional Stokes errors was wrong. It assumed that the error in Stokes I was the channel error, but we're actually using a model of Stokes I, which presumably has much smaller errors.
I've changed it to ignore errors in Stokes I, and just normalize the Stokes Q/U errors by the Stokes I model. Quantifying the errors in the model will be difficult, and should be a very minor contribution for all sources bright enough to have potential detections in polarization. So the fractional Stokes errors will be slightly smaller than before.
I've also changed how the errors propagate through. Before, the actual Stokes Q/U errors were used, but the fractional Stokes were used in the RM synthesis, so there was a contradiction. I've changed it over to using the fractional errors to determine the weights and in the estimation of dFDF. This will down-weight channels with fainter Stokes I, as should be the case.
The changes in the error calculations ripple down through all the output values. For flat spectra nothing will change, but for steep spectra the output values will be a bit different, probably for the better.
from rm-tools.
Related Issues (20)
- Noise in Q and U, and MAD vs theoretical HOT 3
- Progressbar slows progress HOT 1
- Speedup loops with JIT HOT 1
- error message of "TypeError: Object of type float32 is not JSON serializable" on RMtools_1D/do_RMsynth_1D.py HOT 5
- A minor problem on RMtools_3D/mk_test_cube_data.py HOT 1
- RM-Clean does too many iterations
- FDF cubes should have properly set (degenerate) Stokes axis HOT 1
- QU-fitting Stokes I model failing for many double-component RM models (m11 and m4) HOT 4
- Suggest deleting swap file HOT 1
- Odd behaviour in peak outputs from RMsynth3d HOT 3
- Different formulas for rmtf FWHM HOT 1
- Peak fitting maps should output de-rotated Q,U
- peakfitcube tool fails with index error HOT 4
- rmtools_peakfitcube fails with all-NaN maps as output HOT 2
- error message 'ValueError: not enough values to unpack (ecpected 4, got 3)' in rmclean3d HOT 7
- QU-fitting PA uncertainty does not take angle wrapping into account HOT 5
- Release Numpy restriction? HOT 2
- Progress bar for do_fitIcube reverted to old style HOT 3
- Convert Stokes I model to same reference frequency as polarization HOT 1
- Deprecation and versioning notes HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from rm-tools.