Comments (4)
This is not a bug, exactly. The error occurs when the code attempts to evaluate the likelihood of the GP model and the covariance matrix is ill-conditioned (and hence cannot be inverted without large numerical errors). One easy way to fix that problem is to artificially increase the formal uncertainties in your input file. Can you try doing that and see if it fixes the problem?
from k2sc.
I changed the errors in flux values to 10*np.sqrt(flux) but still I can see the same problem.
Here are some lines of my ASCII file:
'''
Dates Cadences Flux1 Uncert1 Xpos Ypos Quality
3072.72870042 145045 595.76412697 244.08279885 610.79639 775.54480 524288
3072.74913296 145046 600.84801035 245.12201255 610.70644 775.56774 0
3072.76956550 145047 588.15472058 242.51901381 610.64158 775.58440 0
3072.78999815 145048 592.47314689 243.40771288 610.58304 775.62371 0
3072.81043059 145049 596.33847980 244.20042584 610.50611 775.64645 0
3072.83086313 145050 595.35398818 243.99876807 610.44960 775.67590 0
3072.85129567 145051 599.34298764 244.81482546 610.39479 775.70505 0
3072.87172821 145052 597.35251292 244.40796078 610.33183 775.74779 0
3072.89216075 145053 598.82663811 244.70934557 610.27831 775.79942 0
3072.91259329 145054 595.56503989 244.04201275 610.22121 775.84784 0
'''
I am not sure if I did the same thing as Suzanne has suggested above. Would you please suggest what else I can check.
from k2sc.
I've added a quick patch that should fix the issue for now. The neglnposterior method in Detrender now returns inf if the covariance matrix inversion fails.
from k2sc.
Hi,
I'm having the same linear algebra error when running k2sc on K2 C16 light curves. It has happened 8 times on the 585 light curves I've processed (so not a huge problem) but was wondering if there was a solution? Thanks!
from k2sc.
Related Issues (20)
- add default splits for campaign 8 HOT 1
- error: file './k2sc/bin/k2plotmast' does not exist HOT 3
- error while detrending the 1 light curve without MPI HOT 4
- k2plot not working with latest version
- .dat input file not recognized HOT 2
- FITSwriter fails on files input with AMCReader
- Add docstrings to standalone.py HOT 4
- bug in training set selection
- Pass on keyword arguments for outlier clipping in detrend()
- standalone.py longer finds the keplerid HOT 2
- Trouble finding george module
- Issue with numpy interface for k2sc HOT 1
- issue regarding running k2sc -- having trouble reproducing the k2sc team's results HOT 23
- sigma_clip returns three arrays HOT 1
- k2sc example fails to run HOT 4
- Text files not reading in under AMCReader with Python 3 HOT 1
- Text files not reading in under AMCReader with Python 3
- Two values called from sigma clipping, but the function outputs three HOT 1
- Linear algebra error crashes detrender HOT 2
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 k2sc.