phyloacc / PhyloAcc

PhyloAcc a software to detect the changes of conservation of a genomic region
GNU General Public License v3.0
26 stars 12 forks source link

PhyloAcc_gBGC output misses logBF1 and logBF2 #24

Open Tkastylevsky opened 3 years ago

Tkastylevsky commented 3 years ago

Dear author, i've been trying to use phyloAcc version corrected for gBGC, and as I was attempting to reproduce your graphs using your R codes, I noticed that the simulated data (and mine) do not output the same columns as the elem_lik.txt files you have for mammals and ratites results (ex ), and that in fact, the columns did not correspond to the ones you describe in your Readme_output.md file. I did not see if there was a specific parameter to add to the param_file in order to get this output. Do you know what I should do ? apart from the species names, this is the parameter file I use :

PHYTREE_FILE input/test_V1.mod SEG_FILE input/test_V1.bed ALIGN_FILE input/test_V1.fasta RESULT_FOLDER output PREFIX test_V1 BURNIN 1000 MCMC 3000 CHAIN 1 TARGETSPECIES OUTGROUP CONSERVE GAPCHAR - NUM_THREAD 4 VERBOSE 0

Best, Timothee