LuChenLab / SCAPE

Single Cell Alternative Polyadenylation using Expectation-maximization
MIT License
6 stars 2 forks source link

Error preparing #5

Closed cnk113 closed 2 years ago

cnk113 commented 2 years ago

Hello,

I'm trying to prepare my annotation, but there is an error:

pybedtools.helpers.BEDToolsError:
Command was:

    bedtools sort -i /tmp/pybedtools.acvu1r7m.tmp

Error message was:
Error: malformed BED entry at line 119. Start was greater than end. Exiting.

Any ideas on how to get through this?

zhou-ran commented 2 years ago

Hi Chang,

Try to replace bedtools to version 2.26.0 in current shell environment

Hope this helps

Ran

cnk113 commented 2 years ago

Hey Ran,

I changed the bedtools version and it stills gives the same error.

zhou-ran commented 2 years ago

Hi,

That's weird. SCAPE with bedtools (2.26.0) and pybedtools(0.8.1) work fine in different linux server.

Counld you send the detial of current environment to my email? Ranzhou1005@gmail.com

cnk113 commented 2 years ago

Thanks for all the help! One more thing, how is the intron output used in the recent update?

zhou-ran commented 2 years ago

Hi,

The intron file was used to call the intronic pA sites, and the old intronic bed was six columns format. I updated the code to output four columns format for keeping consistent with UTR output.

cnk113 commented 2 years ago

Sorry to reopen this issue, but I'm getting a lot divide by zero and invalid error in multiply errors during apamix.

/home/chang/SCAPE/apamix/EM.py:473: RuntimeWarning: divide by zero encountered in log
  log_zmat[:, k] = np.log(ws[k]) + \
/home/chang/SCAPE/apamix/apamix.py:339: RuntimeWarning: divide by zero encountered in log
  all_win_log_lik_mat[:, i] = np.log(
/home/chang/SCAPE/apamix/apamix.py:351: RuntimeWarning: divide by zero encountered in log
  all_win_log_lik_mat[:, i] = np.log(
/home/chang/SCAPE/apamix/apamix.py:368: RuntimeWarning: divide by zero encountered in log
  all_win_log_lik_mat[:, i] = np.log(

It's still continuing through these warnings, it's ok to ignore?

zhou-ran commented 2 years ago

Yes, it's OK

cnk113 commented 2 years ago

Thanks!