Open YingziZhang-github opened 1 year ago
Hi,
I have the same issue also, and I think I can not do the html report because of that. Did you succeed in solving the problem?
Aitor Zabala
Apologies for the late response. I hadn't noticed that GitHub wasn't set up to send me notifications for issues, so this issue remained unseen for far too long.
This is appears to be a bug with pycoMeth handling requests to compare on chromosomes for which there is no data in any one of the MetH5 files. I have just pushed a fix into the dev branch.
Until a new release can be published, a workaround would be to make sure you only run MethComp for chromosomes for which your M5 files actually have methylation data.
Use meth5 list_chunks -i INPUT_FILE.m5
to check what chromosomes are actually in your MetH5 files. Then:
-a
parameter). If you don't, MethComp will perform windowed testing on all chromosomes.Hope this helps. I'll keep this open until the fix is in a release.
Dear pycoMeth developer(s),
I am trying using pycometh Meth_Comp to do methylation differential analysis. My command was
pycometh Meth_Comp -i 1.m5 2.m5 3.m5 4.m5 5.m5 -f reference.fa -b output.bed
.It reported error as the following:
I thought it was because of no valid p-values can be used, so I adjusted
--pvalue_threshold
but it seemed not help. Would you please let me know what should I do to get it work?Thank you very much in advance.
Yingzi