ANGSD: Analysis of next generation Sequencing Data
Latest tar.gz version is (0.938/0.939 on github), see Change_log for changes, and download it here.
User:Thorfinn: Difference between revisions
No edit summary |
No edit summary |
||
Line 48: | Line 48: | ||
# better screen output if unsorted filtering | # better screen output if unsorted filtering | ||
=Anders changes to version 583= | |||
1) | |||
Default minMaf set to -1 (instead of 0) | |||
2) | |||
added | |||
freqs->freq = new double[pars->numSites]; | |||
to phatLoop | |||
added the if part of | |||
if(freq->freq==NULL) | |||
freq->freq = new double[pars->numSites]; | |||
to postFreq and likeFreq funcitons | |||
to stop leak + be able to run phat without postFreq/likeFreq |
Revision as of 17:13, 26 February 2014
- Thetas,Tajima,Neutrality test needs documentation and examples files
- doMajorMinor 4 fjerner sites som ikke har major/minor som ref. total maerkeligt og fucked up
- user custom class missing
- -samglfclean not documented
- soap usage is not documented
- -sites does not work without all chromsomes being in the *sites file
- -sites gives wrong chromomosome in error
- phat does not work alone (-doMaf 8)
- make haploid 1dsfs
- fix wiki doc with -sites should contain -rf with curs from -sites
- describe -trim! does trim 3 remove the first and last 3 bases? or just 2?
- from morten shitoutput
h1 h2 h3 = 7 10 5 has less than 3 blocks. skipping h1 h2 h3 = 8 10 5 has less than 3 blocks. skipping h1 h2 h3 = 10 11 5 has less than 3 blocks. skipping h1 h2 h3 = 7 12 5 has less than 3 blocks. skipping h1 h2 h3 = 10 12 5 has less than 3 blocks. skipping h1 h2 h3 = 10 13 5 has less than 3 blocks. skipping h1 h2 h3 = 10 15 5 has less than 3 blocks. skipping h1 h2 h3 = 7 10 6 has less than 3 blocks. skipping h1 h2 h3 = 8 10 6 has less than 3 blocks. skipping h1 h2 h3 = 9 10 6 has less than 3 blocks. skipping h1…
- fix -doFasta for single chromosomes. put gzwrite && bufstr.l =0 whenever we change cur
- -doMaf 4 is not working if post is not external
- Regarding an open issue of printing %e instead o numbers, using the bombyx/lyca dataset:
- fix }else{
fprintf(stderr,"./a.out [make_bed|do_stat|validate_bed|print]\n"); } in thetaStat
- simplefy filereading. from glf files
- check 1offseted in -sites file
- add more info to -sites file output
- better screen output if unsorted filtering
Anders changes to version 583
1) Default minMaf set to -1 (instead of 0)
2) added
freqs->freq = new double[pars->numSites];
to phatLoop
added the if part of if(freq->freq==NULL)
freq->freq = new double[pars->numSites];
to postFreq and likeFreq funcitons to stop leak + be able to run phat without postFreq/likeFreq